杭电 1007 Quoit Design

题目地址：http://acm.hdu.edu.cn/showproblem.php?pid=1007

意思是，求N个点中距离最短的2个点的距离

算法导论中的几何章节提到，分治法。倒腾了好一阵子，写写记录一下。

#include<cstdio>
#include<cstring>
#include<queue>
#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;

#define MAX 1000001
struct Node
{
	double x,y;
	bool operator <(const struct Node &t)
	{
		if(x==t.x)
			return y<t.y;
		return x<t.x;
	}
};
Node ivec[MAX],tvec[MAX],rvec[MAX];
double getDis(const Node &a,const Node &b)
{
	return sqrt(pow(a.x-b.x,2)+pow(a.y-b.y,2));
}
double getMin(int start,int end)
{
	if(start==end)
		return 0;
	else if(start+1==end)//2个点
		return getDis(ivec[start],ivec[end]);
	else if(start+2==end)//3个点
	{
		double d1=getDis(ivec[start],ivec[start+1]);
		double d2=getDis(ivec[start],ivec[end]);
		double d3=getDis(ivec[start+1],ivec[end]);
		return min(min(d1,d2),d3);
	}
	else
	{
		int mid=(start+end)/2;//按中线划分
		double d=min(getMin(start,mid),getMin(mid+1,end));//两边的最短距离
		int i,k=0;
		//处理X，就是导论中那个& *2&的矩阵
		for(i=start,k=0;i<=end;i++)
		{
			//两点分布在两边或者一左一右中小于这个最小距离的都需要记录下来
			if(fabs(ivec[i].x-ivec[mid].x)<d)
				tvec[k++]=ivec[i];
		}
		int len=k;
		mid=k/2;
		//此时不再需要进行排序。tvec里已然排好序了，这里是处理Y
		for(i=0,k=0;i<len;i++)
			if(fabs(tvec[i].y-tvec[mid].y)<d)
				rvec[k++]=tvec[i];
		len=k;
		//求所有点中的最短距离
		for(int i=0;i<len-1;i++)
			for(int j=i+1;j<len;j++)
				d=min(d,getDis(rvec[i],rvec[j]));
		return d;
	}
}

int main()
{
    unsigned N;
    while((scanf("%d",&N) != EOF) && N)
    {
        for( unsigned i = 0;i < N; ++ i )
        {
            scanf("%lf%lf",&ivec[i].x,&ivec[i].y);
        }
        sort(ivec,ivec + N);
        printf("%.2lf\n",( getMin(0,N - 1)) / 2);
    }
    return 0;
}