1.问题
在一个笛卡尔平面上有若干个点,找出最近的一对点。
2.解析
数据预处理:令P为笛卡尔平面上n>1个点构成的集合,假设集合中的每个都不一样,我们还假设这些点是按照x轴坐标升序排列的。(如果不是这样,可以事先用类似合并排序这样的高效算法对其进行排序。)为了更加方便,我们还可以按照点的y轴坐标在另一个列表中进行升序排列。
问题分类讨论:
当点的数量n满足:n=2时,两点之间的距离就可以直接得出。
当点的数量n满足:n>=3时,我们可以利用点集在x轴方向上的中位数m,在该处作一条垂线,将点集分成大小分别为ceil(n/2)和floor(n/2)的两个子集Pleft和Pright。即使得其中ceil(n/2)个点位于线的左边或线上,floor(n/2)个点位于线的右边或线上。然后就可以通过递归求解子问题Pleft和Pright来得到最近点对问题的解。
注意点,以上分区域比较所得的d不一定是所有点对的最小距离,因为距离最近的两个点可能分别位于分界线的两侧。因此,在合并较小子问题的解时,需要检查是否存在这样的点。显然,我们可以只关注以分隔带为对称的、宽度为2d的垂直带中的点,因为任何其他的点对的距离都至少为d(见图a)。
3.设计
[核心伪代码]
Pair§{
if n<=3
蛮力算法求得的最小距离(距离公式)
Else
分成左右两段
将S的前⌈n/2⌉个点复制到S, 将S中余下的⌊n/2⌋个点复制到S2;
S1=Pair(lf)
S2=Pair(rt)
m=S[⌈n/2⌉-1].x //求分割线的值
将S中所有|x-m|<d的点复制到数组S[]
dmin=d2
for i=0 to num-2
k=i+1
while k<=num-1 and (s[k].y-s[i].y)2 < dmin
dmin=min((s[k].x-s[i].x)2+(s[k].y-s[i].y)2,dminsq)
k++
end if
return sqrt(dmin)
}
4.分析
[算法复杂度推导]
蛮力: O(n^2)
分治: O(nlog2n)
5.源码
https://github.com/zbadbx/suanfa/blob/main/%E6%9C%80%E8%BF%91%E7%82%B9%E5%AF%B9%E9%97%AE%E9%A2%98