Submission #51098

#	Submission time	Handle	Problem	Language	Result	Execution time	Memory
51098	2018-06-16T04:06:55 Z	boook	Circle selection (APIO18_circle_selection)	C++14	57 / 100	3000 ms	63228 KB

circle_selection

/*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()){

Subtask #1

7.0 / 7.0

#	Verdict	Execution time	Memory	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

Subtask #2

12.0 / 12.0

#	Verdict	Execution time	Memory	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

Subtask #3

15.0 / 15.0

#	Verdict	Execution time	Memory	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

Subtask #4

23.0 / 23.0

#	Verdict	Execution time	Memory	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

Subtask #5

0 / 30.0

#	Verdict	Execution time	Memory	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	-

Subtask #6

0 / 13.0

#	Verdict	Execution time	Memory	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	-