2017.3.25 SJY摆石子 思考记录

        终于找到了一个kd树的题、、来看看

     k-d树原理还是不难的，就是代码技巧比较多可能可读性比较差

      它就是把一个平面分成：

找最接近的点原理是利用极值矩形：

这个题就是k-d树板子题了、

下面的注释很详细

码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 1000005
#define inf 1000000000 
using namespace std;
int xycmp,i,n,m,ans,rt,op;

struct dian
{                              //注意：下标范围为[2]的  0指的是x坐标  1指的是y坐标 
int xiao[2],da[2],xy[2];
}d[N];                   //对于d集的：只有xy有用，指每个点读入的坐标   

bool cmp(dian x,dian y) 
{   //cout<<xycmp<<" ";
	return (x.xy[xycmp]<y.xy[xycmp]);                            //比较的是这一次比较的坐标的大小 
}

struct kdtree
{
  	int xia[N][2];   //指向这个节点的两个儿子 
  	dian s[N],q;       //对于s集的：   xiao  指的是矩形左下角的坐标集    da指的是矩形右上角的点集      xy指的是这个点的坐标 
	void upd(int now)  //维护矩形 
	{
		int l=xia[now][0],r=xia[now][1]; 
		for(int i=0;i<2;i++)//依次循环 x、y坐标 
		{
		if(l){s[now].da[i]=max(s[l].da[i],s[now].da[i]);s[now].xiao[i]=min(s[now].xiao[i],s[l].xiao[i]);} //如果有左儿子，就取一遍左儿子x、y  矩形  的   坐标最值  再取最值 
		if(r){s[now].da[i]=max(s[r].da[i],s[now].da[i]);s[now].xiao[i]=min(s[now].xiao[i],s[r].xiao[i]);}//如果有右儿子，就取一遍右儿子x、y  矩形  的   坐标最值  再取最值 
    	}
	}


	void add(int now,dian t)//添加 
	{
		for(int i=0;i<2;i++)s[now].xiao[i]=s[now].da[i]=s[now].xy[i]=t.xy[i];	//初始化最大、最小、这个s集（kd树上的） 值都是 当前点的x、y坐标 
	}
	
		int dis(int now)
	{
		int lin=0;    
		for(int i=0;i<2;i++)lin+=max(s[now].xiao[i]-q.xy[i],0);
		for(int i=0;i<2;i++)lin+=max(q.xy[i]-s[now].da[i],0);
		return lin;
		 //这块比较难理解，如果第一个循环>0第二个循环就一定为0  反之亦然   所以它实际上算的是点q到 now点对应的矩形的 曼哈顿距离 
	}
	
	

	void build(int &k,int l,int r,int xydo)//建kd树   k值要返回，因为没有初值   xydo表示当前比较的什么坐标  0是x 1是y 
	{
		k=(l+r)>>1;   xycmp=xydo;   //用k作为编号方便（依次标也没人管你，这不是线段树，可以乱），   xycmp是比较的坐标 
		nth_element(d+l,d+k,d+r+1,cmp);  //这几乎是为kd树量身打造的函数，详情自行百度或见我 stl 算法总结 
		add(k,d[k]);   //加入这个中间点 
		if(l<k)build(xia[k][0],l,k-1,xydo^1); //建左子树 
		if(r>k)build(xia[k][1],k+1,r,xydo^1); //建右子数 
		upd(k);			//维护矩形 
	}   
	
	void ins(int k,int xydo)	  //插入点   k是当前已有点的编号     
	    {
	    	if(q.xy[xydo]<s[k].xy[xydo])   //因为小的在左子树  大的在右子树 
	    	{
	    		if(xia[k][0])ins(xia[k][0],xydo^1);
	    		else xia[k][0]=++n,add(n,q);//没有的话就找到了归宿   加一个 
	    		
			}else
			{
					if(xia[k][1])ins(xia[k][1],xydo^1);     //
	    		else xia[k][1]=++n,add(n,q);                // 同上 
				
			}
	       upd(k);	//不要忘了更新 
		}
		
	void zhao(int k)
	{
		ans=min(ans,(abs(q.xy[0]-s[k].xy[0])+abs(q.xy[1]-s[k].xy[1])));//答案取这个点和要找点的最小距离 
		int ld=(xia[k][0])?dis(xia[k][0]):inf;//如果在矩形内部就是0，在外部就是最小距离，没有这个子节点就是inf 
		int rd=(xia[k][1])?dis(xia[k][1]):inf;
		if(ld<rd)
		{
			if(ld<ans)zhao(xia[k][0]);if(rd<ans)zhao(xia[k][1]);//用访问顺序来优化（优先访问距离短的来更新答案，这样长的就有可能被卡掉） 
		}else
		{	if(rd<ans)zhao(xia[k][1]);if(ld<ans)zhao(xia[k][0]);
		}
	}
}kd;
int main()
{
	scanf("%d%d",&n,&m);
	
	for(i=1;i<=n;i++)
	{
		scanf("%d%d",&d[i].xy[0],&d[i].xy[1]);		
	}
	kd.build(rt,1,n,0);
	
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&op,&kd.q.xy[0],&kd.q.xy[1]);
		if(op==1)
		{
			kd.ins(rt,0);
			
		}else
		{
			ans=inf;kd.zhao(rt);printf("%d\n",ans);  
		}
	}
	return 0;
 }