#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> vi;
const int MX=300010, inf=2e9;
struct CIRCLE {
int x, y, r, idx, kill;
} C[MX];
ll sq(ll x){ return x*x; }
ll dist(int a, int b){
return sq(C[a].x-C[b].x)+sq(C[a].y-C[b].y);
}
bool hit(int a, int b){
return dist(a,b)<=sq(C[a].r+C[b].r);
}
int n;
void input(){
cin>>n;
for(int i=1; i<=n; i++){
int x, y, r;
cin>>x>>y>>r;
C[i]={x,y,r,i,0};
}
}
vector<vector<vi>> G;
ll t=1<<30;
// [-t, 0), [0, t)
int sect(int idx, int x0, int y0){
// 1 3
// 0 2
int x=C[idx].x, y=C[idx].y;
if(x<x0)
if(y<y0) return 0;
else return 1;
else
if(y<y0) return 2;
else return 3;
}
int range(ll x){
if(x>0) return x/t;
return -((-x+t-1)/t);
}
int xrange(const vi &cell){
return range(C[cell.front()].x);
}
int yrange(const vi &cell){
return range(C[cell.front()].y);
}
int xrange(const vector<vi> &X){
return xrange(X.front());
}
void init_grid(){
vi V[4];
for(int i=1; i<=n; i++){
V[sect(i,0,0)].push_back(i);
}
vector<vi> X[2];
for(int k=0; k<4; k++)
if(!V[k].empty())
X[k/2].push_back(V[k]);
if(!X[0].empty()) G.push_back(X[0]);
if(!X[1].empty()) G.push_back(X[1]);
}
void resize_grid(){
vector<vector<vi>> tmp;
for(auto &X:G){
vector<vi> newX[2];
for(vi &cell:X){
vi V[4];
int x=xrange(cell)*t + t/2;
int y=yrange(cell)*t + t/2;
for(int idx:cell)
V[sect(idx,x,y)].push_back(idx);
for(int k=0; k<4; k++)
if(!V[k].empty())
newX[k/2].push_back(V[k]);
}
if(!newX[0].empty()) tmp.push_back(newX[0]);
if(!newX[1].empty()) tmp.push_back(newX[1]);
}
G=tmp;
t/=2;
}
void process(const vi &cell, int killer){
// cell을 업데이트 안해도 되나?
for(int killed:cell)
if(C[killed].kill==0 && hit(killed, killer))
C[killed].kill=killer;
}
void kill(int idx){
while(C[idx].r*2<t){
resize_grid();
}
int x=C[idx].x, y=C[idx].y;
ll lx=range(x)-2, hx=range(x)+3;
ll ly=range(y)-2, hy=range(y)+3;
int n1=G.size();
int s1=0, e1=n1-1;
while(s1<e1){
int m1=(s1+e1)/2;
if(xrange(G[m1])<lx) s1=m1+1;
else e1=m1;
}
while(s1<n1 && lx<=xrange(G[s1]) && xrange(G[s1])<=hx){
vector<vi> &X=G[s1];
int n2=X.size();
int s2=0, e2=n2-1;
while(s2<e2){
int m2=(s2+e2)/2;
if(yrange(X[m2])<ly) s2=m2+1;
else e2=m2;
}
while(s2<n2 && ly<=yrange(X[s2]) && yrange(X[s2])<=hy){
process(X[s2], idx);
s2++;
}
s1++;
}
}
void solve(){
CIRCLE D[MX];
for(int i=1; i<=n; i++) D[i]=C[i];
sort(D+1, D+n+1, [](CIRCLE &a, CIRCLE &b){
if(a.r==b.r) return a.idx<b.idx;
return a.r>b.r;
});
for(int i=1; i<=n; i++){
int idx=D[i].idx;
if(C[idx].kill!=0) continue;
kill(idx);
}
}
int main(){
ios::sync_with_stdio(0); cin.tie(0);
input();
init_grid();
solve();
sort(C+1, C+n+1, [](CIRCLE &a, CIRCLE &b){ return a.idx<b.idx; });
for(int i=1; i<=n; i++) cout<<C[i].kill<<' ';
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
376 KB |
Output is correct |
2 |
Correct |
2 ms |
488 KB |
Output is correct |
3 |
Correct |
2 ms |
488 KB |
Output is correct |
4 |
Correct |
2 ms |
488 KB |
Output is correct |
5 |
Correct |
3 ms |
540 KB |
Output is correct |
6 |
Correct |
2 ms |
564 KB |
Output is correct |
7 |
Correct |
2 ms |
564 KB |
Output is correct |
8 |
Correct |
2 ms |
564 KB |
Output is correct |
9 |
Correct |
2 ms |
580 KB |
Output is correct |
10 |
Correct |
3 ms |
580 KB |
Output is correct |
11 |
Correct |
4 ms |
612 KB |
Output is correct |
12 |
Correct |
5 ms |
676 KB |
Output is correct |
13 |
Correct |
4 ms |
676 KB |
Output is correct |
14 |
Correct |
3 ms |
676 KB |
Output is correct |
15 |
Correct |
3 ms |
676 KB |
Output is correct |
16 |
Correct |
5 ms |
740 KB |
Output is correct |
17 |
Correct |
4 ms |
740 KB |
Output is correct |
18 |
Correct |
3 ms |
740 KB |
Output is correct |
19 |
Correct |
6 ms |
868 KB |
Output is correct |
20 |
Correct |
9 ms |
868 KB |
Output is correct |
21 |
Correct |
11 ms |
1268 KB |
Output is correct |
22 |
Correct |
48 ms |
2180 KB |
Output is correct |
23 |
Correct |
55 ms |
2180 KB |
Output is correct |
24 |
Correct |
60 ms |
2224 KB |
Output is correct |
25 |
Correct |
61 ms |
2224 KB |
Output is correct |
26 |
Correct |
56 ms |
2224 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
229 ms |
18456 KB |
Output is correct |
2 |
Correct |
239 ms |
18456 KB |
Output is correct |
3 |
Correct |
228 ms |
18456 KB |
Output is correct |
4 |
Correct |
210 ms |
18456 KB |
Output is correct |
5 |
Correct |
711 ms |
44664 KB |
Output is correct |
6 |
Correct |
2166 ms |
95356 KB |
Output is correct |
7 |
Correct |
1474 ms |
95356 KB |
Output is correct |
8 |
Correct |
1556 ms |
95356 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
95356 KB |
Output is correct |
2 |
Correct |
2364 ms |
95356 KB |
Output is correct |
3 |
Execution timed out |
3062 ms |
95356 KB |
Time limit exceeded |
4 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
3047 ms |
95356 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
376 KB |
Output is correct |
2 |
Correct |
2 ms |
488 KB |
Output is correct |
3 |
Correct |
2 ms |
488 KB |
Output is correct |
4 |
Correct |
2 ms |
488 KB |
Output is correct |
5 |
Correct |
3 ms |
540 KB |
Output is correct |
6 |
Correct |
2 ms |
564 KB |
Output is correct |
7 |
Correct |
2 ms |
564 KB |
Output is correct |
8 |
Correct |
2 ms |
564 KB |
Output is correct |
9 |
Correct |
2 ms |
580 KB |
Output is correct |
10 |
Correct |
3 ms |
580 KB |
Output is correct |
11 |
Correct |
4 ms |
612 KB |
Output is correct |
12 |
Correct |
5 ms |
676 KB |
Output is correct |
13 |
Correct |
4 ms |
676 KB |
Output is correct |
14 |
Correct |
3 ms |
676 KB |
Output is correct |
15 |
Correct |
3 ms |
676 KB |
Output is correct |
16 |
Correct |
5 ms |
740 KB |
Output is correct |
17 |
Correct |
4 ms |
740 KB |
Output is correct |
18 |
Correct |
3 ms |
740 KB |
Output is correct |
19 |
Correct |
6 ms |
868 KB |
Output is correct |
20 |
Correct |
9 ms |
868 KB |
Output is correct |
21 |
Correct |
11 ms |
1268 KB |
Output is correct |
22 |
Correct |
48 ms |
2180 KB |
Output is correct |
23 |
Correct |
55 ms |
2180 KB |
Output is correct |
24 |
Correct |
60 ms |
2224 KB |
Output is correct |
25 |
Correct |
61 ms |
2224 KB |
Output is correct |
26 |
Correct |
56 ms |
2224 KB |
Output is correct |
27 |
Correct |
29 ms |
95356 KB |
Output is correct |
28 |
Correct |
9 ms |
95356 KB |
Output is correct |
29 |
Correct |
9 ms |
95356 KB |
Output is correct |
30 |
Correct |
113 ms |
95356 KB |
Output is correct |
31 |
Correct |
131 ms |
95356 KB |
Output is correct |
32 |
Correct |
111 ms |
95356 KB |
Output is correct |
33 |
Correct |
73 ms |
95356 KB |
Output is correct |
34 |
Correct |
76 ms |
95356 KB |
Output is correct |
35 |
Correct |
560 ms |
95356 KB |
Output is correct |
36 |
Correct |
1299 ms |
95356 KB |
Output is correct |
37 |
Correct |
2171 ms |
95356 KB |
Output is correct |
38 |
Correct |
1696 ms |
95356 KB |
Output is correct |
39 |
Correct |
443 ms |
95356 KB |
Output is correct |
40 |
Correct |
418 ms |
95356 KB |
Output is correct |
41 |
Correct |
410 ms |
95356 KB |
Output is correct |
42 |
Correct |
341 ms |
95356 KB |
Output is correct |
43 |
Correct |
700 ms |
95356 KB |
Output is correct |
44 |
Correct |
965 ms |
95356 KB |
Output is correct |
45 |
Correct |
729 ms |
95356 KB |
Output is correct |
46 |
Correct |
733 ms |
95356 KB |
Output is correct |
47 |
Correct |
756 ms |
95356 KB |
Output is correct |
48 |
Correct |
722 ms |
95356 KB |
Output is correct |
49 |
Correct |
726 ms |
95356 KB |
Output is correct |
50 |
Correct |
833 ms |
95356 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
376 KB |
Output is correct |
2 |
Correct |
2 ms |
488 KB |
Output is correct |
3 |
Correct |
2 ms |
488 KB |
Output is correct |
4 |
Correct |
2 ms |
488 KB |
Output is correct |
5 |
Correct |
3 ms |
540 KB |
Output is correct |
6 |
Correct |
2 ms |
564 KB |
Output is correct |
7 |
Correct |
2 ms |
564 KB |
Output is correct |
8 |
Correct |
2 ms |
564 KB |
Output is correct |
9 |
Correct |
2 ms |
580 KB |
Output is correct |
10 |
Correct |
3 ms |
580 KB |
Output is correct |
11 |
Correct |
4 ms |
612 KB |
Output is correct |
12 |
Correct |
5 ms |
676 KB |
Output is correct |
13 |
Correct |
4 ms |
676 KB |
Output is correct |
14 |
Correct |
3 ms |
676 KB |
Output is correct |
15 |
Correct |
3 ms |
676 KB |
Output is correct |
16 |
Correct |
5 ms |
740 KB |
Output is correct |
17 |
Correct |
4 ms |
740 KB |
Output is correct |
18 |
Correct |
3 ms |
740 KB |
Output is correct |
19 |
Correct |
6 ms |
868 KB |
Output is correct |
20 |
Correct |
9 ms |
868 KB |
Output is correct |
21 |
Correct |
11 ms |
1268 KB |
Output is correct |
22 |
Correct |
48 ms |
2180 KB |
Output is correct |
23 |
Correct |
55 ms |
2180 KB |
Output is correct |
24 |
Correct |
60 ms |
2224 KB |
Output is correct |
25 |
Correct |
61 ms |
2224 KB |
Output is correct |
26 |
Correct |
56 ms |
2224 KB |
Output is correct |
27 |
Correct |
229 ms |
18456 KB |
Output is correct |
28 |
Correct |
239 ms |
18456 KB |
Output is correct |
29 |
Correct |
228 ms |
18456 KB |
Output is correct |
30 |
Correct |
210 ms |
18456 KB |
Output is correct |
31 |
Correct |
711 ms |
44664 KB |
Output is correct |
32 |
Correct |
2166 ms |
95356 KB |
Output is correct |
33 |
Correct |
1474 ms |
95356 KB |
Output is correct |
34 |
Correct |
1556 ms |
95356 KB |
Output is correct |
35 |
Correct |
3 ms |
95356 KB |
Output is correct |
36 |
Correct |
2364 ms |
95356 KB |
Output is correct |
37 |
Execution timed out |
3062 ms |
95356 KB |
Time limit exceeded |
38 |
Halted |
0 ms |
0 KB |
- |