杭电OJ 1007：Quoit Design

这个一个经典的利用分治算法解决的问题，这个题目其实是要求平面上很多点中距离最近的两个点的距离的一半（半径）。

解决这个问题的思路是：（1）对所有的点按照x坐标从小到大排序。（2）按照x坐标将所有的点一分为二，分为左边的一半和右边一半。

（2）分别求出左边一半和右边的一半的点的最小距离，假设最小距离分别为minl和minr。取minValue=min(minl,minr)。

（3）对从low到high之间的所有的点判断其是否在x0的minValue范围之内，将在这之间的点加入到一个临时数组T中。

（4）对临时数组T中的点按照y坐标从小到大排序，排完序之后对临时数组T中的每一个点判断它和其他点的最短距离，如果这个距离比minValue小，则更新minValue的值。在判断T中的点和其他的点的距离的过程中，需要剪纸，剪纸的条件是：两个点的y坐标的差值大于minValue。

下面说一说自己对这道题目的看法：在上面的第四步中其实不需要每一次递归都要对数组进行排序，其实没有必要做这么多次排序。可以在输入的这些点的时候另外在建一个数组，对这个数组中的点按照y坐标排序，这样就可以减少排序的次数，降低时间复杂度。但是实际上如果这样做的话运行就会超时，估计是数据的问题吧，要不就是先排序的算法常数参数太大，从而导致算法超时。

附上AC代码：

#include <stdio.h>
#include <math.h>
#include <algorithm>
using namespace std;
const int N = 100010;
struct point{
	double x;
	double y;
}arrx[N],/*arry[N],*/T[N];
bool cmp(point p1,point p2){
	return p1.x<p2.x;
}
bool cmp2(point p1,point p2){
	return p1.y<p2.y;
}
double dist(point p1,point p2){
	return sqrt((p1.x - p2.x)*(p1.x - p2.x) + (p1.y - p2.y)*(p1.y - p2.y));
}
double myabs(double x){
	if(x<0)
		x=-x;
	return x;
}
double myMin(double a,double b){
	if(a<b)
		return a;
	return b;
}
double closetPair(int low,int high,int n){
	if(high-low<=2){
		if(high-low==1){
			return dist(arrx[low],arrx[high]);
		}
		else{
			double min1=dist(arrx[low],arrx[low+1]);
			double min2=dist(arrx[low+1],arrx[low+2]);
			double min3=dist(arrx[low],arrx[low+2]);
			double min=min1<min2?min1:min2;
			return min<min3?min:min3;
		}
	}
	else{
		int mid=(low + high)/2;
		double x0=arrx[mid].x;
		double minl=closetPair(low,mid,n);
		double minr = closetPair(mid +1,high,n);
		double minValue = myMin(minl,minr);
		int index = 0;
		for(int i=low;i<=high;i++){
			/*if(myabs(arry[i].x-x0)<=minValue){
			T[index].x=arry[i].x;
			T[index].y=arry[i].y;
			index++;
			}*/
			if(myabs(arrx[i].x - x0)<=minValue){
				T[index].x = arrx[i].x;
				T[index].y = arrx[i].y;
				index++;
			}
		}
		std::sort(T,T+index,cmp2);
		for(int i=0;i<index;i++){
			//int tmpCount=myMin(i+7,index-1);
			for(int j=i+1;j<index;j++){
				if(T[j].y - T[i].y >minValue)
					break;
				double tmpDist =dist(T[i],T[j]);
				if(tmpDist<minValue)
					minValue = tmpDist;
			}
		}
		return minValue;
	}
}
int main(){
	int n;
	//freopen("1.txt","r",stdin);
	while(scanf("%d",&n)!=EOF){
		if(n<=0)
			break;
		for(int i=0;i<n;i++){
			scanf("%lf%lf",&arrx[i].x,&arrx[i].y);
			//arry[i].x=arrx[i].x;
			//arry[i].y=arrx[i].y;
		}
		sort(arrx,arrx+n,cmp);
		//sort(arry,arry+n,cmp2);
		printf("%.2lf\n",closetPair(0,n-1,n)/2);
	}
	return 0;
}