平面最近点对问题

平面最近点对问题正如其名，给定平面上的$n$个点，找出其中的一对点，使得这对点的距离在所有点对中最小。

 

首先显而易见地我们可以得到这个问题的$O(n^2)$算法，枚举所有点对即可。但是很显然我们可以注意到，这里面有很多点对显然不是最优的，那么我们可以想到一种剪枝方法，就是将只对x坐标差值小于当前已知最小值的点对进行判断（否则必然不是最优解），从而减少判断量。

 

我们考虑使用分治来实现这种剪枝，先将平面上的点分为两部分，分治求出两部分内部的最近点对距离。之后我们要做的就是枚举两个集合之间的点对，并与两部分内部的最近点对距离比较来得到最近点对距离。这里我们是不需要枚举所有点对的，因为我们已经得到了一个两部分各自内部最小的点对距离，因而我们可以结合上面的根据$x$坐标的剪枝方法，只枚举分别属于两部分的$x$坐标差小于已知最小距离的点对。

 

这样做的复杂度近似于$O(n\log^2n)$，至于怎么得到的……我也不知道。_(:зゝ∠)_

 

例题：

1. Vijos 1012 清帝之惑之雍正

链接：https://vijos.org/p/1011

2. 平面最近点对（加强版）

链接：https://www.luogu.org/problem/show?pid=1429

 

另外附上模板：

注意，本模板保留六位小数，不能直接用于提交上面的例题，若要提交请修改输出精度。

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <string>
#include <sstream>
#include <cctype>
#include <cmath>
#include <algorithm>
#define THE_BEST_PONY "Rainbow Dash"

using namespace std;
const int maxn=220000,INF=~0u>>1;

struct Point{
    double x,y;
    bool operator < (const Point &a) const {
        if(x<a.x) return true;
        if(x>a.x) return false;
        return y<a.y;
    }
}p[maxn];

int n;
double ans;

double DisP(int a,int b){
    return sqrt((p[a].x-p[b].x)*(p[a].x-p[b].x)+(p[a].y-p[b].y)*(p[a].y-p[b].y)); 
}

double GetAns(int l,int r){
    int mid=(l+r)>>1;
    if(l==r) return INF;
    if(l==r-1) return DisP(l,r);
    double len=min(GetAns(l,mid),GetAns(mid+1,r));
    for(int i=mid;i>=l;i--){
        if(p[i].x+len<p[mid].x) break;
        for(int j=mid+1;j<=r;j++){
            if(p[mid].x+len<p[j].x) break;
            len=min(len,DisP(i,j));
        }
    }
    return len;
}

int main(){
    scanf("%d",&n);
    for(int i=0;i<n;i++)
        scanf("%lf%lf",&p[i].x,&p[i].y);
    sort(p,p+n);
    ans=GetAns(0,n-1);
    printf("%.6lf",ans);
    return 0;
}