最近点对问题（分治）

首先感谢伯乐的点评，已修改~~

问题描述：

给定n个点的横纵坐标，求它们的最近距离。

例

2

1 0

0 1

输出1.4142

问题分析：

若是直接暴力求取，复杂度为O（n^2)，从n个点选出2个点,n*(n-1)/2嘛

而如果用分治法，复杂度可以降到O(nlogn)，方法是先根据横纵坐标从小到大排序

然后寻找left~mid（点全部在左边）,mid+1~right（点全部在右边）的最小点对，再将它们两个中的最小值和一个点在左边，一个点在右边的情况比较（有点像分治法求最大子序列和），先不管中间的处理，递归到什么情况为止呢，如果left==right则说明只有一个点，返回无穷大，如果left+1==right说明有两个点，直接返回距离

再回过头来考虑中间的情况（一个点在左边，一个点在右边），考虑到之前已经找到最短距离为d,所以中间的点对已经是在距离2d*2d的正方形内（这点很重要！）

严格来说应该是圆内不过圆计算太麻烦且有精度问题。

找出这些点，将他们一一组合可以得到他们的最短距离，如果比d小就更新d.

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cmath>
 using namespace std;
 int n;
 struct node
 {
     double x;
     double y;
 };
int cmp(node a,node b)//两个点按照横坐标从小到大排序，相同则按照纵坐标从小到大
{
     if (a.x==b.x)
     return a.y<b.y;
     return a.x<b.x;
}
node point[10005];
int cmp2(int a,int  b)
 {
     return point[a].y<point[b].y;
 }

 double distance(double c1,double c2,double d1,double d2)
 {
    return sqrt((c1-d1)*(c1-d1)+(c2-d2)*(c2-d2));
 }
 double calculate(int left,int right)//left是最左边的坐标，right是最右边的坐标
 {
     int i,j,mid,tempnumber;
     double d,d1,d2;
     int temp[10005];
     if (left==right)//如果是一个点，那么距离自然是0
        return 999999;//返回一个无穷大
     if (left+1==right)//如果是两个点，直接返回两点的距离
        return  distance(point[left].x,point[left].y,point[right].x,point[right].y);
     mid=(left+right)>>1;
     d1=calculate(left,mid);//左部分的最近点对的距离
     d2=calculate(mid+1,right);//右部分的最近点对的距离
     d=d1<d2?d1:d2;//取较小的，接下来要考虑一个点在左部分，一个点在右部分的情况
     i=mid;
     tempnumber=0;
     while(i>=0&&point[mid].x-point[i].x<d)//从中点往左边找
     {
        temp[tempnumber++]=i;
        i--;
     }
     i=mid+1;
     while(i<n&&point[i].x-point[mid].x<d)//从中点往右边找（注意这次不含mid是从mid+1,避免重复mid)
     {
         temp[tempnumber++]=i;
         i++;
     }
     sort(temp,temp+tempnumber,cmp2);//将得到的temp按照纵坐标从小到大排序，这样一旦某一个j的纵坐标大于i，那么接下来的j也不用考虑（其实个人感觉效率差距不大，毕竟你还要排序，不过网上都这样就随波逐流了...）
     for (i=0;i<tempnumber;i++)
        for (j=i+1;j<tempnumber;j++)
     {
         if (point[temp[j]].y-point[temp[i]].y>d)
            break;//如果已经大于d直接跳出循环
         else
         {
           d1=distance(point[temp[i]].x,point[temp[i]].y,point[temp[j]].x,point[temp[j]].y);//否则看范围内是否还有小于d的点对
           d=d<d1?d:d1;  //更新
         }
     }
     return d;
 }
int main()
{
     int i,j;
     scanf("%d",&n);
     for (i=0;i<n;i++)
     scanf("%lf%lf",&point[i].x,&point[i].y);
     sort(point,point+n,cmp);
     printf("%.4f\n",calculate(0,n-1));
     return 0;
}