最近对问题

问题：设p1=(x1,y1),p2=(x2,y2),…,pn=(xn,yn)是平面上n个点构成的集合S,找出集合S中距离最近的点对。

思路：(1)划分:将集合S分成两个子集S1和S2.
1.pi属于S1，pj属于S2
2.pi属于S1，pj属于S1
3.pi属于S2，pj属于S2
（2）求解子问题：递归求解子问题
（3）合并：取各划分情况下的最近点对的距离最小者
情况3
选取垂直线x=m作为分割线，由此S1={p属于S|xp<=m},S2={q属于S|xq>m}。递归求解得到S1中的最近距离d1,S2中的最大距离d2,d=min(d1,d2)。
p(x,y)是集合P1和集合P2中y坐标最小的点，点p可能在P1或P2中，找出和点P距离小于d的点，y坐标位于[y,y+d]。
将P1和P2中的点按照y坐标升序排列，在[y,y+d]内最多选出8个候选点，计算它们和点p之间的距离。

T(n)=1,k=2
T(n)=2*T(n/2)+n.k>2
T(n)=O(nlog2n)

struct point{
	int x,y;
};
double Closest(point S[],int low,int high){
	double d1,d2,d3,d;
	int mid,i,j,index;
	point P[n];
	if(high-low==1)
		return Distance(S[low],S[high]);
	if(high-low==2){
	d1=Distance(S[low],S[low+1]);
	d1=Distance(S[low+1],S[high]);
	d3=Distance(S[low],S[high]);
	if((d1<d2)&&(d1<d3))
		return d1;
	else if(d2<d3) return d2;
	else return d3;
	}
	mid=(low+high)/2;
	d1=Closest(S,low,mid);
	d2=Closest(S,mid+1,high);
	if(d1<=d2) d=d1;
	else d=d2;
	index=0;
	for(i=mid;(i>=low)&&(s[mid].x-S[i].x<d);i--)
		P[index++]=S[i];
	for(i=mid+1;(i<=high)&&(S[i].x-s[mid].x<d);i--)
		P[index++]=S[i];
	QuickSort(P,0,index-1);
	for(i=0;i<index;i++){
		for(j=i+1;j<index;j++){
			if(P[j].y-P[i].y>=d)
				break;
				else
				{
				d3=Distance(P[i],P[j]);
				id(d3<d) d=d3;
				}
			}
		}
		return d;
	}
	double Distance(point a,point b)
	{
		return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
	}