hdoj 1007 临近点对算法模板

本文介绍了一种利用分治策略求解二维平面上点集中最近点对问题的方法。通过将点集沿X轴排序并递归地划分为两部分,计算每部分内的最短距离以及跨越中线的最短距离,最终合并得到整个点集的最近点对。该算法使用了C++实现,并考虑了边界情况。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

思想 分治
所有点按x坐标排序
每次把整体分成两边求出两边所有点中的最短距离然后求出中间一下片点的最短距离最后合并得总体最短距离

#include<cstdio>
#include<vector>
#include<algorithm>
#include<iostream>
#include<cmath>
#define LL long long
#define pi acos(-1.0)
#define maxn 100005
using namespace std;
struct node{
    double x,y;
}no[maxn],no2[maxn];
int n;
double min(double a,double b)
{
    return a<b?a:b;
}
bool cmpX(node a,node b)
{
     return a.x<b.x;
}
bool cmpY(node a,node b)
{
     return a.y<b.y;
}
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 closet(int s,int e)
{
    if(e-s==1) return dis(no[s],no[e]);
    if(e-s==2) return min(min(dis(no[s],no[s+1]),dis(no[s+1],no[e])),dis(no[s],no[e]));
    int mid = (s+e)>>1;
    double ans = min(closet(s,mid),closet(mid+1,e));

    int l = 0; //合并
    for(int i=s;i<=e;i++) if(no[i].x>no[mid].x-ans&&no[i].x<no[mid].x+ans)no2[l++] = no[i];
    sort(no2,no2+l,cmpY);
    for(int i=0;i<l;i++)
        for(int j=i+1;j<l;j++)
        {
            if(no2[j].y-no2[i].y>=ans)break;
            ans = min(ans,dis(no2[i],no2[j]));
        }
    return ans;
}
int main()
{
    while(scanf("%d",&n)&&n)
    {
        for(int i=0;i<n;i++)scanf("%lf %lf",&no[i].x,&no[i].y);
        sort(no,no+n,cmpX);
        printf("%.2f\n",closet(0,n-1)/2.0);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值