分治法解最近对问题

1、问题

在数组P中存储平面上的n>=2个点，并且按照这些点的x轴坐标升序排列，数组Q中存储与P相同的点，只是它按照这些点的y 轴坐标升序排序，得出最近点之间的欧几里得距离。

2、解析

1、分解
对所有的点按照x坐标（或者y）从小到大排序（排序方法时间复杂度O(nlogn) O(nlogn)O(nlogn)）。
根据下标进行分割，使得点集分为两个集合。
2、解决
递归的寻找两个集合中的最近点对。
取两个集合最近点对中的最小值
3、合并
最近距离不一定存在于两个集合中，可能一个点在集合A，一个点在集合B，而这两点间距离小于两个集合中的最小值。

3、设计

double closestPoint(point s[], int left, int right) {
	if (right - left == 1) {
		return distance(s[right], s[left]);
	}
	if (right - left == 2) {
		double d1 = distance(s[right], s[left]);
		double d2 = distance(s[right], s[left + 1]);
		double d3 = distance(s[left + 1], s[left]);
		d2 = min(d1, d2);
		d3 = min(d2, d3);
		return d3;
	}
	int m = (right + left) / 2;
	int i, j;
	double d1 = closestPoint(s, left, m);
	double d2 = closestPoint(s, m + 1, right);
	double d = min(d1, d2);
	int l = left, r = right;
	while (s[l].x < s[m].x - d && l <= right);
	l++;
	while (s[r].x > s[m].x + d && r >= left)
	r--;
	double d3;
	for (i = l; i <= r; i++) {
		for (j = i + 1; j <= r; j++) {
			if (s[j].y - s[i].y >= d) {
				break;
			}
			else {
				d3 = distance(s[i], s[j]);
				if (d3 < d)
					d = d3;
			}
		}
	}
	return d;
}

4、分析

时间复杂度：O（nlogn）