poj 3714 Raid 分治

先按横坐标排序，然后分治，合并的时候只需要处理两边x坐标差在d范围内的点即可。把这些点找出来之后，再按照y值排序，然后再检查y坐标差在d以内的点，那么可以保证每个点最多检查6个点(否则这些点之间必定存在距离小于d的点对)，那么可以保证o(n)的合并。总的复杂度nlognlogn。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn=1e5+9;
const double inf=1e11;
struct D
{
    int x,y,id;
    bool operator <(const D & xx) const
    {
        return x<xx.x;
    }
}data[maxn<<1],now[maxn<<1];

struct cmp
{
    bool operator()(D a,D b)
    {
        return a.y<b.y;
    }
};

double cal(int t,int s)
{
    if(now[t].id==now[s].id) return inf;
    double a=now[t].x-now[s].x;
    double b=now[t].y-now[s].y;
    a*=a;
    b*=b;
    return sqrt(a+b);
}

double solve(int l,int r)
{
    if(l==r) return inf;
    int mid=l+r>>1;
    double ret=solve(l,mid);
    ret=min(ret,solve(mid+1,r));
    int lon=0;
    for(int i=mid;i>=l;i--)
    {
        if(data[mid+1].x-data[i].x>ret) break;
        now[++lon]=data[i];
    }
    for(int j=mid+1;j<=r;j++)
    {
        if(data[j].x-data[mid].x>ret) break;
        now[++lon]=data[j];
    }
    sort(now+1,now+1+lon);
    for(int i=1;i<=lon;i++)
    for(int j=i+1;j<=lon;j++)
    {
        if(now[j].y-now[i].y>ret) break;
        ret=min(ret,cal(i,j));
    }
    return ret;
}

int main()
{
    int T,n;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(int i=1;i<=n+n;i++)
        {
            scanf("%d %d",&data[i].x,&data[i].y);
            data[i].id=i>n;
        }
        sort(data+1,data+1+n+n);
        double ans=solve(1,n+n);
        printf("%.3f\n",ans);
    }
    return 0;
}