一、题目描述
在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、 采用分治法寻找最近点对时,相较于寻找一维的最近点来说,二维的最近点对寻找要困难许多,难点在于分界线周围的点的