@ysner
2018-08-21T20:46:15.000000Z
字数 2618
阅读 2019
kd-tree
在平面上,有个圆,记为。我们尝试对这些圆运行这个算法:
当被删除时,若循环中第步选择的圆是,我们说被删除。对于每个圆,求出它是被哪一个圆删除的。
题目太火,正解无从想象,准备暴力骗分。
什么东西能够快速维护平面上的东西?
二维树状数组、二维线段树?蒟蒻不会。。。
于是就只剩了。
本质上就是一个经过优化的暴力(优化在缩小搜索范围,只对近邻进行搜索),是靠剪枝吃饭的。
我们可以用它维护每个圆在轴的范围,即。
这样相当于维护一个矩形,但是并不会漏掉答案,并且维护起来很方便。
于是建完树后,从前往后对每个圆通过暴力统计答案即可。
然而这样只有,我们需要更多剪枝。
可以注意到,如果搜到一个圆,在当前统计答案的圆的(矩形)范围之外,这个圆显然不会对答案有贡献,可以跳过。
于是就可以通过洛谷数据了?然而上跑不过最后一档。
把每个点转一个角度就可以。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#define re register
#define il inline
#define ll long long
#define db double
#define ls t[k].l
#define rs t[k].r
#define eps 1e-6
#define pf(x) ((db)(1.0)*(x)*(x))
#define max(a,b) (((a)>(b)?(a):(b)))
#define min(a,b) (((a)<(b)?(a):(b)))
#define fp(i,a,b) for(re int i=a;i<=b;i++)
#define fq(i,a,b) for(re int i=a;i>=b;i--)
using namespace std;
const int N=3e5+100;
int now,n,rt,ans[N];
struct dat
{
int d[2],r,id;
bool operator < (const dat &o) const {return d[now]<o.d[now];}
}a[N];
struct node
{
dat a;
int l,r,mn[2],mx[2];
}t[N];
il int gi()
{
re int x=0,t=1;
re char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
if(ch=='-') t=-1,ch=getchar();
while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
return x*t;
}
il void cmax(re int &x,re int y){x=max(x,y);}
il void cmin(re int &x,re int y){x=min(x,y);}
il bool cmp(dat x,dat y){return (x.r==y.r)?(x.id<y.id):(x.r>y.r);}
struct kd_tree
{
il void upd(re int k,re int p)
{
cmin(t[k].mn[0],t[p].mn[0]);cmax(t[k].mx[0],t[p].mx[0]);
cmin(t[k].mn[1],t[p].mn[1]);cmax(t[k].mx[1],t[p].mx[1]);
}
il void pushup(re int k)
{
re int x=t[k].a.d[0],y=t[k].a.d[1],r=t[k].a.r;
t[k].mn[0]=x-r;t[k].mx[0]=x+r;
t[k].mn[1]=y-r;t[k].mx[1]=y+r;
if(ls) upd(k,ls);if(rs) upd(k,rs);
}
il void Build(re int &k,re int l,re int r,re int tag)
{
re int mid=l+r>>1;now=tag;
nth_element(a+l,a+mid,a+r+1);k=mid;
t[k].a=a[mid];
if(l<mid) Build(ls,l,mid-1,tag^1);else ls=0;
if(mid<r) Build(rs,mid+1,r,tag^1);else rs=0;
pushup(k);
}
il int check(re int k,re dat A)
{
re int x=A.d[0],y=A.d[1],r=A.r+t[k].a.r,xx=t[k].a.d[0],yy=t[k].a.d[1];
return pf(1ll*x-xx)+pf(1ll*y-yy)-eps<=pf(r);
}
il int far(re int k,re dat A)
{
re int x=A.d[0],y=A.d[1],r=A.r;
if(x+r<t[k].mn[0]) return 1;
if(y+r<t[k].mn[1]) return 1;
if(x-r>t[k].mx[0]) return 1;
if(y-r>t[k].mx[1]) return 1;
return 0;
}
il void Query(re int k,re dat A)
{
if(far(k,A)) return;
if(!ans[t[k].a.id]&&check(k,A)) ans[t[k].a.id]=A.id;
if(ls) Query(ls,A);if(rs) Query(rs,A);
}
}kd;
int main()
{
n=gi();
fp(i,1,n)
{
a[i].d[0]=gi();a[i].d[1]=gi();a[i].r=gi();a[i].id=i;
}
kd.Build(rt,1,n,0);
sort(a+1,a+1+n,cmp);
fp(i,1,n) if(!ans[a[i].id]) kd.Query(rt,a[i]);
fp(i,1,n) printf("%d ",ans[i]);puts("");
return 0;
}