AcWing 119袭击【平面点间对问题】

思路：
把所有点按横坐标排序；
开始分治（l,r）;
递归到（L,R）上只有一个点，或者有俩个点；
找到中间点的横坐标，找到递归时最小间间距res；
找到横坐标属于min_x-res 到 min_x+res 的点 放到Temp数组中；
将temp数组按找y值排序；
用暴力找到既属于以边长为res的正方形内 的俩个点，取最小值

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
using namespace std;

//把所有点按横坐标排序，分治 ，，把属于范围的点找出来,按纵坐标排序
int n;
const int INF=1e9+1999;
const int N=100100;
struct node{
    double x,y;
    int type;
}Point[N],Temp[N];
int cmp1(node Px,node Py){
    return Px.x<Py.x;
}
int cmp2(node Px,node Py){
    return Px.y<Py.y;
}

double dis(int a,int b){
    if(Point[a].type==Point[b].type) return INF;
    double Pax=Point[a].x;
    double Pay=Point[a].y;
    double Pbx=Point[b].x;
    double Pby=Point[b].y;
    return sqrt((Pax-Pbx)*(Pax-Pbx)+(Pay-Pby)*(Pay-Pby));}
double dfs(int l,int r){
    if(l==r) return INF;
    if(r-1==l) return dis(l,r);
    int mid=(l+r)>>1;
    double res=min(dfs(l,mid),dfs(mid+1,r));
    double mid_x=Point[mid].x;
    int cnt=0;
    for(int i=l;i<=r;i++){
        if(Point[i].x>=mid_x-res&& Point[i].x<=mid_x+res){
            Temp[++cnt]=Point[i];
        }
    }
    sort(Temp+1,Temp+1+cnt,cmp2);
    for(int i=1;i<cnt;i++){
        for(int j=i+1;j<cnt;j++){
            if(Point[j].y>=Point[i].y+res) break;
            res=min(res,dis(i,j));
        }
    }
    return res;
}
int main(){
    int T;
    cin>>T;
    while(T--){
        cin>>n;
        for(int i=1;i<=n;i++){
            cin>>Point[i].x>>Point[i].y;
            Point[i].type=1;
        }
        for(int i=n+1;i<=2*n;i++){
            cin>>Point[i].x>>Point[i].y;
            Point[i].type=0;
        }
        n<<=1;
        sort(Point+1,Point+1+n,cmp1);
        printf("%.3f\n",dfs(1,n));
    }
    return 0;
}