bzoj 2716 [Violet 3]天使玩偶 【CDQ分治】

 KD-tree可做，但是我不会暂时不考虑

大意：在二维平面内，给定n个点，m个操作。操作A：加入一个点；操作B：询问一个点与平面上加入的点的最近距离

不封装会T不封装会T不封装会T不封装会T不封装会T不封装会T不封装会T不封装会T不封装会T不封装会T不封装会T

把初始存在的点也看成加点操作

首先，曼哈顿距离取绝对值很烦，所以我们可以通过转坐标，把左上 右上 左下 右下通过转坐标都变成左下,最后取个min即可。于是对于(x,y)左下的点(x1,y1)，dis=x-x1+y-y1=(x+y)-(x1+y1)，用树状数组维护前缀最小值x1+y1。

按时间序排列，时间序 x y 三维偏序做CDQ分治即可

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 using namespace std;
 const int N=,inf=1e9;
 int n,m,maxx,ans[N];
 struct BIT
 {
     int c[N];
     int lb(int x)
     {
         return x&(-x);
     }
     void update(int x,int d)
     {
         for(int i=x;i<=maxx;i+=lb(i))
             c[i]=max(c[i],d);
     }
     int ques(int x)
     {
         int ret=;
         for(int i=x;i>=;i-=lb(i))
             ret=max(ret,c[i]);
         return ret;
     }
     void clear(int x)
     {
         for(int i=x;i<=maxx;i+=lb(i))
             c[i] = ;
     }
 }bit;
 struct qwe
 {
     int x,y,k,id;
     bool operator < (const qwe& rhs) const
     {
         if(x != rhs.x)
             return x < rhs.x;
         return id < rhs.id;
     }
 }a[N];
 struct CDQ
 {
     int n;
     qwe T[N];
     void init(int n)
     {
         this->n=n;
         sort(a+,a+n+);
     }
     void solve(int L, int R)
     {
         if(L>=R)
             return;
         int M=(L+R)>>;
         int p=L,q=M+;
         for(int i=L;i<=R;i++)
             if(a[i].id<=M)
                 T[p++]=a[i];
             else
                 T[q++]=a[i];
         for(int i=L;i<=R;i++)
             a[i]=T[i];
         solve(L,M);
         solve(M+,R);
         int i=M+,j=L;
         for(;i<=R;i++)
             if(a[i].k==)
             {
                 for(;j<=M&&a[j].x<=a[i].x;j++)
                     if(a[j].k==)
                         bit.update(a[j].y,a[j].x+a[j].y);
                 int t=bit.ques(a[i].y);
                 if(t)
                     ans[a[i].id]=min(ans[a[i].id],a[i].x+a[i].y-t);
             }
         for(int i=L;i<j;i++)
             if(a[i].k==)
                 bit.clear(a[i].y);
         merge(a+L,a+M+,a+M+,a+R+,T+L);
         for(int i=L;i<=R;i++)
             a[i]=T[i];
     }
 }cdq;
 int read()
 {
     int r=;
     char p=getchar();
     while(p>''||p<'')
         p=getchar();
     while(p>=''&&p<='')
     {
         r=r*+p-;
         p=getchar();
     }
     return r;
 }
 int main()
 {
     n=read(),m=read();
     for(int i=;i<=n;i++)
     {
         a[i].x=read()+;
         a[i].y=read()+;
         a[i].id=i;
         a[i].k=;
         maxx=max(maxx,max(a[i].x,a[i].y));
     }
     for(int i=n+;i<=n+m;i++)
     {
         a[i].k=read();
         a[i].x=read()+;
         a[i].y=read()+;
         a[i].id=i;
         maxx=max(maxx,max(a[i].x,a[i].y));
     }
     maxx++;
     n+=m;
     for(int i=;i<=n;i++)
         ans[i]=inf;
     cdq.init(n);
     cdq.solve(,n);
     for(int i=;i<=n;i++)
         a[i].x=maxx-a[i].x;
     cdq.init(n);
     cdq.solve(,n);
     for(int i=;i<=n;i++)
         a[i].y=maxx-a[i].y;
     cdq.init(n);
     cdq.solve(,n);
     for(int i=;i<=n;i++)
         a[i].x=maxx-a[i].x;
     cdq.init(n);
     cdq.solve(,n);
     for(int i=;i<=n;i++)
         if(ans[i]!=inf)
             printf("%d\n",ans[i]);
     return ;
 }