分治——Luogu1429 平面最近点对（加强版）

题面：Luogu1429
关于求平面最近点对的方法有很多很多。。。
像是那种随机旋转角度然后暴力算的这种方法就很鬼畜QAQ
靠谱一点的算法呢就是分治了
我们计算区间l~r的最近距离
首先计算l~mid和mid~r的最近距离，然后考虑一下mid左边和mid右边对答案的贡献
首先记录下dist=min（l~mid答案，mid~r答案），然后把左边x坐标与mid的差小于等于dist的点统计出来，再把右边x坐标与mid的差小于等于dist的点统计出来
然后左边右边暴力判断统计答案就行了

#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#include<iostream>#include<cstdlib>#include<set>#include<map>#include<string>#include<ctime>#include<queue>#include<vector>using namespace std;inline double sqr(double x){return x*x;}struct ppap{double x,y;}a[200001];ppap p1[200001],p2[200001];inline bool cmp(ppap a,ppap b){return a.x==b.x?a.y<b.y:a.x<b.x;}int n,pp1,pp2;inline double work(int l,int r){    double ans=1e9;    if(l==r)return ans;    if(l+1==r)return sqr(a[l].x-a[r].x)+sqr(a[l].y-a[r].y);    int mid=l+r>>1;    ans=min(work(l,mid),work(mid,r));    pp1=pp2=0;    for(int i=l;i<mid;i++)if(sqr(a[mid].x-a[i].x)<=ans)p1[++pp1]=a[i];    for(int i=mid+1;i<=r;i++)if(sqr(a[i].x-a[mid].x)<=ans)p2[++pp2]=a[i];    for(int i=1;i<=pp1;i++)        for(int j=1;j<=pp2;j++){            double d=sqr(p1[i].x-p2[j].x)+sqr(p1[i].y-p2[j].y);            ans=min(ans,d);        }    return ans;}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+n+1,cmp);    printf("%.4lf",sqrt(work(1,n)));    return 0;}