Submission #109058

# Submission time Handle Problem Language Result Execution time Memory
109058 2019-05-04T07:17:45 Z ryansee Circle selection (APIO18_circle_selection) C++14
60 / 100
3000 ms 97984 KB
#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, ll> 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;

ll magic = 5; // odd number pls
bool 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 <ll, ll> mp;
vector <ll> d;
vector<pi>ind[MAXN],tmp[MAXN];
ll latest = 0;
void split() { 
	FOR(i,0,MAXN) 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) { mp[C[ind[i][0].s].f/new_R] = co; tmp[co++].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,MAXN) ind[i]=tmp[i];
  return;
	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 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)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 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: In function 'void st4()':
circle_selection.cpp:140: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 15 ms 14464 KB Output is correct
2 Correct 17 ms 14464 KB Output is correct
3 Correct 17 ms 14464 KB Output is correct
4 Correct 18 ms 14464 KB Output is correct
5 Correct 14 ms 14464 KB Output is correct
6 Correct 15 ms 14464 KB Output is correct
7 Correct 17 ms 14540 KB Output is correct
8 Correct 15 ms 14592 KB Output is correct
9 Correct 19 ms 14464 KB Output is correct
10 Correct 18 ms 14464 KB Output is correct
11 Correct 31 ms 14464 KB Output is correct
12 Correct 29 ms 14464 KB Output is correct
13 Correct 25 ms 14464 KB Output is correct
14 Correct 36 ms 14532 KB Output is correct
15 Correct 35 ms 14464 KB Output is correct
16 Correct 28 ms 14592 KB Output is correct
17 Correct 15 ms 14592 KB Output is correct
18 Correct 19 ms 14592 KB Output is correct
19 Correct 20 ms 14848 KB Output is correct
20 Correct 22 ms 14848 KB Output is correct
21 Correct 39 ms 15096 KB Output is correct
22 Correct 60 ms 15480 KB Output is correct
23 Correct 79 ms 15864 KB Output is correct
24 Correct 85 ms 15352 KB Output is correct
25 Correct 68 ms 15204 KB Output is correct
26 Correct 89 ms 15352 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 250 ms 44024 KB Output is correct
2 Correct 312 ms 44304 KB Output is correct
3 Correct 243 ms 39392 KB Output is correct
4 Correct 308 ms 43972 KB Output is correct
5 Correct 678 ms 71308 KB Output is correct
6 Incorrect 1666 ms 70056 KB Output isn't correct
7 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 25 ms 14464 KB Output is correct
2 Correct 1356 ms 38676 KB Output is correct
3 Execution timed out 3035 ms 85220 KB Time limit exceeded
4 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1014 ms 42192 KB Output is correct
2 Correct 1104 ms 53244 KB Output is correct
3 Correct 1673 ms 37092 KB Output is correct
4 Correct 895 ms 52192 KB Output is correct
5 Correct 927 ms 52704 KB Output is correct
6 Correct 1224 ms 35808 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 15 ms 14464 KB Output is correct
2 Correct 17 ms 14464 KB Output is correct
3 Correct 17 ms 14464 KB Output is correct
4 Correct 18 ms 14464 KB Output is correct
5 Correct 14 ms 14464 KB Output is correct
6 Correct 15 ms 14464 KB Output is correct
7 Correct 17 ms 14540 KB Output is correct
8 Correct 15 ms 14592 KB Output is correct
9 Correct 19 ms 14464 KB Output is correct
10 Correct 18 ms 14464 KB Output is correct
11 Correct 31 ms 14464 KB Output is correct
12 Correct 29 ms 14464 KB Output is correct
13 Correct 25 ms 14464 KB Output is correct
14 Correct 36 ms 14532 KB Output is correct
15 Correct 35 ms 14464 KB Output is correct
16 Correct 28 ms 14592 KB Output is correct
17 Correct 15 ms 14592 KB Output is correct
18 Correct 19 ms 14592 KB Output is correct
19 Correct 20 ms 14848 KB Output is correct
20 Correct 22 ms 14848 KB Output is correct
21 Correct 39 ms 15096 KB Output is correct
22 Correct 60 ms 15480 KB Output is correct
23 Correct 79 ms 15864 KB Output is correct
24 Correct 85 ms 15352 KB Output is correct
25 Correct 68 ms 15204 KB Output is correct
26 Correct 89 ms 15352 KB Output is correct
27 Correct 57 ms 16376 KB Output is correct
28 Correct 22 ms 15360 KB Output is correct
29 Correct 23 ms 15360 KB Output is correct
30 Correct 117 ms 16732 KB Output is correct
31 Correct 109 ms 16244 KB Output is correct
32 Correct 112 ms 16248 KB Output is correct
33 Correct 135 ms 21576 KB Output is correct
34 Correct 84 ms 22356 KB Output is correct
35 Correct 340 ms 48268 KB Output is correct
36 Correct 1078 ms 33604 KB Output is correct
37 Correct 1208 ms 33988 KB Output is correct
38 Correct 1276 ms 35616 KB Output is correct
39 Correct 1231 ms 97984 KB Output is correct
40 Correct 1156 ms 91448 KB Output is correct
41 Correct 1296 ms 93592 KB Output is correct
42 Correct 419 ms 24044 KB Output is correct
43 Correct 318 ms 29308 KB Output is correct
44 Correct 288 ms 29416 KB Output is correct
45 Correct 250 ms 29444 KB Output is correct
46 Correct 240 ms 29548 KB Output is correct
47 Correct 244 ms 29432 KB Output is correct
48 Correct 243 ms 29408 KB Output is correct
49 Correct 294 ms 29552 KB Output is correct
50 Correct 263 ms 29420 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 15 ms 14464 KB Output is correct
2 Correct 17 ms 14464 KB Output is correct
3 Correct 17 ms 14464 KB Output is correct
4 Correct 18 ms 14464 KB Output is correct
5 Correct 14 ms 14464 KB Output is correct
6 Correct 15 ms 14464 KB Output is correct
7 Correct 17 ms 14540 KB Output is correct
8 Correct 15 ms 14592 KB Output is correct
9 Correct 19 ms 14464 KB Output is correct
10 Correct 18 ms 14464 KB Output is correct
11 Correct 31 ms 14464 KB Output is correct
12 Correct 29 ms 14464 KB Output is correct
13 Correct 25 ms 14464 KB Output is correct
14 Correct 36 ms 14532 KB Output is correct
15 Correct 35 ms 14464 KB Output is correct
16 Correct 28 ms 14592 KB Output is correct
17 Correct 15 ms 14592 KB Output is correct
18 Correct 19 ms 14592 KB Output is correct
19 Correct 20 ms 14848 KB Output is correct
20 Correct 22 ms 14848 KB Output is correct
21 Correct 39 ms 15096 KB Output is correct
22 Correct 60 ms 15480 KB Output is correct
23 Correct 79 ms 15864 KB Output is correct
24 Correct 85 ms 15352 KB Output is correct
25 Correct 68 ms 15204 KB Output is correct
26 Correct 89 ms 15352 KB Output is correct
27 Correct 250 ms 44024 KB Output is correct
28 Correct 312 ms 44304 KB Output is correct
29 Correct 243 ms 39392 KB Output is correct
30 Correct 308 ms 43972 KB Output is correct
31 Correct 678 ms 71308 KB Output is correct
32 Incorrect 1666 ms 70056 KB Output isn't correct
33 Halted 0 ms 0 KB -