最近对问题(C++)

该博客讨论了如何使用C++解决最近点对问题,首先从一维情况简化问题,然后逐步推广到二维情况。博主通过将点按x坐标排序并分割,采用递归方式在左右两侧寻找最近点对,再考虑中间区域的6个关键点进行计算,以降低复杂度。文章还提及了算法设计、复杂度分析和GitHub上的源码链接。
摘要由CSDN通过智能技术生成

最近对问题(C++)

问题

问题描述:平面上有n个点P1,P2,…,Pn,n>1,Pi的直角坐标为(Xi,Yi),i=1,2,…,n.求距离最近的两个点及他们之间的欧几里得距离。

解析

一维情况

为了简单化问题,我们先考虑一维情形,即所有的点都在x轴上(假如顺序是有序的,也可以排好序,为简单,直接输入有序序列),x轴上某个点m将集合划分为s1和s2,并且两个集合的点数都是相同的,递归的求解两个子集最近对的点,(p1,p2)和(q1,q2)。如果,两个集合的最近点对都在子集s1和s2中,则

 distance = min(|p1-p2|,|q1-q2|)

如果,最近点对分别在s1和s2中,则一定是(p3,q3)。|p3-q3|<d,则p3和q3两者与m的距离都不超过d,且在区间(m-d,d]和(d,m+d]各有且仅有一个点。这样,就可以在线性时间内实现合并。
在这里插入图片描述

二维情况

推广一维,来看二维,我们仿照一维的情况先把所有点按照x(横坐标)从左到右升序排列。以X横坐标中间的点作为分界线.将平面的点分成左边和右边。找到左边的点中最近点对的距离d1,和右边最近点对的距离d2。distance = min(|p1-p2|,|q1-q2|)。如下图所示:
在这里插入图片描述
当前的d1,d2都是两边各自的点的最近距离,但是没有考虑到两边的点相互配对的情况,即是点对一个在左边区域,一个在右边区域。可以用上面找到的d来限制配对.即是明显超过d的两点不需要配对,如下:
在这里插入图片描述
以中间线为中心,dmin 为半径划分一个长带,最小点对还有可能存在于s1和s2的交界处,,p点和q点分别位于S1和S1的虚线范围内,只有在这个范围内,p点和q点之间的距离才会小于dmin。

不选取2dmin 带中所有的点进行计算,而是对于S1虚框范围内的p点,只选取S2虚框范围内长2dmin,宽为dmin的中的点进行计算,由dmin的意义可知S2中任何2个S中的点的距离都不小于d。由此可以推出矩形R中最多只有6个S2中的点。
具体为什么只有6个点我们不用知道,前人留下来的结论直接用就行࿰

点对问题(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
发出的红包

打赏作者

『gorgeous』

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值