oj.uz

Submission #109062

# Submission time Handle Problem Language Result Execution time Memory
109062 2019-05-04T07:21:26 Z ryansee Circle selection (APIO18_circle_selection) C++14
60 / 100
 3000 ms 97696 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];
	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:139: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 16 ms 14464 KB Output is correct
2 Correct 17 ms 14464 KB Output is correct
3 Correct 20 ms 14436 KB Output is correct
4 Correct 18 ms 14464 KB Output is correct
5 Correct 17 ms 14464 KB Output is correct
6 Correct 21 ms 14456 KB Output is correct
7 Correct 19 ms 14464 KB Output is correct
8 Correct 16 ms 14592 KB Output is correct
9 Correct 19 ms 14464 KB Output is correct
10 Correct 15 ms 14464 KB Output is correct
11 Correct 30 ms 14464 KB Output is correct
12 Correct 43 ms 14464 KB Output is correct
13 Correct 40 ms 14464 KB Output is correct
14 Correct 38 ms 14436 KB Output is correct
15 Correct 35 ms 14464 KB Output is correct
16 Correct 28 ms 14592 KB Output is correct
17 Correct 16 ms 14592 KB Output is correct
18 Correct 18 ms 14720 KB Output is correct
19 Correct 20 ms 14848 KB Output is correct
20 Correct 20 ms 14848 KB Output is correct
21 Correct 43 ms 15096 KB Output is correct
22 Correct 96 ms 15352 KB Output is correct
23 Correct 72 ms 15864 KB Output is correct
24 Correct 77 ms 15480 KB Output is correct
25 Correct 61 ms 15360 KB Output is correct
26 Correct 64 ms 15480 KB Output is correct

Subtask #2
0 / 12.0

# Verdict Execution time Memory Grader output
1 Correct 275 ms 44084 KB Output is correct
2 Correct 298 ms 44128 KB Output is correct
3 Correct 257 ms 39392 KB Output is correct
4 Correct 258 ms 43872 KB Output is correct
5 Correct 681 ms 71264 KB Output is correct
6 Incorrect 1970 ms 70112 KB Output isn't correct
7 Halted 0 ms 0 KB -

Subtask #3
0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 31 ms 14464 KB Output is correct
2 Correct 1526 ms 38676 KB Output is correct
3 Execution timed out 3012 ms 84564 KB Time limit exceeded
4 Halted 0 ms 0 KB -

Subtask #4
23.0 / 23.0

# Verdict Execution time Memory Grader output
1 Correct 843 ms 41936 KB Output is correct
2 Correct 1010 ms 53112 KB Output is correct
3 Correct 1576 ms 37060 KB Output is correct
4 Correct 921 ms 52192 KB Output is correct
5 Correct 969 ms 52676 KB Output is correct
6 Correct 1240 ms 35684 KB Output is correct

Subtask #5
30.0 / 30.0

# Verdict Execution time Memory Grader output
1 Correct 16 ms 14464 KB Output is correct
2 Correct 17 ms 14464 KB Output is correct
3 Correct 20 ms 14436 KB Output is correct
4 Correct 18 ms 14464 KB Output is correct
5 Correct 17 ms 14464 KB Output is correct
6 Correct 21 ms 14456 KB Output is correct
7 Correct 19 ms 14464 KB Output is correct
8 Correct 16 ms 14592 KB Output is correct
9 Correct 19 ms 14464 KB Output is correct
10 Correct 15 ms 14464 KB Output is correct
11 Correct 30 ms 14464 KB Output is correct
12 Correct 43 ms 14464 KB Output is correct
13 Correct 40 ms 14464 KB Output is correct
14 Correct 38 ms 14436 KB Output is correct
15 Correct 35 ms 14464 KB Output is correct
16 Correct 28 ms 14592 KB Output is correct
17 Correct 16 ms 14592 KB Output is correct
18 Correct 18 ms 14720 KB Output is correct
19 Correct 20 ms 14848 KB Output is correct
20 Correct 20 ms 14848 KB Output is correct
21 Correct 43 ms 15096 KB Output is correct
22 Correct 96 ms 15352 KB Output is correct
23 Correct 72 ms 15864 KB Output is correct
24 Correct 77 ms 15480 KB Output is correct
25 Correct 61 ms 15360 KB Output is correct
26 Correct 64 ms 15480 KB Output is correct
27 Correct 60 ms 16504 KB Output is correct
28 Correct 22 ms 15488 KB Output is correct
29 Correct 27 ms 15608 KB Output is correct
30 Correct 101 ms 16888 KB Output is correct
31 Correct 117 ms 16432 KB Output is correct
32 Correct 148 ms 16376 KB Output is correct
33 Correct 104 ms 22252 KB Output is correct
34 Correct 100 ms 23024 KB Output is correct
35 Correct 366 ms 48836 KB Output is correct
36 Correct 1309 ms 34416 KB Output is correct
37 Correct 1350 ms 33904 KB Output is correct
38 Correct 1186 ms 34928 KB Output is correct
39 Correct 1189 ms 97696 KB Output is correct
40 Correct 1137 ms 90376 KB Output is correct
41 Correct 1144 ms 92604 KB Output is correct
42 Correct 430 ms 22636 KB Output is correct
43 Correct 249 ms 26860 KB Output is correct
44 Correct 281 ms 26864 KB Output is correct
45 Correct 272 ms 26788 KB Output is correct
46 Correct 274 ms 26932 KB Output is correct
47 Correct 256 ms 26864 KB Output is correct
48 Correct 275 ms 26832 KB Output is correct
49 Correct 268 ms 26860 KB Output is correct
50 Correct 279 ms 26916 KB Output is correct

Subtask #6
0 / 13.0

# Verdict Execution time Memory Grader output
1 Correct 16 ms 14464 KB Output is correct
2 Correct 17 ms 14464 KB Output is correct
3 Correct 20 ms 14436 KB Output is correct
4 Correct 18 ms 14464 KB Output is correct
5 Correct 17 ms 14464 KB Output is correct
6 Correct 21 ms 14456 KB Output is correct
7 Correct 19 ms 14464 KB Output is correct
8 Correct 16 ms 14592 KB Output is correct
9 Correct 19 ms 14464 KB Output is correct
10 Correct 15 ms 14464 KB Output is correct
11 Correct 30 ms 14464 KB Output is correct
12 Correct 43 ms 14464 KB Output is correct
13 Correct 40 ms 14464 KB Output is correct
14 Correct 38 ms 14436 KB Output is correct
15 Correct 35 ms 14464 KB Output is correct
16 Correct 28 ms 14592 KB Output is correct
17 Correct 16 ms 14592 KB Output is correct
18 Correct 18 ms 14720 KB Output is correct
19 Correct 20 ms 14848 KB Output is correct
20 Correct 20 ms 14848 KB Output is correct
21 Correct 43 ms 15096 KB Output is correct
22 Correct 96 ms 15352 KB Output is correct
23 Correct 72 ms 15864 KB Output is correct
24 Correct 77 ms 15480 KB Output is correct
25 Correct 61 ms 15360 KB Output is correct
26 Correct 64 ms 15480 KB Output is correct
27 Correct 275 ms 44084 KB Output is correct
28 Correct 298 ms 44128 KB Output is correct
29 Correct 257 ms 39392 KB Output is correct
30 Correct 258 ms 43872 KB Output is correct
31 Correct 681 ms 71264 KB Output is correct
32 Incorrect 1970 ms 70112 KB Output isn't correct
33 Halted 0 ms 0 KB -