首先感谢伯乐的点评,已修改~~
问题描述:
给定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;
}