最近点对问题 c++

最近点对的算法还是比较清楚的,在实际代码阶段要让效率达到最高,需要先对x排好序,递归过程中保证不再更改(因为x只用于分割这一步)。之后对y的排序要融合归并排序进去,即把元素分割到1个,这样保证排好序。之后在合并步骤只要合并两个排序好的数组,O(N)即可解决。

#include <bits/stdc++.h>
typedef long long ll;
using namespace std;

typedef tuple<double, double> P;
int INF = 1 << 30;

bool cmp_y(P a, P b)
{
   
    return get<1>(a) < get<1>(b);
}

double calculate_dist(P a, P b)
{
   
    double x = get<0>(a) - get<0>(b);
    double y = get<1>(a) - get<1>(b);
    return sqrt(x * x + y * y);
}

// 调用时,横坐标都是排好序的
double min_dist(vector<P>& v, int l, int r)
{
   
    if (l == r)return INF;
    int mid = l + (r - l) / 2;
    double d = min(min_dist(v, l, mid), min_dist(v, mid + 1, r));
    double mid_x = get<0>(v[mid]
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
问题(closest pair problem)是指在一个平面上给定n个,找到其中距离最近的两个的距离。这个问题可以用分治法来解决。具体的实现方法可以参考以下步骤: 1. 将所有按照x轴坐标排序 2. 分治,将集分为左右两部分,递归查找左右两部分的最小距离 3. 在中间区域内(横跨左右两部分的区域),找到距离中线最近的两个的距离 4. 比较以上三个距离,取最小值作为最终结果 下面是一个简单的C++实现代码: ```c++ #include <iostream> #include <algorithm> #include <cmath> using namespace std; struct Point { int x, y; }; bool cmpX(Point a, Point b) { return a.x < b.x; } bool cmpY(Point a, Point b) { return a.y < b.y; } double getDistance(Point a, Point b) { return sqrt(pow(a.x - b.x, 2) + pow(a.y - b.y, 2)); } double merge(Point* points, int left, int mid, int right, double d) { Point* sortedPoints = new Point[right - left + 1]; int index = 0; int i = left, j = mid + 1; while (i <= mid && j <= right) { if (points[i].y < points[j].y) { sortedPoints[index++] = points[i++]; } else { sortedPoints[index++] = points[j++]; } } while (i <= mid) sortedPoints[index++] = points[i++]; while (j <= right) sortedPoints[index++] = points[j++]; for (int i = 0; i < index; i++) { for (int j = i + 1; j < index && sortedPoints[j].y - sortedPoints[i].y < d; j++) { double distance = getDistance(sortedPoints[i], sortedPoints[j]); if (distance < d) { d = distance; } } } return d; } double closestPair(Point* points, int left, int right) { if (left >= right) { return INT_MAX; } int mid = (left + right) / 2; double d1 = closestPair(points, left, mid); double d2 = closestPair(points, mid + 1, right); double d = min(d1, d2); return merge(points, left, mid, right, d); } int main() { int n; cin >> n; Point* points = new Point[n]; for (int i = 0; i < n; i++) { cin >> points[i].x >> points[i].y; } sort(points, points + n, cmpX); double ans = closestPair(points, 0, n - 1); cout << ans << endl; return 0; } ``` 其中,cmpX和cmpY是用于排序的比较函数,getDistance用于计算两之间的距离,merge用于计算跨越左右两部分的区域内的最小距离,closestPair用于递归查找左右两部分的最小距离并返回最终结果。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值