HDU 1007 Quoit Design

http://acm.hdu.edu.cn/showproblem.php?pid=1007

求最近两个点的距离
《算法导论》p591的分治算法 O(nlogn)
再改造下《算法竞赛入门经典》p143 归并排序的代码

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

const double eps=1e-8;
int dcmp(double x)
{
    return x<-eps?-1:x>eps;
}
struct point
{
    double x,y;
    void input()
    {
        scanf("%lf%lf",&x,&y);
    }
    bool operator < (const point &p)const
    {
        return y==p.y?x<p.x:y<p.y;
    }
};
bool cmpx(const point &a,const point b)
{
    return a.x==b.x?a.y<b.y:a.x<b.x;
}
double dis(const point &a,const point &b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

const int N=100010;
const double INF=1e30;
point pp[N],qq[N];
double mergesort(point *a,int l,int r,point *b)
{
    if(r-l<=1) return INF;
    int m=(l+r)/2;
    double x=a[m].x;
    int p=l, q=m, i=l;
    double t1=mergesort(a,l,m,b);
    double t2=mergesort(a,m,r,b);
    double t=t1<t2?t1:t2,d=t;
    while(p<m || q<r)
    {
        if(q>=r || (p<m && a[p]<a[q])) b[i++]=a[p++];
        else b[i++]=a[q++];
    }
    for(int i=l;i<r;i++) a[i]=b[i];
    for(int i=l;i<r;i++) if(dcmp(fabs(a[i].x-x)-d)<0)
    {
        int cnt=0;
        for(int j=i+1;j<r && cnt<7;j++)
            if(dcmp(fabs(a[j].x-x)-d)<0)
            {
                double c=dis(a[i],a[j]);
                if(c<t) t=c;
                cnt++;
            }
    }
    return t;
}
int main()
{
    int n;
    while (scanf("%d",&n),n)
    {
        for(int i=0;i<n;i++) pp[i].input();
        sort(pp,pp+n,cmpx);
        double ans=mergesort(pp,0,n,qq);
        printf("%.2lf\n",ans/2);
    }
    return 0;
}