《算法导论》实验五：最近点对算法（C++）

一、题目描述

在n>=2个点的集合Q中寻找最近点对。
“最近”是指通常意义下的欧几里得距离：即点p1（x1，y1）和p2（x2，y2）之间的距离为：sqrt（（x1-x2）2 +（y1-y2）2）。

二、算法设计与分析

算法主要思想就是分治。
情况1：点数小于等于二时：直接计算，求该两点之间的距离。
情况2：集合中有三个点：两两比较，求三个点中的最近的两个点距离。
情况3：点数大于三时：首先划分集合S为SL和SR，使得SL中的每一个点位于SR中每一个点的左边，并且SL和SR中点数相同。分别在SL和SR中解决最近点对问题，得到DL和DR，分别表示SL和SR中的最近点对的距离。令d=min(DL,DR)。如果S中的最近点对(P1,P2)。P1、P2两点一个在SL和一个在SR中，那么P1和P2一定在以L为中心的间隙内，以L-d和L+d为界。
如下图-1所示：

图-1 点数大于三时时分治

如果在SL中的点P和在SR中的点Q成为最近点对，那么P和Q的距离必定小于d。因此对间隙中的每一个点，在合并步骤中，只需要检验yp+d和yp-d内的点即可。
步骤1：根据点的y值和x值对S中的点排序。
步骤2：找出中线L将S划分为SL和SR
步骤3：将步骤2递归的应用解决SL和SR的最近点对问题，并令d=min(dL,dR)。
步骤4：将L-d~L+d内的点以y值排序，对于每一个点(x1,y1)找出y值在y1-d~y1+d内的接下来的7个点，计算距离为d’。如果d’小于d，令d=d’，最后的d值就是答案。

三、实验结果与分析

本实验的所用到的测试用例是用随机函数生成的，所以每次的测试用例都有所不同。但可通过控制台设定需要测试点的个数N。
(一)、当点个数N=1时（输入不合理）。见下图-2：

图-2 N=1时（输入不合理）

(二)、当点个数N=2时，见下图-3：

图-3 N=2的输出结果

(三)、当点个数N=3时，见下图-4：

图-4 N=3的输出结果

(四)、当点个数N>3时，见下图-5：

图-5 N=15的输出结果

四、实验总结

1、 采用分治法寻找最近点对时，相较于寻找一维的最近点来说，二维的最近点对寻找要困难许多，难点在于分界线周围的点的