Submission #109068

# Submission time Handle Problem Language Result Execution time Memory
109068 2019-05-04T07:32:45 Z ryansee Circle selection (APIO18_circle_selection) C++14
72 / 100
3000 ms 97284 KB
#pragma GCC target ("avx2")
#pragma GCC optimization ("Ofast")
#pragma GCC optimization ("unroll-loops")

#include "bits/stdc++.h"
using namespace std;
#define FAST ios_base::sync_with_stdio(false); cin.tie(0);
#define LLINF ((long long) 1e18)//1234567890987654321
#define INF 1234567890ll
#define pb push_back
#define ins insert
#define f first
#define s second	
#define db 0
#define EPS (1e-7)    //0.0000001 the value
#define PI (acos(-1))
#define MAXN (300006)
#define MAXK 26
#define MAXX 15000006
#define ll long long int 
#define ld long double
mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());    //can be used by calling rng() or shuffle(A, A+n, rng)
#define FOR(ii, ss, ee) for(ll ii = ss; ii < ee; ++ii)
#define space " "
#define cbr cerr << "hi\n"
#define mmst(x, v) memset((x), v, sizeof ((x)))
#define siz(x) ((ll)x.size())
#define ph push
#define btinpct(x) __builtin_popcountll(x)
#define p2(x) (1LL<<(x))
#define all(x) (x).begin(), (x).end()
#define lbd(x, y) lower_bound(all(x), y)
#define ubd(x, y) upper_bound(all(x), y)
typedef pair <ll, int> pi;
typedef pair <ll, pi> spi;
typedef pair <pi, pi> dpi;
inline ll rand(ll x, ll y) { ++y; return (rng() % (y-x)) + x; } //inclusivesss
ll n; bool s2 = 1;
pi A[MAXN], C[MAXN];
bool done[MAXN];
int ans[MAXN];
#define sq(x) ((x)*(x))
ld sqr(ll x) { return (x ? sqrtl(x) : 0); }
ld mdist(ll a, ll b, ll c, ll d) { return sqr( sq(llabs(a-c)) + sq(llabs(b-d)) ); }
vector<pi>in;
bool s4 = 1; ll R = -1;
	#define sdb __attribute((optimize("-Ofast"), target("arch=sandybridge")))
ll magic = 5; // odd number pls
inline bool sdb inter(ll i, ll j) {
	if(done[j]) return 1;
	if(mdist(C[i].f, C[i].s, C[j].f, C[j].s) <= (ld)A[i].f + (ld) A[j].f) {
// 		assert(i < j);
		done[j] = 1;
		ans[j] = i;
		return 1;
	}
	else return 0;
}

unordered_map <int, int> mp;
vector <int> d;
vector<pi>ind[MAXN],tmp[MAXN];
int latest = 0;
void sdb split() { 
	FOR(i,0,latest+1) tmp[i].clear();
	ll new_R = R/2; // if(new_R%2)new_R++;
	ll co = 1;
	mp.clear();
	FOR(i,1,latest+1) {
		if(ind[i].empty()) { assert(0); continue; } assert(co < MAXN);
		if(ind[i].size()==1) { if(mp[C[ind[i][0].s].f/new_R] == 0) mp[C[ind[i][0].s].f/new_R] = co++; tmp[mp[C[ind[i][0].s].f/new_R]].pb(ind[i][0]); }
		else {
			ll one = LLINF;
			ll two = 0;
			for(auto j : ind[i]) two = max(two, C[j.s].f/new_R);
			for(auto j : ind[i]) one = min(one, C[j.s].f/new_R);
			ll three = -LLINF;
			for(auto j : ind[i]) three = ((C[j.s].f/new_R) != one && (C[j.s].f/new_R) != two) ? C[j.s].f/new_R : three;
			// if(three != -LLINF) {
				// for(auto j : ind[i]) {
					// cer
					// ll cmp = C[j.s].f/new_R;
					// if(one == cmp || two == cmp || three == cmp);
					// else {
						// cout << "BRUHHH\n";						
						// exit(0);
					// }
				// }
			// }
			// cerr << "1,2,3: " << one << ' ' << two << ' ' << three << '\n';
			if(one == two) { // 1
				assert(three==-LLINF);
				if(mp[one]==0)mp[one] = co ++;
			} else if(three == -LLINF) { // 2
				if(mp[one]==0)mp[one] = co ++;
				if(mp[two]==0)mp[two] = co++;
			} else { // 3
				swap(two,three);
				if(mp[one]==0)mp[one] = co ++;
				if(mp[two]==0)mp[two] = co ++;
				if(mp[three]==0)mp[three] = co ++;
			}
			assert(co+3 < MAXN);
			for(auto j : ind[i]) {
				// cerr << C[j.s].f << ' ' << new_R << ' ' << C[j.s].f/new_R << ' ' << C[j.s].f/R << '\n';
					tmp[mp[C[j.s].f/new_R]].pb(j);
			}
			// cerr << '\n' << '\n';
		}
	}
	R=new_R; latest = co-1;
			assert(co < MAXN);
	FOR(i,0,latest+1) ind[i]=tmp[i];
	// ll p = -LLINF;
	// FOR(i,1,latest+1) {
		// assert(ind[i].size());
		// ll hmm = C[ind[i][0].s].f/R;
		// for(auto j : ind[i]) if(hmm != C[j.s].f/R) assert(0);
		// assert(hmm >= p);
		// p = hmm;
	// }
}
void sdb st4() {
	FOR(i,0,n) {
		d.pb(C[i].f/R);
	}
	sort(all(d)); d.resize(unique(all(d))-d.begin());
	ll co = 1; for(auto i : d) mp[i] = co++; 
	FOR(i,0,n) {
		ind[mp[C[i].f/R]].pb(pi(C[i].s,i));
		latest = max(latest, mp[C[i].f/R]);
	}
	// FOR(i,1,latest+1) sort(all(ind[i]));
	int p[n+5];
	for(ll i=0;i<n;i++) p[i]=i;
	sort(p,p+n,[](ll x, ll y){if(A[x].f==A[y].f)return A[x].s<A[y].s; else return A[x].f>A[y].f;});
// 	assert(ind[0].empty());
	FOR(iii,0,n) {
		ll i = p[iii]; // cerr << i+1 << ' ';
		if(done[i]) continue;
		done[i] = 1;
		ans[i] = i; if(A[i].f < R/2) { split(); }
		ll x = mp[C[i].f/R], y = C[i].s;
		ll ii = i;
		FOR(i,x-magic/2,x+magic/2+1) {
			if(i <= 0||i>latest)continue;
// 			ll d = llabs(x-i) * R + R;
// 			auto lower = ind[i].lower_bound(pi(sqr((4ll*sq(R))-sq(d))-y, 0));
// 			auto upper = (ind[i].upper_bound(pi(sqr((4ll*sq(R))-sq(d))+y, LLINF)));
// 			if(upper != ind[i].end() && pi(lower->f,lower->s) > pi(upper->f,upper->s)) continue;
// 			if(lower == ind[i].end()) continue;
			for(auto j = ind[i].begin(); j != ind[i].end();j++) {
				if(inter(ii,j->s)) {}
			}
		}
	}
	FOR(i,0,n) assert(done[i]);
	FOR(i,0,n) cout << ans[i] + 1 << ' ';
}
int sdb main()
{
	// freopen("int","r",stdin); // freopen("out","w",stdout);
	FAST
	cin >> n;
	FOR(i,0,n) {
		cin >> C[i].f >> C[i].s >> A[i].f; if(0)in.pb(pi(C[i].f,A[i].f));
		C[i].f+=1e9;
		if(C[i].s) s2 = 0;
		A[i].s = i;
		if(R == -1) R = A[i].f;
		if(R != A[i].f) s4 = 0;
		R=max(R,A[i].f);
	} // if(R%2)R++;
	// if(st1() != st2()) {
		
		// cerr << n << '\n';
		// for(auto i : in) cerr << i.f << ' ' << i.s << '\n';
		// assert(0);
	// }
	// assert(st1() == st2());
	st4();
}
// 1 10 1 4 5 6 7 8 4 10 6
// 1 2 1 4 5 6 7 8 4 2 6
/*


3
5 0 3
10 0 2
20 0 8

4
1 0 3
10 0 5
15 0 1
20 0 10


8
1 4
14 12
6 10
0 6
14 0
9 6
3 2
0 0
* 
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
*/

Compilation message

circle_selection.cpp:2:0: warning: ignoring #pragma GCC optimization [-Wunknown-pragmas]
 #pragma GCC optimization ("Ofast")
 
circle_selection.cpp:3:0: warning: ignoring #pragma GCC optimization [-Wunknown-pragmas]
 #pragma GCC optimization ("unroll-loops")
 
circle_selection.cpp: In function 'void st4()':
circle_selection.cpp:143:24: warning: unused variable 'y' [-Wunused-variable]
   ll x = mp[C[i].f/R], y = C[i].s;
                        ^
# Verdict Execution time Memory Grader output
1 Correct 18 ms 14464 KB Output is correct
2 Correct 19 ms 14464 KB Output is correct
3 Correct 16 ms 14464 KB Output is correct
4 Correct 17 ms 14464 KB Output is correct
5 Correct 16 ms 14464 KB Output is correct
6 Correct 18 ms 14464 KB Output is correct
7 Correct 19 ms 14464 KB Output is correct
8 Correct 16 ms 14464 KB Output is correct
9 Correct 15 ms 14464 KB Output is correct
10 Correct 17 ms 14464 KB Output is correct
11 Correct 15 ms 14464 KB Output is correct
12 Correct 16 ms 14464 KB Output is correct
13 Correct 17 ms 14464 KB Output is correct
14 Correct 19 ms 14464 KB Output is correct
15 Correct 20 ms 14464 KB Output is correct
16 Correct 18 ms 14592 KB Output is correct
17 Correct 17 ms 14592 KB Output is correct
18 Correct 17 ms 14592 KB Output is correct
19 Correct 19 ms 14948 KB Output is correct
20 Correct 20 ms 14976 KB Output is correct
21 Correct 27 ms 15228 KB Output is correct
22 Correct 46 ms 15360 KB Output is correct
23 Correct 40 ms 15812 KB Output is correct
24 Correct 35 ms 15352 KB Output is correct
25 Correct 42 ms 15232 KB Output is correct
26 Correct 37 ms 15232 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 323 ms 43000 KB Output is correct
2 Correct 340 ms 42236 KB Output is correct
3 Correct 271 ms 38120 KB Output is correct
4 Correct 329 ms 42468 KB Output is correct
5 Correct 816 ms 70176 KB Output is correct
6 Correct 2087 ms 69208 KB Output is correct
7 Correct 1280 ms 83696 KB Output is correct
8 Correct 1364 ms 76724 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 15 ms 14464 KB Output is correct
2 Correct 1388 ms 38340 KB Output is correct
3 Execution timed out 3030 ms 81004 KB Time limit exceeded
4 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1206 ms 40756 KB Output is correct
2 Correct 966 ms 51920 KB Output is correct
3 Correct 2473 ms 35844 KB Output is correct
4 Correct 1018 ms 51048 KB Output is correct
5 Correct 990 ms 51636 KB Output is correct
6 Correct 2059 ms 34636 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 18 ms 14464 KB Output is correct
2 Correct 19 ms 14464 KB Output is correct
3 Correct 16 ms 14464 KB Output is correct
4 Correct 17 ms 14464 KB Output is correct
5 Correct 16 ms 14464 KB Output is correct
6 Correct 18 ms 14464 KB Output is correct
7 Correct 19 ms 14464 KB Output is correct
8 Correct 16 ms 14464 KB Output is correct
9 Correct 15 ms 14464 KB Output is correct
10 Correct 17 ms 14464 KB Output is correct
11 Correct 15 ms 14464 KB Output is correct
12 Correct 16 ms 14464 KB Output is correct
13 Correct 17 ms 14464 KB Output is correct
14 Correct 19 ms 14464 KB Output is correct
15 Correct 20 ms 14464 KB Output is correct
16 Correct 18 ms 14592 KB Output is correct
17 Correct 17 ms 14592 KB Output is correct
18 Correct 17 ms 14592 KB Output is correct
19 Correct 19 ms 14948 KB Output is correct
20 Correct 20 ms 14976 KB Output is correct
21 Correct 27 ms 15228 KB Output is correct
22 Correct 46 ms 15360 KB Output is correct
23 Correct 40 ms 15812 KB Output is correct
24 Correct 35 ms 15352 KB Output is correct
25 Correct 42 ms 15232 KB Output is correct
26 Correct 37 ms 15232 KB Output is correct
27 Correct 29 ms 16248 KB Output is correct
28 Correct 23 ms 15440 KB Output is correct
29 Correct 22 ms 15352 KB Output is correct
30 Correct 77 ms 16632 KB Output is correct
31 Correct 52 ms 16220 KB Output is correct
32 Correct 60 ms 16248 KB Output is correct
33 Correct 104 ms 21316 KB Output is correct
34 Correct 124 ms 21388 KB Output is correct
35 Correct 387 ms 47904 KB Output is correct
36 Correct 1241 ms 33312 KB Output is correct
37 Correct 1168 ms 33568 KB Output is correct
38 Correct 1137 ms 34772 KB Output is correct
39 Correct 1369 ms 97284 KB Output is correct
40 Correct 1507 ms 90308 KB Output is correct
41 Correct 1392 ms 92276 KB Output is correct
42 Correct 788 ms 22132 KB Output is correct
43 Correct 303 ms 26452 KB Output is correct
44 Correct 313 ms 26452 KB Output is correct
45 Correct 303 ms 26452 KB Output is correct
46 Correct 297 ms 26504 KB Output is correct
47 Correct 287 ms 26572 KB Output is correct
48 Correct 323 ms 26400 KB Output is correct
49 Correct 336 ms 26580 KB Output is correct
50 Correct 345 ms 26508 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 18 ms 14464 KB Output is correct
2 Correct 19 ms 14464 KB Output is correct
3 Correct 16 ms 14464 KB Output is correct
4 Correct 17 ms 14464 KB Output is correct
5 Correct 16 ms 14464 KB Output is correct
6 Correct 18 ms 14464 KB Output is correct
7 Correct 19 ms 14464 KB Output is correct
8 Correct 16 ms 14464 KB Output is correct
9 Correct 15 ms 14464 KB Output is correct
10 Correct 17 ms 14464 KB Output is correct
11 Correct 15 ms 14464 KB Output is correct
12 Correct 16 ms 14464 KB Output is correct
13 Correct 17 ms 14464 KB Output is correct
14 Correct 19 ms 14464 KB Output is correct
15 Correct 20 ms 14464 KB Output is correct
16 Correct 18 ms 14592 KB Output is correct
17 Correct 17 ms 14592 KB Output is correct
18 Correct 17 ms 14592 KB Output is correct
19 Correct 19 ms 14948 KB Output is correct
20 Correct 20 ms 14976 KB Output is correct
21 Correct 27 ms 15228 KB Output is correct
22 Correct 46 ms 15360 KB Output is correct
23 Correct 40 ms 15812 KB Output is correct
24 Correct 35 ms 15352 KB Output is correct
25 Correct 42 ms 15232 KB Output is correct
26 Correct 37 ms 15232 KB Output is correct
27 Correct 323 ms 43000 KB Output is correct
28 Correct 340 ms 42236 KB Output is correct
29 Correct 271 ms 38120 KB Output is correct
30 Correct 329 ms 42468 KB Output is correct
31 Correct 816 ms 70176 KB Output is correct
32 Correct 2087 ms 69208 KB Output is correct
33 Correct 1280 ms 83696 KB Output is correct
34 Correct 1364 ms 76724 KB Output is correct
35 Correct 15 ms 14464 KB Output is correct
36 Correct 1388 ms 38340 KB Output is correct
37 Execution timed out 3030 ms 81004 KB Time limit exceeded
38 Halted 0 ms 0 KB -