分治法实现最近点对问题(C++)

最近点对问题

要求

  • 随机生成30个浮点数(1-10)之间
  • 求出最近距离的两点
  • C语言实现
#include <iostream>
#include <math.h>

using namespace std;

struct dot{  //点结构
    float x;
    float y;
};
void sortByX(dot dots[6],int num){  //将点按照X坐标从小到大排列好(冒泡)
    for (int i = num-1; i >= 0; --i) {
        for (int j = 0; j < i; ++j) {
            if (dots[j].x > dots[j + 1].x) {
                dot flag = dots[j];
                dots[j] = dots[j + 1];
                dots[j + 1] = flag;

            }
        }
    }
}
void sortByY(dot dots[6],int num){  //将点按照Y坐标从小到大排列好(冒泡)
    for (int i = num-1; i >= 0; --i) {
        for (int j = 0; j < i; ++j) {
            if (dots[j].y > dots[j + 1].y) {
                dot flag = dots[j];
                dots[j] = dots[j + 1];
              
  • 1
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
最近点对问题是指在平面上给定n个点,找出其中距离最近的两个点的距离。分治是求解最近点对问题的经典算之一。 1. 将所有点按照横坐标排序,然后分成左右两个子问题; 2. 对于左右两个子问题,分别递归求解最近点对; 3. 在左右两个子问题最近点对中,选取距离更近的一对作为候选解; 4. 计算横跨左右两个子问题的点对的距离,并将它们按照纵坐标排序; 5. 在纵向排序后的点集中,依次计算每个点与后面最多7个点之间的距离,并选取距离更近的一对作为候选解; 6. 在步骤3和步骤5中选取距离更近的一对作为最终解。 下面是C++代码实现: #include <iostream> #include <algorithm> #include <cmath> #include <vector> using namespace std; struct Point { int x; int y; }; bool cmp_x(Point a, Point b) { return a.x < b.x; } bool cmp_y(Point a, Point b) { return a.y < b.y; } double dist(Point a, Point b) { return sqrt(pow(a.x - b.x, 2) + pow(a.y - b.y, 2)); } double brute_force(vector<Point>& points, int lo, int hi) { double min_dist = 1e9; for (int i = lo; i <= hi; i++) { for (int j = i + 1; j <= hi; j++) { min_dist = min(min_dist, dist(points[i], points[j])); } } return min_dist; } double closest_pair(vector<Point>& points, int lo, int hi) { if (lo == hi) { return 1e9; } if (hi - lo == 1) { return dist(points[lo], points[hi]); } int mid = (lo + hi) / 2; double d1 = closest_pair(points, lo, mid); double d2 = closest_pair(points, mid + 1, hi); double d = min(d1, d2); vector<Point> strip; for (int i = lo; i <= hi; i++) { if (abs(points[i].x - points[mid].x) < d) { strip.push_back(points[i]); } } sort(strip.begin(), strip.end(), cmp_y); int n = strip.size(); for (int i = 0; i < n; i++) { for (int j = i + 1; j < n && strip[j].y - strip[i].y < d; j++) { d = min(d, dist(strip[i], strip[j])); } } return d; } int main() { int n; cin >> n; vector<Point> points(n); for (int i = 0; i < n; i++) { cin >> points[i].x >> points[i].y; } sort(points.begin(), points.end(), cmp_x); double ans = closest_pair(points, 0, n - 1); if (ans < 1e9) { cout << ans << endl; } else { cout << "No solution" << endl; } return 0; }

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值