답안 #51098

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
51098 2018-06-16T04:06:55 Z boook 원 고르기 (APIO18_circle_selection) C++14
57 / 100
3000 ms 63228 KB
/*input
11
9 9 2
13 2 1
11 8 2
3 3 2
3 12 1
12 14 1
9 8 5
2 8 2
5 2 1
14 4 2
14 14 1
*/
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
using namespace std;
using namespace __gnu_pbds;
#define REP(i,j,k)     for(int i = j ; i < k ; ++i)
#define RREP(i,j,k)    for(int i = j ; i >=k ; --i)
#define A    first
#define B    second
#define mp   make_pair
#define pb   emplace_back
#define PII pair<int , int>
#define MEM(i,j)   memset(i , j , sizeof i)
#define ALL(i)     i.begin() , i.end()
#define DBGG(i,j)     cout << i << " " << j << endl
#define DB4(i,j,k,l)  cout << i << " " << j << " " << k << " " << l << endl
#define IOS cin.tie(0) , cout.sync_with_stdio(0)
#define endl "\n"
///------------------------------------------------------------
#define MAX 300090
#define getchar getchar_unlocked
#define INF 2e9

gp_hash_table<int_fast64_t , int> cc;

vector<PII> v[MAX];
vector<int> las;
int_fast64_t R = 1e14;
int n , pos;
int x[MAX] , y[MAX] , r[MAX] , sol[MAX] , ans[MAX];
inline int_fast64_t trs(int a , int b){ return 1LL * a * 1000000009 + b; }
inline int cmp(int a , int b){ return mp(r[a] , -a) > mp(r[b] , -b); }
inline int tch(int a , int b){
	int_fast64_t val = 1LL * (x[a] - x[b]) * (x[a] - x[b]) + 1LL * (y[a] - y[b]) * (y[a] - y[b]);
	return (val <= 1LL * (r[a] + r[b]) * (r[a] + r[b]));
}
inline int rit(){
	int val = 0 , f = 1;
	char c = getchar();
	while(!isdigit(c)){
		if(c == '-') f = -1;
		c = getchar();
	}
	while(isdigit(c)){
		val = (val << 1) + (val << 3) + (c & 15);
		c = getchar();
	}
	return val * f;
}
int32_t main(){
	n = rit();
	REP(i , 1 , n + 1) x[i] = rit() , y[i] = rit() , r[i] = rit();
	// REP(i , 1 , n + 1) cin >> x[i] >> y[i] >> r[i];

	REP(i , 1 , n + 1) sol[i] = i;
	sort(sol + 1 , sol + 1 + n , cmp);
	// REP(idx , 1 , n + 1){ int i = sol[idx];
	// 	if(ans[i]) continue;
	// 	REP(j , 1 , n + 1) if(ans[j] == 0 && tch(i , j)) ans[j] = i;
	// }
	REP(i , 1 , n + 1) las.pb(i);
	int del = 1000000;
	REP(idx , 1 , n + 1){ int i = sol[idx];
		if(ans[i]) continue;
		// DBGG("solve = " , i);
		if(R / 5 > r[i] && del > 1000){
			del = 0;
			cc.clear();
			R = r[sol[1]];
			vector<int> nw;
			for(auto j : las){
				if(ans[j]) continue;
				nw.pb(j);
				int xx = x[j] / R , yy = y[j] / R;
				if(cc.find(trs(xx , yy)) == cc.end()){
					v[pos].clear();
					v[pos].pb(0 , 1);
					cc[trs(xx , yy)] = pos ++;
					if(pos == MAX) pos = 0;
				}
				v[cc[trs(xx , yy)]].pb(j , v[cc[trs(xx , yy)]].size() + 1);
			}
			swap(nw , las);
		}
		int xx = x[i] / R , yy = y[i] / R;
		REP(ii , xx - 2 , xx + 3) REP(jj , yy - 2 , yy + 3){
			if(cc.find(trs(ii , jj)) != cc.end()){
				int id = cc[trs(ii , jj)] , pre = 0;
				while(v[id][pre].B < v[id].size()){
					int now = v[id][pre].B;
					if(tch(i , v[id][now].A)){
						ans[v[id][now].A] = i;
						del ++;
						v[id][pre].B = v[id][now].B;
					}
					else pre = now;
				}
				// for(auto to : v[id]){
				// 	if(ans[to.A] == 0 && tch(i , to.A)) ans[to.A] = i;
				// }
			}
		}
	}
	REP(i , 1 , n + 1) printf("%d " , ans[i]);
	puts("");
    return 0;
}

// 7 2 7 4 5 6 7 7 4 7 6

Compilation message

circle_selection.cpp: In function 'int32_t main()':
circle_selection.cpp:102:24: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
     while(v[id][pre].B < v[id].size()){
# 결과 실행 시간 메모리 Grader output
1 Correct 8 ms 7416 KB Output is correct
2 Correct 8 ms 7416 KB Output is correct
3 Correct 7 ms 7432 KB Output is correct
4 Correct 8 ms 7508 KB Output is correct
5 Correct 8 ms 7528 KB Output is correct
6 Correct 8 ms 7612 KB Output is correct
7 Correct 7 ms 7612 KB Output is correct
8 Correct 9 ms 7612 KB Output is correct
9 Correct 8 ms 7612 KB Output is correct
10 Correct 7 ms 7612 KB Output is correct
11 Correct 8 ms 7612 KB Output is correct
12 Correct 8 ms 7612 KB Output is correct
13 Correct 8 ms 7612 KB Output is correct
14 Correct 7 ms 7612 KB Output is correct
15 Correct 8 ms 7612 KB Output is correct
16 Correct 10 ms 7728 KB Output is correct
17 Correct 9 ms 7728 KB Output is correct
18 Correct 10 ms 7812 KB Output is correct
19 Correct 10 ms 7932 KB Output is correct
20 Correct 11 ms 7932 KB Output is correct
21 Correct 11 ms 7948 KB Output is correct
22 Correct 14 ms 8604 KB Output is correct
23 Correct 13 ms 8608 KB Output is correct
24 Correct 13 ms 8620 KB Output is correct
25 Correct 15 ms 8620 KB Output is correct
26 Correct 16 ms 8636 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 158 ms 21128 KB Output is correct
2 Correct 153 ms 21128 KB Output is correct
3 Correct 152 ms 21128 KB Output is correct
4 Correct 182 ms 21544 KB Output is correct
5 Correct 214 ms 21544 KB Output is correct
6 Correct 488 ms 40212 KB Output is correct
7 Correct 202 ms 40212 KB Output is correct
8 Correct 277 ms 40212 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 8 ms 40212 KB Output is correct
2 Correct 243 ms 40212 KB Output is correct
3 Correct 966 ms 61160 KB Output is correct
4 Correct 928 ms 61160 KB Output is correct
5 Correct 1210 ms 61160 KB Output is correct
6 Correct 232 ms 61160 KB Output is correct
7 Correct 105 ms 61160 KB Output is correct
8 Correct 24 ms 61160 KB Output is correct
9 Correct 2000 ms 61160 KB Output is correct
10 Correct 1029 ms 61160 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 598 ms 63104 KB Output is correct
2 Correct 420 ms 63228 KB Output is correct
3 Correct 379 ms 63228 KB Output is correct
4 Correct 447 ms 63228 KB Output is correct
5 Correct 446 ms 63228 KB Output is correct
6 Correct 145 ms 63228 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 8 ms 7416 KB Output is correct
2 Correct 8 ms 7416 KB Output is correct
3 Correct 7 ms 7432 KB Output is correct
4 Correct 8 ms 7508 KB Output is correct
5 Correct 8 ms 7528 KB Output is correct
6 Correct 8 ms 7612 KB Output is correct
7 Correct 7 ms 7612 KB Output is correct
8 Correct 9 ms 7612 KB Output is correct
9 Correct 8 ms 7612 KB Output is correct
10 Correct 7 ms 7612 KB Output is correct
11 Correct 8 ms 7612 KB Output is correct
12 Correct 8 ms 7612 KB Output is correct
13 Correct 8 ms 7612 KB Output is correct
14 Correct 7 ms 7612 KB Output is correct
15 Correct 8 ms 7612 KB Output is correct
16 Correct 10 ms 7728 KB Output is correct
17 Correct 9 ms 7728 KB Output is correct
18 Correct 10 ms 7812 KB Output is correct
19 Correct 10 ms 7932 KB Output is correct
20 Correct 11 ms 7932 KB Output is correct
21 Correct 11 ms 7948 KB Output is correct
22 Correct 14 ms 8604 KB Output is correct
23 Correct 13 ms 8608 KB Output is correct
24 Correct 13 ms 8620 KB Output is correct
25 Correct 15 ms 8620 KB Output is correct
26 Correct 16 ms 8636 KB Output is correct
27 Correct 11 ms 63228 KB Output is correct
28 Correct 12 ms 63228 KB Output is correct
29 Correct 15 ms 63228 KB Output is correct
30 Correct 21 ms 63228 KB Output is correct
31 Correct 19 ms 63228 KB Output is correct
32 Correct 19 ms 63228 KB Output is correct
33 Correct 72 ms 63228 KB Output is correct
34 Correct 73 ms 63228 KB Output is correct
35 Correct 53 ms 63228 KB Output is correct
36 Correct 217 ms 63228 KB Output is correct
37 Correct 208 ms 63228 KB Output is correct
38 Correct 209 ms 63228 KB Output is correct
39 Execution timed out 3051 ms 63228 KB Time limit exceeded
40 Halted 0 ms 0 KB -
# 결과 실행 시간 메모리 Grader output
1 Correct 8 ms 7416 KB Output is correct
2 Correct 8 ms 7416 KB Output is correct
3 Correct 7 ms 7432 KB Output is correct
4 Correct 8 ms 7508 KB Output is correct
5 Correct 8 ms 7528 KB Output is correct
6 Correct 8 ms 7612 KB Output is correct
7 Correct 7 ms 7612 KB Output is correct
8 Correct 9 ms 7612 KB Output is correct
9 Correct 8 ms 7612 KB Output is correct
10 Correct 7 ms 7612 KB Output is correct
11 Correct 8 ms 7612 KB Output is correct
12 Correct 8 ms 7612 KB Output is correct
13 Correct 8 ms 7612 KB Output is correct
14 Correct 7 ms 7612 KB Output is correct
15 Correct 8 ms 7612 KB Output is correct
16 Correct 10 ms 7728 KB Output is correct
17 Correct 9 ms 7728 KB Output is correct
18 Correct 10 ms 7812 KB Output is correct
19 Correct 10 ms 7932 KB Output is correct
20 Correct 11 ms 7932 KB Output is correct
21 Correct 11 ms 7948 KB Output is correct
22 Correct 14 ms 8604 KB Output is correct
23 Correct 13 ms 8608 KB Output is correct
24 Correct 13 ms 8620 KB Output is correct
25 Correct 15 ms 8620 KB Output is correct
26 Correct 16 ms 8636 KB Output is correct
27 Correct 158 ms 21128 KB Output is correct
28 Correct 153 ms 21128 KB Output is correct
29 Correct 152 ms 21128 KB Output is correct
30 Correct 182 ms 21544 KB Output is correct
31 Correct 214 ms 21544 KB Output is correct
32 Correct 488 ms 40212 KB Output is correct
33 Correct 202 ms 40212 KB Output is correct
34 Correct 277 ms 40212 KB Output is correct
35 Correct 8 ms 40212 KB Output is correct
36 Correct 243 ms 40212 KB Output is correct
37 Correct 966 ms 61160 KB Output is correct
38 Correct 928 ms 61160 KB Output is correct
39 Correct 1210 ms 61160 KB Output is correct
40 Correct 232 ms 61160 KB Output is correct
41 Correct 105 ms 61160 KB Output is correct
42 Correct 24 ms 61160 KB Output is correct
43 Correct 2000 ms 61160 KB Output is correct
44 Correct 1029 ms 61160 KB Output is correct
45 Correct 598 ms 63104 KB Output is correct
46 Correct 420 ms 63228 KB Output is correct
47 Correct 379 ms 63228 KB Output is correct
48 Correct 447 ms 63228 KB Output is correct
49 Correct 446 ms 63228 KB Output is correct
50 Correct 145 ms 63228 KB Output is correct
51 Correct 11 ms 63228 KB Output is correct
52 Correct 12 ms 63228 KB Output is correct
53 Correct 15 ms 63228 KB Output is correct
54 Correct 21 ms 63228 KB Output is correct
55 Correct 19 ms 63228 KB Output is correct
56 Correct 19 ms 63228 KB Output is correct
57 Correct 72 ms 63228 KB Output is correct
58 Correct 73 ms 63228 KB Output is correct
59 Correct 53 ms 63228 KB Output is correct
60 Correct 217 ms 63228 KB Output is correct
61 Correct 208 ms 63228 KB Output is correct
62 Correct 209 ms 63228 KB Output is correct
63 Execution timed out 3051 ms 63228 KB Time limit exceeded
64 Halted 0 ms 0 KB -