oj.uz

Submission #109078

# Submission time Handle Problem Language Result Execution time Memory
109078 2019-05-04T07:51:36 Z ryansee Circle selection (APIO18_circle_selection) C++14
72 / 100
 3000 ms 97196 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, 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;
ll magic = 5; // odd number pls
inline 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 <int, int> mp;
vector <int> d;
vector<pi>ind[MAXN],tmp[MAXN];
int latest = 0;
void 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 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 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:138:24: warning: unused variable 'y' [-Wunused-variable]
   ll x = mp[C[i].f/R], y = C[i].s;
                        ^

Subtask #1
7.0 / 7.0

# 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 16 ms 14396 KB Output is correct
4 Correct 18 ms 14464 KB Output is correct
5 Correct 17 ms 14464 KB Output is correct
6 Correct 16 ms 14464 KB Output is correct
7 Correct 17 ms 14584 KB Output is correct
8 Correct 17 ms 14464 KB Output is correct
9 Correct 15 ms 14464 KB Output is correct
10 Correct 16 ms 14464 KB Output is correct
11 Correct 14 ms 14464 KB Output is correct
12 Correct 17 ms 14464 KB Output is correct
13 Correct 16 ms 14692 KB Output is correct
14 Correct 15 ms 14464 KB Output is correct
15 Correct 16 ms 14464 KB Output is correct
16 Correct 16 ms 14592 KB Output is correct
17 Correct 19 ms 14524 KB Output is correct
18 Correct 15 ms 14592 KB Output is correct
19 Correct 22 ms 14976 KB Output is correct
20 Correct 25 ms 14976 KB Output is correct
21 Correct 26 ms 15104 KB Output is correct
22 Correct 37 ms 15480 KB Output is correct
23 Correct 39 ms 15736 KB Output is correct
24 Correct 36 ms 15356 KB Output is correct
25 Correct 33 ms 15360 KB Output is correct
26 Correct 33 ms 15232 KB Output is correct

Subtask #2
12.0 / 12.0

# Verdict Execution time Memory Grader output
1 Correct 278 ms 42908 KB Output is correct
2 Correct 294 ms 42012 KB Output is correct
3 Correct 296 ms 38216 KB Output is correct
4 Correct 301 ms 42488 KB Output is correct
5 Correct 666 ms 70096 KB Output is correct
6 Correct 1810 ms 69044 KB Output is correct
7 Correct 1172 ms 83768 KB Output is correct
8 Correct 1391 ms 76644 KB Output is correct

Subtask #3
0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 16 ms 14464 KB Output is correct
2 Correct 1274 ms 38276 KB Output is correct
3 Execution timed out 3032 ms 83964 KB Time limit exceeded
4 Halted 0 ms 0 KB -

Subtask #4
23.0 / 23.0

# Verdict Execution time Memory Grader output
1 Correct 734 ms 40932 KB Output is correct
2 Correct 853 ms 51892 KB Output is correct
3 Correct 968 ms 36072 KB Output is correct
4 Correct 904 ms 50920 KB Output is correct
5 Correct 868 ms 51396 KB Output is correct
6 Correct 1038 ms 34668 KB Output is correct

Subtask #5
30.0 / 30.0

# 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 16 ms 14396 KB Output is correct
4 Correct 18 ms 14464 KB Output is correct
5 Correct 17 ms 14464 KB Output is correct
6 Correct 16 ms 14464 KB Output is correct
7 Correct 17 ms 14584 KB Output is correct
8 Correct 17 ms 14464 KB Output is correct
9 Correct 15 ms 14464 KB Output is correct
10 Correct 16 ms 14464 KB Output is correct
11 Correct 14 ms 14464 KB Output is correct
12 Correct 17 ms 14464 KB Output is correct
13 Correct 16 ms 14692 KB Output is correct
14 Correct 15 ms 14464 KB Output is correct
15 Correct 16 ms 14464 KB Output is correct
16 Correct 16 ms 14592 KB Output is correct
17 Correct 19 ms 14524 KB Output is correct
18 Correct 15 ms 14592 KB Output is correct
19 Correct 22 ms 14976 KB Output is correct
20 Correct 25 ms 14976 KB Output is correct
21 Correct 26 ms 15104 KB Output is correct
22 Correct 37 ms 15480 KB Output is correct
23 Correct 39 ms 15736 KB Output is correct
24 Correct 36 ms 15356 KB Output is correct
25 Correct 33 ms 15360 KB Output is correct
26 Correct 33 ms 15232 KB Output is correct
27 Correct 32 ms 16252 KB Output is correct
28 Correct 22 ms 15332 KB Output is correct
29 Correct 23 ms 15360 KB Output is correct
30 Correct 63 ms 16632 KB Output is correct
31 Correct 59 ms 16296 KB Output is correct
32 Correct 63 ms 16220 KB Output is correct
33 Correct 90 ms 21228 KB Output is correct
34 Correct 88 ms 21436 KB Output is correct
35 Correct 332 ms 47936 KB Output is correct
36 Correct 1209 ms 33336 KB Output is correct
37 Correct 1163 ms 33568 KB Output is correct
38 Correct 1241 ms 34712 KB Output is correct
39 Correct 1008 ms 97196 KB Output is correct
40 Correct 1000 ms 90064 KB Output is correct
41 Correct 925 ms 92232 KB Output is correct
42 Correct 451 ms 22128 KB Output is correct
43 Correct 270 ms 26456 KB Output is correct
44 Correct 309 ms 26448 KB Output is correct
45 Correct 245 ms 26452 KB Output is correct
46 Correct 239 ms 26420 KB Output is correct
47 Correct 275 ms 26452 KB Output is correct
48 Correct 276 ms 26612 KB Output is correct
49 Correct 265 ms 26424 KB Output is correct
50 Correct 273 ms 26448 KB Output is correct

Subtask #6
0 / 13.0

# 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 16 ms 14396 KB Output is correct
4 Correct 18 ms 14464 KB Output is correct
5 Correct 17 ms 14464 KB Output is correct
6 Correct 16 ms 14464 KB Output is correct
7 Correct 17 ms 14584 KB Output is correct
8 Correct 17 ms 14464 KB Output is correct
9 Correct 15 ms 14464 KB Output is correct
10 Correct 16 ms 14464 KB Output is correct
11 Correct 14 ms 14464 KB Output is correct
12 Correct 17 ms 14464 KB Output is correct
13 Correct 16 ms 14692 KB Output is correct
14 Correct 15 ms 14464 KB Output is correct
15 Correct 16 ms 14464 KB Output is correct
16 Correct 16 ms 14592 KB Output is correct
17 Correct 19 ms 14524 KB Output is correct
18 Correct 15 ms 14592 KB Output is correct
19 Correct 22 ms 14976 KB Output is correct
20 Correct 25 ms 14976 KB Output is correct
21 Correct 26 ms 15104 KB Output is correct
22 Correct 37 ms 15480 KB Output is correct
23 Correct 39 ms 15736 KB Output is correct
24 Correct 36 ms 15356 KB Output is correct
25 Correct 33 ms 15360 KB Output is correct
26 Correct 33 ms 15232 KB Output is correct
27 Correct 278 ms 42908 KB Output is correct
28 Correct 294 ms 42012 KB Output is correct
29 Correct 296 ms 38216 KB Output is correct
30 Correct 301 ms 42488 KB Output is correct
31 Correct 666 ms 70096 KB Output is correct
32 Correct 1810 ms 69044 KB Output is correct
33 Correct 1172 ms 83768 KB Output is correct
34 Correct 1391 ms 76644 KB Output is correct
35 Correct 16 ms 14464 KB Output is correct
36 Correct 1274 ms 38276 KB Output is correct
37 Execution timed out 3032 ms 83964 KB Time limit exceeded
38 Halted 0 ms 0 KB -