#include <bits/stdc++.h>
using namespace std;
class SegmentTree {
public:
int sz;
vector<set<pair<int, int>>> tree;
SegmentTree() {}
SegmentTree(int sz) : sz(sz), tree(2 * sz) {}
void Add(int l, int r, int v, int id = -1) {
for (l += sz, r += sz + 1; l < r; l /= 2, r /= 2) {
if (l & 1) tree[l++].emplace(v, id);
if (r & 1) tree[--r].emplace(v, id);
}
}
int FindMax(int p, int x) { // Find maximum <= x
int res = 0;
for (p += sz; p > 0; p /= 2) {
if (!tree[p].empty()) {
auto it = tree[p].upper_bound({x, sz + 1});
if (it != begin(tree[p])) {
res = max(res, prev(it)->first);
}
}
}
return res;
}
int FindMin(int p, int x) { // Find minimum >= x
int res = sz - 1;
for (p += sz; p > 0; p /= 2) {
if (!tree[p].empty()) {
auto it = tree[p].lower_bound({x, -1});
if (it != end(tree[p])) {
res = min(res, it->first);
}
}
}
return res;
}
vector<int> FindIntersection(int p, int x, int y) {
vector<int> res;
for (p += sz; p > 0; p /= 2) {
if (tree[p].empty()) {
continue;
}
auto &s = tree[p];
auto it = s.lower_bound({x, -1});
while (it != end(s) && it->first <= y) {
res.emplace_back(it->second);
it = s.erase(it);
}
}
sort(begin(res), end(res));
res.resize(unique(begin(res), end(res)) - begin(res));
return res;
}
};
void Read(int &S, int &T, int &U, int &V, int &N,
vector<int> &A, vector<int> &B, vector<int> &C, vector<int> &D) {
cin >> S >> T >> U >> V >> N;
A.resize(N), B.resize(N), C.resize(N), D.resize(N);
vector<int> xcoord = {S, U};
vector<int> ycoord = {T, V};
for (int i = 0; i < N; i++) {
cin >> A[i] >> B[i] >> C[i] >> D[i];
xcoord.emplace_back(A[i]);
xcoord.emplace_back(B[i]);
ycoord.emplace_back(C[i]);
ycoord.emplace_back(D[i]);
}
sort(begin(xcoord), end(xcoord));
sort(begin(ycoord), end(ycoord));
xcoord.resize(unique(begin(xcoord), end(xcoord)) - begin(xcoord));
ycoord.resize(unique(begin(ycoord), end(ycoord)) - begin(ycoord));
const auto GetPosition = [&](const vector<int> &a, int x) {
return lower_bound(begin(a), end(a), x) - begin(a);
};
S = GetPosition(xcoord, S);
T = GetPosition(ycoord, T);
U = GetPosition(xcoord, U);
V = GetPosition(ycoord, V);
for (int i = 0; i < N; i++) {
A[i] = GetPosition(xcoord, A[i]);
B[i] = GetPosition(xcoord, B[i]);
C[i] = GetPosition(ycoord, C[i]);
D[i] = GetPosition(ycoord, D[i]);
}
}
int main() {
ios::sync_with_stdio(0);
cin.tie(0), cout.tie(0);
int S, T, U, V, N;
vector<int> A, B, C, D;
Read(S, T, U, V, N, A, B, C, D);
vector<SegmentTree> boundary(2, SegmentTree(2 * N + 2));
for (int i = 0; i < N; i++) {
boundary[0].Add(C[i] + 1, D[i] - 1, A[i]);
boundary[0].Add(C[i] + 1, D[i] - 1, B[i]);
boundary[1].Add(A[i] + 1, B[i] - 1, C[i]);
boundary[1].Add(A[i] + 1, B[i] - 1, D[i]);
}
vector<array<int, 4>> lines;
for (int i = 0; i < N; i++) {
lines.push_back({A[i], C[i], D[i], 0});
lines.push_back({B[i], C[i], D[i], 0});
lines.push_back({C[i], A[i], B[i], 1});
lines.push_back({D[i], A[i], B[i], 1});
}
lines.push_back({S, T, T, 0});
lines.push_back({T, S, S, 1});
lines.push_back({U, V, V, 0});
lines.push_back({V, U, U, 1});
for (auto &l : lines) {
l[1] = boundary[l[3] ^ 1].FindMax(l[0], l[1] - 1);
l[2] = boundary[l[3] ^ 1].FindMin(l[0], l[2] + 1);
}
vector<SegmentTree> intersect(2, SegmentTree(2 * N + 2));
for (int i = 0; i < (int) lines.size(); i++) {
intersect[lines[i][3]].Add(lines[i][1], lines[i][2], lines[i][0], i);
}
queue<int> q;
vector<int> dist(4 * N + 4, -1);
for (int i = 0; i < (int) lines.size(); i++) {
if (lines[i][3] == 0) { // vertical
if (lines[i][0] == S && lines[i][1] <= T && T <= lines[i][2]) {
q.emplace(i);
dist[i] = 1;
}
} else { // horizontal
if (lines[i][0] == T && lines[i][1] <= S && S <= lines[i][2]) {
q.emplace(i);
dist[i] = 1;
}
}
}
while (!q.empty()) {
int u = q.front();
q.pop();
auto l = lines[u];
auto adj = intersect[l[3] ^ 1].FindIntersection(l[0], l[1], l[2]);
for (auto v : adj) {
if (dist[v] == -1) {
dist[v] = dist[u] + 1;
q.emplace(v);
}
}
}
int ans = 1e9;
for (int i = 0; i < (int) lines.size(); i++) {
if (lines[i][3] == 0) { // vertical
if (lines[i][0] == U && lines[i][1] <= V && V <= lines[i][2]) {
ans = min(ans, dist[i]);
}
} else { // horizontal
if (lines[i][0] == V && lines[i][1] <= U && U <= lines[i][2]) {
ans = min(ans, dist[i]);
}
}
}
cout << ans << "\n";
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
384 KB |
Output is correct |
2 |
Correct |
0 ms |
384 KB |
Output is correct |
3 |
Correct |
1 ms |
384 KB |
Output is correct |
4 |
Correct |
1 ms |
512 KB |
Output is correct |
5 |
Correct |
13 ms |
2760 KB |
Output is correct |
6 |
Correct |
13 ms |
2756 KB |
Output is correct |
7 |
Correct |
12 ms |
2760 KB |
Output is correct |
8 |
Correct |
13 ms |
2884 KB |
Output is correct |
9 |
Correct |
14 ms |
2884 KB |
Output is correct |
10 |
Correct |
13 ms |
2884 KB |
Output is correct |
11 |
Correct |
13 ms |
2884 KB |
Output is correct |
12 |
Correct |
16 ms |
2756 KB |
Output is correct |
13 |
Correct |
13 ms |
2756 KB |
Output is correct |
14 |
Correct |
13 ms |
2884 KB |
Output is correct |
15 |
Correct |
5 ms |
1536 KB |
Output is correct |
16 |
Correct |
13 ms |
2920 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
384 KB |
Output is correct |
2 |
Correct |
0 ms |
384 KB |
Output is correct |
3 |
Correct |
1 ms |
384 KB |
Output is correct |
4 |
Correct |
1 ms |
512 KB |
Output is correct |
5 |
Correct |
13 ms |
2760 KB |
Output is correct |
6 |
Correct |
13 ms |
2756 KB |
Output is correct |
7 |
Correct |
12 ms |
2760 KB |
Output is correct |
8 |
Correct |
13 ms |
2884 KB |
Output is correct |
9 |
Correct |
14 ms |
2884 KB |
Output is correct |
10 |
Correct |
13 ms |
2884 KB |
Output is correct |
11 |
Correct |
13 ms |
2884 KB |
Output is correct |
12 |
Correct |
16 ms |
2756 KB |
Output is correct |
13 |
Correct |
13 ms |
2756 KB |
Output is correct |
14 |
Correct |
13 ms |
2884 KB |
Output is correct |
15 |
Correct |
5 ms |
1536 KB |
Output is correct |
16 |
Correct |
13 ms |
2920 KB |
Output is correct |
17 |
Correct |
16 ms |
3272 KB |
Output is correct |
18 |
Correct |
19 ms |
3176 KB |
Output is correct |
19 |
Correct |
15 ms |
3268 KB |
Output is correct |
20 |
Correct |
15 ms |
3268 KB |
Output is correct |
21 |
Correct |
18 ms |
3268 KB |
Output is correct |
22 |
Correct |
16 ms |
3396 KB |
Output is correct |
23 |
Correct |
15 ms |
3268 KB |
Output is correct |
24 |
Correct |
15 ms |
3268 KB |
Output is correct |
25 |
Correct |
15 ms |
3268 KB |
Output is correct |
26 |
Correct |
15 ms |
3268 KB |
Output is correct |
27 |
Correct |
5 ms |
1792 KB |
Output is correct |
28 |
Correct |
11 ms |
3160 KB |
Output is correct |
29 |
Correct |
11 ms |
3160 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
384 KB |
Output is correct |
2 |
Correct |
0 ms |
384 KB |
Output is correct |
3 |
Correct |
1 ms |
384 KB |
Output is correct |
4 |
Correct |
1 ms |
512 KB |
Output is correct |
5 |
Correct |
13 ms |
2760 KB |
Output is correct |
6 |
Correct |
13 ms |
2756 KB |
Output is correct |
7 |
Correct |
12 ms |
2760 KB |
Output is correct |
8 |
Correct |
13 ms |
2884 KB |
Output is correct |
9 |
Correct |
14 ms |
2884 KB |
Output is correct |
10 |
Correct |
13 ms |
2884 KB |
Output is correct |
11 |
Correct |
13 ms |
2884 KB |
Output is correct |
12 |
Correct |
16 ms |
2756 KB |
Output is correct |
13 |
Correct |
13 ms |
2756 KB |
Output is correct |
14 |
Correct |
13 ms |
2884 KB |
Output is correct |
15 |
Correct |
5 ms |
1536 KB |
Output is correct |
16 |
Correct |
13 ms |
2920 KB |
Output is correct |
17 |
Correct |
16 ms |
3272 KB |
Output is correct |
18 |
Correct |
19 ms |
3176 KB |
Output is correct |
19 |
Correct |
15 ms |
3268 KB |
Output is correct |
20 |
Correct |
15 ms |
3268 KB |
Output is correct |
21 |
Correct |
18 ms |
3268 KB |
Output is correct |
22 |
Correct |
16 ms |
3396 KB |
Output is correct |
23 |
Correct |
15 ms |
3268 KB |
Output is correct |
24 |
Correct |
15 ms |
3268 KB |
Output is correct |
25 |
Correct |
15 ms |
3268 KB |
Output is correct |
26 |
Correct |
15 ms |
3268 KB |
Output is correct |
27 |
Correct |
5 ms |
1792 KB |
Output is correct |
28 |
Correct |
11 ms |
3160 KB |
Output is correct |
29 |
Correct |
11 ms |
3160 KB |
Output is correct |
30 |
Correct |
5815 ms |
413968 KB |
Output is correct |
31 |
Correct |
6003 ms |
417272 KB |
Output is correct |
32 |
Correct |
6483 ms |
421208 KB |
Output is correct |
33 |
Correct |
6103 ms |
421576 KB |
Output is correct |
34 |
Correct |
5887 ms |
426312 KB |
Output is correct |
35 |
Correct |
6371 ms |
425820 KB |
Output is correct |
36 |
Correct |
6773 ms |
425428 KB |
Output is correct |
37 |
Correct |
6283 ms |
424836 KB |
Output is correct |
38 |
Correct |
6173 ms |
425904 KB |
Output is correct |
39 |
Correct |
6233 ms |
425012 KB |
Output is correct |
40 |
Correct |
4345 ms |
405180 KB |
Output is correct |
41 |
Correct |
4428 ms |
404904 KB |
Output is correct |
42 |
Correct |
4470 ms |
409100 KB |
Output is correct |
43 |
Correct |
4483 ms |
409116 KB |
Output is correct |
44 |
Correct |
4490 ms |
413260 KB |
Output is correct |
45 |
Correct |
4550 ms |
413244 KB |
Output is correct |
46 |
Correct |
4484 ms |
413232 KB |
Output is correct |
47 |
Correct |
4508 ms |
413024 KB |
Output is correct |
48 |
Correct |
4366 ms |
413024 KB |
Output is correct |
49 |
Correct |
4419 ms |
413204 KB |
Output is correct |
50 |
Correct |
11 ms |
3164 KB |
Output is correct |
51 |
Correct |
11 ms |
3160 KB |
Output is correct |
52 |
Correct |
11 ms |
3160 KB |
Output is correct |