#pragma GCC optimize("O3,unroll-loops")
#pragma GCC target("avx2,bmi,bmi2,lzcnt,popcnt")
#include <bits/stdc++.h>
using namespace std;
typedef pair<int, int> pii;
const int K = 2;
int N;
int X[303030];
int Y[303030];
int R[303030];
int XL[303030];
int XR[303030];
int ans[303030];
int cur;
char rr;
void rf(int &x) {
x = 0;
int sgn = 0;
while (rr < 48 && rr != '-') rr = getchar();
if (rr == '-') { sgn = 1; rr = getchar(); }
while (47 < rr) { x = (x << 3) + (x << 1) + (rr & 15); rr = getchar(); }
if (sgn) x = -x;
}
inline bool intersect(int i, int j) {
long long dx = X[i] - X[j], dy = Y[i] - Y[j], r = R[i] + R[j];
return dx * dx + dy * dy <= r * r;
}
inline void chk(int cand) {
int &t = ans[cur];
if (R[cand] < R[t] || (R[cand] == R[t] && cand > t)) return;
if (intersect(cur, cand)) t = cand;
}
struct SegTree {
struct Node {
set<pii> s;
};
int n, base;
vector<Node> T;
vector<int> L;
SegTree(int _n) : n(_n) {
for (base = 1; base < n; base <<= 1);
T.resize(base + base);
L.resize(base + base);
for (int i = base; i < base + base; i++) L[i] = i - base + 1;
for (int i = base - 1; i >= 1; i--) L[i] = L[i << 1];
}
void chk_node(int idx, int y, int r, int li) {
auto it = T[idx].s.lower_bound({y, 0});
auto tmp = it;
if (it != T[idx].s.end() && XL[it->second] > li) {
chk(it->second);
if (it->first < (long long)y + r + r) {
it++;
if (it != T[idx].s.end()) chk(it->second);
}
}
// for (int i = 0; i < K; i++) {
// if (it == T[idx].end()) break;
// if (XL[it->second] <= li) break;
// chk(it->second);
// if (it->first >= (long long)y + r + r) break;
// if (i < K - 1) it++;
// }
it = tmp;
if (it != T[idx].s.begin()) {
it--;
if (XL[it->second] > li) {
chk(it->second);
if (it->first > (long long)y - r - r && it != T[idx].s.begin()) {
chk(prev(it)->second);
}
}
}
// for (int i = 0; i < K; i++) {
// if (it == T[idx].begin()) break;
// it--;
// if (XL[it->second] <= li) break;
// chk(it->second);
// if (it->first <= (long long)y - r - r) break;
// }
}
void upd(int p, int q, pii x) {
p += base; q += base;
p--; q--;
while (p <= q) {
if (p & 1) T[p++].s.insert(x);
if (~q & 1) T[q--].s.insert(x);
p >>= 1; q >>= 1;
}
}
void get(int p, int y, int r, int li) {
p += base; p--;
while (p >= 1) {
if (L[p] <= li) break;
chk_node(p, y, r, li);
p >>= 1;
}
}
};
int main() {
rf(N);
vector<int> v;
for (int i = 1; i <= N; i++) {
rf(X[i]); rf(Y[i]); rf(R[i]);
XL[i] = X[i] - R[i], XR[i] = X[i] + R[i];
v.push_back(XL[i]);
v.push_back(XR[i]);
}
sort(v.begin(), v.end());
v.erase(unique(v.begin(), v.end()), v.end());
for (int i = 1; i <= N; i++) {
XL[i] = lower_bound(v.begin(), v.end(), XL[i]) - v.begin() + 1;
XR[i] = lower_bound(v.begin(), v.end(), XR[i]) - v.begin() + 1;
}
int M = v.size();
SegTree seg(M);
vector<int> idx(N);
iota(idx.begin(), idx.end(), 1);
sort(idx.begin(), idx.end(), [&](int a, int b) {
return pii(R[a], -a) > pii(R[b], -b);
});
for (int i : idx) {
cur = i;
ans[i] = i;
seg.get(XL[i], Y[i], R[i], 0);
seg.get(XR[i], Y[i], R[i], XL[i]);
if (ans[i] == i) seg.upd(XL[i], XR[i], {Y[i], i});
}
for (int i = 1; i <= N; i++) printf("%d ", ans[i]); puts("");
return 0;
}
Compilation message
circle_selection.cpp: In function 'int main()':
circle_selection.cpp:146:2: warning: this 'for' clause does not guard... [-Wmisleading-indentation]
146 | for (int i = 1; i <= N; i++) printf("%d ", ans[i]); puts("");
| ^~~
circle_selection.cpp:146:54: note: ...this statement, but the latter is misleadingly indented as if it were guarded by the 'for'
146 | for (int i = 1; i <= N; i++) printf("%d ", ans[i]); puts("");
| ^~~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
6488 KB |
Output is correct |
2 |
Correct |
1 ms |
6492 KB |
Output is correct |
3 |
Correct |
1 ms |
6492 KB |
Output is correct |
4 |
Correct |
1 ms |
6492 KB |
Output is correct |
5 |
Correct |
1 ms |
6492 KB |
Output is correct |
6 |
Correct |
1 ms |
6492 KB |
Output is correct |
7 |
Correct |
1 ms |
6492 KB |
Output is correct |
8 |
Correct |
1 ms |
6492 KB |
Output is correct |
9 |
Correct |
1 ms |
6492 KB |
Output is correct |
10 |
Correct |
1 ms |
6492 KB |
Output is correct |
11 |
Correct |
1 ms |
6492 KB |
Output is correct |
12 |
Correct |
1 ms |
6492 KB |
Output is correct |
13 |
Correct |
1 ms |
6492 KB |
Output is correct |
14 |
Correct |
1 ms |
6492 KB |
Output is correct |
15 |
Correct |
1 ms |
6488 KB |
Output is correct |
16 |
Correct |
2 ms |
6748 KB |
Output is correct |
17 |
Correct |
2 ms |
6748 KB |
Output is correct |
18 |
Correct |
1 ms |
6748 KB |
Output is correct |
19 |
Correct |
4 ms |
8536 KB |
Output is correct |
20 |
Correct |
4 ms |
8284 KB |
Output is correct |
21 |
Correct |
5 ms |
8284 KB |
Output is correct |
22 |
Correct |
6 ms |
8796 KB |
Output is correct |
23 |
Correct |
9 ms |
9548 KB |
Output is correct |
24 |
Correct |
5 ms |
8796 KB |
Output is correct |
25 |
Correct |
5 ms |
8540 KB |
Output is correct |
26 |
Correct |
5 ms |
8796 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
343 ms |
119640 KB |
Output is correct |
2 |
Correct |
356 ms |
120128 KB |
Output is correct |
3 |
Correct |
345 ms |
119992 KB |
Output is correct |
4 |
Correct |
351 ms |
121264 KB |
Output is correct |
5 |
Correct |
286 ms |
68724 KB |
Output is correct |
6 |
Correct |
346 ms |
131040 KB |
Output is correct |
7 |
Correct |
332 ms |
122940 KB |
Output is correct |
8 |
Correct |
339 ms |
125352 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
6492 KB |
Output is correct |
2 |
Correct |
236 ms |
57436 KB |
Output is correct |
3 |
Correct |
912 ms |
186036 KB |
Output is correct |
4 |
Correct |
884 ms |
185232 KB |
Output is correct |
5 |
Correct |
1141 ms |
183092 KB |
Output is correct |
6 |
Correct |
392 ms |
85952 KB |
Output is correct |
7 |
Correct |
124 ms |
44232 KB |
Output is correct |
8 |
Correct |
15 ms |
14552 KB |
Output is correct |
9 |
Correct |
993 ms |
188244 KB |
Output is correct |
10 |
Correct |
1446 ms |
179040 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
662 ms |
176580 KB |
Output is correct |
2 |
Correct |
296 ms |
142780 KB |
Output is correct |
3 |
Correct |
647 ms |
127540 KB |
Output is correct |
4 |
Correct |
344 ms |
145588 KB |
Output is correct |
5 |
Correct |
320 ms |
144056 KB |
Output is correct |
6 |
Correct |
556 ms |
121336 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
6488 KB |
Output is correct |
2 |
Correct |
1 ms |
6492 KB |
Output is correct |
3 |
Correct |
1 ms |
6492 KB |
Output is correct |
4 |
Correct |
1 ms |
6492 KB |
Output is correct |
5 |
Correct |
1 ms |
6492 KB |
Output is correct |
6 |
Correct |
1 ms |
6492 KB |
Output is correct |
7 |
Correct |
1 ms |
6492 KB |
Output is correct |
8 |
Correct |
1 ms |
6492 KB |
Output is correct |
9 |
Correct |
1 ms |
6492 KB |
Output is correct |
10 |
Correct |
1 ms |
6492 KB |
Output is correct |
11 |
Correct |
1 ms |
6492 KB |
Output is correct |
12 |
Correct |
1 ms |
6492 KB |
Output is correct |
13 |
Correct |
1 ms |
6492 KB |
Output is correct |
14 |
Correct |
1 ms |
6492 KB |
Output is correct |
15 |
Correct |
1 ms |
6488 KB |
Output is correct |
16 |
Correct |
2 ms |
6748 KB |
Output is correct |
17 |
Correct |
2 ms |
6748 KB |
Output is correct |
18 |
Correct |
1 ms |
6748 KB |
Output is correct |
19 |
Correct |
4 ms |
8536 KB |
Output is correct |
20 |
Correct |
4 ms |
8284 KB |
Output is correct |
21 |
Correct |
5 ms |
8284 KB |
Output is correct |
22 |
Correct |
6 ms |
8796 KB |
Output is correct |
23 |
Correct |
9 ms |
9548 KB |
Output is correct |
24 |
Correct |
5 ms |
8796 KB |
Output is correct |
25 |
Correct |
5 ms |
8540 KB |
Output is correct |
26 |
Correct |
5 ms |
8796 KB |
Output is correct |
27 |
Correct |
8 ms |
10072 KB |
Output is correct |
28 |
Correct |
8 ms |
10076 KB |
Output is correct |
29 |
Correct |
8 ms |
10076 KB |
Output is correct |
30 |
Correct |
12 ms |
11948 KB |
Output is correct |
31 |
Correct |
10 ms |
11100 KB |
Output is correct |
32 |
Correct |
10 ms |
11100 KB |
Output is correct |
33 |
Correct |
87 ms |
35172 KB |
Output is correct |
34 |
Correct |
94 ms |
35012 KB |
Output is correct |
35 |
Correct |
101 ms |
35224 KB |
Output is correct |
36 |
Correct |
136 ms |
47812 KB |
Output is correct |
37 |
Correct |
138 ms |
48152 KB |
Output is correct |
38 |
Correct |
162 ms |
50380 KB |
Output is correct |
39 |
Correct |
106 ms |
11460 KB |
Output is correct |
40 |
Correct |
110 ms |
11424 KB |
Output is correct |
41 |
Correct |
110 ms |
11456 KB |
Output is correct |
42 |
Correct |
115 ms |
11716 KB |
Output is correct |
43 |
Correct |
90 ms |
42080 KB |
Output is correct |
44 |
Correct |
88 ms |
41316 KB |
Output is correct |
45 |
Correct |
90 ms |
41828 KB |
Output is correct |
46 |
Correct |
90 ms |
41928 KB |
Output is correct |
47 |
Correct |
89 ms |
41152 KB |
Output is correct |
48 |
Correct |
90 ms |
42060 KB |
Output is correct |
49 |
Correct |
85 ms |
40892 KB |
Output is correct |
50 |
Correct |
89 ms |
41412 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
6488 KB |
Output is correct |
2 |
Correct |
1 ms |
6492 KB |
Output is correct |
3 |
Correct |
1 ms |
6492 KB |
Output is correct |
4 |
Correct |
1 ms |
6492 KB |
Output is correct |
5 |
Correct |
1 ms |
6492 KB |
Output is correct |
6 |
Correct |
1 ms |
6492 KB |
Output is correct |
7 |
Correct |
1 ms |
6492 KB |
Output is correct |
8 |
Correct |
1 ms |
6492 KB |
Output is correct |
9 |
Correct |
1 ms |
6492 KB |
Output is correct |
10 |
Correct |
1 ms |
6492 KB |
Output is correct |
11 |
Correct |
1 ms |
6492 KB |
Output is correct |
12 |
Correct |
1 ms |
6492 KB |
Output is correct |
13 |
Correct |
1 ms |
6492 KB |
Output is correct |
14 |
Correct |
1 ms |
6492 KB |
Output is correct |
15 |
Correct |
1 ms |
6488 KB |
Output is correct |
16 |
Correct |
2 ms |
6748 KB |
Output is correct |
17 |
Correct |
2 ms |
6748 KB |
Output is correct |
18 |
Correct |
1 ms |
6748 KB |
Output is correct |
19 |
Correct |
4 ms |
8536 KB |
Output is correct |
20 |
Correct |
4 ms |
8284 KB |
Output is correct |
21 |
Correct |
5 ms |
8284 KB |
Output is correct |
22 |
Correct |
6 ms |
8796 KB |
Output is correct |
23 |
Correct |
9 ms |
9548 KB |
Output is correct |
24 |
Correct |
5 ms |
8796 KB |
Output is correct |
25 |
Correct |
5 ms |
8540 KB |
Output is correct |
26 |
Correct |
5 ms |
8796 KB |
Output is correct |
27 |
Correct |
343 ms |
119640 KB |
Output is correct |
28 |
Correct |
356 ms |
120128 KB |
Output is correct |
29 |
Correct |
345 ms |
119992 KB |
Output is correct |
30 |
Correct |
351 ms |
121264 KB |
Output is correct |
31 |
Correct |
286 ms |
68724 KB |
Output is correct |
32 |
Correct |
346 ms |
131040 KB |
Output is correct |
33 |
Correct |
332 ms |
122940 KB |
Output is correct |
34 |
Correct |
339 ms |
125352 KB |
Output is correct |
35 |
Correct |
1 ms |
6492 KB |
Output is correct |
36 |
Correct |
236 ms |
57436 KB |
Output is correct |
37 |
Correct |
912 ms |
186036 KB |
Output is correct |
38 |
Correct |
884 ms |
185232 KB |
Output is correct |
39 |
Correct |
1141 ms |
183092 KB |
Output is correct |
40 |
Correct |
392 ms |
85952 KB |
Output is correct |
41 |
Correct |
124 ms |
44232 KB |
Output is correct |
42 |
Correct |
15 ms |
14552 KB |
Output is correct |
43 |
Correct |
993 ms |
188244 KB |
Output is correct |
44 |
Correct |
1446 ms |
179040 KB |
Output is correct |
45 |
Correct |
662 ms |
176580 KB |
Output is correct |
46 |
Correct |
296 ms |
142780 KB |
Output is correct |
47 |
Correct |
647 ms |
127540 KB |
Output is correct |
48 |
Correct |
344 ms |
145588 KB |
Output is correct |
49 |
Correct |
320 ms |
144056 KB |
Output is correct |
50 |
Correct |
556 ms |
121336 KB |
Output is correct |
51 |
Correct |
8 ms |
10072 KB |
Output is correct |
52 |
Correct |
8 ms |
10076 KB |
Output is correct |
53 |
Correct |
8 ms |
10076 KB |
Output is correct |
54 |
Correct |
12 ms |
11948 KB |
Output is correct |
55 |
Correct |
10 ms |
11100 KB |
Output is correct |
56 |
Correct |
10 ms |
11100 KB |
Output is correct |
57 |
Correct |
87 ms |
35172 KB |
Output is correct |
58 |
Correct |
94 ms |
35012 KB |
Output is correct |
59 |
Correct |
101 ms |
35224 KB |
Output is correct |
60 |
Correct |
136 ms |
47812 KB |
Output is correct |
61 |
Correct |
138 ms |
48152 KB |
Output is correct |
62 |
Correct |
162 ms |
50380 KB |
Output is correct |
63 |
Correct |
106 ms |
11460 KB |
Output is correct |
64 |
Correct |
110 ms |
11424 KB |
Output is correct |
65 |
Correct |
110 ms |
11456 KB |
Output is correct |
66 |
Correct |
115 ms |
11716 KB |
Output is correct |
67 |
Correct |
90 ms |
42080 KB |
Output is correct |
68 |
Correct |
88 ms |
41316 KB |
Output is correct |
69 |
Correct |
90 ms |
41828 KB |
Output is correct |
70 |
Correct |
90 ms |
41928 KB |
Output is correct |
71 |
Correct |
89 ms |
41152 KB |
Output is correct |
72 |
Correct |
90 ms |
42060 KB |
Output is correct |
73 |
Correct |
85 ms |
40892 KB |
Output is correct |
74 |
Correct |
89 ms |
41412 KB |
Output is correct |
75 |
Correct |
360 ms |
120400 KB |
Output is correct |
76 |
Correct |
358 ms |
121580 KB |
Output is correct |
77 |
Correct |
320 ms |
120704 KB |
Output is correct |
78 |
Correct |
340 ms |
120504 KB |
Output is correct |
79 |
Correct |
417 ms |
121112 KB |
Output is correct |
80 |
Correct |
347 ms |
120268 KB |
Output is correct |
81 |
Correct |
801 ms |
181428 KB |
Output is correct |
82 |
Correct |
664 ms |
172892 KB |
Output is correct |
83 |
Correct |
688 ms |
172980 KB |
Output is correct |
84 |
Correct |
1102 ms |
196364 KB |
Output is correct |
85 |
Correct |
761 ms |
177336 KB |
Output is correct |
86 |
Correct |
642 ms |
169396 KB |
Output is correct |
87 |
Correct |
2575 ms |
240372 KB |
Output is correct |
88 |
Correct |
388 ms |
22968 KB |
Output is correct |
89 |
Correct |
377 ms |
23656 KB |
Output is correct |
90 |
Correct |
380 ms |
22968 KB |
Output is correct |
91 |
Correct |
371 ms |
22068 KB |
Output is correct |
92 |
Correct |
389 ms |
22184 KB |
Output is correct |
93 |
Correct |
340 ms |
142264 KB |
Output is correct |
94 |
Correct |
2756 ms |
90196 KB |
Output is correct |
95 |
Correct |
337 ms |
140928 KB |
Output is correct |
96 |
Correct |
339 ms |
141492 KB |
Output is correct |
97 |
Correct |
2056 ms |
67148 KB |
Output is correct |
98 |
Correct |
365 ms |
132532 KB |
Output is correct |
99 |
Correct |
339 ms |
141080 KB |
Output is correct |
100 |
Correct |
312 ms |
142704 KB |
Output is correct |
101 |
Correct |
481 ms |
145588 KB |
Output is correct |
102 |
Correct |
347 ms |
142088 KB |
Output is correct |
103 |
Execution timed out |
3012 ms |
87576 KB |
Time limit exceeded |
104 |
Halted |
0 ms |
0 KB |
- |