HDU1007 平面最近点对

只判六个点的抽屉原理好像不太懂，，打了个类似暴力的东西但是听说复杂度可证？反正也过去了。。

先将区间分成左右两个子区间，递归的计算这两个子区间内的点对
的最近距离，记这两个最小值中较小的一个为res。

然后计算跨越这两个子区间分界线的点对的最近距离，显然，枚举的
时候不需要枚举坐标之差的绝对值大于res的点对，同理，y坐标之差的绝
对值大于res的点对也不需要进行枚举。

因此，我们在递归的同时用归并排序的方法将区间[l,r]内的点按y坐标
从小到大排序，将x坐标不在[mid−res,mid+ees]范围内的点剔除(mid为x坐

标分界线)，枚举剩下的y坐标差的绝对值不大于res的点对，更新答案。

以上为st的讲课题解

这道题简单代码粗暴写得很快

//QWsin
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const double INF=1e10;
const int maxn=100000+10;

struct point{
    double x,y;int t;
    inline void input(){scanf("%lf%lf",&x,&y);} 
    inline bool operator < (const point &rhs)const{
        return y<rhs.y||(y==rhs.y&&x<rhs.x);    
    }
}a[maxn],t[maxn];

inline int cmp(point a,point b){return a.x<b.x;}

inline double dis(point a,point b)
{
    double dx=a.x-b.x,dy=a.y-b.y;
    return sqrt(dx*dx+dy*dy);
}

int n;

double solve(int l,int r)
{
    if(l==r) return INF;
    if(l+1==r) return dis(a[l],a[r]);

    int mid=(l+r)>>1;
    double d=min(solve(l,mid),solve(mid+1,r));

    int k=0;
    for(int i=l;i<=r;++i) 
        if(fabs(a[i].x-a[mid].x) <= d) t[++k]=a[i];

    sort(t+1,t+k+1);

    for(int i=1;i<=k;++i)
        for(int j=i+1;j<=k&&t[j].y-t[i].y <= d;++j)
            d=min(d,dis(t[i],t[j]));

    return d;
}

inline void solve()
{
    for(int i=1;i<=n;++i) a[i].input();
    sort(a+1,a+n+1,cmp);
    printf("%.2f\n",solve(1,n)/2);
}

int main()
{
    while(scanf("%d",&n)==1&&n) solve();
    return 0;
}