Submission #51090

# Submission time Handle Problem Language Result Execution time Memory
51090 2018-06-16T03:47:56 Z boook Circle selection (APIO18_circle_selection) C++14
57 / 100
3000 ms 63180 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 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 * 1000000007 + 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 - '0';
		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;
			pos = 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 ++;
				}
				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:100:24: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
     while(v[id][pre].B < v[id].size()){
# Verdict Execution time Memory Grader output
1 Correct 7 ms 7416 KB Output is correct
2 Correct 8 ms 7512 KB Output is correct
3 Correct 8 ms 7512 KB Output is correct
4 Correct 7 ms 7512 KB Output is correct
5 Correct 8 ms 7520 KB Output is correct
6 Correct 8 ms 7520 KB Output is correct
7 Correct 8 ms 7520 KB Output is correct
8 Correct 7 ms 7568 KB Output is correct
9 Correct 8 ms 7612 KB Output is correct
10 Correct 7 ms 7616 KB Output is correct
11 Correct 8 ms 7616 KB Output is correct
12 Correct 8 ms 7616 KB Output is correct
13 Correct 8 ms 7616 KB Output is correct
14 Correct 8 ms 7616 KB Output is correct
15 Correct 10 ms 7628 KB Output is correct
16 Correct 10 ms 7628 KB Output is correct
17 Correct 10 ms 7628 KB Output is correct
18 Correct 10 ms 7632 KB Output is correct
19 Correct 13 ms 7888 KB Output is correct
20 Correct 11 ms 7888 KB Output is correct
21 Correct 11 ms 7904 KB Output is correct
22 Correct 16 ms 8544 KB Output is correct
23 Correct 19 ms 8544 KB Output is correct
24 Correct 19 ms 8616 KB Output is correct
25 Correct 15 ms 8616 KB Output is correct
26 Correct 19 ms 8616 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 233 ms 19732 KB Output is correct
2 Correct 233 ms 20024 KB Output is correct
3 Correct 233 ms 20216 KB Output is correct
4 Correct 214 ms 21468 KB Output is correct
5 Correct 203 ms 21468 KB Output is correct
6 Correct 577 ms 40284 KB Output is correct
7 Correct 246 ms 40284 KB Output is correct
8 Correct 330 ms 40284 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 9 ms 40284 KB Output is correct
2 Correct 250 ms 40284 KB Output is correct
3 Correct 1069 ms 61024 KB Output is correct
4 Correct 987 ms 61084 KB Output is correct
5 Correct 1379 ms 61084 KB Output is correct
6 Correct 288 ms 61084 KB Output is correct
7 Correct 151 ms 61084 KB Output is correct
8 Correct 40 ms 61084 KB Output is correct
9 Correct 2103 ms 61084 KB Output is correct
10 Correct 994 ms 61084 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 559 ms 63056 KB Output is correct
2 Correct 507 ms 63056 KB Output is correct
3 Correct 327 ms 63056 KB Output is correct
4 Correct 515 ms 63056 KB Output is correct
5 Correct 524 ms 63180 KB Output is correct
6 Correct 211 ms 63180 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 7 ms 7416 KB Output is correct
2 Correct 8 ms 7512 KB Output is correct
3 Correct 8 ms 7512 KB Output is correct
4 Correct 7 ms 7512 KB Output is correct
5 Correct 8 ms 7520 KB Output is correct
6 Correct 8 ms 7520 KB Output is correct
7 Correct 8 ms 7520 KB Output is correct
8 Correct 7 ms 7568 KB Output is correct
9 Correct 8 ms 7612 KB Output is correct
10 Correct 7 ms 7616 KB Output is correct
11 Correct 8 ms 7616 KB Output is correct
12 Correct 8 ms 7616 KB Output is correct
13 Correct 8 ms 7616 KB Output is correct
14 Correct 8 ms 7616 KB Output is correct
15 Correct 10 ms 7628 KB Output is correct
16 Correct 10 ms 7628 KB Output is correct
17 Correct 10 ms 7628 KB Output is correct
18 Correct 10 ms 7632 KB Output is correct
19 Correct 13 ms 7888 KB Output is correct
20 Correct 11 ms 7888 KB Output is correct
21 Correct 11 ms 7904 KB Output is correct
22 Correct 16 ms 8544 KB Output is correct
23 Correct 19 ms 8544 KB Output is correct
24 Correct 19 ms 8616 KB Output is correct
25 Correct 15 ms 8616 KB Output is correct
26 Correct 19 ms 8616 KB Output is correct
27 Correct 13 ms 63180 KB Output is correct
28 Correct 15 ms 63180 KB Output is correct
29 Correct 14 ms 63180 KB Output is correct
30 Correct 20 ms 63180 KB Output is correct
31 Correct 25 ms 63180 KB Output is correct
32 Correct 23 ms 63180 KB Output is correct
33 Correct 69 ms 63180 KB Output is correct
34 Correct 110 ms 63180 KB Output is correct
35 Correct 109 ms 63180 KB Output is correct
36 Correct 344 ms 63180 KB Output is correct
37 Correct 246 ms 63180 KB Output is correct
38 Correct 231 ms 63180 KB Output is correct
39 Execution timed out 3063 ms 63180 KB Time limit exceeded
40 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 7 ms 7416 KB Output is correct
2 Correct 8 ms 7512 KB Output is correct
3 Correct 8 ms 7512 KB Output is correct
4 Correct 7 ms 7512 KB Output is correct
5 Correct 8 ms 7520 KB Output is correct
6 Correct 8 ms 7520 KB Output is correct
7 Correct 8 ms 7520 KB Output is correct
8 Correct 7 ms 7568 KB Output is correct
9 Correct 8 ms 7612 KB Output is correct
10 Correct 7 ms 7616 KB Output is correct
11 Correct 8 ms 7616 KB Output is correct
12 Correct 8 ms 7616 KB Output is correct
13 Correct 8 ms 7616 KB Output is correct
14 Correct 8 ms 7616 KB Output is correct
15 Correct 10 ms 7628 KB Output is correct
16 Correct 10 ms 7628 KB Output is correct
17 Correct 10 ms 7628 KB Output is correct
18 Correct 10 ms 7632 KB Output is correct
19 Correct 13 ms 7888 KB Output is correct
20 Correct 11 ms 7888 KB Output is correct
21 Correct 11 ms 7904 KB Output is correct
22 Correct 16 ms 8544 KB Output is correct
23 Correct 19 ms 8544 KB Output is correct
24 Correct 19 ms 8616 KB Output is correct
25 Correct 15 ms 8616 KB Output is correct
26 Correct 19 ms 8616 KB Output is correct
27 Correct 233 ms 19732 KB Output is correct
28 Correct 233 ms 20024 KB Output is correct
29 Correct 233 ms 20216 KB Output is correct
30 Correct 214 ms 21468 KB Output is correct
31 Correct 203 ms 21468 KB Output is correct
32 Correct 577 ms 40284 KB Output is correct
33 Correct 246 ms 40284 KB Output is correct
34 Correct 330 ms 40284 KB Output is correct
35 Correct 9 ms 40284 KB Output is correct
36 Correct 250 ms 40284 KB Output is correct
37 Correct 1069 ms 61024 KB Output is correct
38 Correct 987 ms 61084 KB Output is correct
39 Correct 1379 ms 61084 KB Output is correct
40 Correct 288 ms 61084 KB Output is correct
41 Correct 151 ms 61084 KB Output is correct
42 Correct 40 ms 61084 KB Output is correct
43 Correct 2103 ms 61084 KB Output is correct
44 Correct 994 ms 61084 KB Output is correct
45 Correct 559 ms 63056 KB Output is correct
46 Correct 507 ms 63056 KB Output is correct
47 Correct 327 ms 63056 KB Output is correct
48 Correct 515 ms 63056 KB Output is correct
49 Correct 524 ms 63180 KB Output is correct
50 Correct 211 ms 63180 KB Output is correct
51 Correct 13 ms 63180 KB Output is correct
52 Correct 15 ms 63180 KB Output is correct
53 Correct 14 ms 63180 KB Output is correct
54 Correct 20 ms 63180 KB Output is correct
55 Correct 25 ms 63180 KB Output is correct
56 Correct 23 ms 63180 KB Output is correct
57 Correct 69 ms 63180 KB Output is correct
58 Correct 110 ms 63180 KB Output is correct
59 Correct 109 ms 63180 KB Output is correct
60 Correct 344 ms 63180 KB Output is correct
61 Correct 246 ms 63180 KB Output is correct
62 Correct 231 ms 63180 KB Output is correct
63 Execution timed out 3063 ms 63180 KB Time limit exceeded
64 Halted 0 ms 0 KB -