平面最近点对算法 ClosestPair

1. 问题

• 平面最近点对问题，给定一个平面内所有点的坐标，找出这些点中最近的两个点的距离。

2. 解析

• 朴素暴力算法：

  • 求解平面最近点对问题，需要枚举两个点对，然后记录最小的点对距离，故总时间复杂度为 $O(N^2)$

• 分治法求解最近点对：

  • $Divide$：

    • 将当前的集合问题转化为两个子集合问题，即将当前 $[l,r]$ 的集合分割成 $[l,mid]$ 和 $[mid+1,r]$ 两个集合求解。

    • 需要先对整个集合根据点的横坐标进行排序

  • $Merge$：

    • 我们定义集合 $[l,r]$ 中最近点对的距离为 $D$

    • 那么 $d$ 的归属有三种情况：1. 属于 $[l,mid]$ 集合 2. 属于 $[mid+1,r]$ 集合 3. 两个点分别属于左右两个集合。

    • 我们可以通过递归分治求解当前的1和2情况，得到左右区间的最小点对距离为 $d$ ，所以我们只需要考虑3情况即可。

    • 对于3的情况，由于左右集合的最近点对距离是 $d$。在一维的情形下，我们可以根据中位数的横坐标X，划分出一个 $[X-d,X+d]$ 的集合，并将这个范围内的点纳入我们的考虑集合 $S$ 中，对于范围之外的点对的最近距离只可能大于等于 $d$ ，故直接排除。

    • 接着我们考虑二维的情形，即集合中点的纵坐标的影响，我们可以对集合 $S$ 中的点根据纵坐标升序排序，然后枚举点对然后计算点对的距离 $dis$, 并记录下最近点对距离，即 $D=min(D,dis)$，如果两个点的纵坐标差距大于 $D$ 可以直接结束当前循环，最后记录下来的 $D$ 即是最小点对距离。根据鸽笼原理可以证明某个点最多只需要检查与其相邻的6个点；

3. 设计

平面最近点对算法：

#include<bits/stdc++.h>
using namespace std;
const int N = 1e5+10;

int n,m,k;

struct point {
    double x,y;
}e[N];

int cmpX(point p1,point p2) {
    return p1.x < p2.x;
}

int cmpY(point p1,point p2) {
    return p1.y < p2.y;
}


double dis(point p1,point p2) {
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}

double getMinDistance(point e[],int l,int r) {
    //如果刚好有两个点
    if (r-l==1) {
        return dis(e[l],e[r]);
    }

    //如果刚好有三个点
    if (r-l==2) {
        return min(dis(e[l],e[l+1]),min(dis(e[l],e[r]),dis(e[l+1],e[r])));
    }

    //如果多于三个点需要分治
    //先根据一维的x坐标进行排序和分治，根据中位数分成左右两部分点
    int mid = l+r>>1;
    
    //得到左右集合最小的距离
    double d = min(getMinDistance(e,l,mid),getMinDistance(e,mid+1,r));

    //有可能两个最短点对分别在左右两个集合还需要进一步讨论
    point p[r-l+1];
    int cnt=0;

    //将中位数旁边[x-d,x+d]范围之内的点都容纳进来
    for (int i=l;i<=r;i++) {
        if (fabs(e[i].x-e[mid].x)<d) {
            p[++cnt] = e[i];
        }
    }
    
    //根据y坐标进行升序排序
    sort(p+1,p+cnt+1,cmpY);

    for (int i=1;i<=cnt;i++) {
        for (int j=i+1;j<=cnt;j++) {
            
            //由于是升序排序,如果当前的点对y坐标差距比d还要大
            //说明接下来的点的距离都要大于d,故break
            if (p[j].y - p[i].y >= d) break;

            //否则可能是最小点对,需要更新一下最小点对距离
            d = min(d,dis(p[i],p[j]));
        }
    }
    return d;
}


void run() {
    while (scanf("%d",&n),n) {
        for (int i=1;i<=n;i++) {
            scanf("%lf %lf",&e[i].x,&e[i].y);
        }
        sort(e+1,e+n+1,cmpX);
        printf("The minimum distance of point pair is = %lf",getMinDistance(e,1,n));
    }
}

int main() {
    run();
    return 0;
}

4. 分析

• 平面点对算法需要多次分治和进行子问题合并操作

• 合并的复杂度为 $O(n)$，分治的子问题复杂度为 $O(\frac{n}{2})$

• 故总时间复杂度约为 $O(\frac{n}{2})+O(\frac{n}{2})+O(n)=O(nlogn)$

5. 源码

https://github.com/a894985555/Algorithm/tree/main/%E6%9C%80%E8%BF%91%E7%82%B9%E5%AF%B9