P1257 平面上的最接近点对 题解

没被 hack 的一种做法。

将点集按 $x$ 坐标排序后递归划分，最近点对要么在左半、要么在右半、要么跨中线。

**主要步骤**：

1. 递归求出左右半的最近距离 $d = \min(d_l, d_r)$。
2. 跨中线情况只需考虑 $|x_i - midX| < d$ 的点
3. 将这些点按 $y$ 排序后，每个点只需与后面 $y$ 差 $< d$ 的常数个点比较（理论保证 $\le 6$ 个）

**时间复杂度**：

$T(n) = 2T(n/2) + O(n \log n) = O(n \log^2 n)$。

若合并时用归并按 $y$ 排序，可优化至 $O(n \log n)$。

```cpp
#include <bits/stdc++.h>
using namespace std;
const int N=10005;
const long double INF=1e20L;
struct P
{
    long double x,y;
}p[N],t[N];
int n;
bool cmpx(const P& a,const P& b)
{
    return a.x<b.x;
}
bool cmpy(const P& a,const P& b)
{
    return a.y<b.y;
}
long double d(const P& a,const P& b)
{
    long double dx=a.x-b.x,dy=a.y-b.y;
    return sqrtl(dx*dx+dy*dy);
}
long double f(int l,int r)
{
    if(l==r) return INF;
    if(l+1==r) return d(p[l],p[r]);
    int m=(l+r)/2;
    long double mx=p[m].x;
    long double ans=min(f(l,m),f(m+1,r));
    int k=0;
    for(int i=l;i<=r;i++)
    {
        if(fabsl(p[i].x-mx)<ans) t[k++]=p[i];
    }
    sort(t,t+k,cmpy);
    for(int i=0;i<k;i++)
    {
        for(int j=i+1;j<k&&t[j].y-t[i].y<ans;j++)
        {
            ans=min(ans,d(t[i],t[j]));
        }
    }
    return ans;
}
int main()
{
    scanf("%d",&n);
    for(int i=0;i<n;i++)
    {
        double x,y;
        scanf("%lf%lf",&x,&y);
        p[i].x=x;p[i].y=y;
    }
    sort(p,p+n,cmpx);
    printf("%.4Lf\n",f(0,n-1));
    return 0;
}
```