最近点对

问题描述

给定一组随机点，求任意两点间距离的最小值。

暴力方法

将任意两点间距离全部算出来，然后找到其中最小值。
但显然其效率为O(n^2)。

一般方法

使用分治的思想，将点划分成更小的集合来进行处理。

先从一维的点开始考虑

在x轴上分布着n个点，对点的下标进行二分。
其中出现的问题： 最小点对恰好跨过了二分的界限，必须要由当前程序完成对此区域的搜索。
1. 左递归，找左半区间最近点对dl 。
2. 右递归，找右半区间最近点对dr 。
3. 得到d=min(dr，dl)；
4. 再以分割点左右各扩展d，找扩展区间最近点对，并返回。

二维中最近点对

同上思想，只不过是把划分的线变成了划分为一块区域（中位点左右各一块）
1. 以x为序排序
2. 递归二分
1 . 以区间中位数（以x为准）来分割区间
2. 左右递归，分别求各区间最近点对dr，dl，并由此得到d。
3. 再以d左右扩展分割线，在此区间找最近点对，来更新d。
4. 对此区间的点另存一数组，以y轴为准排序，排序后相邻两点的y值之差大于d无需计算。

附上代码

#include"cstdio"
#include"iostream"
#include"cmath"
#include"cstdlib"
using namespace std;

struct point
{
    int x,y;
};

point a[100001],b[100001];

double dist(point sp,point op)
{
    return sqrt((sp.x-op.x)*(sp.x-op.x)+(sp.y-op.y)*(sp.y-op.y));
}

int cmp1(const void *a,const void *b)
{
    point *p=(point *)a;
    point *q=(point *)b;
    if(p->x > q->x )return 1;
    return -1; 
}

int cmp2(const void *a,const void *b)
{
    point *p=(point *)a;
    point *q=(point *)b;
    if(p->y > q->y )return 1;
    return -1; 
}

double erfen(int m,int n)
{
    if(n-m<1)return 1<<29;
    if(n-m==1)return dist(a[m],a[n]);
    int z=(n-m)/2+m;
    double dl=erfen(m,z),dr=erfen(z,n);
    double d=min(dr,dl);
    int i,j,l=0;
    for(i=z;i>=m;i--)
    {
        if(a[z].x-a[i].x>d)break;
        b[l++]=a[i];
    }
    for(i=z+1;i<=n;i++)
    {
        if(a[i].x-a[z].x>d)break;
        b[l++]=a[i];
    }
    l--;
    qsort(b,l,sizeof(point),cmp2);

/*  cout<<endl;
    for(i=0;i<l;i++)
        printf("%d %d\n",b[i].x,b[i].y);
    cout<<endl;    
*/  
    for(i=0;i<l;i++)
        for(j=i+1;j<l;j++)
        {
            if(a[j].y-a[i].y>d)break;
            double distence=dist(b[i],b[j]);
            if(distence<d)d=distence;
        }
// printf("~~~~~~~~~~~~~~~~~~~~%lf~~~~~~~~~~~~~~~~~~~~\n",d);
    return d;
}

int main()
{
    freopen("nearest.in","r",stdin);
    freopen("nearest.out","w",stdout);
    int n,i;
    scanf("%d",&n);
    for(i=0;i<n;i++)
        scanf("%d %d",&a[i].x,&a[i].y);
    qsort(a,n,sizeof(point),cmp1);
//  for(i=0;i<n;i++)
//      printf("%d %d\n",a[i].x,a[i].y);
    printf("%.2lf",erfen(0,n-1));
    return 0;
}