求最短距离

使用分治法求点之间的最短距离

伪代码：只有一个点时，输出0
			两个点的时候，直接输出距离
			三个点的时候，分成左右两个子集，继续进行递归。对于递归的结果，先求出左右半边的两个点的距离d，然后进行一次查找，找到小于该距离的点。把这些点按照y排序，可以得到一个序列，然后将其中小于d的点赋值给d，直至找到最小的距离。
		
（2）代码：
#include<iostream>
#include <cstdio>  
#include <cstdlib>  
#include <cstring>  
#include <algorithm>  
#include <cmath>  
using namespace std;  
struct Point  
{  
    double x,y;  
}p[100],temp[100];  
const double eps = 1e-8;  
const int INF = 0x7fffffff;  
int n;    
bool cmpy(Point a, Point b)//比较y直的大小
{  
    return a.y < b.y;  
}  
bool cmpx(Point a,Point b)//先把n个点从左到右排序一波 
{   
    if(a.x!=b.x)  
        return a.x<b.x;  
    return a.y<b.y;  
}  
double Dis(Point a, Point b)//求距离函数 
{  
    return sqrt((a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y));  
}  
  
double Closest_Pair(int left, int right)  
{  
    double d = INF;  
    if(left == right)//一个点的情况 
        return d;  
    if(left +1 == right) //判断到2个点之间的时候，直接判断距离值  
        return Dis(p[left],p[right]);  
    int mid = (left+right)>>1;    
    double d1 = Closest_Pair(left,mid);// 利用分而治之求左半段点的的最小距离
    double d2 = Closest_Pair(mid,right);// 求右半段点的最小距离
    d = min(d1,d2);//d为2点均为左右两边的点的距离最小值  
    int k = 0;  
    for(int i = left; i <= right; i++)// temp存的是这个区域与中线距离小于d的点（因为已经有距离为d的点了）  
    {  
        if(fabs(p[mid].x - p[i].x) <= d)  
            temp[k++] = p[i];  
    }  
    sort(temp,temp+k,cmpy); //调用sort函数把这些点按照y排序  
    for(int i = 0; i < k; i++)  
    {  
        for(int j = i+1; j < k && temp[j].y - temp[i].y < d; j++)  
        {  
            double d3 = Dis(temp[i],temp[j]);  
            d = min(d,d3);  
        }  
    }  
    return d;  
}  
int main()  
{   
    cin>>n;  
    for(int i=0; i<n; i++)  
    {  
        double a,b;   
		cin>>a>>b;
        p[i].x=a;  
        p[i].y=b;  
    }  
    sort(p,p+n,cmpx);
	cout<<Closest_Pair(0,n-1)<<endl; 
}