【bzoj2648】【sjy摆棋子】【kd树】

Description

这天，SJY显得无聊。在家自己玩。在一个棋盘上，有N个黑色棋子。他每次要么放到棋盘上一个黑色棋子，要么放上一个白色棋子，如果是白色棋子，他会找出距离这个白色棋子最近的黑色棋子。此处的距离是 曼哈顿距离 即(|x1-x2|+|y1-y2|) 。现在给出N<=500000个初始棋子。和M<=500000个操作。对于每个白色棋子，输出距离这个白色棋子最近的黑色棋子的距离。同一个格子可能有多个棋子。

 

Input

第一行两个数 N M

以后M行，每行3个数 t x y

如果t=1 那么放下一个黑色棋子

如果t=2 那么放下一个白色棋子

Output

对于每个T=2 输出一个最小距离

 

Sample Input

2 3
1 1
2 3
2 1 2
1 3 3
2 4 2

Sample Output

1
2

HINT

 

kdtree可以过

题解：比较裸的kd树.

           插入点的时候暴力走到底即可.

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define N 600010
#define inf 2100000000
using namespace std;
int n,x,y,f,opt,rt,m;
struct use{
  int mn[2],mx[2],d[2],l,r;
  int& operator[](int x){return d[x];}
  friend bool operator<(use a,use b){return a[f]<b[f];}
  friend int dis(use a,use b){return abs(a[0]-b[0])+abs(a[1]-b[1]);}
}p[N];
struct kdtree{
  use t[1000010],T;int ans;
  void update(int x){
  	int l=t[x].l,r=t[x].r;
	for (int i=0;i<=1;i++){
  	  t[x].mn[i]=t[x].mx[i]=t[x][i];
  	  if (l) t[x].mn[i]=min(t[x].mn[i],t[l].mn[i]);
  	  if (l) t[x].mx[i]=max(t[x].mx[i],t[l].mx[i]);
  	  if (r) t[x].mn[i]=min(t[x].mn[i],t[r].mn[i]);
  	  if (r) t[x].mx[i]=max(t[x].mx[i],t[r].mx[i]);
  	} 
  }
  int build(int l,int r,int k){
  	f=k;int mid=(l+r)>>1;
  	nth_element(p+l,p+mid,p+r+1);
    t[mid]=p[mid];
    for (int i=0;i<=1;i++)
      t[mid].mx[i]=t[mid].mn[i]=t[mid][i];
    if (l<mid) t[mid].l=build(l,mid-1,k^1);
    if (r>mid) t[mid].r=build(mid+1,r,k^1);
    update(mid);return mid;
  }
  int getmn(use a){
    int ans(0);
	for (int i=0;i<=1;i++){
	   ans+=max(T[i]-a.mx[i],0);
	   ans+=max(a.mn[i]-T[i],0);
	} 
	return ans;
  }
  void qmn(int x){
    ans=min(ans,dis(T,t[x]));
	int l=t[x].l,r=t[x].r,ll=inf,rr=inf;
    if (l) ll=getmn(t[l]);if (r) rr=getmn(t[r]);
	if (ll<rr){
	  if (ll<ans) qmn(l);
	  if (rr<ans) qmn(r);
	} 
	else{
	  if (rr<ans) qmn(r);
	  if (ll<ans) qmn(l);
	}
  }
  void insert(int k,int now){
    int l=t[k].l,r=t[k].r;
	if (T[now]>=t[k][now]){
      if (r) insert(r,now^1);
      else{
        t[k].r=++n;t[n]=T;
        for (int i=0;i<=1;i++)
          t[n].mx[i]=t[n].mn[i]=t[n][i];
	  }
    } 
    else{
      if (l) insert(l,now^1);
	  else{
	    t[k].l=++n;t[n]=T;
	    for (int i=0;i<=1;i++)
	      t[n].mx[i]=t[n].mn[i]=t[n][i];
	  } 
    }
    update(k);
  }
  void ins(int x,int y){
    T[0]=x;T[1]=y;
    insert(rt,0);
  }
  int query(int x,int y){
    T[0]=x;T[1]=y;ans=inf;
    qmn(rt);return ans;
  }
}kd;
int main(){
   scanf("%d%d",&n,&m);
   for (int i=1;i<=n;i++){
     scanf("%d%d",&x,&y);
     p[i][0]=x;p[i][1]=y;
   }
   rt=kd.build(1,n,0);
   for (int i=1;i<=m;i++){
     scanf("%d",&opt);
     if (opt==1){
       scanf("%d%d",&x,&y);
       kd.ins(x,y);
     }
     else{
       scanf("%d%d",&x,&y);
       printf("%d\n",kd.query(x,y));
     } 
   }
}