BZOJ 2648 SJY摆棋子

题目传送门
分析：
就一KD树板题。。。

这里其实可以先建树在直接暴力加点，然后查询就好了。。。
但是这样的做法BZOJ能过，洛谷上会T（应该是第11号点是后来加的hack数据（难受））

所以考虑离线，然后在每个黑子上打是否已添加的标记。。。
很有道理2333
但是BZOJ上T了，洛谷上过了

呵呵呵。。。呵呵呵呵。。。。哈哈哈哈

然后就有神仙说沿用一下替罪羊树的思想，过于畸形的子树直接拍扁重构
那这样左旋右旋也可以呗。。。。

然后还有神仙会CDQ分治，貌似也挺有道理的。。。

我前两种都写了一下。。。其实差不多吧，都很好写。。。

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>

#define maxn 1000005
#define INF 0x3f3f3f3f

using namespace std;

inline int getint()
{
    int num=0,flag=1;char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=num*10+c-48,c=getchar();
    return num*flag;
}

int n,m,rt;
struct node{
    int x[2];
}P[maxn];
struct pt{
    node p;
    int lc,rc;
    int mx[2],mn[2];
}a[maxn];
int ans,cur;
inline bool cmp0(node x,node y){return x.x[0]<y.x[0];}
inline bool cmp1(node x,node y){return x.x[1]<y.x[1];}

inline void pushup(int x,int y)
{
    a[x].mx[0]=max(a[x].mx[0],a[y].mx[0]);
    a[x].mx[1]=max(a[x].mx[1],a[y].mx[1]);
    a[x].mn[0]=min(a[x].mn[0],a[y].mn[0]);
    a[x].mn[1]=min(a[x].mn[1],a[y].mn[1]);
}

inline void build(int &now,int l,int r,int k)
{
    int mid=now=(l+r)>>1;
    if(k)nth_element(P+l+1,P+mid+1,P+r+1,cmp1);
    else nth_element(P+l+1,P+mid+1,P+r+1,cmp0);
    a[now].p=P[mid];
    a[now].mn[0]=a[now].mx[0]=P[mid].x[0];
    a[now].mn[1]=a[now].mx[1]=P[mid].x[1];
    if(l<now)build(a[now].lc,l,mid-1,k^1),pushup(now,a[now].lc);
    if(r>now)build(a[now].rc,mid+1,r,k^1),pushup(now,a[now].rc);
}

inline int getdis(node x,node y)
{return abs(x.x[0]-y.x[0])+abs(x.x[1]-y.x[1]);}

inline int msr(int now,node p)
{
    if(!now)return INF;
    int tmp1=(a[now].mn[0]<=p.x[0]&&a[now].mx[0]>=p.x[0])?0:min(abs(p.x[0]-a[now].mn[0]),abs(p.x[0]-a[now].mx[0]));
    int tmp2=(a[now].mn[1]<=p.x[1]&&a[now].mx[1]>=p.x[1])?0:min(abs(p.x[1]-a[now].mn[1]),abs(p.x[1]-a[now].mx[1]));
    return tmp1+tmp2;
}

inline void insert(int &now,node p,int k)
{
    if(!now)
    {
        now=++cur,a[now].p=p;
        a[now].mn[0]=a[now].mx[0]=p.x[0];
        a[now].mn[1]=a[now].mx[1]=p.x[1];
        return;
    }
    if(p.x[k]<=a[now].p.x[k])insert(a[now].lc,p,k^1),pushup(now,a[now].lc);
    else insert(a[now].rc,p,k^1),pushup(now,a[now].rc);
}

inline void query(int rt,node p)
{
    int tmp=getdis(p,a[rt].p);
    if(tmp<ans)ans=tmp;
    int tmp1=msr(a[rt].lc,p),tmp2=msr(a[rt].rc,p);
    if(tmp1<tmp2)
    {
        if(a[rt].lc&&tmp1<ans)query(a[rt].lc,p);
        if(a[rt].rc&&tmp2<ans)query(a[rt].rc,p);
    }
    else
    {
        if(a[rt].rc&&tmp2<ans)query(a[rt].rc,p);
        if(a[rt].lc&&tmp1<ans)query(a[rt].lc,p);
    }
}

int main()
{
    n=cur=getint(),m=getint();
    for(int i=1;i<=n;i++)
        P[i].x[0]=getint(),P[i].x[1]=getint();
    build(rt,1,n,0);
    for(int i=1;i<=m;i++)
    {
        int op=getint();
        node tmp;tmp.x[0]=getint(),tmp.x[1]=getint();
        if(op==1)insert(rt,tmp,0);
        else ans=INF,query(rt,tmp),printf("%d\n",ans);
    }
}