bzoj2648SJY摆棋子&&bzoj2716[Violet 3]天使玩偶*

bzoj2648SJY摆棋子

bzoj2716[Violet 3]天使玩偶

题意：

棋盘上有n个棋子，现在有m个操作，一种是加棋子，一种是查询离某个点最近的棋子。n，m≤500000。

题解：

先将已有的棋子建kd树，然后加棋子就直接向kd树插入节点。因为本题数据弱，所以直接插节点不会T，如果是一些数据比较强的题目，需要在插入一定量节点后重构整棵树。

代码：

 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #define inc(i,j,k) for(int i=j;i<=k;i++)
 #define maxn 500010
 #define INF 0x3fffffff
 using namespace std;

 inline int read(){
     char ch=getchar(); int f=,x=;
     while(ch<''||ch>''){if(ch=='-')f=-; ch=getchar();}
     while(ch>=''&&ch<='')x=x*+ch-'',ch=getchar();
     return f*x;
 }
 int n,m,f,rt,ans;
 struct p{int pos[]; bool operator < (const p &a)const{return pos[f]<a.pos[f];}}ps[maxn];
 struct nd{p pos; int mx[],mn[],lc,rc;}nds[maxn*];
 int dis(p a,p b){return abs(a.pos[]-b.pos[])+abs(a.pos[]-b.pos[]);}
 void update(int x){
     inc(i,,){
         if(nds[x].lc)
             nds[x].mx[i]=max(nds[x].mx[i],nds[nds[x].lc].mx[i]),
             nds[x].mn[i]=min(nds[x].mn[i],nds[nds[x].lc].mn[i]);
         if(nds[x].rc)
             nds[x].mx[i]=max(nds[x].mx[i],nds[nds[x].rc].mx[i]),
             nds[x].mn[i]=min(nds[x].mn[i],nds[nds[x].rc].mn[i]);
     }
 }
 int build(int l,int r,int now){
     f=now; int mid=(l+r)>>; nth_element(ps+l,ps+mid,ps+r+);
     inc(i,,)nds[mid].mx[i]=nds[mid].mn[i]=ps[mid].pos[i]; nds[mid].pos=ps[mid];
     if(l<mid)nds[mid].lc=build(l,mid-,now^); if(mid<r)nds[mid].rc=build(mid+,r,now^);
     update(mid); return mid;
 }
 void insert(int &x,p a,int now){
     if(!x){x=++n; inc(i,,)nds[x].mx[i]=nds[x].mn[i]=a.pos[i]; nds[x].pos=a; return;}
     f=now; if(a<nds[x].pos)insert(nds[x].lc,a,now^);else insert(nds[x].rc,a,now^); update(x);
 }
 int get(int x,p a){
     int q=; inc(i,,)q+=max(a.pos[i]-nds[x].mx[i],),q+=max(nds[x].mn[i]-a.pos[i],); return q;
 }
 void query(int x,p a){
     ans=min(ans,dis(a,nds[x].pos)); int dl=INF-,dr=INF-;
     if(nds[x].lc)dl=get(nds[x].lc,a); if(nds[x].rc)dr=get(nds[x].rc,a);
     if(dl<dr){
         if(ans>dl)query(nds[x].lc,a); if(ans>dr)query(nds[x].rc,a);
     }else{
         if(ans>dr)query(nds[x].rc,a); if(ans>dl)query(nds[x].lc,a);
     }
 }
 int main(){
     n=read(); m=read(); inc(i,,n)ps[i].pos[]=read(),ps[i].pos[]=read(); rt=build(,n,);
     inc(i,,m){
         int t=read(),x=read(),y=read();
         if(t==)insert(rt,(p){x,y},);else ans=INF,query(rt,(p){x,y}),printf("%d\n",ans);
     }
     return ;
 }

20160906