bzoj2716: [Violet 3]天使玩偶

传送门
这题好像有CDQ分治做法
不过我太菜了没听懂
于是我打了一发kdtree
度娘：
k-d树（k-dimensional树的简称），是一种分割k维数据空间的数据结构。主要应用于多维空间关键数据的搜索（如：范围搜索和最近邻搜索）。K-D树是二进制空间分割树的特殊的情况。
大概思路就是每次确定一个轴，将点集平均分成两份，再递归建树。
至于怎么划分：
我们可以按照方差最大的方案来做。
但是会产生这样的情况：

所以我们强制每一次划分坐标必须和上一次不同。
然后长成这样：

这就好看多了。
每一次插入我们按照划分顺序找到一个新的节点，并且插入进去。
好像可以类似于替罪羊树的思路来保证时间复杂度。
对于每一询问最近节点，我们可以大力求出边界离要求的点P最近的部分先递归，再递归较远的部分。
但是会出现这种情况：

这样时间复杂度会退化到O（N^2)
但是你要相信出题人是不会来卡你的kdtree的
做这题就是kdtree的简单应用。
看着代码理解一下，你会发现我是选择看hzwer的、

#include<cmath>#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>#include<algorithm>using namespace std;int n,m,fl,x,y,rt,D;struct lzh{    int d[2],mn[2],mx[2],l,r;    int operator [](int x){        return d[x];    }    bool operator <(const lzh &x)const{        return d[D]<x.d[D];    }    lzh(int x=0,int y=0){        d[0]=x; d[1]=y;        l=r=0;    }}p[500005];int dis(lzh a,lzh b){    return abs(a.d[0]-b.d[0])+abs(a.d[1]-b.d[1]);}struct Lzh{    int ans;    lzh t[1000005],T;    void update(int k){        lzh l=t[t[k].l],r=t[t[k].r];        for (int i=0;i<2;i++){            if (t[k].l){                t[k].mn[i]=min(t[k].mn[i],l.mn[i]);                t[k].mx[i]=max(t[k].mx[i],l.mx[i]);            }            if (t[k].r){                t[k].mn[i]=min(t[k].mn[i],r.mn[i]);                t[k].mx[i]=max(t[k].mx[i],r.mx[i]);            }        }    }    int build(int l,int r,int now){        D=now;        int mid=(l+r)/2;        nth_element(p+l,p+mid,p+r+1);        t[mid]=p[mid];        for (int i=0;i<2;i++)            t[mid].mn[i]=t[mid].mx[i]=t[mid][i];        if (l<mid) t[mid].l=build(l,mid-1,now^1);        if (r>mid) t[mid].r=build(mid+1,r,now^1);        update(mid);        return mid;    }    int get(int k,lzh p){        int s=0;        for (int i=0;i<2;i++)            s+=max(0,t[k].mn[i]-p[i])+max(0,p[i]-t[k].mx[i]);        return s;    }    void insert(int k,int now){        for (int i=1;i<=2000000000;i++);        if (T[now]>=t[k][now]){            if (t[k].r) insert(t[k].r,now^1);            else{                t[k].r=++n; t[n]=T;                for (int i=0;i<2;i++)                    t[n].mn[i]=t[n].mx[i]=t[n][i];            }        }        else{            if (t[k].l) insert(t[k].l,now^1);            else{                t[k].l=++n; t[n]=T;                for (int i=0;i<2;i++)                    t[n].mn[i]=t[n].mx[i]=t[n][i];            }        }        update(k);    }    void ask(int k,int now){        int d=dis(t[k],T),dl=2e9,dr=2e9;        ans=min(ans,d);        if (t[k].l) dl=get(t[k].l,T);        if (t[k].r) dr=get(t[k].r,T);        if (dl<dr){            if (dl<ans) ask(t[k].l,now^1);            if (dr<ans) ask(t[k].r,now^1);        }        else{            if (dr<ans) ask(t[k].r,now^1);            if (dl<ans) ask(t[k].l,now^1);        }    }    int ask(lzh p){        ans=2e9; T=p;        ask(rt,0);        return ans;    }    void insert(lzh p){        T=p;        insert(rt,0);    }}Lzh;int main(){    scanf("%d%d",&n,&m);    for (int i=1;i<=n;i++)        scanf("%d%d",&p[i].d[0],&p[i].d[1]);    rt=Lzh.build(1,n,0);    while (m--){        scanf("%d%d%d",&fl,&x,&y);        if (fl==1) Lzh.insert(lzh(x,y));        else printf("%d\n",Lzh.ask(lzh(x,y)));    } }