作业5：分治法求最近对

1.问题

在包含有n个点的集合S中，找出距离最近的两个点。设 p1(x1,y1)，p2(x2,y2)，……，pn(xn,yn)是平面的n个点。 严格地将，最近点对可能不止一对，此例输出一对即可。

2.解析

分治法求解：我们先根据x坐标排序，进行递归算出每一部分的距离，在根据y坐标排序计算每一部分之间的最近距离，最后得到答案。 同时我们在进行第二部分计算的时候，通过鸽舍原理易知，最多需要计算距离中心点y坐标最近的6个点来得到答案。 在利用分治法思想解决此问题时，首先考虑将最近对问题进行分治，设计其分治策略。将集合S分成两个子集S1和S2，根据平衡子问题原则，每个子集中的点数大致都为n/2。这样分治后，最近点对将会出现三种情况：在S1中，在S2中或者最近点对分别在集合S1和S2中。利用递归分析法分别计算前两种情况，第三种方法另外分析。求解出三类子情况后，再合并三类情况，比较分析后输出三者中最小的距离。

3.设计

EfficientClosestPair(P, Q)
//使用分治算法来求解最近点对问题
//输入: 数组P中存储了平面上的n≥2个点，并且按照这些点的x轴坐标升序排列
数组Q存储了与P相同的点，只是它是按照这点的y轴坐标升序排列
//输出: 最近点对之间的欧几里得距离
If n≤3
返回由蛮力算法求出的最小距离
else
将P的前[n/2]个点复制到P1
将Q的前[n/2]个点复制到Q1
将P中余下的[n/2]个点复制到Pr
将Q中余下的[n/2 ]个点复制到Qr
d1←EfficientClosestPair(P1 Q1)
dr←EfficientClosestPair(Pr Qr)
d←min{d1 dr}
m←P[「n/2」-1].x
将Q中所有|x- m| <d的点复制到数组s[0到num-1]
dminsq←d^2
for i←0 to mum—2 do
k←i+ 1
while k≤num-1 and (S[k].y- S[i].y)^2 < dminsq .
dminsq←min( (S[k].x- S[i].x)^2 + (S[k].y- S[i].y)^2，dminsq)
k←-k+1
return sqrt(dminsq)

4.分析

在分治算法中，当求解n个点的集合的最近点对时，对于上述三类情况中的前两者可由递归算得， 而分析可得第三类情况的时间代价 ，合并后问题求解的总的时间按复杂度可由下列公式来T（n）=1 （n=2）T（n）=2T（n/2）+O（n） 当n>2时 T（n）=O（nlog2n）

5.源码

https://github.com/vin-zqw/test5