#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
#define ALL(v) (v).begin(), (v).end()
#define MASK(i) (1LL << (i))
#define GETBIT(mask, i) (((mask) >> (i)) & 1)
// mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
mt19937_64 rng(1);
ll rngesus(ll l, ll r){return ((ull) rng()) % (r - l + 1) + l;}
ll max(ll a, ll b){return (a > b) ? a : b;}
ll min(ll a, ll b){return (a < b) ? a : b;}
ll LASTBIT(ll mask){return mask & (-mask);}
ll pop_cnt(ll mask){return __builtin_popcountll(mask);}
ll ctz(ll mask){return __builtin_ctzll(mask);}
ll clz(ll mask){return __builtin_clzll(mask);}
ll logOf(ll mask){return 63 - clz(mask);}
template <class T1, class T2>
bool minimize(T1 &a, T2 b){
if (a > b){a = b; return true;}
return false;
}
template <class T1, class T2>
bool maximize(T1 &a, T2 b){
if (a < b){a = b; return true;}
return false;
}
template <class T>
void printArr(T& a, string separator = " ", string finish = "\n", ostream& out = cout){
for(auto i: a) out << i << separator;
out << finish;
}
template <class T>
void remove_dup(vector<T> &a){
sort(ALL(a));
a.resize(unique(ALL(a)) - a.begin());
}
const int N = 3e5 + 69;
const ll INF = 1e9;
int n;
int ans[N];
vector<array<ll, 4>> a;
ll sqr(ll x){return x * x;}
ll block_sz;
map<pair<ll, ll>, vector<int>> grid;
pair<ll, ll> dumb_down(pair<ll, ll> a){
return {a.first / block_sz, a.second / block_sz};
}
void set_up(ll arg){
block_sz = arg;
grid.clear();
for(int i= 1; i<=n; ++i) if (ans[i] == 0){
grid[dumb_down(make_pair(a[i][0], a[i][1]))].push_back(i);
}
}
ll dis(pair<ll, ll> a, pair<ll, ll> b){return sqr(a.first - b.first) + sqr(a.second - b.second);}
bool check(pair<ll, ll> point, ll r, pair<ll, ll> des){
if (des.first < 0 || des.second < 0) return false;
for(int x = 0; x <= 1; ++x) for(int y = 0; y<=1; ++y){
pair<ll, ll> ligma = {des.first * block_sz, des.second * block_sz};
if (x) ligma.first += block_sz - 1;
if (y) ligma.second += block_sz - 1;
if (dis(point, ligma) <= sqr(r)) return true;
}
return false;
}
bool query(int i, int j){
return sqr(a[i][0] - a[j][0]) + sqr(a[i][1] - a[j][1]) <= sqr(a[i][2] + a[j][2]);
}
int main(void){
ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
cin >> n;
a.resize(n+1);
for(int i = 1; i<=n; ++i) {
for(int j = 0; j < 3; ++j) cin >> a[i][j];
a[i][0] += INF; a[i][1] += INF;
a[i][3] = i;
}
vector<array<ll, 4>> b = a;
sort(1 + ALL(b), [](array<ll, 4> x, array<ll, 4> y){
if (x[2] != y[2]) return x[2] > y[2];
return x[3] < y[3];
});
set_up((b[1][2] + 1) / 2);
for(int i = 1; i<=n; ++i) if (ans[b[i][3]] == 0){
int _i = b[i][3];
ans[_i] = _i;
if (b[i][2] < block_sz)
set_up((b[i][2] + 1) / 2);
pair<ll, ll> cur = make_pair(b[i][0], b[i][1]);
pair<ll, ll> magnifico = dumb_down(cur);
const int EXTEND = 10;
for(int x = -EXTEND; x <= EXTEND; ++x) for(int y = -EXTEND; y<=EXTEND; ++y) {
pair<ll, ll> cu = make_pair(magnifico.first + x, magnifico.second + y);
if (max(abs(x), abs(y)) == EXTEND + 1 && check(cur, b[i][2] * 2, cu)) exit(1);
if (!check(cur, b[i][2] * 2, cu) || !grid.count(cu)) continue;
for(int j: grid[cu]){
if (ans[j] == 0 && query(_i, j)) ans[j] = _i;
}
}
}
for(int i = 1; i<=n; ++i)
cout << ans[i] << " ";
cout << "\n";
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
1 ms |
348 KB |
Output is correct |
4 |
Correct |
1 ms |
348 KB |
Output is correct |
5 |
Correct |
1 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
1 ms |
348 KB |
Output is correct |
9 |
Correct |
0 ms |
348 KB |
Output is correct |
10 |
Correct |
0 ms |
348 KB |
Output is correct |
11 |
Correct |
1 ms |
344 KB |
Output is correct |
12 |
Correct |
1 ms |
348 KB |
Output is correct |
13 |
Correct |
1 ms |
348 KB |
Output is correct |
14 |
Correct |
1 ms |
600 KB |
Output is correct |
15 |
Correct |
2 ms |
348 KB |
Output is correct |
16 |
Correct |
1 ms |
348 KB |
Output is correct |
17 |
Correct |
1 ms |
344 KB |
Output is correct |
18 |
Correct |
1 ms |
348 KB |
Output is correct |
19 |
Correct |
3 ms |
788 KB |
Output is correct |
20 |
Correct |
3 ms |
1048 KB |
Output is correct |
21 |
Correct |
2 ms |
860 KB |
Output is correct |
22 |
Correct |
27 ms |
1500 KB |
Output is correct |
23 |
Correct |
37 ms |
1624 KB |
Output is correct |
24 |
Correct |
30 ms |
1748 KB |
Output is correct |
25 |
Correct |
28 ms |
1372 KB |
Output is correct |
26 |
Correct |
31 ms |
1372 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
113 ms |
23032 KB |
Output is correct |
2 |
Correct |
128 ms |
22708 KB |
Output is correct |
3 |
Correct |
121 ms |
23468 KB |
Output is correct |
4 |
Correct |
127 ms |
22524 KB |
Output is correct |
5 |
Incorrect |
195 ms |
30340 KB |
Output isn't correct |
6 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
344 KB |
Output is correct |
2 |
Correct |
776 ms |
18512 KB |
Output is correct |
3 |
Correct |
2953 ms |
55076 KB |
Output is correct |
4 |
Correct |
2728 ms |
55076 KB |
Output is correct |
5 |
Correct |
2433 ms |
51700 KB |
Output is correct |
6 |
Correct |
796 ms |
26648 KB |
Output is correct |
7 |
Correct |
344 ms |
14392 KB |
Output is correct |
8 |
Correct |
54 ms |
3332 KB |
Output is correct |
9 |
Correct |
2917 ms |
54784 KB |
Output is correct |
10 |
Correct |
2062 ms |
51364 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
2842 ms |
55116 KB |
Output is correct |
2 |
Correct |
1670 ms |
55128 KB |
Output is correct |
3 |
Incorrect |
583 ms |
44152 KB |
Output isn't correct |
4 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
1 ms |
348 KB |
Output is correct |
4 |
Correct |
1 ms |
348 KB |
Output is correct |
5 |
Correct |
1 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
1 ms |
348 KB |
Output is correct |
9 |
Correct |
0 ms |
348 KB |
Output is correct |
10 |
Correct |
0 ms |
348 KB |
Output is correct |
11 |
Correct |
1 ms |
344 KB |
Output is correct |
12 |
Correct |
1 ms |
348 KB |
Output is correct |
13 |
Correct |
1 ms |
348 KB |
Output is correct |
14 |
Correct |
1 ms |
600 KB |
Output is correct |
15 |
Correct |
2 ms |
348 KB |
Output is correct |
16 |
Correct |
1 ms |
348 KB |
Output is correct |
17 |
Correct |
1 ms |
344 KB |
Output is correct |
18 |
Correct |
1 ms |
348 KB |
Output is correct |
19 |
Correct |
3 ms |
788 KB |
Output is correct |
20 |
Correct |
3 ms |
1048 KB |
Output is correct |
21 |
Correct |
2 ms |
860 KB |
Output is correct |
22 |
Correct |
27 ms |
1500 KB |
Output is correct |
23 |
Correct |
37 ms |
1624 KB |
Output is correct |
24 |
Correct |
30 ms |
1748 KB |
Output is correct |
25 |
Correct |
28 ms |
1372 KB |
Output is correct |
26 |
Correct |
31 ms |
1372 KB |
Output is correct |
27 |
Correct |
5 ms |
1360 KB |
Output is correct |
28 |
Correct |
5 ms |
1368 KB |
Output is correct |
29 |
Correct |
6 ms |
1368 KB |
Output is correct |
30 |
Correct |
60 ms |
2484 KB |
Output is correct |
31 |
Correct |
76 ms |
2392 KB |
Output is correct |
32 |
Correct |
55 ms |
2440 KB |
Output is correct |
33 |
Correct |
49 ms |
11708 KB |
Output is correct |
34 |
Correct |
50 ms |
11472 KB |
Output is correct |
35 |
Correct |
56 ms |
10716 KB |
Output is correct |
36 |
Correct |
644 ms |
21116 KB |
Output is correct |
37 |
Correct |
759 ms |
21216 KB |
Output is correct |
38 |
Correct |
720 ms |
21260 KB |
Output is correct |
39 |
Correct |
305 ms |
18128 KB |
Output is correct |
40 |
Correct |
328 ms |
17920 KB |
Output is correct |
41 |
Correct |
308 ms |
17972 KB |
Output is correct |
42 |
Correct |
309 ms |
20076 KB |
Output is correct |
43 |
Correct |
626 ms |
21076 KB |
Output is correct |
44 |
Correct |
632 ms |
21084 KB |
Output is correct |
45 |
Correct |
635 ms |
21060 KB |
Output is correct |
46 |
Correct |
640 ms |
21076 KB |
Output is correct |
47 |
Correct |
620 ms |
21072 KB |
Output is correct |
48 |
Correct |
677 ms |
21088 KB |
Output is correct |
49 |
Correct |
637 ms |
21072 KB |
Output is correct |
50 |
Correct |
608 ms |
21076 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
1 ms |
348 KB |
Output is correct |
4 |
Correct |
1 ms |
348 KB |
Output is correct |
5 |
Correct |
1 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
1 ms |
348 KB |
Output is correct |
9 |
Correct |
0 ms |
348 KB |
Output is correct |
10 |
Correct |
0 ms |
348 KB |
Output is correct |
11 |
Correct |
1 ms |
344 KB |
Output is correct |
12 |
Correct |
1 ms |
348 KB |
Output is correct |
13 |
Correct |
1 ms |
348 KB |
Output is correct |
14 |
Correct |
1 ms |
600 KB |
Output is correct |
15 |
Correct |
2 ms |
348 KB |
Output is correct |
16 |
Correct |
1 ms |
348 KB |
Output is correct |
17 |
Correct |
1 ms |
344 KB |
Output is correct |
18 |
Correct |
1 ms |
348 KB |
Output is correct |
19 |
Correct |
3 ms |
788 KB |
Output is correct |
20 |
Correct |
3 ms |
1048 KB |
Output is correct |
21 |
Correct |
2 ms |
860 KB |
Output is correct |
22 |
Correct |
27 ms |
1500 KB |
Output is correct |
23 |
Correct |
37 ms |
1624 KB |
Output is correct |
24 |
Correct |
30 ms |
1748 KB |
Output is correct |
25 |
Correct |
28 ms |
1372 KB |
Output is correct |
26 |
Correct |
31 ms |
1372 KB |
Output is correct |
27 |
Correct |
113 ms |
23032 KB |
Output is correct |
28 |
Correct |
128 ms |
22708 KB |
Output is correct |
29 |
Correct |
121 ms |
23468 KB |
Output is correct |
30 |
Correct |
127 ms |
22524 KB |
Output is correct |
31 |
Incorrect |
195 ms |
30340 KB |
Output isn't correct |
32 |
Halted |
0 ms |
0 KB |
- |