• HDU 1007 Quoit Design最近点对( 分治法)


    题意:

      给出平面上的n个点,问任意点对之间的最短距离是多少? 

    思路:

      先将所有点按照x坐标排序,用二分法将n个点一分为二个部分,递归下去直到剩下两或一个点。对于一个部分,左右部分的答案分别都知道,那么可能的答案就是min(left_ans,right_ans) 。注意更小的点对可能一个在左,一个在右。所以还得处理两个边内的紧靠着的部分,如果左边的一个点到达中线的距离已经超过当前最短距离,那么这个点到达右边任意一个点也不会是最短距离了。同时,若一左一右的两个点的y距离已经超过目前最短距离,那么也不可能是最短距离,可能为答案的点对并不太多,大概O(n)对。因此,还得在x距离满足的情况下,y也得满足,所以y也得排序。

     1 #include <bits/stdc++.h>
     2 using namespace std;
     3 
     4 struct node
     5 { 
     6     double x,y;
     7 }a[100005];  
     8 
     9 int c[100005];  
    10 
    11 double cmpy(int t1,int t2) { return a[t1].y<a[t2].y;}
    12 bool cmp(node t1,node t2) { return t1.x<t2.x; }  
    13 double dis(node t1,node t2) {return sqrt((t1.x-t2.x)*(t1.x-t2.x)+(t1.y-t2.y)*(t1.y-t2.y));}  
    14 
    15 double find(int left,int right)  
    16 {  
    17     if(left+1==right)    return dis(a[left],a[right]);  
    18     if(left==right)        return 2147483.999;  
    19 
    20     int mid=(left+right)>>1;
    21     double aa=find(left,mid);
    22     double bb=find(mid+1,right);
    23     double ans=min(aa,bb); //当前最小
    24 
    25     int cnt=0;
    26     double x=a[mid].x;
    27     for(int i=left; i<=right; i++)
    28         if(fabs(a[i].x-x)<ans)    c[cnt++]=i;
    29     sort(c,c+cnt,cmpy);        //按y来排序
    30 
    31 
    32     for(int i=0;i<cnt;i++)        //计算x=[mid-ans,mid+ans]
    33         for(int j=i+1;j<cnt;j++)    
    34         {    
    35             if(a[c[j]].y-a[c[i]].y>ans)    break; //两点的距离已超过ans
    36             ans=min(ans,dis(a[c[i]],a[c[j]]));    //求最小距离
    37         }   
    38     return ans;
    39 }  
    40 
    41 int main()  
    42 {  
    43     int n,i;  
    44     while(cin>>n,n)  
    45     {  
    46         for(i=0;i<n;i++)     
    47             scanf("%lf%lf",&a[i].x,&a[i].y);
    48         sort(a,a+n,cmp);      //按X排序
    49         printf("%.2lf
    ",find(0,n-1)/2);
    50     }  
    51     return 0;  
    52 }
    AC代码
  • 相关阅读:
    写在vue总结之前(一)
    前端应该掌握的web基础和网络知识
    sass之为什么要使用预处理器
    ThinkPHP简单的验证码实现
    ajax接收php返回得到一堆html代码
    Bootstrap 4,“未捕获错误:Bootstrap工具提示需要Tether(http://github.hubspot.com/tether/)”
    百度AI开放平台- API实战调用
    最短路径算法—Dijkstra(迪杰斯特拉)算法分析与实现(C/C++)
    C#避免踩坑之如何添加paint事件
    php插入mysql中文数据出现乱码
  • 原文地址:https://www.cnblogs.com/xcw0754/p/4881625.html
Copyright © 2020-2023  润新知