bzoj 2648: SJY摆棋子（kd tree）

kd tree，裸的模板题。初学，看着别人模板敲了一遍，一直wa，但还找不到原因。。
学习资料：
kd 树算法之思路篇：https://zhuanlan.zhihu.com/p/22557068
kd 树算法之详细篇：https://www.joinquant.com/post/2843
上边两篇和下边这篇还是有那么一点地方不同的。
K-D tree 数据结构：http://blog.csdn.net/zhjchengfeng5/article/details/7855241

代码参考：http://www.cnblogs.com/ljh2000-jump/p/5513620.html

#include <stdio.h>#include <algorithm>using namespace std;typedef long long LL;const int MAXN = 1000011;const int inf = 0x7fffffff;int n,m;int nowD;int root;int ans;int ql,qr;struct node{    int Min[2],Max[2];    int d[2];    int l,r;} t[MAXN*2];int getint(){    int w=0,q=0;    char c=getchar();    while((c<'0' || c>'9') && c!='-') c=getchar();    if (c=='-')  q=1, c=getchar();    while (c>='0' && c<='9') w=w*10+c-'0', c=getchar();    return q ? -w : w;}bool cmp(node q,node qq){    if(q.d[nowD]==qq.d[nowD]) return q.d[!nowD]<qq.d[!nowD];    return q.d[nowD]<qq.d[nowD];}void kd_updata(int now){    if(t[now].l)    {        if(t[t[now].l].Max[0]>t[now].Max[0]) t[now].Max[0]=t[t[now].l].Max[0];        if(t[t[now].l].Max[1]>t[now].Max[1]) t[now].Max[1]=t[t[now].l].Max[1];        if(t[t[now].l].Min[0]<t[now].Min[0]) t[now].Min[0]=t[t[now].l].Min[0];        if(t[t[now].l].Min[1]<t[now].Min[1]) t[now].Min[1]=t[t[now].l].Min[1];    }    if(t[now].r)    {        if(t[t[now].r].Max[0]>t[now].Max[0]) t[now].Max[0]=t[t[now].r].Max[0];        if(t[t[now].r].Max[1]>t[now].Max[1]) t[now].Max[1]=t[t[now].r].Max[1];        if(t[t[now].r].Min[0]<t[now].Min[0]) t[now].Min[0]=t[t[now].r].Min[0];        if(t[t[now].r].Min[1]<t[now].Min[1]) t[now].Min[1]=t[t[now].r].Min[1];    }}int kd_build(int l,int r,int D){    int mid=(l+r) >> 1;    nowD=D;    nth_element(t+l+1,t+mid+1,t+r+1,cmp);    if(l!=mid)        t[mid].l=kd_build(l,mid-1,!D);    if(r!=mid)        t[mid].r=kd_build(mid+1,r,!D);    t[mid].Max[0]=t[mid].Min[0]=t[mid].d[0];    t[mid].Max[1]=t[mid].Min[1]=t[mid].d[1];    kd_updata(mid);    return mid;}int dist(int p){    int dis=0;    if(ql < t[p].Min[0])        dis+=t[p].Min[0]-ql;    if(ql > t[p].Max[0])        dis+=ql-t[p].Max[0];    if(qr < t[p].Min[1])        dis+=t[p].Min[1]-qr;    if(qr > t[p].Max[1])        dis+=qr-t[p].Max[1];    return dis;}void kd_query(int p){    int dl,dr,d0;    d0=abs(t[p].d[0]-ql)+abs(t[p].d[1]-qr);    if(d0<ans) ans=d0;    if(t[p].l) dl=dist(t[p].l);    else dl=inf;    if(t[p].r) dr=dist(t[p].r);    else dr=inf;    if(dl<dr)    {        if(dl<ans) kd_query(t[p].l);        if(dr<ans) kd_query(t[p].r);    }    else    {        if(dr<ans) kd_query(t[p].r);        if(dl<ans) kd_query(t[p].l);    }}inline void kd_insert(int now){    int p=root,D=0;    while(true)    {        if(t[now].Max[0]>t[p].Max[0])            t[p].Max[0]=t[now].Max[0];        if(t[now].Max[1]>t[p].Max[1])            t[p].Max[1]=t[now].Max[1];        if(t[now].Min[0]<t[p].Min[0])            t[p].Min[0]=t[now].Min[0];        if(t[now].Min[1]<t[p].Min[1])            t[p].Min[1]=t[now].Min[1];        if(t[now].d[D]>=t[p].d[D])        {            if(!t[p].r)            {                t[p].r=now;                return;            }            else p=t[p].r;        }        else        {            if(!t[p].l)            {                t[p].l=now;                return ;            }            else p=t[p].l;        }        D=!D;    }}int main(){    scanf("%d %d",&n,&m);    for(int i=1; i<=n; i++)        scanf("%d %d",&t[i].d[0],&t[i].d[1]);    root=kd_build(1,n,0);    int x,y,z;    for(int i=1; i<=m; i++)    {        scanf("%d %d %d",&x,&y,&z);        if(x==1)        {            ++n;            t[n].Max[0]=t[n].Min[0]=t[n].d[0]=y;            t[n].Max[1]=t[n].Min[1]=t[n].d[1]=z;            kd_insert(n);        }        else        {            ans=inf;            ql=y,qr=z;            kd_query(root);            printf("%d\n",ans);        }    }    return 0;}