HDU 1007 Quoit Design（平面最小点对）

题记：首先将这n个点按照x从小到大排序，用分治的思想把这些点拆分，使得两边点的数目都是当前点的个数的一半。

S集合分成左集合S1和右集合S2
d记录当前状态下的最小距离

分情况讨论合并时的问题：

1、这个集合只有一个点，返回一个很大的数。（取最小值时会去掉这个数）

2、这个集合有两个点，返回这两个点的距离。然后与d取最小值。

3、这个集合有三个或三个以上的点：
(1)所有点都在左集合或者是右集合，那么直接合并即可。

(2)有些点在左集合，有些点在右集合（难点）。这时我们要在S1和S2的中间点p[mid]附近找到左右离它小于d的点。然后对这些点的y从小到大排序，之后暴力法遍历这些存下来的点更新d的值。

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<cmath>
using namespace std;
const double eps=1e-8;
const int MAXN=1e5+10;
const double INF=1e20;

struct point{
    double x,y;//结构体存坐标
};

point p[MAXN];
int a[MAXN];

bool cmpx(const point &a,const point b){
    return a.x<b.x;//按照x从小到大排序
}

bool cmpy(const int &a,const int &b){
    return p[a].y<p[b].y;
}

double dis(int a,int b){
    return sqrt((p[a].x-p[b].x)*(p[a].x-p[b].x)+(p[a].y-p[b].y)*(p[a].y-p[b].y));
    //求两坐标距离(x2-x1)^2+(y2-y1)^2
}


double cloest(int left,int right){
    if(left==right)//只有一个点的时候返回一个很大的值
        return 1000000;

    if(left+1==right)//当只有两个点的时候计算两个点的距离
        return dis(left,right);

    int mid=(left+right)>>1;//取中间值
    double d1=cloest(left,mid);//递归左边部分
    double d2=cloest(mid+1,right);//递归右边部分

    double d=min(d1,d2);//当前记录距离最小取最小值

    int i,j,k=0;
    for(i=left;i<=right;i++){
        if(fabs(p[mid].x-p[i].x)<d)//找出和mid位置点小于d的点的位置
            a[k++]=i;//a数组存的是位置
    }

    sort(a,a+k,cmpy);//按y排序

    for( i = 0; i < k - 1; i++ )
	{
		for( j = i+1; j < i + 7 && j < k; j++ )
		{
			if(p[a[j]].y - p[a[i]].y >= d)
				break;
			d = min(d,dis(a[i],a[j]));//每次都跟新一下最小距离d
		}
	}
    return d;
}




int main(){
    int n;
    while(~scanf("%d",&n)&&n){
        for(int i=0;i<n;i++) scanf("%lf %lf",&p[i].x,&p[i].y);

        sort(p,p+n,cmpx);//按x排序

        printf("%.2f\n",cloest(0,n-1)/2);
    }
    return 0;
}