bzoj4533 [BeiJing2014 WinterCamp] 数据（K-D tree）

模板题改了两天，wc要gg了。
曼哈顿距离最小值注意分类讨论。

#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
using namespace std;
#define ll long long
#define inf 0x7fffffff
#define N 200010
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,D,rt,ans;
struct point{
    int d[2];
    int& operator[](int x){return d[x];}
    friend bool operator<(point a,point b){return a[D]<b[D];}
}P[N],T;
struct node{
    point x;int lc,rc,mx[2],mn[2];
}tr[N];
inline void pushup(int p){
    int l=tr[p].lc,r=tr[p].rc;
    for(int i=0;i<2;++i){
        tr[p].mx[i]=max(max(tr[l].mx[i],tr[r].mx[i]),tr[p].x[i]);
        tr[p].mn[i]=min(min(tr[l].mn[i],tr[r].mn[i]),tr[p].x[i]);
    }
}
inline void build(int &p,int l,int r,int op){
    int mid=l+r>>1;p=mid;D=op;tr[p].lc=tr[p].rc=0;
    nth_element(P+l,P+mid,P+r+1);tr[p].x=P[mid];
    for(int i=0;i<2;++i) tr[p].mx[i]=tr[p].mn[i]=tr[p].x[i];
    if(l<mid) build(tr[p].lc,l,mid-1,op^1);
    if(r>mid) build(tr[p].rc,mid+1,r,op^1);
    pushup(p);
}
inline void ins(int &p,int op){
    if(!p){
        p=n;tr[p].x=P[n];
        for(int i=0;i<2;++i) tr[p].mx[i]=tr[p].mn[i]=tr[p].x[i];return;
    }if(P[n][op]<=tr[p].x[op]) ins(tr[p].lc,op^1);
    else ins(tr[p].rc,op^1);pushup(p); 
}
inline int abs(int x){return x<0?-x:x;}
inline int dis(point a,point b){return abs(a[0]-b[0])+abs(a[1]-b[1]);}
inline int calmn(int p){
    if(!p) return inf;int res=0;
    for(int i=0;i<2;++i){
        if(T[i]<tr[p].mn[i]) res+=tr[p].mn[i]-T[i];
        if(T[i]>tr[p].mx[i]) res+=T[i]-tr[p].mx[i];
    }return res;
}
inline int calmx(int p){
    if(!p) return 0;int res=0;
    for(int i=0;i<2;++i) res+=max(abs(tr[p].mn[i]-T[i]),abs(tr[p].mx[i]-T[i]));
    return res;
}
inline void qmin(int p){
    ans=min(ans,dis(tr[p].x,T));
    int dl=calmn(tr[p].lc),dr=calmn(tr[p].rc);
    if(dl<dr){
        if(dl<ans) qmin(tr[p].lc);
        if(dr<ans) qmin(tr[p].rc);
    }else{
        if(dr<ans) qmin(tr[p].rc);
        if(dl<ans) qmin(tr[p].lc);
    }
}
inline void qmax(int p){
    ans=max(ans,dis(tr[p].x,T));
    int dl=calmx(tr[p].lc),dr=calmx(tr[p].rc);
    if(dl>dr){
        if(dl>ans) qmax(tr[p].lc);
        if(dr>ans) qmax(tr[p].rc);
    }else{
        if(dr>ans) qmax(tr[p].rc);
        if(dl>ans) qmax(tr[p].lc);
    }
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();
    for(int i=0;i<2;++i) tr[0].mn[i]=inf,tr[0].mx[i]=0;
    for(int i=1;i<=n;++i) P[i][0]=read(),P[i][1]=read();
    build(rt,1,n,0);m=read();
    while(m--){
        int op=read(),x=read(),y=read();
        if(op==0){
            P[++n][0]=x;P[n][1]=y;
            if(n%30000) ins(rt,0);
            else build(rt,1,n,0);continue;
        }T[0]=x;T[1]=y;
        if(op==1){ans=inf;qmin(rt);printf("%d\n",ans);}
        else {ans=0;qmax(rt);printf("%d\n",ans);}
    }return 0;
}