算法分析与实践-作业5-最近对问题

1. 问题

要求在一个包含n个点的集合中，找出距离最近的两个点

2. 解析

分解

• 对所有的点按照x坐标（或者y）从小到大排序
• 根据下标进行分割，使得点集分为两个集合

解决

• 递归的寻找两个集合中的最近点对。
• 取两个集合最近点对中的最小值$min(distanc{e}_{left}，distanc{e}_{right})$

合并

• 最近距离不一定存在于两个集合中，可能一个点在集合A，一个点在集合B，而这两点间距离小于dis

3. 设计

EfficientClosestPair(point p[], point q[], int n)
使用分治算法求解最近对问题
输入：数组p中存储了平面上的$n\geq2$个点，并且按照这些点的x轴坐标升序排列
	   数组q存储了与p相同的点，只是它是按照这点的y轴坐标升序排列
输出：最近点对的距离
if (n == 2) {
		return distance(p[0], p[1]);
}
else if(n==3) {
	return min(distance(p[0], p[1]), distance(p[1], p[2]));
}
else{
	将p的前[n/2]个点复制到p1
	将q的前[n/2]个点复制到q1
	将p中余下的[n/2]个点复制到pt
	将q中余下的[n/2]个点复制到qt
	d1 = EfficientClosestPair(p1, q1, n / 2);
	dt = EfficientClosestPair(pt, qt, n - n / 2);
	d = min(d1, dt);
	m = p[index].x;
	将q中所有|x-m|<d的点复制到数组s
	dminsq = d * d;
	for (i = 0; i < num - 2; i++) {
		k = i + 1;
		while (k <= num - 1 && (s[k].y - s[i].y)*(s[k].y - s[i].y) < dminsq) {
			dminsq = min((s[k].x - s[i].x)*(s[k].x - s[i].x) + (s[k].y - s[i].y)*(s[k].y - s[i].y), dminsq);
			k++;		
		}
	}
}
	return sqrt(dminsq);
	

4. 分析

$$T(n)=2T(\frac{n}{2})+f(n)$$

5. 源码

https://github.com/QAQnoname/homework/blob/master/%E4%BD%9C%E4%B8%9A5/%E6%9C%80%E8%BF%91%E5%AF%B9%E9%97%AE%E9%A2%98.cpp