算法分析与设计 实验5 平面最近点对

1、问题:

求平面中距离最近的两个点

2、解析:

可以通过枚举两两顶点来得到结果,但是时间复杂度来到了 O ( n 2 ) O(n^2) On2,显然是很慢的,所以我们采用分治的策略来分析这个问题。
这里的分治算法的主要思路是将平面上的n个点分为两个子集S1,S2。每个子集中有n/2个点,然后递归的在每个子集中求解最近点对。
算法步骤:

  • 把所有的点按照横坐标排序
  • 将所有点分成左右两部分
  • 递归分别算出左右部分的最近点对,设为d1和d2
  • 求出一个点在左部分,一个点在右部分的最近点对,设为d3
  • 结果就是min(d1, d2, d3)

3、设计(核心伪代码):

void Closest_Pair(int left, int right, int flag) {       //最近点对距离  先对x排序,flag表示禁用的点
    if(left == right) return ;    //只剩1个点
    if(left + 1 == right){
        if(P[left].id != flag && P[right].id != flag){
            if(dis > Dist(P[left], P[right])){
                dis = Dist(P[left], P[right]);
                if(flag == -1){
                    pos[0] = P[left].id;
                    pos[1] = P[right].id;
                }
            }
        }
        return ;
    }
    int mid = (left + right) / 2;           //分治
    Closest_Pair(left, mid, flag);    // 求s1(集合)中的最近点对
    Closest_Pair(mid + 1, right, flag);   //求s2中的最近点对
    int k = 0;
    for(int i = left; i <= right; i++){        //在s1和s2中间附近找可能的最近点对
        if(abs(P[mid].x - P[i].x) <= dis && P[i].id != flag)      //按x坐标来找
            tmp_P[k++] = P[i];
    }
    sort(tmp_P, tmp_P + k, cmpy);     //对y排序 用于剪枝
    for(int i = 0; i < k; i++){
        for(int j = i + 1; j < k && tmp_P[j].y - tmp_P[i].y < dis; j++){
            if(dis > Dist(tmp_P[i], tmp_P[j])){
                if(flag == -1){
                    pos[0] = tmp_P[i].id;
                    pos[1] = tmp_P[j].id;
                }
                dis = Dist(tmp_P[i], tmp_P[j]);
            }
        }
    }
}

4、分析:

时间复杂度为 O ( n l o g n ) O(nlogn) Onlogn

5、源码:

平面最近点对源码

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值