HDU 1007 Quoit Design(分治法求最近点对问题)

本题的大致题意,求出最近两个点的距离的一半。
首先先看数据范围,如果暴力的话肯定会超时的,所以要想一个优化点的方法,就是对其进行分治。

大致思路:先把一个大区间分为两个小区间,然后分别对两个区间进行排序,筛选。然后再把筛选后的点放入一个新的数组里面进行求值。然后不断的更新最小值。具体看代码吧,代码中我也给出了大量的注释。
AC代码:

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 100005
struct node
{
    double x,y;
} point[N];

int point1[N];//存储符合标准的点的坐标

bool cmpx(node a,node b)//按x坐标从小到大排序
{
    return a.x<b.x;
}

bool cmpy(int a,int b)//按y坐标从小到大进行排序
{
    return point[a].y<point[b].y;
}

double dist(node a,node b)//求出两点间的距离
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

double near(int l,int r)//利用分治法找出最近的两个点的距离
{
    if(l>=r)
        return 1e10;
    int mid=(l+r)/2;));
    int n=0;
    double distan=min(near(l,mid),near(mid+1,r));//利用二分法分为两个区间
    for(int i=mid; i>=l; --i)//把前一半符合条件的点的位置存入数组q中
    {
        if(point[mid].x-point[i].x<distan)
        {
            point1[n++]=i;
        }
        else break;
    }
    for(int i=mid+1; i<=r; ++i)//把后一半符合条件的点的位置存入数组q中
    {
        if(point[i].x-point[mid].x<distan)
        {
            point1[n++]=i;
        }
        else break;
    }
    sort(point1,point1+n,cmpy);//对这些点按纵坐标进行升序排序
    for(int i=0; i<n; i++)
    {
        for(int j=i+1; j<n; j++)
        {
            if(point[point1[j]].y-point[point1[i]].y<distan)//找出这些点中距离的最小值
            {
                distan=min(dist(point[point1[i]],point[point1[j]]),distan);
            }
            else break;
        }
    }
    return distan;
}

int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        if(n==0)
            break;
        for(int i=0; i<n; i++)
        {
            scanf("%lf %lf",&point[i].x,&point[i].y);
        }
        sort(point,point+n,cmpx);
        printf("%.2lf\n",near(0,n-1)/2);
    }
    return 0;
}

你可能感兴趣的:(HDU)