Submission #51096

#	Submission time	Handle	Problem	Language	Result	Execution time	Memory
51096	2018-06-16T04:04:34 Z	boook	Circle selection (APIO18_circle_selection)	C++14	57 / 100	3000 ms	63168 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 - '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;
			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	9 ms	7416 KB	Output is correct
2	Correct	8 ms	7528 KB	Output is correct
3	Correct	8 ms	7528 KB	Output is correct
4	Correct	7 ms	7528 KB	Output is correct
5	Correct	7 ms	7528 KB	Output is correct
6	Correct	7 ms	7560 KB	Output is correct
7	Correct	7 ms	7560 KB	Output is correct
8	Correct	8 ms	7576 KB	Output is correct
9	Correct	8 ms	7576 KB	Output is correct
10	Correct	8 ms	7576 KB	Output is correct
11	Correct	7 ms	7624 KB	Output is correct
12	Correct	8 ms	7628 KB	Output is correct
13	Correct	7 ms	7712 KB	Output is correct
14	Correct	7 ms	7712 KB	Output is correct
15	Correct	7 ms	7712 KB	Output is correct
16	Correct	10 ms	7712 KB	Output is correct
17	Correct	10 ms	7712 KB	Output is correct
18	Correct	9 ms	7712 KB	Output is correct
19	Correct	9 ms	7884 KB	Output is correct
20	Correct	9 ms	7884 KB	Output is correct
21	Correct	10 ms	7900 KB	Output is correct
22	Correct	13 ms	8556 KB	Output is correct
23	Correct	14 ms	8692 KB	Output is correct
24	Correct	14 ms	8692 KB	Output is correct
25	Correct	13 ms	8692 KB	Output is correct
26	Correct	14 ms	8692 KB	Output is correct

Subtask #2

12.0 / 12.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	155 ms	21208 KB	Output is correct
2	Correct	173 ms	21208 KB	Output is correct
3	Correct	146 ms	21208 KB	Output is correct
4	Correct	160 ms	21468 KB	Output is correct
5	Correct	189 ms	21468 KB	Output is correct
6	Correct	557 ms	40288 KB	Output is correct
7	Correct	247 ms	40288 KB	Output is correct
8	Correct	287 ms	40288 KB	Output is correct

Subtask #3

15.0 / 15.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	9 ms	40288 KB	Output is correct
2	Correct	216 ms	40288 KB	Output is correct
3	Correct	933 ms	61120 KB	Output is correct
4	Correct	916 ms	61120 KB	Output is correct
5	Correct	1181 ms	61120 KB	Output is correct
6	Correct	244 ms	61120 KB	Output is correct
7	Correct	108 ms	61120 KB	Output is correct
8	Correct	23 ms	61120 KB	Output is correct
9	Correct	1831 ms	61120 KB	Output is correct
10	Correct	1045 ms	61120 KB	Output is correct

Subtask #4

23.0 / 23.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	526 ms	62952 KB	Output is correct
2	Correct	440 ms	63040 KB	Output is correct
3	Correct	469 ms	63040 KB	Output is correct
4	Correct	513 ms	63040 KB	Output is correct
5	Correct	433 ms	63168 KB	Output is correct
6	Correct	168 ms	63168 KB	Output is correct

Subtask #5

0 / 30.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	9 ms	7416 KB	Output is correct
2	Correct	8 ms	7528 KB	Output is correct
3	Correct	8 ms	7528 KB	Output is correct
4	Correct	7 ms	7528 KB	Output is correct
5	Correct	7 ms	7528 KB	Output is correct
6	Correct	7 ms	7560 KB	Output is correct
7	Correct	7 ms	7560 KB	Output is correct
8	Correct	8 ms	7576 KB	Output is correct
9	Correct	8 ms	7576 KB	Output is correct
10	Correct	8 ms	7576 KB	Output is correct
11	Correct	7 ms	7624 KB	Output is correct
12	Correct	8 ms	7628 KB	Output is correct
13	Correct	7 ms	7712 KB	Output is correct
14	Correct	7 ms	7712 KB	Output is correct
15	Correct	7 ms	7712 KB	Output is correct
16	Correct	10 ms	7712 KB	Output is correct
17	Correct	10 ms	7712 KB	Output is correct
18	Correct	9 ms	7712 KB	Output is correct
19	Correct	9 ms	7884 KB	Output is correct
20	Correct	9 ms	7884 KB	Output is correct
21	Correct	10 ms	7900 KB	Output is correct
22	Correct	13 ms	8556 KB	Output is correct
23	Correct	14 ms	8692 KB	Output is correct
24	Correct	14 ms	8692 KB	Output is correct
25	Correct	13 ms	8692 KB	Output is correct
26	Correct	14 ms	8692 KB	Output is correct
27	Correct	12 ms	63168 KB	Output is correct
28	Correct	11 ms	63168 KB	Output is correct
29	Correct	11 ms	63168 KB	Output is correct
30	Correct	21 ms	63168 KB	Output is correct
31	Correct	18 ms	63168 KB	Output is correct
32	Correct	18 ms	63168 KB	Output is correct
33	Correct	50 ms	63168 KB	Output is correct
34	Correct	60 ms	63168 KB	Output is correct
35	Correct	55 ms	63168 KB	Output is correct
36	Correct	235 ms	63168 KB	Output is correct
37	Correct	229 ms	63168 KB	Output is correct
38	Correct	251 ms	63168 KB	Output is correct
39	Execution timed out	3059 ms	63168 KB	Time limit exceeded
40	Halted	0 ms	0 KB	-

Subtask #6

0 / 13.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	9 ms	7416 KB	Output is correct
2	Correct	8 ms	7528 KB	Output is correct
3	Correct	8 ms	7528 KB	Output is correct
4	Correct	7 ms	7528 KB	Output is correct
5	Correct	7 ms	7528 KB	Output is correct
6	Correct	7 ms	7560 KB	Output is correct
7	Correct	7 ms	7560 KB	Output is correct
8	Correct	8 ms	7576 KB	Output is correct
9	Correct	8 ms	7576 KB	Output is correct
10	Correct	8 ms	7576 KB	Output is correct
11	Correct	7 ms	7624 KB	Output is correct
12	Correct	8 ms	7628 KB	Output is correct
13	Correct	7 ms	7712 KB	Output is correct
14	Correct	7 ms	7712 KB	Output is correct
15	Correct	7 ms	7712 KB	Output is correct
16	Correct	10 ms	7712 KB	Output is correct
17	Correct	10 ms	7712 KB	Output is correct
18	Correct	9 ms	7712 KB	Output is correct
19	Correct	9 ms	7884 KB	Output is correct
20	Correct	9 ms	7884 KB	Output is correct
21	Correct	10 ms	7900 KB	Output is correct
22	Correct	13 ms	8556 KB	Output is correct
23	Correct	14 ms	8692 KB	Output is correct
24	Correct	14 ms	8692 KB	Output is correct
25	Correct	13 ms	8692 KB	Output is correct
26	Correct	14 ms	8692 KB	Output is correct
27	Correct	155 ms	21208 KB	Output is correct
28	Correct	173 ms	21208 KB	Output is correct
29	Correct	146 ms	21208 KB	Output is correct
30	Correct	160 ms	21468 KB	Output is correct
31	Correct	189 ms	21468 KB	Output is correct
32	Correct	557 ms	40288 KB	Output is correct
33	Correct	247 ms	40288 KB	Output is correct
34	Correct	287 ms	40288 KB	Output is correct
35	Correct	9 ms	40288 KB	Output is correct
36	Correct	216 ms	40288 KB	Output is correct
37	Correct	933 ms	61120 KB	Output is correct
38	Correct	916 ms	61120 KB	Output is correct
39	Correct	1181 ms	61120 KB	Output is correct
40	Correct	244 ms	61120 KB	Output is correct
41	Correct	108 ms	61120 KB	Output is correct
42	Correct	23 ms	61120 KB	Output is correct
43	Correct	1831 ms	61120 KB	Output is correct
44	Correct	1045 ms	61120 KB	Output is correct
45	Correct	526 ms	62952 KB	Output is correct
46	Correct	440 ms	63040 KB	Output is correct
47	Correct	469 ms	63040 KB	Output is correct
48	Correct	513 ms	63040 KB	Output is correct
49	Correct	433 ms	63168 KB	Output is correct
50	Correct	168 ms	63168 KB	Output is correct
51	Correct	12 ms	63168 KB	Output is correct
52	Correct	11 ms	63168 KB	Output is correct
53	Correct	11 ms	63168 KB	Output is correct
54	Correct	21 ms	63168 KB	Output is correct
55	Correct	18 ms	63168 KB	Output is correct
56	Correct	18 ms	63168 KB	Output is correct
57	Correct	50 ms	63168 KB	Output is correct
58	Correct	60 ms	63168 KB	Output is correct
59	Correct	55 ms	63168 KB	Output is correct
60	Correct	235 ms	63168 KB	Output is correct
61	Correct	229 ms	63168 KB	Output is correct
62	Correct	251 ms	63168 KB	Output is correct
63	Execution timed out	3059 ms	63168 KB	Time limit exceeded
64	Halted	0 ms	0 KB	-