[Bzoj2716/2648]天使玩偶/SJY摆棋子

题意:平面上有一些点,之后还会在平面上插入一些点,还会询问某一个点到平面中最近的点的距离

$K-Dtree$,最近的点像$[SDOI2010]Hide$ $and$ $Seek$这样求就好了

然后这题是带插入的,和平衡树一样,可能会被卡掉

然而$K-Dtree$又不能旋转,所以只能像替罪羊树那样重构了

然后不会$K-Dtree$的看这里,再看这里,这个太长了,感觉将就着看吧

#include<bits/stdc++.h>
#define fp(i,a,b) for(register int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(register int i=a,I=b-1;i>I;--i)
#define go(u) for(register int i=fi[u],v=e[i].to;i;v=e[i=e[i].nx].to)
#define file(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}
using namespace std;
char ss[1<<17],*A=ss,*B=ss;
inline char gc(){return A==B&&(B=(A=ss)+fread(ss,1,1<<17,stdin),A==B)?-1:*A++;}
template<class T>inline void sd(T&x){
    char c;T y=1;while(c=gc(),(c<48||57<c)&&c!=-1)if(c==45)y=-1;x=c-48;
    while(c=gc(),47<c&&c<58)x=x*10+c-48;x*=y;
}
char sr[1<<21],z[20];int C=-1,Z;
inline void Ot(){fwrite(sr,1,C+1,stdout),C=-1;}
template<class T>inline void we(T x){
    if(C>1<<20)Ot();if(x<0)sr[++C]=45,x=-x;
    while(z[++Z]=x%10+48,x/=10);
    while(sr[++C]=z[Z],--Z);sr[++C]='\n';
}
const int N=5e5+5,M=2*N,inf=-1u>>1;
const double alpha=.77;
typedef int arr[N];
int n,m,rt,D,ans,cb,cn,tot,bin[M];
struct po{
    int d[2],mi[2],mx[2],l,r,sz;
    inline int&operator[](int x){return d[x];}
    inline void clr(){fp(i,0,1)mi[i]=mx[i]=d[i];l=r=0,sz=1;}
    inline bool operator<(po b)const{return d[D]<b[D];}
}a[N],tr[M],t; 
inline int dis(po a,po b){int s=0;fp(i,0,1)s+=abs(a[i]-b[i]);return s;}
struct KD_Tree{
    #define lc tr[p].l
    #define rc tr[p].r
    inline void up(int p){
        tr[p].sz=tr[lc].sz+tr[rc].sz+1;
        fp(i,0,1)
            tr[p].mi[i]=min(tr[p][i],min(tr[lc].mi[i],tr[rc].mi[i])),
            tr[p].mx[i]=max(tr[p][i],max(tr[lc].mx[i],tr[rc].mx[i]));
    }
    inline int New(){return cb?bin[cb--]:++cn;}
    int build(int L,int R,int k){
        int mid=(L+R)>>1,p=New();D=k;
        nth_element(a+L,a+mid,a+R+1);tr[p]=a[mid];
        lc=L<mid?build(L,mid-1,k^1):0;
        rc=mid<R?build(mid+1,R,k^1):0;
        return up(p),p;
    }
    void cle(int&p){if(!p)return;bin[++cb]=p;cle(lc),cle(rc),a[++tot]=tr[p],a[tot].clr(),p=0;}
    void chk(int&p,int k){
        if(!p)return;
        if(tr[lc].sz>tr[p].sz*alpha||tr[rc].sz>tr[p].sz*alpha)
            return tot=0,cle(p),p=build(1,tot,k),void();
        if(t[k]<tr[p][k])chk(lc,k^1);
        else chk(rc,k^1);
    }
    void ins(int&p,int k){
        if(!p)return tr[p=New()]=t,void();
        if(t[k]<tr[p][k])ins(lc,k^1);
        else ins(rc,k^1);up(p);
    }
    inline int h(int p){return!p?inf:(max(tr[p].mi[0]-t[0],0)+max(t[0]-tr[p].mx[0],0)+max(tr[p].mi[1]-t[1],0)+max(t[1]-tr[p].mx[1],0));}
    void qmin(int p){
        int dl=h(lc),dr=h(rc),d0=dis(tr[p],t);
        cmin(ans,d0);
        if(dl<dr){
            if(dl<ans)qmin(lc);
            if(dr<ans)qmin(rc);
        }else{
            if(dr<ans)qmin(rc);
            if(dl<ans)qmin(lc);
        }
    }
}kd;
int main(){
    #ifndef ONLINE_JUDGE
        file("s");
    #endif
    fp(i,0,1)a[0].mx[i]=-inf,a[0].mi[i]=inf;tr[0]=a[0];
    sd(n);sd(m);int op;
    fp(i,1,n){fp(j,0,1)sd(a[i][j]);a[i].clr();}rt=kd.build(1,n,0);
    while(m--){
        sd(op);sd(t[0]),sd(t[1]);t.clr();
        if(op&1)kd.ins(rt,0),kd.chk(rt,0);
        else ans=inf,kd.qmin(rt),we(ans);
    }
return Ot(),0;
}