平面上的最接近点对

平面上最接近点对

一维最接近点对问题：

#include <cstdio>
#include <algorithm>
using namespace std;
const int inf=0x3f3f3f3f;
double s[100];
double mn;

int main(){
    int n;
    while(~scanf("%d",&n)){
        for(int i=0;i<n;i++){
            scanf("%lf",s+i);
        }
        sort(s,s+n);
        mn=inf;
        for(int i=0;i<n-1;i++){
            mn=s[i+1]-s[i]>mn?mn:s[i+1]-s[i];
        }
        printf("%lf\n",mn);
    }
}

使用分治求解：

S中的n个点为x轴上的n个实数x1，x2，…，xn。最接近点对即为这n个实数中相差最小的两个实数。显然可以先将点排好序，然后线性扫描就可以了（上述程序实现）。但我们为了便于推广到二维的情形，为下面二维，尝试用分治法解决这个问题。

假设我们用m点将S分为S1和S2两个集合，这样一来，对于所有的p(S1中的点)和q(S2中的点)，有p<q。

递归地在S1和S2上找出其最接近点对{p1,p2}和{q1,q2}，并设

　　　　　　　　　　　　　　　　　　d = min{ |p1-p2| , |q1-q2| } 　　　　　　　　　　　　　　　　　　

由此易知，S中最接近点对或者是{p1,p2}，或者是{q1,q2}，或者是某个{q3,p3}，如下图所示。

    如果最接近点对是{q3,p3}，即|p3-q3|<d，则p3和q3两者与m的距离都不超过d，且在区间(m-d,d]和(d,m+d]各有且仅有一个点。这样，就可以在线性时间内实现合并。

    此时，一维情形下的最近点对时间复杂度为O(nlogn)。

#include <cstdio>
#include <cstdlib>
#include <algorithm>
using namespace std;
const int inf=0x3f3f3f3f;
double s[100];
double mn;
double closest(int low,int high){
    if(low+1==high) return s[high]-s[low];
    if(low+2==high) return min(s[low+1]-s[low],s[high]-s[low+1]);
    int mid=(low+high)>>1;
    double ans=min(closest(low,mid),closest(mid+1,high));
    if(s[mid+1]-s[mid]<ans) ans=s[mid+1]-s[mid];
    return ans;
}

double Random()
{
    double result=rand()%10000;
     return result*0.01;
}
int main(){
    int n,m;
    //随机生成种子值
    srand(1);
    while(~scanf("%d",&n)){

        for(int i=0;i<n;i++){
            s[i]=rand()%10000;
            printf("%lf ",s[i]);
        }
        printf("\n");
        sort(s,s+n);
        for(int i=0;i<n;i++)
            printf("%lf ",s[i]);
        printf("\n");
        printf("%.4lf\n",closest(0,n-1));
    }
}

二维最接近点对问题：

分治：

1) 把它分成两个或多个更小的问题；

 2) 分别解决每个小问题； 

3) 把各小问题的解答组合起来，即可得到原问题的解答。小问题通常与原问题相似，可以递归地使用分而治之策略来解决。

其实，这里 用到了分治的思想。将所给平面上n个点的集合S分成两个子集S1和S2，每个子集中约有n/2个点。然后在每个子集中递归地求最接近的点对。在这里，一个 关键的问题是如何实现分治法中的合并步骤，即由S1和S2的最接近点对，如何求得原集合S中的最接近点对。如果这两个点分别在S1和S2中，问题就变得复 杂了。

 

在二维的情况下:

我们仿照一维的情况先把所有点按照x(横坐标)从左到右升序排列.

以X横坐标中间的点作为分界线.将平面的点分成左边和右边,以上图为例,分为左边5个右边5个.

然后找到左边的点中最近点对的距离d1,和右边最近点对的距离d2。

令d=min{d1,d2}.那么d是当前整个平面的最近点对的距离吗?

答案不是!!

当前的d1,d2都是两边各自的点的最近距离,但是没有考虑到两边的点相互配对的情况,即是点对一个在左边区域,一个在右边区域.

如果我们这个时候把两边相互配对的所有情况全部罗列出来那么时间复杂度就变为第一种方法的O(n²).

这个时候我们便可以用上面找到的d来限制配对.即是明显超过d的两点不需要配对.如下：

那么在范围内的只有三个点,也就是说只有这个三个点相互配对的距离才可能比d小.那么再按照方法一一样在这些点中找到最短距离再跟d去比较.小的就是当前整个平面中所考虑的所有点中最近点对的距离。

然而在以上问题处理上,有一个问题尚未解决就是如何找到两边区域中的最近点对的距离d1,d2 ?

我们可以发现在处理上面这个问题的时候，和最开始处理所有平面的点最近点距离的问题类似，只是点的数目减半了.

那么我们可以递归以上问题.直到划分的区域中只有二个或者三个点.这样,该区域的最近点对距离就是对应的距离。这也是递归终止的条件。

 

 

平面上的最接近点对

题目描述

给定平面上n个点，找出其中的一对点的距离，使得在这n个点的所有点对中，该距离为所有点对中最小的。

输入

第一行：n；2≤n≤60000
接下来n行：每行两个实数：x y，表示一个点的行坐标和列坐标，中间用一个空格隔开。

输出

仅一行，一个实数，表示最短距离，精确到小数点后面4位。

样例输入

3
1 1
1 2
2 2

样例输出

1.0000

平面上的最短距离。

#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;

平面上的最短距离。
struct point{
    double x,y;
} p[60005];
int a[60000];

int cmpx( point a,point b){
    //if(a.x==b.x)return a.y<b.y;
    return a.x<b.x;
}

int cmpy(int a,int b){
    return p[a].y<p[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(int low,int high){
    if(low+1==high) return dis(p[low],p[high]);
    if(low+2==high)return min(dis(p[low],p[high]),min(dis(p[low],p[low+1]),dis(p[low+1],p

[high])));
    int mid = (low+high)>>1;
    double ans= min(closest(low,mid),closest(mid+1,high));

   // printf("%d %d %lf\n",low,high,ans);
    int i,j,cnt=0;
    for(i=low;i<=high;i++){//最小值对在左边有一个在右边有一个情况，找出在[mid-ans,mid]和

[mid,mid+ans]的所有值
        if(p[i].x>=p[mid].x-ans&&p[i].x<=p[mid].x+ans)  a[cnt++]=i;
    }
    //printf("%d \n",cnt);
    sort(a,a+cnt,cmpy);
    for(int i=0;i<cnt;i++){ //在[mid-ans,mid]和[mid,mid+ans]区间找到比ans小的值。
        for(j=i+1;j<cnt;j++){
            if(p[a[j]].y-p[a[i]].y>=ans)break;
            ans=min(ans,dis(p[a[i]],p[a[j]]));
        }
    }
    //printf("%lf\n",ans);
    return ans;
}

int main()
{
    freopen("data.in","r",stdin);
    int i,n;
    while(~scanf("%d",&n)){
        for(i=0;i<n;i++){
            scanf("%lf %lf",&p[i].x,&p[i].y);
        }
        sort(p,p+n,cmpx);
        for(int i=0;i<n;i++)printf("%lf %lf\n",p[i].x,p[i].y);
        printf("%.4lf\n",closest(0,n-1));
    }
    return 0;
}

 

参考：http://blog.csdn.net/hackbuteer1/article/details/7482232
　　　http://www.mamicode.com/info-detail-875129.html
　　　http://my.oschina.net/bgbfbsdchenzheng/blog/498823