#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
struct point
{
ll x, y;
};
struct circle
{
point centre;
ll radius;
};
const int maxn = 5010;
int n;
circle c[maxn];
void input()
{
cin >> n;
for (int i = 1; i <= n; i ++)
cin >> c[i].centre.x >> c[i].centre.y >> c[i].radius;
}
int used[maxn], par[maxn];
int find_biggest()
{
int idx = -1;
for (int i = 1; i <= n; i ++)
{
if (used[i])
continue;
if (idx == -1 || c[i].radius > c[idx].radius)
idx = i;
}
return idx;
}
ll distance(point A, point B)
{
ll dx = (A.x - B.x);
ll dy = (A.y - B.y);
return dx * dx + dy * dy;
}
void simulate()
{
while(true)
{
int idx = find_biggest();
if (idx == -1)
break;
///cout << "current " << idx << endl;
for (int i = 1; i <= n; i ++)
{
if (used[i] == 1)
continue;
if (distance(c[idx].centre, c[i].centre) <=
(c[idx].radius + c[i].radius) * (c[idx].radius + c[i].radius))
used[i] = 1, par[i] = idx;
}
}
}
void print()
{
for (int i = 1; i <= n; i ++)
cout << par[i] << " ";
cout << endl;
}
void solve_subtask_2()
{
set < pair < int, int > > lf, rf;
set < pair < int, int > > act;
for (int i = 1; i <= n; i ++)
{
lf.insert(make_pair(c[i].centre.x - c[i].radius, i));
rf.insert(make_pair(c[i].centre.x + c[i].radius, i));
act.insert(make_pair(c[i].radius, i));
}
while(!act.empty())
{
int idx = act.begin() -> second;
int left = c[idx].centre.x - c[idx].radius;
int right = c[idx].centre.x + c[idx].radius;
while(true)
{
set < pair < int, int > > :: iterator it = lf.lower_bound({left, -1});
if (it -> first <= right)
{
int i = it -> second;
par[i] = idx;
lf.erase(make_pair(c[i].centre.x - c[i].radius, i));
rf.erase(make_pair(c[i].centre.x + c[i].radius, i));
act.erase(make_pair(c[i].radius, i));
}
else
{
break;
}
}
while(true)
{
set < pair < int, int > > :: iterator it = rf.lower_bound({left, -1});
if (it -> first <= right)
{
int i = it -> second;
par[i] = idx;
lf.erase(make_pair(c[i].centre.x - c[i].radius, i));
rf.erase(make_pair(c[i].centre.x + c[i].radius, i));
act.erase(make_pair(c[i].radius, i));
}
else
{
break;
}
}
}
}
void solve()
{
input();
bool all_on_x = true;
for (int i = 1; i <= n; i ++)
{
if (c[i].centre.y != 0)
all_on_x = false;
}
if (all_on_x)
solve_subtask_2();
else
simulate();
print();
}
int main()
{
solve();
return 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
*/
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
344 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 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 |
0 ms |
348 KB |
Output is correct |
12 |
Correct |
0 ms |
348 KB |
Output is correct |
13 |
Correct |
0 ms |
348 KB |
Output is correct |
14 |
Correct |
0 ms |
348 KB |
Output is correct |
15 |
Correct |
0 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 |
5 ms |
448 KB |
Output is correct |
20 |
Correct |
5 ms |
600 KB |
Output is correct |
21 |
Correct |
5 ms |
600 KB |
Output is correct |
22 |
Correct |
110 ms |
348 KB |
Output is correct |
23 |
Correct |
112 ms |
600 KB |
Output is correct |
24 |
Correct |
110 ms |
576 KB |
Output is correct |
25 |
Correct |
110 ms |
344 KB |
Output is correct |
26 |
Correct |
109 ms |
564 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
5 ms |
600 KB |
Execution killed with signal 11 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Runtime error |
5 ms |
780 KB |
Execution killed with signal 11 |
3 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
5 ms |
600 KB |
Execution killed with signal 11 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
344 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 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 |
0 ms |
348 KB |
Output is correct |
12 |
Correct |
0 ms |
348 KB |
Output is correct |
13 |
Correct |
0 ms |
348 KB |
Output is correct |
14 |
Correct |
0 ms |
348 KB |
Output is correct |
15 |
Correct |
0 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 |
5 ms |
448 KB |
Output is correct |
20 |
Correct |
5 ms |
600 KB |
Output is correct |
21 |
Correct |
5 ms |
600 KB |
Output is correct |
22 |
Correct |
110 ms |
348 KB |
Output is correct |
23 |
Correct |
112 ms |
600 KB |
Output is correct |
24 |
Correct |
110 ms |
576 KB |
Output is correct |
25 |
Correct |
110 ms |
344 KB |
Output is correct |
26 |
Correct |
109 ms |
564 KB |
Output is correct |
27 |
Runtime error |
5 ms |
600 KB |
Execution killed with signal 11 |
28 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
344 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 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 |
0 ms |
348 KB |
Output is correct |
12 |
Correct |
0 ms |
348 KB |
Output is correct |
13 |
Correct |
0 ms |
348 KB |
Output is correct |
14 |
Correct |
0 ms |
348 KB |
Output is correct |
15 |
Correct |
0 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 |
5 ms |
448 KB |
Output is correct |
20 |
Correct |
5 ms |
600 KB |
Output is correct |
21 |
Correct |
5 ms |
600 KB |
Output is correct |
22 |
Correct |
110 ms |
348 KB |
Output is correct |
23 |
Correct |
112 ms |
600 KB |
Output is correct |
24 |
Correct |
110 ms |
576 KB |
Output is correct |
25 |
Correct |
110 ms |
344 KB |
Output is correct |
26 |
Correct |
109 ms |
564 KB |
Output is correct |
27 |
Runtime error |
5 ms |
600 KB |
Execution killed with signal 11 |
28 |
Halted |
0 ms |
0 KB |
- |