bzoj2648SJY摆棋子

bzoj2648SJY摆棋子

bzoj2716[Violet 3]天使玩偶

题意：

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

题解：

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

代码：

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <algorithm>
 4 #define inc(i,j,k) for(int i=j;i<=k;i++)
 5 #define maxn 500010
 6 #define INF 0x3fffffff
 7 using namespace std;
 8
 9 inline int read(){
10     char ch=getchar(); int f=1,x=0;
11     while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1; ch=getchar();}
12     while(ch>=‘0‘&&ch<=‘9‘)x=x*10+ch-‘0‘,ch=getchar();
13     return f*x;
14 }
15 int n,m,f,rt,ans;
16 struct p{int pos[2]; bool operator < (const p &a)const{return pos[f]<a.pos[f];}}ps[maxn];
17 struct nd{p pos; int mx[2],mn[2],lc,rc;}nds[maxn*2];
18 int dis(p a,p b){return abs(a.pos[0]-b.pos[0])+abs(a.pos[1]-b.pos[1]);}
19 void update(int x){
20     inc(i,0,1){
21         if(nds[x].lc)
22             nds[x].mx[i]=max(nds[x].mx[i],nds[nds[x].lc].mx[i]),
23             nds[x].mn[i]=min(nds[x].mn[i],nds[nds[x].lc].mn[i]);
24         if(nds[x].rc)
25             nds[x].mx[i]=max(nds[x].mx[i],nds[nds[x].rc].mx[i]),
26             nds[x].mn[i]=min(nds[x].mn[i],nds[nds[x].rc].mn[i]);
27     }
28 }
29 int build(int l,int r,int now){
30     f=now; int mid=(l+r)>>1; nth_element(ps+l,ps+mid,ps+r+1);
31     inc(i,0,1)nds[mid].mx[i]=nds[mid].mn[i]=ps[mid].pos[i]; nds[mid].pos=ps[mid];
32     if(l<mid)nds[mid].lc=build(l,mid-1,now^1); if(mid<r)nds[mid].rc=build(mid+1,r,now^1);
33     update(mid); return mid;
34 }
35 void insert(int &x,p a,int now){
36     if(!x){x=++n; inc(i,0,1)nds[x].mx[i]=nds[x].mn[i]=a.pos[i]; nds[x].pos=a; return;}
37     f=now; if(a<nds[x].pos)insert(nds[x].lc,a,now^1);else insert(nds[x].rc,a,now^1); update(x);
38 }
39 int get(int x,p a){
40     int q=0; inc(i,0,1)q+=max(a.pos[i]-nds[x].mx[i],0),q+=max(nds[x].mn[i]-a.pos[i],0); return q;
41 }
42 void query(int x,p a){
43     ans=min(ans,dis(a,nds[x].pos)); int dl=INF-1,dr=INF-1;
44     if(nds[x].lc)dl=get(nds[x].lc,a); if(nds[x].rc)dr=get(nds[x].rc,a);
45     if(dl<dr){
46         if(ans>dl)query(nds[x].lc,a); if(ans>dr)query(nds[x].rc,a);
47     }else{
48         if(ans>dr)query(nds[x].rc,a); if(ans>dl)query(nds[x].lc,a);
49     }
50 }
51 int main(){
52     n=read(); m=read(); inc(i,1,n)ps[i].pos[0]=read(),ps[i].pos[1]=read(); rt=build(1,n,0);
53     inc(i,1,m){
54         int t=read(),x=read(),y=read();
55         if(t==1)insert(rt,(p){x,y},0);else ans=INF,query(rt,(p){x,y}),printf("%d\n",ans);
56     }
57     return 0;
58 }

20160906