[APIO2018]Circle selection

https://www.zybuluo.com/ysner/note/1257597

题面

在平面上，有\(n\)个圆，记为\(c_1,c_2,...,c_n\)。我们尝试对这些圆运行这个算法：

• 找到这些圆中半径最大的。如果有多个半径最大的圆，选择编号最小的。记为\(c_i\)。
• 删除\(c_i\)及与其相交的所有圆。
• 重复上面两个步骤直到所有的圆都被删除。

当\(c_i\)被删除时，若循环中第\(1\)步选择的圆是\(c_j\)，我们说\(c_i\)被\(c_j\)删除。对于每个圆，求出它是被哪一个圆删除的。

解析

题目太火，正解无从想象，准备暴力骗分。

什么东西能够快速维护平面上的东西？

二维树状数组、二维线段树？蒟蒻不会。。。

于是就只剩\(kd-tree\)了。

\(kd-tree\)本质上就是一个经过优化的暴力（优化在缩小搜索范围，只对近邻进行搜索），是靠剪枝吃饭的。

我们可以用它维护每个圆在\(x,y\)轴的范围，即\([x-r,x+r],[y-r,y+r]\)。

这样相当于维护一个矩形，但是并不会漏掉答案，并且维护起来很方便。

于是建完树后，从前往后对每个圆通过\(kd-tree\)暴力统计答案即可。

然而这样只有\(42pts\)，我们需要更多剪枝。

可以注意到，如果搜到一个圆，在当前统计答案的圆的（矩形）范围之外，这个圆显然不会对答案有贡献，可以跳过。

于是就可以通过洛谷数据了？然而\(loj\)上跑不过最后一档。

把每个点转一个角度就可以\(AC\)。

#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;
}