参考B站讲解:最近点对问题
1.最近点对问题描述:
给定n个二维空间中的点 ( x,1 )()...(n),寻找出距离最近的一对点。严格来说可能不止一对,这里只求一对就行了。
例如:正漫步于校园的你,突然需要打一个电话,在遍布于校园之中的各个公用电话中,你当然想找离自己最近的那部。W而分治策略则能设计有效算法来找出最近点对。
2.解题思路:
(1)将二维X轴上点和Y轴上的点排序存到相应数组X【】和Y【】
(2)找到x轴中间值进行划分,得到左半区域和右半区域 求出两区域中最小距离Dx和Dy,取它们最小值D1
(3)找出一个点位于左边,另外一个点位于右边的最小值D2,比较D1和D2,递归得出最小距离D
☆重点注意:鸽笼原理(一如果有 n+1 个对象放入 n 个容器中,那么至少有一个容器中会包含两个或更多的对象。)
(1)推理:根据左右区域最小距离推理出右边区域方块(长为& 高为2&)中最多有六个点
(2)得出上下4个的比较点:由下图可得,左遍区域的点与在右边区域中可能小于&的点有右边区域上下4个点(共8个点)范围内。
(3)得出结果:历遍左边区域的点和右边相应8个点的距离,得出左右两边的最小距离。
3.伪代码:
代码:
#include<iostream>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
struct Point{
double x;
double y;
};
// 计算两个点之间的欧几里得距离
double calculateDistance(Point a, Point b) {
double dx = a.x - b.x;
double dy = a.y - b.y;
return sqrt(dx * dx + dy * dy);
}
//让X进行升序排序
bool comparePointsX(const Point& p1, const Point& p2) {
return p1.x < p2.x;
}
double zuiJingDiandui( vector<Point> points)
{
int n = points.size();
double minDistance;
if(n <= 1)
{
return 0;
}
if(n == 2)
{
double d1 = calculateDistance(points[0],points[1]);
return d1;
}
if(n == 3)
{
double d1 = calculateDistance(points[0],points[1]);
double d2 = calculateDistance(points[1],points[2]);
double d3 = calculateDistance(points[2],points[3]);
return min(min(d1,d2),min(d2,d3));
}
if(n > 3)
{
// 按照x坐标排序
sort(points.begin(), points.end(),comparePointsX);
// 分割点集
int mid = n / 2; // 找到中间点索引
vector<Point> leftPoints(points.begin(), points.begin() + mid); // 分割出左侧的点集
vector<Point> rightPoints(points.begin() + mid, points.end()); // 分割出右侧的点集
// 递归计算左右子集的最近点距离
double leftDistance = zuiJingDiandui(leftPoints);
double rightDistance = zuiJingDiandui(rightPoints);
// 找到左右子集的最小距离
minDistance = min(leftDistance, rightDistance);
// 找到可能跨越子集边界的最近点对
vector<Point> stripPoints;
double midLine = points[mid].x; // 中线的x坐标
for (int i = 0; i < n; i++) {
if (fabs(points[i].x - midLine) < minDistance) { // 将横坐标与中线横坐标距离小于最小距离的点放入stripPoints
stripPoints.push_back(points[i]);
}
}
// 在stripPoints中寻找更小的距离
int stripSize = stripPoints.size();
for (int i = 0; i < stripSize; i++) {
for (int j = i + 1; j < stripSize && (stripPoints[j].y - stripPoints[i].y) < minDistance; j++) {
double distance = calculateDistance(stripPoints[i], stripPoints[j]);
minDistance = min(minDistance, distance);
}
}
}
return minDistance;
}
int main()
{
int N;
cout<<"请输入点的个数"<<endl;
cin>>N;
vector<Point> points;
double x,y;
for(int i=0 ;i<N;i++)
{
cout<<"请输入第"<<i+1<<"个点的x值和y的值:"<<endl;
cin>>x>>y;
points.push_back(Point{x, y});
}
if(zuiJingDiandui(points)==0)
{
cout<<"点数小于2 不存在最近点"<<endl;
}
else
cout<<"最小点对距离为:"<<zuiJingDiandui(points)<<endl;
return 0;
}
运行结果 :