#pragma GCC optimize("Ofast")
#include <bits/stdc++.h>
using namespace std;
#define cerr cout
#define F first
#define S second
#define FOR(i,a,b) for (auto i = (a); i <= (b); ++i)
#define NFOR(i,a,b) for(auto i = (a); i >= (b); --i)
#define all(x) (x).begin(), (x).end()
#define sz(x) int(x.size())
typedef long long ll; typedef pair <int, int> ii; typedef vector <int> vi; const int inf = 1e9 + 7;
string to_string(string s) { return '"' + s + '"';}
string to_string(char s) { return string(1, s);}
string to_string(const char* s) { return to_string((string) s);}
string to_string(bool b) { return (b ? "true" : "false");}
template <typename A> string to_string(A);
template <typename A, typename B>string to_string(pair<A, B> p) {return "(" + to_string(p.first) + ", " + to_string(p.second) + ")";}
template <typename A> string to_string(A v) {bool f = 1; string r = "{"; for (const auto &x : v) {if (!f)r += ", "; f = 0; r += to_string(x);} return r + "}";}
void debug_out() { cerr << endl; }
template <typename Head, typename... Tail> void debug_out(Head H, Tail... T) {cerr << " " << to_string(H); debug_out(T...);}
#define pr(...) cerr << "[" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__)
const int N = 3e5 + 10;
vi x, y, r;
struct node {
int l, r, sum;
vector<int> a;
node(): l(0), r(0), sum(0) {}
};
vector<node> T;
vi xc, yc;
vi par;
inline bool intersect(int i, int j) {
return (x[i] - x[j]) * 1LL * (x[i] - x[j]) + (y[i] - y[j]) * 1LL * (y[i] - y[j]) <= (r[i] + r[j]) * 1LL * (r[i] + r[j]);
}
inline int make_node() {
T.emplace_back();
return sz(T) - 1;
}
int build(int no, bool dim, int xl, int xr, int yl, int yr, const vi& ids) {
if (ids.empty()) return no;
if (no == 0) no = make_node();
bool g = 1;
FOR (i, 1, sz(ids) - 1) if (x[ids[i]] != x[ids[0]] or y[ids[i]] != y[ids[0]]) {g = 0; break;}
if (g) {
T[no].sum = sz(ids);
T[no].a = ids;
return no;
}
if (dim) {
int mid = (yl + yr) >> 1;
vi lids, rids;
for (int it : ids) {
if (y[it] <= yc[mid]) lids.emplace_back(it);
else rids.emplace_back(it);
}
int L = build(0, 0, xl, xr, yl, mid, lids);
int R = build(0, 0, xl, xr, mid + 1, yr, rids);
T[no].l = L, T[no].r = R;
T[no].sum = T[T[no].l].sum + T[T[no].r].sum;
} else {
int mid = (xl + xr) >> 1;
vi lids, rids;
for (int it : ids) {
if (x[it] <= xc[mid]) lids.emplace_back(it);
else rids.emplace_back(it);
}
int L = build(0, 1, xl, mid, yl, yr, lids);
int R = build(0, 1, mid + 1, xr, yl, yr, rids);
T[no].l = L, T[no].r = R;
T[no].sum = T[T[no].l].sum + T[T[no].r].sum;
}
return no;
}
int xu, xv, yu, yv, bc;
vi rem;
void get(int no, bool dim, int xl, int xr, int yl, int yr) {
if (T[no].sum == 0) return;
if (xr < xu or xl > xv) return;
if (yr < yu or yl > yv) return;
if (T[no].l == 0 and T[no].r == 0) {
for (int it: T[no].a) if (intersect(it, bc)) {
par[it] = bc;
rem.emplace_back(it);
T[no].sum--;
}
reverse(all(rem));
vi& a = T[no].a;
int j = 0;
FOR (i, 0, sz(a) - 1) {
if (rem.empty() or a[i] != rem.back()) {
a[j] = a[i];
++j;
} else {
rem.pop_back();
}
}
a.resize(j);
rem.clear();
return;
}
if (dim) {
int mid = (yl+yr)>>1;
get(T[no].l, 0, xl,xr,yl, mid);
get(T[no].r, 0, xl, xr, mid+1, yr);
T[no].sum = T[T[no].l].sum + T[T[no].r].sum;
} else {
int mid = (xl+xr)>>1;
get(T[no].l, 1, xl, mid, yl, yr);
get(T[no].r, 1, mid+1, xr, yl, yr);
T[no].sum = T[T[no].l].sum + T[T[no].r].sum;
}
}
int main()
{
ios::sync_with_stdio(0); cin.tie(0);
int n; cin >> n;
x.resize(n), y.resize(n), r.resize(n);
par = vi(n, -1);
FOR (i, 0, n - 1) {
cin >> x[i] >> y[i] >> r[i];
xc.emplace_back(x[i]);
yc.emplace_back(y[i]);
}
xc.emplace_back(-inf), xc.emplace_back(inf);
yc.emplace_back(-inf), yc.emplace_back(inf);
sort(all(xc)); xc.erase(unique(all(xc)), xc.end());
sort(all(yc)); yc.erase(unique(all(yc)), yc.end());
vi id(n); iota(all(id), 0);
sort(all(id), [&](int x, int y) {
if (r[x] == r[y]) return x < y;
return r[x] > r[y];
});
T.reserve(3e6);
T.emplace_back();
int root = build(0, 0, 0, sz(xc)-1, 0, sz(yc)-1, id);
for (int i : id) if (par[i] == -1) {
bc = i;
xu = max(x[i] - 2LL * r[i], (ll) - inf);
xu = lower_bound(all(xc), xu) - xc.begin();
xv = min(x[i] + 2LL * r[i], (ll)inf);
xv = upper_bound(all(xc), xv) - xc.begin()-1;
yu = max(y[i] - 2LL * r[i], (ll) - 1e9);
yu = lower_bound(all(yc), yu) - yc.begin();
yv = min(y[i] + 2LL * r[i], (ll)1e9);
yv = upper_bound(all(yc), yv) - yc.begin()-1;
get(root, 0, 0, sz(xc)-1, 0, sz(yc)-1);
}
FOR (i, 0, n - 1) cout << par[i] + 1 << " ";
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
2 ms |
376 KB |
Output is correct |
2 |
Correct |
2 ms |
376 KB |
Output is correct |
3 |
Correct |
2 ms |
376 KB |
Output is correct |
4 |
Correct |
2 ms |
376 KB |
Output is correct |
5 |
Correct |
2 ms |
380 KB |
Output is correct |
6 |
Correct |
2 ms |
376 KB |
Output is correct |
7 |
Correct |
2 ms |
348 KB |
Output is correct |
8 |
Correct |
2 ms |
376 KB |
Output is correct |
9 |
Correct |
2 ms |
376 KB |
Output is correct |
10 |
Correct |
2 ms |
376 KB |
Output is correct |
11 |
Correct |
2 ms |
376 KB |
Output is correct |
12 |
Correct |
2 ms |
376 KB |
Output is correct |
13 |
Correct |
1 ms |
376 KB |
Output is correct |
14 |
Correct |
2 ms |
376 KB |
Output is correct |
15 |
Correct |
2 ms |
376 KB |
Output is correct |
16 |
Correct |
3 ms |
504 KB |
Output is correct |
17 |
Correct |
3 ms |
504 KB |
Output is correct |
18 |
Correct |
3 ms |
504 KB |
Output is correct |
19 |
Correct |
8 ms |
1400 KB |
Output is correct |
20 |
Correct |
8 ms |
1400 KB |
Output is correct |
21 |
Correct |
8 ms |
1400 KB |
Output is correct |
22 |
Correct |
11 ms |
1272 KB |
Output is correct |
23 |
Correct |
12 ms |
1400 KB |
Output is correct |
24 |
Correct |
11 ms |
1272 KB |
Output is correct |
25 |
Correct |
11 ms |
1400 KB |
Output is correct |
26 |
Correct |
11 ms |
1400 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
522 ms |
61052 KB |
Output is correct |
2 |
Correct |
536 ms |
60500 KB |
Output is correct |
3 |
Correct |
788 ms |
60272 KB |
Output is correct |
4 |
Correct |
513 ms |
60776 KB |
Output is correct |
5 |
Correct |
488 ms |
45804 KB |
Output is correct |
6 |
Correct |
709 ms |
59952 KB |
Output is correct |
7 |
Correct |
577 ms |
58316 KB |
Output is correct |
8 |
Correct |
581 ms |
59592 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
2 ms |
376 KB |
Output is correct |
2 |
Correct |
293 ms |
17964 KB |
Output is correct |
3 |
Correct |
1134 ms |
52876 KB |
Output is correct |
4 |
Correct |
1102 ms |
52568 KB |
Output is correct |
5 |
Correct |
1138 ms |
63452 KB |
Output is correct |
6 |
Correct |
496 ms |
41700 KB |
Output is correct |
7 |
Correct |
250 ms |
24556 KB |
Output is correct |
8 |
Correct |
45 ms |
6392 KB |
Output is correct |
9 |
Correct |
1137 ms |
56900 KB |
Output is correct |
10 |
Correct |
915 ms |
65884 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1064 ms |
52732 KB |
Output is correct |
2 |
Correct |
929 ms |
52800 KB |
Output is correct |
3 |
Correct |
537 ms |
52136 KB |
Output is correct |
4 |
Correct |
962 ms |
52324 KB |
Output is correct |
5 |
Correct |
932 ms |
52596 KB |
Output is correct |
6 |
Correct |
466 ms |
51900 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
2 ms |
376 KB |
Output is correct |
2 |
Correct |
2 ms |
376 KB |
Output is correct |
3 |
Correct |
2 ms |
376 KB |
Output is correct |
4 |
Correct |
2 ms |
376 KB |
Output is correct |
5 |
Correct |
2 ms |
380 KB |
Output is correct |
6 |
Correct |
2 ms |
376 KB |
Output is correct |
7 |
Correct |
2 ms |
348 KB |
Output is correct |
8 |
Correct |
2 ms |
376 KB |
Output is correct |
9 |
Correct |
2 ms |
376 KB |
Output is correct |
10 |
Correct |
2 ms |
376 KB |
Output is correct |
11 |
Correct |
2 ms |
376 KB |
Output is correct |
12 |
Correct |
2 ms |
376 KB |
Output is correct |
13 |
Correct |
1 ms |
376 KB |
Output is correct |
14 |
Correct |
2 ms |
376 KB |
Output is correct |
15 |
Correct |
2 ms |
376 KB |
Output is correct |
16 |
Correct |
3 ms |
504 KB |
Output is correct |
17 |
Correct |
3 ms |
504 KB |
Output is correct |
18 |
Correct |
3 ms |
504 KB |
Output is correct |
19 |
Correct |
8 ms |
1400 KB |
Output is correct |
20 |
Correct |
8 ms |
1400 KB |
Output is correct |
21 |
Correct |
8 ms |
1400 KB |
Output is correct |
22 |
Correct |
11 ms |
1272 KB |
Output is correct |
23 |
Correct |
12 ms |
1400 KB |
Output is correct |
24 |
Correct |
11 ms |
1272 KB |
Output is correct |
25 |
Correct |
11 ms |
1400 KB |
Output is correct |
26 |
Correct |
11 ms |
1400 KB |
Output is correct |
27 |
Correct |
14 ms |
2296 KB |
Output is correct |
28 |
Correct |
14 ms |
2424 KB |
Output is correct |
29 |
Correct |
15 ms |
2424 KB |
Output is correct |
30 |
Correct |
22 ms |
2296 KB |
Output is correct |
31 |
Correct |
22 ms |
2300 KB |
Output is correct |
32 |
Correct |
21 ms |
2296 KB |
Output is correct |
33 |
Correct |
186 ms |
17960 KB |
Output is correct |
34 |
Correct |
158 ms |
18156 KB |
Output is correct |
35 |
Correct |
151 ms |
17980 KB |
Output is correct |
36 |
Correct |
304 ms |
18012 KB |
Output is correct |
37 |
Correct |
297 ms |
17972 KB |
Output is correct |
38 |
Correct |
296 ms |
18028 KB |
Output is correct |
39 |
Correct |
928 ms |
15568 KB |
Output is correct |
40 |
Correct |
988 ms |
15600 KB |
Output is correct |
41 |
Correct |
700 ms |
15600 KB |
Output is correct |
42 |
Correct |
157 ms |
16976 KB |
Output is correct |
43 |
Correct |
346 ms |
18284 KB |
Output is correct |
44 |
Correct |
334 ms |
19644 KB |
Output is correct |
45 |
Correct |
294 ms |
18608 KB |
Output is correct |
46 |
Correct |
300 ms |
19948 KB |
Output is correct |
47 |
Correct |
299 ms |
19180 KB |
Output is correct |
48 |
Correct |
282 ms |
20204 KB |
Output is correct |
49 |
Correct |
370 ms |
20076 KB |
Output is correct |
50 |
Correct |
313 ms |
20336 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
2 ms |
376 KB |
Output is correct |
2 |
Correct |
2 ms |
376 KB |
Output is correct |
3 |
Correct |
2 ms |
376 KB |
Output is correct |
4 |
Correct |
2 ms |
376 KB |
Output is correct |
5 |
Correct |
2 ms |
380 KB |
Output is correct |
6 |
Correct |
2 ms |
376 KB |
Output is correct |
7 |
Correct |
2 ms |
348 KB |
Output is correct |
8 |
Correct |
2 ms |
376 KB |
Output is correct |
9 |
Correct |
2 ms |
376 KB |
Output is correct |
10 |
Correct |
2 ms |
376 KB |
Output is correct |
11 |
Correct |
2 ms |
376 KB |
Output is correct |
12 |
Correct |
2 ms |
376 KB |
Output is correct |
13 |
Correct |
1 ms |
376 KB |
Output is correct |
14 |
Correct |
2 ms |
376 KB |
Output is correct |
15 |
Correct |
2 ms |
376 KB |
Output is correct |
16 |
Correct |
3 ms |
504 KB |
Output is correct |
17 |
Correct |
3 ms |
504 KB |
Output is correct |
18 |
Correct |
3 ms |
504 KB |
Output is correct |
19 |
Correct |
8 ms |
1400 KB |
Output is correct |
20 |
Correct |
8 ms |
1400 KB |
Output is correct |
21 |
Correct |
8 ms |
1400 KB |
Output is correct |
22 |
Correct |
11 ms |
1272 KB |
Output is correct |
23 |
Correct |
12 ms |
1400 KB |
Output is correct |
24 |
Correct |
11 ms |
1272 KB |
Output is correct |
25 |
Correct |
11 ms |
1400 KB |
Output is correct |
26 |
Correct |
11 ms |
1400 KB |
Output is correct |
27 |
Correct |
522 ms |
61052 KB |
Output is correct |
28 |
Correct |
536 ms |
60500 KB |
Output is correct |
29 |
Correct |
788 ms |
60272 KB |
Output is correct |
30 |
Correct |
513 ms |
60776 KB |
Output is correct |
31 |
Correct |
488 ms |
45804 KB |
Output is correct |
32 |
Correct |
709 ms |
59952 KB |
Output is correct |
33 |
Correct |
577 ms |
58316 KB |
Output is correct |
34 |
Correct |
581 ms |
59592 KB |
Output is correct |
35 |
Correct |
2 ms |
376 KB |
Output is correct |
36 |
Correct |
293 ms |
17964 KB |
Output is correct |
37 |
Correct |
1134 ms |
52876 KB |
Output is correct |
38 |
Correct |
1102 ms |
52568 KB |
Output is correct |
39 |
Correct |
1138 ms |
63452 KB |
Output is correct |
40 |
Correct |
496 ms |
41700 KB |
Output is correct |
41 |
Correct |
250 ms |
24556 KB |
Output is correct |
42 |
Correct |
45 ms |
6392 KB |
Output is correct |
43 |
Correct |
1137 ms |
56900 KB |
Output is correct |
44 |
Correct |
915 ms |
65884 KB |
Output is correct |
45 |
Correct |
1064 ms |
52732 KB |
Output is correct |
46 |
Correct |
929 ms |
52800 KB |
Output is correct |
47 |
Correct |
537 ms |
52136 KB |
Output is correct |
48 |
Correct |
962 ms |
52324 KB |
Output is correct |
49 |
Correct |
932 ms |
52596 KB |
Output is correct |
50 |
Correct |
466 ms |
51900 KB |
Output is correct |
51 |
Correct |
14 ms |
2296 KB |
Output is correct |
52 |
Correct |
14 ms |
2424 KB |
Output is correct |
53 |
Correct |
15 ms |
2424 KB |
Output is correct |
54 |
Correct |
22 ms |
2296 KB |
Output is correct |
55 |
Correct |
22 ms |
2300 KB |
Output is correct |
56 |
Correct |
21 ms |
2296 KB |
Output is correct |
57 |
Correct |
186 ms |
17960 KB |
Output is correct |
58 |
Correct |
158 ms |
18156 KB |
Output is correct |
59 |
Correct |
151 ms |
17980 KB |
Output is correct |
60 |
Correct |
304 ms |
18012 KB |
Output is correct |
61 |
Correct |
297 ms |
17972 KB |
Output is correct |
62 |
Correct |
296 ms |
18028 KB |
Output is correct |
63 |
Correct |
928 ms |
15568 KB |
Output is correct |
64 |
Correct |
988 ms |
15600 KB |
Output is correct |
65 |
Correct |
700 ms |
15600 KB |
Output is correct |
66 |
Correct |
157 ms |
16976 KB |
Output is correct |
67 |
Correct |
346 ms |
18284 KB |
Output is correct |
68 |
Correct |
334 ms |
19644 KB |
Output is correct |
69 |
Correct |
294 ms |
18608 KB |
Output is correct |
70 |
Correct |
300 ms |
19948 KB |
Output is correct |
71 |
Correct |
299 ms |
19180 KB |
Output is correct |
72 |
Correct |
282 ms |
20204 KB |
Output is correct |
73 |
Correct |
370 ms |
20076 KB |
Output is correct |
74 |
Correct |
313 ms |
20336 KB |
Output is correct |
75 |
Correct |
509 ms |
52200 KB |
Output is correct |
76 |
Correct |
504 ms |
52828 KB |
Output is correct |
77 |
Correct |
519 ms |
52884 KB |
Output is correct |
78 |
Correct |
580 ms |
52316 KB |
Output is correct |
79 |
Correct |
572 ms |
52956 KB |
Output is correct |
80 |
Correct |
511 ms |
52080 KB |
Output is correct |
81 |
Correct |
1115 ms |
52796 KB |
Output is correct |
82 |
Correct |
1099 ms |
52828 KB |
Output is correct |
83 |
Correct |
1077 ms |
52572 KB |
Output is correct |
84 |
Correct |
1082 ms |
52444 KB |
Output is correct |
85 |
Correct |
1088 ms |
53212 KB |
Output is correct |
86 |
Correct |
1106 ms |
52668 KB |
Output is correct |
87 |
Correct |
1195 ms |
52672 KB |
Output is correct |
88 |
Execution timed out |
3032 ms |
42592 KB |
Time limit exceeded |
89 |
Halted |
0 ms |
0 KB |
- |