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;  
}