[Luogu1429]平面最近点对（加强版）

题目大意：
　　平面最近点对。

思路：
　　分治。
　　首先将所有点排序
　　每次把当前区间分为两半，递归求解两个区间内部的情况，然后枚举区间两边的点。

#include<cmath>
#include<cstdio>
#include<cctype>
#include<algorithm>
inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
typedef std::pair<double,double> Point;
const int N=200000;
Point p[N];
int num[N];
double ans=1e10;
inline double sqr(const double &x) {
    return x*x;
}
inline double dis(const Point &a,const Point &b) {
    return sqrt(sqr(std::abs(a.first-b.first))+sqr(std::abs(a.second-b.second)));
}
void solve(const int &l,const int &r) {
    if(l==r) return;
    const int mid=(l+r)/2;
    solve(l,mid);1`6
    solve(mid+1,r);
    num[0]=0;
    for(register int i=mid+1;i<=r;i++) {
        if(p[i].first-p[mid].first<ans) {
            num[++num[0]]=i;
        }
    }
    for(register int i=l;i<=mid;i++) {
        if(p[mid].first-p[i].first<ans) {
            for(register int j=1;j<=num[0];j++) {
                ans=std::min(ans,dis(p[i],p[num[j]]));
            }
        }
    }
}
int main() {
    const int n=getint();
    for(register int i=0;i<n;i++) {
        scanf("%lf%lf",&p[i].first,&p[i].second);
    }
    std::sort(&p[0],&p[n]);
    solve(0,n-1);
    printf("%.4f\n",ans);
    return 0;
}