分治递归--二维最近点对问题

   参考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;

}

运行结果 ：