【POJ3714】Raid

题意

在两个点集中，找出距离最近的两个点，求出这个距离

分析

又是学习之后才会做，最小点对。

思想是二分，将点以x坐标排序，然后分为两部分，在这两部分中分别求出最近的距离，然后可以得到一个ans=min(ansl,ansr)。

但是还需要处理中间部分，可知能比ans这个距离还小的，横坐标一定是在(mid-ans,mid+ans)之间，否则光横坐标就已经GG了

还有一个纵坐标的问题。将刚刚筛选出来可能使答案更小的点以y坐标排序，同理，当y坐标的跨度>ans后就GG了

看起来枚举这些x,y时间复杂度很大，然而这个枚举是6次，根据抽屉原理可以证明，我研究了挺久也不会，直面自己的弱......

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cmath>
using  namespace std;
#define N 300100
#define oo 1e15  
struct email
{
    double x,y;
    int c;
}a[N],tmp[N];
double ans; 
int t,n,mid,cnt;
bool cmpx(email a,email b)
{
    return a.x<b.x;
}
bool cmpy(email a,email b)
{
    return a.y<b.y;
}
double d(email a,email b)
{
    if(a.c==b.c)
        return oo;
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

double check(int l,int r)
{
    if(l==r) return oo;
    mid=l+r>>1,cnt=0;
    ans=min(check(l,mid),check(mid+1,r));
    for(int i=l;i<=r;i++)
        if(fabs(a[i].x-a[mid].x)<ans)
             tmp[++cnt]=a[i];
    sort(tmp+1,tmp+1+cnt,cmpy);
    for(int i=1;i<=cnt;i++)
        for(int j=i+1;j<=cnt;j++)
        {
            if(tmp[j].y-tmp[i].y>ans)break;
            ans=min(ans,d(tmp[i],tmp[j]));
        }
    return ans;        
}

int main()
{
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
            scanf("%lf%lf",&a[i].x,&a[i].y),a[i].c=0;
        for(int i=1;i<=n;i++)
            scanf("%lf%lf",&a[i+n].x,&a[i+n].y),a[i+n].c=1;
        sort(a+1,a+n*2+1,cmpx);
        printf("%.3lf
",check(1,2*n));
    }
    return 0;
}