学习笔记#平面最近点对

引入

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

思考

嗯，n方枚举求解，好像海星，但是如果n有20w呢
这时候我们需要用nlogn的算法来求解

分治法

用到分治的思想，将点按照x坐标排序，画一幅图

我们以中间的点为mid来二分查找，将左边的点分为一组，右边的点分为一组，求解左右两组内部点与点之间的最小距离
设左右两组内部最小距离为u，那么我们只需要找到一个区间[l,r]，使这个区间里的点距离mid点不超过u，如图

为什么呢？因为如果我们要更新答案，一定是左右区间的点来更新，而距离mid点超过u，一定不会更新答案
这时候暴力枚举就好了，经过计算之后，可以推出，每一次暴力枚举最多只需要计算6对点对直接的距离

ac代码

#include<bits/stdc++.h>
using namespace std;
struct p{double x,y;}a[200010];
int n;
bool cmp(p i,p j){return i.x<j.x;}
double f(int i,int j){return sqrt((a[i].x-a[j].x)*(a[i].x-a[j].x)+(a[i].y-a[j].y)*(a[i].y-a[j].y));}
double erfen(int l,int r)
{
    if(r-l==1)return f(l,r);
    if(r-l==2)return min(min(f(l,l+1),f(l+1,l+2)),f(l,l+2));
    int mid=(l+r)>>1;
    double u=min(erfen(l,mid),erfen(mid+1,r));
    if(u==0)return 0;
    int ll=mid,rr=mid;
    while(ll>l&&a[ll-1].x+u>a[mid].x)ll--;
    while(rr<r&&a[rr+1].x-u<a[mid].x)rr++;
    for(int i=ll;i<=mid;i++)for(int j=mid+1;j<=rr;j++)u=min(u,f(i,j));
    return u;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%lf%lf",&a[i].x,&a[i].y);
    sort(a+1,a+1+n,cmp),printf("%.4lf\n",erfen(1,n));
    return 0;
}