Hdu 1007Quoit Design 归并法

  题目大意:
 

  给出二维平面上的n个点，求其中最近的两个点的距离的一半。

  

 采用分治方法，先把n个点按照x坐标排序，分成左右两个部分，再对每个部分求最近点的距离，然后合并

  合并过程应该检查宽为2d的带状区间是否有两个点分别属于两个集合且距离小于d，最多可能有n个点，

合并最坏情况下O(n^2) 

 对于左边的点，右边中的点必定落在一个d*2d的矩形中，且最多检查六个点(鸽巢原理)

 这样，先将带状区间的点按照y坐标排序，然后线性扫描。

#include<iostream>
#include<cmath>
#include<algorithm>
#include<stdio.h>
#include<iomanip>
using namespace std;
double MAX = 1e10;
int a,b;  //a,b返回的最近点的编号从0开始 
struct Node{
	double x,y;
	int key;
}
ar[100005],br[100005];
bool cmpx(Node a,Node b){
	return a.x<b.x;
}
bool cmpy(Node a,Node b){
	return a.y<b.y;
}
double min(double a,double b){
	return a<b?a:b;
}
double dis(Node a,Node b){
	return sqrt(pow(a.x-b.x,2)+pow(a.y-b.y,2));
}
double cal(int s,int e){
	int mid,i,j,tail=0;
	double d;
	if(s==e){
		return MAX;
	} 
	mid=(s+e)/2;
	d=min(cal(s,mid),cal(mid+1,e));
	for(i=mid;i>=s&&ar[mid].x-ar[i].x<d;i--){
		br[tail++]=ar[i];
	}
	for(i=mid+1;i<e&&ar[i].x-ar[mid].x<d;i++){
		br[tail++]=ar[i];
	}
	sort(br,br+tail,cmpy);
	for(i=0;i<tail;i++){
		for(j=i+1;j<tail&&br[j].y-br[i].y<d;j++){
			if(d>dis(br[i],br[j])){
				a=min(br[i].key,br[j].key);
				b=br[i].key+br[j].key;
				d=min(d,dis(br[i],br[j]));
			}
		}
	}
	return d;
}
int main(){
    int n;
	while(cin>>n&&n){
		for(int i=0;i<n;i++){
			ar[i].key=i+1;
			scanf("%lf%lf",&ar[i].x,&ar[i].y);
		}
		sort(ar,ar+n,cmpx);
		double d = cal(0,n);
		printf("%.2lf\n",d/2.0);
	}	
	return 0;
}