// And you curse yourself for things you never done
#include<bits/stdc++.h>
#define F first
#define S second
#define PB push_back
#define sz(s) int((s).size())
#define bit(n,k) (((n)>>(k))&1)
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
typedef pair<ll, ll> pll;
const int maxn = 3e5 + 10, mod = 1e9 + 7, inf = 1e9 + 10;
pair<pll, ll> p[maxn];
int arr[maxn], ans[maxn];
pll operator - (pll a, pll b){
return {a.F - b.F, a.S - b.S};
}
pll operator + (pll a, pll b){
return {a.F + b.F, a.S + b.S};
}
ll operator ^ (pll a, pll b){
return 1ll * a.F * b.F + 1ll * a.S * b.S;
}
bool ok(int i, int j){
return 1ll * (p[j].S + p[i].S) * (p[j].S + p[i].S) >= ((p[i].F - p[j].F) ^ (p[i].F - p[j].F));
}
map<pii, vector<int> > mp, mp2;
int Log = 32;
int n;
int cnt[maxn];
vector<int> vvv[2][2];
void reSize(){
--Log;
mp2.clear();
for(auto &it : mp){
vvv[0][0].clear();
vvv[0][1].clear();
vvv[1][0].clear();
vvv[1][1].clear();
for(auto x : it.S){
if(ans[x] == -1)
vvv[bit(p[x].F.F, Log)][bit(p[x].F.S, Log)].PB(x);
}
if(!vvv[0][0].empty())
mp2[{2 * it.F.F, 2 * it.F.S}] = vvv[0][0];
if(!vvv[0][1].empty())
mp2[{2 * it.F.F, 2 * it.F.S + 1}] = vvv[0][1];
if(!vvv[1][0].empty())
mp2[{2 * it.F.F + 1, 2 * it.F.S}] = vvv[1][0];
if(!vvv[1][1].empty())
mp2[{2 * it.F.F + 1, 2 * it.F.S + 1}] = vvv[1][1];
}
mp = mp2;
}
int main(){
ios_base::sync_with_stdio(false); cin.tie(0); cout.tie();
memset(ans, -1, sizeof ans);
cin >> n;
for(int i = 0; i < n; i++){
cin >> p[i].F.F >> p[i].F.S >> p[i].S;
p[i].F.F+= inf, p[i].F.S+= inf;
}
auto cmp = [](int i, int j){
return (p[i].S == p[j].S ? i > j : p[i].S < p[j].S);
};
iota(arr, arr + n, 0);
sort(arr, arr + n, cmp);
vector<int> chert;
for(int i = 0; i < n; i++)
chert.PB(i);
mp[{0,0}] = chert;
for(int I = n-1; I >= 0; I--){
int i = arr[I];
while(Log > 0 && (1<<(Log-1)) >= p[i].S)
reSize();
if(ans[i] == -1){
ans[i] = i;
int x = p[i].F.F >> Log, y = p[i].F.S >> Log;
for(int X = x-2; X <= x+2; X++){
for(int Y = y-2; Y <= y+2; Y++){
if(mp.count({X, Y})){
vector<int> vec;
for(int id : mp[{X, Y}]){
if(ans[id] == -1 && ok(id, i))
ans[id] = i;
if(ans[id] == -1)
vec.PB(id);//, assert(++cnt[id] <= 10);
}
if(!vec.empty())
mp[{X, Y}] = vec;
}
}
}
}
}
for(int i = 0; i < n; i++){
cout << ans[i] + 1 << " ";
}
return cout << endl, 0;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
1536 KB |
Output is correct |
2 |
Correct |
1 ms |
1536 KB |
Output is correct |
3 |
Correct |
1 ms |
1536 KB |
Output is correct |
4 |
Correct |
1 ms |
1536 KB |
Output is correct |
5 |
Correct |
1 ms |
1536 KB |
Output is correct |
6 |
Correct |
1 ms |
1536 KB |
Output is correct |
7 |
Correct |
1 ms |
1536 KB |
Output is correct |
8 |
Correct |
1 ms |
1536 KB |
Output is correct |
9 |
Correct |
1 ms |
1536 KB |
Output is correct |
10 |
Correct |
1 ms |
1536 KB |
Output is correct |
11 |
Correct |
1 ms |
1536 KB |
Output is correct |
12 |
Correct |
1 ms |
1536 KB |
Output is correct |
13 |
Correct |
1 ms |
1536 KB |
Output is correct |
14 |
Correct |
1 ms |
1536 KB |
Output is correct |
15 |
Correct |
1 ms |
1536 KB |
Output is correct |
16 |
Correct |
2 ms |
1664 KB |
Output is correct |
17 |
Correct |
2 ms |
1536 KB |
Output is correct |
18 |
Correct |
2 ms |
1536 KB |
Output is correct |
19 |
Correct |
4 ms |
1792 KB |
Output is correct |
20 |
Correct |
6 ms |
1792 KB |
Output is correct |
21 |
Correct |
4 ms |
1792 KB |
Output is correct |
22 |
Correct |
55 ms |
1792 KB |
Output is correct |
23 |
Correct |
56 ms |
1792 KB |
Output is correct |
24 |
Correct |
57 ms |
1792 KB |
Output is correct |
25 |
Correct |
55 ms |
1792 KB |
Output is correct |
26 |
Correct |
54 ms |
1792 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
187 ms |
14376 KB |
Output is correct |
2 |
Correct |
189 ms |
14356 KB |
Output is correct |
3 |
Correct |
191 ms |
14440 KB |
Output is correct |
4 |
Correct |
197 ms |
14440 KB |
Output is correct |
5 |
Correct |
2774 ms |
16188 KB |
Output is correct |
6 |
Execution timed out |
3090 ms |
15760 KB |
Time limit exceeded |
7 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
1536 KB |
Output is correct |
2 |
Execution timed out |
3091 ms |
6460 KB |
Time limit exceeded |
3 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
3093 ms |
15608 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
1536 KB |
Output is correct |
2 |
Correct |
1 ms |
1536 KB |
Output is correct |
3 |
Correct |
1 ms |
1536 KB |
Output is correct |
4 |
Correct |
1 ms |
1536 KB |
Output is correct |
5 |
Correct |
1 ms |
1536 KB |
Output is correct |
6 |
Correct |
1 ms |
1536 KB |
Output is correct |
7 |
Correct |
1 ms |
1536 KB |
Output is correct |
8 |
Correct |
1 ms |
1536 KB |
Output is correct |
9 |
Correct |
1 ms |
1536 KB |
Output is correct |
10 |
Correct |
1 ms |
1536 KB |
Output is correct |
11 |
Correct |
1 ms |
1536 KB |
Output is correct |
12 |
Correct |
1 ms |
1536 KB |
Output is correct |
13 |
Correct |
1 ms |
1536 KB |
Output is correct |
14 |
Correct |
1 ms |
1536 KB |
Output is correct |
15 |
Correct |
1 ms |
1536 KB |
Output is correct |
16 |
Correct |
2 ms |
1664 KB |
Output is correct |
17 |
Correct |
2 ms |
1536 KB |
Output is correct |
18 |
Correct |
2 ms |
1536 KB |
Output is correct |
19 |
Correct |
4 ms |
1792 KB |
Output is correct |
20 |
Correct |
6 ms |
1792 KB |
Output is correct |
21 |
Correct |
4 ms |
1792 KB |
Output is correct |
22 |
Correct |
55 ms |
1792 KB |
Output is correct |
23 |
Correct |
56 ms |
1792 KB |
Output is correct |
24 |
Correct |
57 ms |
1792 KB |
Output is correct |
25 |
Correct |
55 ms |
1792 KB |
Output is correct |
26 |
Correct |
54 ms |
1792 KB |
Output is correct |
27 |
Correct |
12 ms |
2048 KB |
Output is correct |
28 |
Correct |
7 ms |
2048 KB |
Output is correct |
29 |
Correct |
7 ms |
2048 KB |
Output is correct |
30 |
Correct |
196 ms |
2168 KB |
Output is correct |
31 |
Correct |
204 ms |
2176 KB |
Output is correct |
32 |
Correct |
201 ms |
2176 KB |
Output is correct |
33 |
Correct |
67 ms |
5876 KB |
Output is correct |
34 |
Correct |
72 ms |
5876 KB |
Output is correct |
35 |
Correct |
65 ms |
6052 KB |
Output is correct |
36 |
Execution timed out |
3093 ms |
6564 KB |
Time limit exceeded |
37 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
1536 KB |
Output is correct |
2 |
Correct |
1 ms |
1536 KB |
Output is correct |
3 |
Correct |
1 ms |
1536 KB |
Output is correct |
4 |
Correct |
1 ms |
1536 KB |
Output is correct |
5 |
Correct |
1 ms |
1536 KB |
Output is correct |
6 |
Correct |
1 ms |
1536 KB |
Output is correct |
7 |
Correct |
1 ms |
1536 KB |
Output is correct |
8 |
Correct |
1 ms |
1536 KB |
Output is correct |
9 |
Correct |
1 ms |
1536 KB |
Output is correct |
10 |
Correct |
1 ms |
1536 KB |
Output is correct |
11 |
Correct |
1 ms |
1536 KB |
Output is correct |
12 |
Correct |
1 ms |
1536 KB |
Output is correct |
13 |
Correct |
1 ms |
1536 KB |
Output is correct |
14 |
Correct |
1 ms |
1536 KB |
Output is correct |
15 |
Correct |
1 ms |
1536 KB |
Output is correct |
16 |
Correct |
2 ms |
1664 KB |
Output is correct |
17 |
Correct |
2 ms |
1536 KB |
Output is correct |
18 |
Correct |
2 ms |
1536 KB |
Output is correct |
19 |
Correct |
4 ms |
1792 KB |
Output is correct |
20 |
Correct |
6 ms |
1792 KB |
Output is correct |
21 |
Correct |
4 ms |
1792 KB |
Output is correct |
22 |
Correct |
55 ms |
1792 KB |
Output is correct |
23 |
Correct |
56 ms |
1792 KB |
Output is correct |
24 |
Correct |
57 ms |
1792 KB |
Output is correct |
25 |
Correct |
55 ms |
1792 KB |
Output is correct |
26 |
Correct |
54 ms |
1792 KB |
Output is correct |
27 |
Correct |
187 ms |
14376 KB |
Output is correct |
28 |
Correct |
189 ms |
14356 KB |
Output is correct |
29 |
Correct |
191 ms |
14440 KB |
Output is correct |
30 |
Correct |
197 ms |
14440 KB |
Output is correct |
31 |
Correct |
2774 ms |
16188 KB |
Output is correct |
32 |
Execution timed out |
3090 ms |
15760 KB |
Time limit exceeded |
33 |
Halted |
0 ms |
0 KB |
- |