#include <bits/stdc++.h>
#pragma GCC optimize ("Ofast")
#define FOR(i,a,b) for (int i=(a),_b=(b); i<=_b; ++i)
#define FORR(i,a,b) for (int i=(a),_b=(b); i>=_b; --i)
#define REP(i,b) for (int i=0,_b=(b); i<_b; ++i)
#define endl '\n'
#define sz(x) (int) x.size()
#define mod %
#define fillchar(x,y,z) memset(x,z,y)
#define pii pair<int,int>
#define fi first
#define se second
#define mp make_pair
typedef long long int64;
typedef unsigned long long qword;
using namespace std;
int ReadInt()
{
char c;
for (c = getchar(); c!='-'&&(c < '0' || c > '9'); c = getchar());
int nega=0;
int ans = c - '0';
if (c=='-') nega=1,ans=0;
for (c = getchar(); c >= '0' && c <= '9'; c = getchar())
ans = ans * 10 + c - '0';
if (nega) return -ans;
return ans;
}
void WriteInt(int x)
{
if (x > 9) WriteInt(x / 10);
putchar(x % 10 + '0');
}
const int maxn=3e5+5;
inline int64 sqr(int64 x)
{
return x*x;
}
struct Circle{
int64 x,y,r;
int id;
bool IsIntersect(const Circle &c)
{
return sqr(x-c.x)+sqr(y-c.y)<=sqr(r+c.r);
}
}cir[maxn];
bool operator <(const Circle &b,const Circle &c)
{
return mp(-b.r,b.id)<mp(-c.r,c.id);
}
vector<set<pii> > Point;
vector<int> X;
int n;
int ans[maxn];
pii tmp[maxn];
int main()
{
#ifdef meomeomeooooo
freopen("input.txt","r",stdin);
//freopen(".txt","w",stdout);
#endif // meomeomeooooo
iostream::sync_with_stdio(false);
cin.tie(0);
n=ReadInt();
set<Circle> myset;
FOR(i,1,n) {
//cin >> cir[i].x >> cir[i].y >> cir[i].r;
cir[i].x=ReadInt(); cir[i].y=ReadInt(); cir[i].r=ReadInt();
X.push_back(cir[i].x);
cir[i].id=i;
myset.insert(cir[i]);
}
sort(X.begin(),X.end());
X.resize(unique(X.begin(),X.end())-X.begin());
#define GetDec(x) (lower_bound(X.begin(),X.end(),x)-X.begin())
Point.resize(X.size());
FOR(i,1,n) {
int x=GetDec(cir[i].x);
Point[x].insert({cir[i].y,cir[i].id});
}
while (myset.size()) {
Circle cur=*myset.begin(); myset.erase(cur);
ans[cur.id]=cur.id;
int x=cur.x,y=cur.y,r=cur.r;
int Lx=max(1ll*x-2*r,-1000000000ll);
int Rx=min(1ll*x+2*r,1000000000ll);
Lx=GetDec(Lx);
int Ly=max(1ll*y-2*r,-1000000000ll);
int Ry=min(1ll*y+2*r,1000000000ll);
for (int i=Lx; i<X.size()&&X[i]<=Rx; ++i) {
int top=0;
int64 meomeo=sqrt(4ll*r*r-sqr(X[i]-x));
// int Ly=max(1ll*y-meomeo,-1000000000ll);
// int Ry=min(1ll*y+meomeo,1000000000ll);
for (auto j=Point[i].lower_bound(mp(Ly,0)); j!=Point[i].end()&&(*j).fi<=Ry; ++j) {
int id=(*j).se;
if (cur.IsIntersect(cir[id])) {
assert(ans[id]==0||ans[id]==cur.id);
ans[id]=cur.id;
myset.erase(cir[id]);
tmp[++top]=*j;
}
}
FOR(j,1,top)
Point[i].erase(tmp[j]);
}
}
FOR(i,1,n) {
WriteInt(ans[i]); putchar(' ');
}
return 0;
}
Compilation message
circle_selection.cpp: In function 'int main()':
circle_selection.cpp:92:25: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (int i=Lx; i<X.size()&&X[i]<=Rx; ++i) {
~^~~~~~~~~
circle_selection.cpp:94:19: warning: unused variable 'meomeo' [-Wunused-variable]
int64 meomeo=sqrt(4ll*r*r-sqr(X[i]-x));
^~~~~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
384 KB |
Output is correct |
2 |
Correct |
2 ms |
384 KB |
Output is correct |
3 |
Correct |
2 ms |
384 KB |
Output is correct |
4 |
Correct |
2 ms |
384 KB |
Output is correct |
5 |
Correct |
2 ms |
384 KB |
Output is correct |
6 |
Correct |
2 ms |
384 KB |
Output is correct |
7 |
Correct |
2 ms |
384 KB |
Output is correct |
8 |
Correct |
2 ms |
384 KB |
Output is correct |
9 |
Correct |
2 ms |
384 KB |
Output is correct |
10 |
Correct |
2 ms |
384 KB |
Output is correct |
11 |
Correct |
2 ms |
384 KB |
Output is correct |
12 |
Correct |
2 ms |
384 KB |
Output is correct |
13 |
Correct |
3 ms |
384 KB |
Output is correct |
14 |
Correct |
2 ms |
380 KB |
Output is correct |
15 |
Correct |
2 ms |
384 KB |
Output is correct |
16 |
Correct |
4 ms |
512 KB |
Output is correct |
17 |
Correct |
3 ms |
612 KB |
Output is correct |
18 |
Correct |
3 ms |
512 KB |
Output is correct |
19 |
Correct |
12 ms |
1536 KB |
Output is correct |
20 |
Correct |
9 ms |
1380 KB |
Output is correct |
21 |
Correct |
8 ms |
1408 KB |
Output is correct |
22 |
Correct |
10 ms |
1408 KB |
Output is correct |
23 |
Correct |
12 ms |
1408 KB |
Output is correct |
24 |
Correct |
8 ms |
1408 KB |
Output is correct |
25 |
Correct |
9 ms |
1408 KB |
Output is correct |
26 |
Correct |
10 ms |
1408 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1247 ms |
66064 KB |
Output is correct |
2 |
Correct |
1304 ms |
66024 KB |
Output is correct |
3 |
Correct |
1261 ms |
65844 KB |
Output is correct |
4 |
Correct |
1231 ms |
66000 KB |
Output is correct |
5 |
Correct |
1149 ms |
61000 KB |
Output is correct |
6 |
Correct |
1062 ms |
65752 KB |
Output is correct |
7 |
Correct |
1184 ms |
64852 KB |
Output is correct |
8 |
Correct |
1068 ms |
65380 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
512 KB |
Output is correct |
2 |
Correct |
437 ms |
22096 KB |
Output is correct |
3 |
Correct |
1610 ms |
65864 KB |
Output is correct |
4 |
Correct |
1662 ms |
65860 KB |
Output is correct |
5 |
Correct |
2307 ms |
64072 KB |
Output is correct |
6 |
Correct |
779 ms |
34412 KB |
Output is correct |
7 |
Correct |
294 ms |
18284 KB |
Output is correct |
8 |
Correct |
39 ms |
4088 KB |
Output is correct |
9 |
Correct |
1848 ms |
65864 KB |
Output is correct |
10 |
Execution timed out |
3026 ms |
63920 KB |
Time limit exceeded |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
806 ms |
65752 KB |
Output is correct |
2 |
Correct |
576 ms |
65740 KB |
Output is correct |
3 |
Correct |
2271 ms |
65344 KB |
Output is correct |
4 |
Correct |
605 ms |
65788 KB |
Output is correct |
5 |
Correct |
599 ms |
65864 KB |
Output is correct |
6 |
Correct |
2268 ms |
65124 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
384 KB |
Output is correct |
2 |
Correct |
2 ms |
384 KB |
Output is correct |
3 |
Correct |
2 ms |
384 KB |
Output is correct |
4 |
Correct |
2 ms |
384 KB |
Output is correct |
5 |
Correct |
2 ms |
384 KB |
Output is correct |
6 |
Correct |
2 ms |
384 KB |
Output is correct |
7 |
Correct |
2 ms |
384 KB |
Output is correct |
8 |
Correct |
2 ms |
384 KB |
Output is correct |
9 |
Correct |
2 ms |
384 KB |
Output is correct |
10 |
Correct |
2 ms |
384 KB |
Output is correct |
11 |
Correct |
2 ms |
384 KB |
Output is correct |
12 |
Correct |
2 ms |
384 KB |
Output is correct |
13 |
Correct |
3 ms |
384 KB |
Output is correct |
14 |
Correct |
2 ms |
380 KB |
Output is correct |
15 |
Correct |
2 ms |
384 KB |
Output is correct |
16 |
Correct |
4 ms |
512 KB |
Output is correct |
17 |
Correct |
3 ms |
612 KB |
Output is correct |
18 |
Correct |
3 ms |
512 KB |
Output is correct |
19 |
Correct |
12 ms |
1536 KB |
Output is correct |
20 |
Correct |
9 ms |
1380 KB |
Output is correct |
21 |
Correct |
8 ms |
1408 KB |
Output is correct |
22 |
Correct |
10 ms |
1408 KB |
Output is correct |
23 |
Correct |
12 ms |
1408 KB |
Output is correct |
24 |
Correct |
8 ms |
1408 KB |
Output is correct |
25 |
Correct |
9 ms |
1408 KB |
Output is correct |
26 |
Correct |
10 ms |
1408 KB |
Output is correct |
27 |
Correct |
19 ms |
2532 KB |
Output is correct |
28 |
Correct |
17 ms |
2560 KB |
Output is correct |
29 |
Correct |
19 ms |
2688 KB |
Output is correct |
30 |
Correct |
17 ms |
2560 KB |
Output is correct |
31 |
Correct |
22 ms |
2560 KB |
Output is correct |
32 |
Correct |
16 ms |
2560 KB |
Output is correct |
33 |
Correct |
325 ms |
22128 KB |
Output is correct |
34 |
Correct |
307 ms |
22128 KB |
Output is correct |
35 |
Correct |
325 ms |
22128 KB |
Output is correct |
36 |
Correct |
392 ms |
22128 KB |
Output is correct |
37 |
Correct |
323 ms |
22304 KB |
Output is correct |
38 |
Correct |
343 ms |
22308 KB |
Output is correct |
39 |
Correct |
1265 ms |
17512 KB |
Output is correct |
40 |
Correct |
1290 ms |
17520 KB |
Output is correct |
41 |
Correct |
1327 ms |
17520 KB |
Output is correct |
42 |
Correct |
309 ms |
17448 KB |
Output is correct |
43 |
Correct |
199 ms |
22128 KB |
Output is correct |
44 |
Correct |
187 ms |
22128 KB |
Output is correct |
45 |
Correct |
182 ms |
22128 KB |
Output is correct |
46 |
Correct |
183 ms |
22080 KB |
Output is correct |
47 |
Correct |
178 ms |
22128 KB |
Output is correct |
48 |
Correct |
175 ms |
22100 KB |
Output is correct |
49 |
Correct |
182 ms |
22128 KB |
Output is correct |
50 |
Correct |
215 ms |
22172 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
384 KB |
Output is correct |
2 |
Correct |
2 ms |
384 KB |
Output is correct |
3 |
Correct |
2 ms |
384 KB |
Output is correct |
4 |
Correct |
2 ms |
384 KB |
Output is correct |
5 |
Correct |
2 ms |
384 KB |
Output is correct |
6 |
Correct |
2 ms |
384 KB |
Output is correct |
7 |
Correct |
2 ms |
384 KB |
Output is correct |
8 |
Correct |
2 ms |
384 KB |
Output is correct |
9 |
Correct |
2 ms |
384 KB |
Output is correct |
10 |
Correct |
2 ms |
384 KB |
Output is correct |
11 |
Correct |
2 ms |
384 KB |
Output is correct |
12 |
Correct |
2 ms |
384 KB |
Output is correct |
13 |
Correct |
3 ms |
384 KB |
Output is correct |
14 |
Correct |
2 ms |
380 KB |
Output is correct |
15 |
Correct |
2 ms |
384 KB |
Output is correct |
16 |
Correct |
4 ms |
512 KB |
Output is correct |
17 |
Correct |
3 ms |
612 KB |
Output is correct |
18 |
Correct |
3 ms |
512 KB |
Output is correct |
19 |
Correct |
12 ms |
1536 KB |
Output is correct |
20 |
Correct |
9 ms |
1380 KB |
Output is correct |
21 |
Correct |
8 ms |
1408 KB |
Output is correct |
22 |
Correct |
10 ms |
1408 KB |
Output is correct |
23 |
Correct |
12 ms |
1408 KB |
Output is correct |
24 |
Correct |
8 ms |
1408 KB |
Output is correct |
25 |
Correct |
9 ms |
1408 KB |
Output is correct |
26 |
Correct |
10 ms |
1408 KB |
Output is correct |
27 |
Correct |
1247 ms |
66064 KB |
Output is correct |
28 |
Correct |
1304 ms |
66024 KB |
Output is correct |
29 |
Correct |
1261 ms |
65844 KB |
Output is correct |
30 |
Correct |
1231 ms |
66000 KB |
Output is correct |
31 |
Correct |
1149 ms |
61000 KB |
Output is correct |
32 |
Correct |
1062 ms |
65752 KB |
Output is correct |
33 |
Correct |
1184 ms |
64852 KB |
Output is correct |
34 |
Correct |
1068 ms |
65380 KB |
Output is correct |
35 |
Correct |
3 ms |
512 KB |
Output is correct |
36 |
Correct |
437 ms |
22096 KB |
Output is correct |
37 |
Correct |
1610 ms |
65864 KB |
Output is correct |
38 |
Correct |
1662 ms |
65860 KB |
Output is correct |
39 |
Correct |
2307 ms |
64072 KB |
Output is correct |
40 |
Correct |
779 ms |
34412 KB |
Output is correct |
41 |
Correct |
294 ms |
18284 KB |
Output is correct |
42 |
Correct |
39 ms |
4088 KB |
Output is correct |
43 |
Correct |
1848 ms |
65864 KB |
Output is correct |
44 |
Execution timed out |
3026 ms |
63920 KB |
Time limit exceeded |