답안 #49663

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
49663 2018-06-01T18:43:54 Z gs14004 원 고르기 (APIO18_circle_selection) C++17
64 / 100
4000 ms 95164 KB
#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];
 
inline ll sq(ll x){ return x*x; }
inline ll dist(int a, int b){
	return sq(C[a].x-C[b].x)+sq(C[a].y-C[b].y);
}
inline 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)
 
inline 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;
}
 
inline int range(ll x){
	if(x>0) return x/t;
	return -((-x+t-1)/t);
}
inline int xrange(const vi &cell){
	return range(C[cell.front()].x);
}
inline int yrange(const vi &cell){
	return range(C[cell.front()].y);
}
inline 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){
				if(C[idx].kill!=0) continue;
				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]);
	}
	// slow!!!
	G=tmp;
	t/=2;
}
void process(vi &cell, int killer){
	vi tmp;
	for(int killed:cell){
		if(C[killed].kill==0 && hit(killed, killer))
			C[killed].kill=killer;
		if(C[killed].kill==0) tmp.push_back(killed);
	}
	cell=tmp;
}
 
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;
	}
	int pos=s1;
	while(pos<min(n1,s1+5) && lx<=xrange(G[pos]) && xrange(G[pos])<=hx){
		vector<vi> &X=G[pos];
		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++;
		}
		pos++;
	}
}
 
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();
	
	for(int i=1; i<=n; i++) cout<<C[i].kill<<' ';
	return 0;
}
# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 376 KB Output is correct
2 Correct 2 ms 472 KB Output is correct
3 Correct 3 ms 548 KB Output is correct
4 Correct 3 ms 552 KB Output is correct
5 Correct 2 ms 552 KB Output is correct
6 Correct 3 ms 552 KB Output is correct
7 Correct 2 ms 552 KB Output is correct
8 Correct 2 ms 552 KB Output is correct
9 Correct 3 ms 552 KB Output is correct
10 Correct 2 ms 572 KB Output is correct
11 Correct 3 ms 592 KB Output is correct
12 Correct 3 ms 640 KB Output is correct
13 Correct 4 ms 640 KB Output is correct
14 Correct 4 ms 640 KB Output is correct
15 Correct 3 ms 640 KB Output is correct
16 Correct 3 ms 700 KB Output is correct
17 Correct 3 ms 700 KB Output is correct
18 Correct 3 ms 700 KB Output is correct
19 Correct 5 ms 892 KB Output is correct
20 Correct 6 ms 892 KB Output is correct
21 Correct 6 ms 908 KB Output is correct
22 Correct 33 ms 1932 KB Output is correct
23 Correct 30 ms 1932 KB Output is correct
24 Correct 30 ms 1964 KB Output is correct
25 Correct 38 ms 2204 KB Output is correct
26 Correct 35 ms 2204 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 208 ms 15956 KB Output is correct
2 Correct 219 ms 17176 KB Output is correct
3 Correct 227 ms 17700 KB Output is correct
4 Correct 190 ms 17700 KB Output is correct
5 Correct 346 ms 18820 KB Output is correct
6 Correct 1169 ms 38340 KB Output is correct
7 Correct 476 ms 38340 KB Output is correct
8 Correct 529 ms 38340 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 38340 KB Output is correct
2 Correct 720 ms 38340 KB Output is correct
3 Correct 2352 ms 63140 KB Output is correct
4 Correct 2134 ms 63188 KB Output is correct
5 Correct 1485 ms 63188 KB Output is correct
6 Correct 616 ms 63188 KB Output is correct
7 Correct 476 ms 63188 KB Output is correct
8 Correct 63 ms 63188 KB Output is correct
9 Correct 2722 ms 63188 KB Output is correct
10 Correct 1440 ms 63188 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 2400 ms 68612 KB Output is correct
2 Execution timed out 4051 ms 95164 KB Time limit exceeded
3 Halted 0 ms 0 KB -
# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 376 KB Output is correct
2 Correct 2 ms 472 KB Output is correct
3 Correct 3 ms 548 KB Output is correct
4 Correct 3 ms 552 KB Output is correct
5 Correct 2 ms 552 KB Output is correct
6 Correct 3 ms 552 KB Output is correct
7 Correct 2 ms 552 KB Output is correct
8 Correct 2 ms 552 KB Output is correct
9 Correct 3 ms 552 KB Output is correct
10 Correct 2 ms 572 KB Output is correct
11 Correct 3 ms 592 KB Output is correct
12 Correct 3 ms 640 KB Output is correct
13 Correct 4 ms 640 KB Output is correct
14 Correct 4 ms 640 KB Output is correct
15 Correct 3 ms 640 KB Output is correct
16 Correct 3 ms 700 KB Output is correct
17 Correct 3 ms 700 KB Output is correct
18 Correct 3 ms 700 KB Output is correct
19 Correct 5 ms 892 KB Output is correct
20 Correct 6 ms 892 KB Output is correct
21 Correct 6 ms 908 KB Output is correct
22 Correct 33 ms 1932 KB Output is correct
23 Correct 30 ms 1932 KB Output is correct
24 Correct 30 ms 1964 KB Output is correct
25 Correct 38 ms 2204 KB Output is correct
26 Correct 35 ms 2204 KB Output is correct
27 Correct 9 ms 95164 KB Output is correct
28 Correct 13 ms 95164 KB Output is correct
29 Correct 9 ms 95164 KB Output is correct
30 Correct 57 ms 95164 KB Output is correct
31 Correct 59 ms 95164 KB Output is correct
32 Correct 56 ms 95164 KB Output is correct
33 Correct 75 ms 95164 KB Output is correct
34 Correct 72 ms 95164 KB Output is correct
35 Correct 69 ms 95164 KB Output is correct
36 Correct 884 ms 95164 KB Output is correct
37 Correct 766 ms 95164 KB Output is correct
38 Correct 643 ms 95164 KB Output is correct
39 Correct 404 ms 95164 KB Output is correct
40 Correct 471 ms 95164 KB Output is correct
41 Correct 425 ms 95164 KB Output is correct
42 Correct 283 ms 95164 KB Output is correct
43 Correct 705 ms 95164 KB Output is correct
44 Correct 909 ms 95164 KB Output is correct
45 Correct 728 ms 95164 KB Output is correct
46 Correct 704 ms 95164 KB Output is correct
47 Correct 568 ms 95164 KB Output is correct
48 Correct 534 ms 95164 KB Output is correct
49 Correct 569 ms 95164 KB Output is correct
50 Correct 593 ms 95164 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 376 KB Output is correct
2 Correct 2 ms 472 KB Output is correct
3 Correct 3 ms 548 KB Output is correct
4 Correct 3 ms 552 KB Output is correct
5 Correct 2 ms 552 KB Output is correct
6 Correct 3 ms 552 KB Output is correct
7 Correct 2 ms 552 KB Output is correct
8 Correct 2 ms 552 KB Output is correct
9 Correct 3 ms 552 KB Output is correct
10 Correct 2 ms 572 KB Output is correct
11 Correct 3 ms 592 KB Output is correct
12 Correct 3 ms 640 KB Output is correct
13 Correct 4 ms 640 KB Output is correct
14 Correct 4 ms 640 KB Output is correct
15 Correct 3 ms 640 KB Output is correct
16 Correct 3 ms 700 KB Output is correct
17 Correct 3 ms 700 KB Output is correct
18 Correct 3 ms 700 KB Output is correct
19 Correct 5 ms 892 KB Output is correct
20 Correct 6 ms 892 KB Output is correct
21 Correct 6 ms 908 KB Output is correct
22 Correct 33 ms 1932 KB Output is correct
23 Correct 30 ms 1932 KB Output is correct
24 Correct 30 ms 1964 KB Output is correct
25 Correct 38 ms 2204 KB Output is correct
26 Correct 35 ms 2204 KB Output is correct
27 Correct 208 ms 15956 KB Output is correct
28 Correct 219 ms 17176 KB Output is correct
29 Correct 227 ms 17700 KB Output is correct
30 Correct 190 ms 17700 KB Output is correct
31 Correct 346 ms 18820 KB Output is correct
32 Correct 1169 ms 38340 KB Output is correct
33 Correct 476 ms 38340 KB Output is correct
34 Correct 529 ms 38340 KB Output is correct
35 Correct 2 ms 38340 KB Output is correct
36 Correct 720 ms 38340 KB Output is correct
37 Correct 2352 ms 63140 KB Output is correct
38 Correct 2134 ms 63188 KB Output is correct
39 Correct 1485 ms 63188 KB Output is correct
40 Correct 616 ms 63188 KB Output is correct
41 Correct 476 ms 63188 KB Output is correct
42 Correct 63 ms 63188 KB Output is correct
43 Correct 2722 ms 63188 KB Output is correct
44 Correct 1440 ms 63188 KB Output is correct
45 Correct 2400 ms 68612 KB Output is correct
46 Execution timed out 4051 ms 95164 KB Time limit exceeded
47 Halted 0 ms 0 KB -