算法分析：求最近点对问题（c++）

##算法分析：求最近点对问题（c++）
1、 准备部分：
对于最近点问题，需要将问题代码化，所以要先创建一个点类，来使得计算使用时候更加方面，

class Points    //点类
{
   
public:
        float x;
        float y;
};

因为会多次计算距离，需要一个计算两点之间距离的函数。

float Distance(Points p1,Points p2)//求两点距离函数
{
       return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
但这里改进将开根号放在外面可以使得计算量减少。
float Distance(Points p1,Points p2)//求两点距离函数
{
       return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}

2、 穷举法
算法思想：
遍历每个个点，将每个点与另外其他所有个点求距离，然后找到最短的那根距离，同时保留点对返回。

实现代码：
遍历每个点，使用两个循环，每次求和，第i个只用和第i+1到n的数去比较，因为前面的数已经比较过，不必重复计算。一开始最小值为无穷大，计算出来的距离一个个比较，如果比最小值小，就将其的值替换，保存点对。

for (int i=0; i<length-1; i++)
    {
          for(int j=i+1; j<length; j++)
          {
               float temp = Distance(p[i],p[j]);
               if (temp < min_dis)
               {
                    min_dis = temp;
                    a = p[i];
                    b = p[j];
                }
          }
    }

3、 分治法
函数名：

float merge_distance(Points* po,int len,Points &a,Points &b)

参数：
po排好序的点集，len为点的数量，a保存第一个点，b保存第二点。

算法思想及代码解释：
把二维先看成一维数组，按照x轴将所有点排序。
代码实现：

    sort(p, p+length, CmpX);                       //排序
min_dis = merge_distance(p, length, a, b);  //调用函数

点数较少情形。

代码实现：

    if(len < 2)   return na_n;//无穷大
    if(len == 2)
    {
        a = po[0];
         b = po[1];
         dis = Distance(po[0],po[1]);
     }

点数|S|>