hdoj 1007 Quoit Design(分治法)

題意:給你n個點，讓你求出n個點裏距離最近的兩點距離。

思路:分治法先把點按照座標x進行排序，劃分爲左右兩部分，

最近點對有三種情況，都在左半部分，都在右半部分，兩點分別在左邊和右邊

遞歸求出都在左邊和都在右邊的情況，選一個最小值m1。

然後求出兩邊都有的情況。分別在兩邊的情況可以對點進行進一步的篩選，將點的下標

保存在t[]，對篩選出的點按照y座標進行排序，方便求出最小值。

代碼如下:

#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<math.h>
using namespace std;
struct point
{
	double x,y;
};
point p[100005];
int t[100005]; 
double min(double a,double b)
{
	return (a>b)?b:a;
}
double dis(point a,point b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
bool cmp(const point &a,const point &b)
{
	return a.x<b.x;
}
bool cmp1(int a,int b)
{
	return p[a].y<p[b].y;
}
double spit(int a,int b)
{
	int i,j,c;
	double m1;
	if(a+1==b)
	return dis(p[a],p[b]);
	if(a+2==b)
	return min(dis(p[a],p[a+1]),min(dis(p[a],p[b]),dis(p[b],p[a+1])));
    c=(a+b)>>1;  
    m1=min(spit(a,c),spit(c+1,b));  
	int k=0;
	for(i=a;i<=b;i++) 
	{
		if(p[i].x>=p[c].x-m1&&p[i].x<=p[c].x+m1)
        t[k++]=i;
	} 
	sort(t,t+k,cmp1);
    for(i=0;i<k;i++)  
        for(j=i+1;j<k;j++)
		{                       
            if(p[t[j]].y-p[t[i]].y>m1) break;  
            m1=min(m1,dis(p[t[i]],p[t[j]]));
        }  
    return m1; 
}
int main()
{  
    int n;  
    int i;  
    double r;  
    while(scanf("%d",&n)){  
        if(n==0) break;  
        for(i=0;i<n;i++) scanf("%lf %lf",&p[i].x,&p[i].y);  
        sort(p,p+n,cmp);  
        r=spit(0,n-1);  
        r=r/2;  
        printf("%.2lf\n",r);  
    }  
    return 0;  
}