| # | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
|---|---|---|---|---|---|---|---|
| 49981 | zetapi | 원 고르기 (APIO18_circle_selection) | C++14 | 0 ms | 0 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
using namespace std;
#define pb push_back
#define mp make_pair
#define int long long
#define itr ::iterator
typedef pair<int,int> pii;
const int MAX=4e5;
const int INF=1e12;
vector< pair<int,pii> > vec;
set<pii> st;
set<pii> itr it,it2;
int N,ind,X[MAX],Y[MAX],R[MAX],par[MAX];
void ok(int F,int S)
{
if((X[F]-X[S])*(X[F]-X[S])+(Y[F]-Y[S])*(Y[F]-Y[S])>(R[F]+R[S])*(R[F]+R[S]))
return ;
if(R[F]==R[S] and F>S)
swap(F,S);
if(R[F]>R[S])
swap(F,S);
par[F]=S;
return ;
}
signed main()
{
ios_base::sync_with_stdio(false);
/*cin.tie(0);
cout.tie(0);*/
cin>>N;
for(int A=1;A<=N;A++)
{
par[A]=A;
cin>>X[A]>>Y[A]>>R[A];
vec.pb(mp(Y[A]-R[A],mp(0,A)));
vec.pb(mp(Y[A]+R[A],mp(1,A)));
}
sort(vec.begin(),vec.end(),cmp);
for(int A=0;A<vec.size();A++)
{
int pre=0;
for(auto A:st)
{
if(pre)
ok(pre,A.second);
pre=A.second;
}
ind=vec[A].second.second;
if(vec[A].second.first==0)
{
st.insert(mp(X[ind],ind));
/*it=st.find(mp(X[ind],ind));
it++;
if(it!=st.end())
ok(ind,it->second);
it--;
if(it!=st.begin())
{
it--;
ok(ind,it->second);
}*/
}
else
{
it=st.find(mp(X[ind],ind));
/*it2=it;
it2++;
if(it2!=st.end() and it!=st.begin())
{
it--;
ok(it->second,it2->second);
}*/
st.erase(mp(X[ind],ind));
}
}
int pre=0;
for(auto A:st)
{
if(pre)
ok(pre,A.second);
pre=A.second;
}
for(int A=1;A<=N;A++)
cout<<par[A]<<" ";
return 0;
}