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)));
    } 
}