#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
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 (pii(R[cand], -cand) < pii(R[t], -t)) return;
if (intersect(cur, cand)) t = cand;
}
struct SegTree {
struct Node {
set<ll> 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];
}
inline ll cvt(pii x) { return (ll)x.first * N + x.second - 1; }
void chk_node(int idx, int y, int r) {
auto it = T[idx].s.lower_bound(cvt({y, 1}));
auto tmp = it;
for (int i = 0; i < K; i++) {
if (it == T[idx].s.end()) break;
chk(*it % N + 1);
if (*it / N >= (long long)y + r + r) break;
if (i < K - 1) it++;
}
it = tmp;
for (int i = 0; i < K; i++) {
if (it == T[idx].s.begin()) break;
it--;
chk(*it % N + 1);
if (*it / N <= (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(cvt(x));
if (~q & 1) T[q--].s.insert(cvt(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);
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:129:2: warning: this 'for' clause does not guard... [-Wmisleading-indentation]
129 | for (int i = 1; i <= N; i++) printf("%d ", ans[i]); puts("");
| ^~~
circle_selection.cpp:129:54: note: ...this statement, but the latter is misleadingly indented as if it were guarded by the 'for'
129 | 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 |
6488 KB |
Output is correct |
3 |
Correct |
1 ms |
6492 KB |
Output is correct |
4 |
Correct |
1 ms |
6488 KB |
Output is correct |
5 |
Correct |
1 ms |
6488 KB |
Output is correct |
6 |
Correct |
1 ms |
6488 KB |
Output is correct |
7 |
Incorrect |
1 ms |
6488 KB |
Output isn't correct |
8 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
388 ms |
121528 KB |
Output is correct |
2 |
Correct |
374 ms |
119828 KB |
Output is correct |
3 |
Correct |
395 ms |
120216 KB |
Output is correct |
4 |
Correct |
355 ms |
120964 KB |
Output is correct |
5 |
Correct |
323 ms |
68632 KB |
Output is correct |
6 |
Correct |
368 ms |
131244 KB |
Output is correct |
7 |
Correct |
371 ms |
122940 KB |
Output is correct |
8 |
Correct |
380 ms |
124600 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
1 ms |
6488 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
667 ms |
176896 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
6488 KB |
Output is correct |
2 |
Correct |
1 ms |
6488 KB |
Output is correct |
3 |
Correct |
1 ms |
6492 KB |
Output is correct |
4 |
Correct |
1 ms |
6488 KB |
Output is correct |
5 |
Correct |
1 ms |
6488 KB |
Output is correct |
6 |
Correct |
1 ms |
6488 KB |
Output is correct |
7 |
Incorrect |
1 ms |
6488 KB |
Output isn't correct |
8 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
6488 KB |
Output is correct |
2 |
Correct |
1 ms |
6488 KB |
Output is correct |
3 |
Correct |
1 ms |
6492 KB |
Output is correct |
4 |
Correct |
1 ms |
6488 KB |
Output is correct |
5 |
Correct |
1 ms |
6488 KB |
Output is correct |
6 |
Correct |
1 ms |
6488 KB |
Output is correct |
7 |
Incorrect |
1 ms |
6488 KB |
Output isn't correct |
8 |
Halted |
0 ms |
0 KB |
- |