hdu 1007

题目链接

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

简单地描述一下题目的

求两个点之间的最小距离的一半

(0,0) (1,1)

之间的距离是 2开根号 再除以2 大小是 0.707140...

#include<cmath>
#include<algorithm>
using namespace std;
int n;
struct node
{
    double x;
    double y;
}p[100005];
int a[100005];
double cmpx(node a,node b)
{
    return a.x<b.x;
}
double cmpy(int a,int b)
{
    return p[a].y<p[b].y;
}
double min(double a,double b)
{
    return a<b?a:b;
}
double dis(node a,node b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double find(int l,int r)
{
    // 相邻的两个点直接距离
    if(r == l+1)
        return dis(p[l],p[r]);
    // 三个点，情况也不多，直接计算
    
    if(l+2 == r)
        return min(dis(p[l],p[r]),min(dis(p[l],p[l+1]),dis(p[l+1],p[r])));
    //  取出中间点
    
    int mid = (l+r)/2;
    //  继续往下递归 ， 有点模仿二分查找  找两个点之间的距离
    double ans=min(find(l,mid),find(mid+1,r));
    int i,j,cnt=0;
    //  找到距离比相邻两个点还小的点（以x 排序）
     
    for(i=l;i<=r;i++)
    {
        if(p[i].x>=p[mid].x-ans && p[i].x<=p[mid].x+ans)
            a[cnt++]=i;
    }
    //  按照Y的坐标进行排序
    
    sort(a,a+cnt,cmpy);
    
     // 从所有比最小的距离还小点中亮亮进行比较找到 激励最近的那个点
    for(i=0;i<cnt;i++)
    {
        for(j=i+1;j<cnt;j++)
        {
            /**
             *  一开始是按照下进行排序的 
             *  现在按照y进行排序，如果y的值大于当前最小值的话，直接就结束了
             *  后边的值肯定是大于当前的最小值的
             *  这个优化很关键
             */
            if(p[a[j]].y-p[a[i]].y>=ans) break;
            ans=min(ans,dis(p[a[i]],p[a[j]]));
        }
    }
    return ans;
}
int main()
{
    int i;
    while(scanf("%d",&n)!=EOF)
    {
        if(!n) break;
        for(i=0;i<n;i++)
            scanf("%lf %lf",&p[i].x,&p[i].y);
        //  先按照 x 坐标进行排序
        sort(p,p+n,cmpx);
        printf("%.2lf\n",find(0,n-1)/2);
    }
    return 0;
}