最近点对问题

分治——最近点对问题

问题描述
n个点在公共空间中，求出所有点对的欧几里得距离最小的点对。

分治法求解

分解
对所有的点按照x坐标（或者y）从小到大排序（排序方法时间复杂度O ( n l o g n ）。
根据下标进行分割，使得点集分为两个集合。
解决
递归的寻找两个集合中的最近点对。
取两个集合最近点对中的最小值m i n ( d i s l e f t ， d i s r i g h t )。
合并
最近距离不一定存在于两个集合中，可能一个点在集合A，一个点在集合B，而这两点间距离小于dis。

有已mid为中心，

// divideandconquer.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"



struct point {
	double x;
	double y;
}P[100];
double distance(point p1, point p2) {
	return sqrt((p1.x - p2.x)*(p1.x - p2.x) + (p1.y - p2.y)*(p1.y - p2.y));
}
bool cmp1(point p1, point p2) {
	return p1.x < p2.x;
}
bool cmp2(point p1, point p2) {
	return p1.y < p2.y;
}
//分治法
double nearest_pair(point S[], int left, int right) {
	cout << left << " " << right << endl;
	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[right + 1]);
		double d3 = distance(S[right + 1], S[left]);
		d2 = min(d1, d2);
		d3 = min(d2, d3);
		return d3;
	}

	int m = (right + left) / 2;
	double d1 = nearest_pair(S, left, m);
	double d2 = nearest_pair(S, m + 1, right);
	//sort(S+right, S+left, cmp2);
	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++;
	sort(S + l, S + r + 1, cmp2);
	double d3;
	for (int i = l; i <= r; i++) {
		for (int 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;
}



int _tmain(int argc, _TCHAR* argv[])
{
	int n;
	cout << "Input n:";
	cin >> n;
	for (int i = 1; i <= n; i++) {
		cout << "Input the " << i << "th number:";
		cin >> P[i].x >> P[i].y;
	}

	sort(P + 1, P + n + 1, cmp1);

	for (int i = 1; i <= n; i++) {
		cout << P[i].x << " " << P[i].y << endl;
	}

	double m2 = nearest_pair(P, 1, n);
	cout << m2 << endl;



	system("pause");
	return 0;
}