ACM-HDU-1007

ACM-HDU-1007

题意大概:求n个点中最近两个点的距离的一半
思路:
用结构体来存x,y值,对结构体的x或y进行排序,从头到尾判断是否有x,y个相等的两个点,有就输出0.00
否则就用分治思想的函数:先按x坐标将进行排序,排好序的(x,y)点,分成两部分（左右）,基点定为node[mid],mid=（l+r）/2,递归的分下去直到某一部分(左或右)只剩下2个或者3个点(可以直接算出距离比较得出min)，
最后min只会在3种情况中出现:①.min在右半部分,②,min在左半部分,③.min是跨越左右两部分的点的距离的最小值.
对于①,②种情况，递归下去l== r+1 || l==r+2就可以得出min,
③种情况需要对每个点与基点进行判断(基点为node[mid]，左右两部分以x坐标排序的点横跨两部分必会经过mid的x坐标)：如果他们的x坐标距离就>=min了,就没什么意义了,再将<min的点存起来(点与点之间的距离可能会小于min),再将这些点的y坐标进行排序,将排好序的点从头到尾相比较如果y坐标的距离>=min就跳出（break，因为按y坐标排序后,不管是大—>小,还是小->大，如果连当前点与后面的点都>=min，后面未遍历更远的点就没意义进行比较了）,对于满足x,y的点进行枚举,min=min(函数)(dist(i,j),min);得到min值
代码:

#include<bits/stdc++.h>
#define maxn 100105
using namespace std;

struct spot{
    double x,y;
}node[maxn];
//按x从小到大进行排序
bool cmp(spot a,spot b){
     return a.x<b.x;
}
//按y从大到小进行排序
bool cmpy(int a,int b){
     return node[a].y>node[b].y;
}
//计算任意两点的距离
double dist(spot a,spot b){
     return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
//寻找最短距离的函数
double solve(int l,int r){
     double ans;
     int i,j,mid,cnt=0;
     //终止条件
     if(l+1==r)
        return dist(node[l],node[r]);
     if(l+2==r)
        return min(dist(node[l],node[l+1]),min(dist(node[l],node[r]),dist(node[l+1],node[r])));
     mid=(r+l)/2;
     //分治
     ans=min(solve(l,mid),solve(mid+1,r));
     int a[maxn];
     for(i=l;i<=r;i++){
        if(fabs(node[i].x-node[mid].x)<ans){
            a[cnt++]=i;
        }
     }
     sort(a,a+cnt,cmpy);
     for(i=0;i<cnt;i++){
        for(j=i+1;j<cnt;j++){
            if(node[a[i]].y-node[a[j]].y>=ans)
                break;
            ans=min(dist(node[a[i]],node[a[j]]),ans);
        }
     }
     return ans;
}
int main()
{
    int n,i,flag;
    scanf("%d",&n);
    while(n!=0){
    flag=0;
    for(i=0;i<n;i++){
        scanf("%lf%lf",&node[i].x,&node[i].y);
    }
    sort(node,node+n,cmp);
    //看是否有相同的两点
    for(i=0;i<n-1;i++){
        if(node[i].x==node[i+1].x&&node[i].y==node[i+1].y){
            flag=1;
            printf("0.00\n");
            break;
        }
    }
    if(!flag)
    printf("%.2f\n",solve(0,n-1)/2);
    scanf("%d",&n);
    }
    return 0;
}

如果对你有帮助的，小伙伴麻烦点个赞,我会持续更新HDU的题目