// M
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 300005;
int n, X[N], Y[N], R[N], M[N];
vector < int > UX;
vector < pair < int , int > > V[N * 2];
set < pair < int , int > > ST[N * 2];
inline int GetId(int a)
{
return (int)(lower_bound(UX.begin(), UX.end(), a) - UX.begin());
}
inline bool Check(int i, int j)
{
return ((ll)(X[i] - X[j]) * (ll)(X[i] - X[j]) + (ll)(Y[i] - Y[j]) * (ll)(Y[i] - Y[j]) <= (ll)(R[i] + R[j]) * (ll)(R[i] + R[j]));
}
inline void Get(int l, int r, int le, int ri, int i)
{
l += (int)UX.size();
r += (int)UX.size();
for (; l < r; l >>= 1, r >>= 1)
{
if (l & 1)
{
auto it = ST[l].lower_bound({le, -1});
while (it != ST[l].end() && it->first < ri)
{
if (M[it->second])
it = ST[l].erase(it);
else if (Check(i, it->second))
M[it->second] = i, it = ST[l].erase(it);
else
it ++;
}
l ++;
}
if (r & 1)
{
r --;
auto it = ST[r].lower_bound({le, -1});
while (it != ST[r].end() && it->first < ri)
{
if (M[it->second])
it = ST[r].erase(it);
else if (Check(i, it->second))
M[it->second] = i, it = ST[r].erase(it);
else
it ++;
}
}
}
}
inline void Do(int i)
{
M[i] = i;
int lx = GetId(X[i] - R[i] * 2);
int rx = GetId(X[i] + R[i] * 2 + 1);
int ly = Y[i] - R[i] * 2;
int ry = Y[i] + R[i] * 2 + 1;
Get(lx, rx, ly, ry, i);
}
int main()
{
scanf("%d", &n);
for (int i = 1; i <= n; i ++)
scanf("%d%d%d", &X[i], &Y[i], &R[i]), UX.push_back(X[i]);
sort(UX.begin(), UX.end());
UX.resize(unique(UX.begin(), UX.end()) - UX.begin());
for (int i = 1; i <= n; i ++)
{
int x = GetId(X[i]);
V[x + (int)UX.size()].push_back({Y[i], i});
}
for (int i = (int)UX.size(); i < (int)UX.size() * 2; i ++)
sort(V[i].begin(), V[i].end());
for (int i = (int)UX.size() - 1; i; i --)
merge(V[i << 1].begin(), V[i << 1].end(), V[i << 1 ^ 1].begin(), V[i << 1 ^ 1].end(), back_inserter(V[i]));
for (int i = 0; i < (int)UX.size() * 2; i ++)
ST[i] = set < pair < int , int > > (V[i].begin(), V[i].end());
set < pair < int , int > > TS;
for (int i = 1; i <= n; i ++)
TS.insert({-R[i], i});
while (TS.size())
{
while (TS.size() && M[TS.begin()->second])
TS.erase(TS.begin());
if (!TS.size())
break;
int i = TS.begin()->second;
TS.erase(TS.begin());
Do(i);
}
for (int i = 1; i <= n; i ++)
printf("%d ", M[i]);
printf("\n");
return 0;
}
Compilation message
circle_selection.cpp: In function 'int main()':
circle_selection.cpp:67:14: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
scanf("%d", &n);
~~~~~^~~~~~~~~~
circle_selection.cpp:69:53: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
scanf("%d%d%d", &X[i], &Y[i], &R[i]), UX.push_back(X[i]);
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^~~~~~~~~~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
26 ms |
42624 KB |
Output is correct |
2 |
Correct |
25 ms |
42616 KB |
Output is correct |
3 |
Correct |
25 ms |
42624 KB |
Output is correct |
4 |
Correct |
25 ms |
42624 KB |
Output is correct |
5 |
Correct |
25 ms |
42624 KB |
Output is correct |
6 |
Incorrect |
26 ms |
42624 KB |
Output isn't correct |
7 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1386 ms |
395804 KB |
Output is correct |
2 |
Correct |
1476 ms |
395584 KB |
Output is correct |
3 |
Correct |
1492 ms |
395960 KB |
Output is correct |
4 |
Incorrect |
1423 ms |
395140 KB |
Output isn't correct |
5 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
25 ms |
42624 KB |
Output is correct |
2 |
Correct |
564 ms |
152024 KB |
Output is correct |
3 |
Correct |
1908 ms |
394880 KB |
Output is correct |
4 |
Correct |
1896 ms |
394724 KB |
Output is correct |
5 |
Correct |
1913 ms |
384908 KB |
Output is correct |
6 |
Correct |
882 ms |
218968 KB |
Output is correct |
7 |
Correct |
441 ms |
131692 KB |
Output is correct |
8 |
Correct |
86 ms |
58996 KB |
Output is correct |
9 |
Correct |
2027 ms |
395364 KB |
Output is correct |
10 |
Correct |
1970 ms |
395256 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1597 ms |
395952 KB |
Output is correct |
2 |
Correct |
1434 ms |
399080 KB |
Output is correct |
3 |
Correct |
1494 ms |
394740 KB |
Output is correct |
4 |
Correct |
1521 ms |
395108 KB |
Output is correct |
5 |
Correct |
1500 ms |
394228 KB |
Output is correct |
6 |
Correct |
1405 ms |
395344 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
26 ms |
42624 KB |
Output is correct |
2 |
Correct |
25 ms |
42616 KB |
Output is correct |
3 |
Correct |
25 ms |
42624 KB |
Output is correct |
4 |
Correct |
25 ms |
42624 KB |
Output is correct |
5 |
Correct |
25 ms |
42624 KB |
Output is correct |
6 |
Incorrect |
26 ms |
42624 KB |
Output isn't correct |
7 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
26 ms |
42624 KB |
Output is correct |
2 |
Correct |
25 ms |
42616 KB |
Output is correct |
3 |
Correct |
25 ms |
42624 KB |
Output is correct |
4 |
Correct |
25 ms |
42624 KB |
Output is correct |
5 |
Correct |
25 ms |
42624 KB |
Output is correct |
6 |
Incorrect |
26 ms |
42624 KB |
Output isn't correct |
7 |
Halted |
0 ms |
0 KB |
- |