hdu 1007 Quoit Design(最近点对)

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

一看题目就是一个很裸的最近点对。。暴力的做法是o(n^2)而对于o(nlogn)算法知识理论上理解了，但是还没实现过，好不容西写出来还是tle最后在gsb的指点下总算是A了。。汗。。

#include <iostream>

#include <cstdio>

#include <algorithm>

#include <cmath>

const int max_s = 100007;

using namespace std;

const int M = 99999999.0;

struct node

{

    double x,y;

}p[max_s];

int a[max_s];

int cmp_x(node a,node b)

{

  return a.x>b.x;

}

int cmp_y(int s,int e)

{

    return p[s].y>p[e].y;

}

double dis(node t1,node t2)

{

    double x=t1.x-t2.x;

    double y=t1.y-t2.y;

    return sqrt(x*x+y*y);

}

double short_s(int l,int r)//分治的思想

{

    int i,j;

    if(l==r) return M;

    else if(l+1==r) return dis(p[l],p[r]);

    int m=(l+r)>>1;

    double Min=min(short_s(l,m),short_s(m+1,r));

    int k=0;

    for(i=l;i<=r;i++)//这里是从L开始的。。

    {

        if(p[i].x>=p[m].x-Min&&p[i].x<=p[m].x+Min)

        a[k++]=i;

    }

    sort(a,a+k,cmp_y);

    for(i=0;i<k-1;i++)

    {

        for(j=i+1;j<k;j++)

        {

            if(p[a[i]].y-p[a[j]].y>=Min)

            break;

            Min=min(Min,dis(p[a[j]],p[a[i]]));

        }

    }

    return Min;

}

int main()

{

    //freopen("d.txt","r",stdin);

    int n,i;

    while(scanf("%d",&n),n)

    {

        for(i=0;i<n;i++)

        scanf("%lf%lf",&p[i].x,&p[i].y);

        sort(p,p+n,cmp_x);//sort的速度要比qsort的快一些貌似。。

        printf("%.2lf\n",short_s(0,n-1)/2);

    }

    return 0;

}