算法作业5

问题：
给定平面S上n个点，找其中的一对点，使得在n个点组成的所有点对中，该点对间的距离最小。

解析：
在n个点中找出距离最近的一对点，我们首先可以想到的方法就是暴力穷举的方法：遍历整个点集，计算每一对点之间的距离d，就能找出最短距离dmin。
显然这样的方法时间复杂度为O(n2),时间复杂度太高，若问题规模较大，穷举的算法并不是一个合适的算法。

在利用分治法解决该问题时，我们先要考虑划分，即将整个点集划分为两各部分S1和S2，每个子点集的点的数量都为n/2,最近点对会出现以下三种情况：

点对在S1中；
点对在S2中；
点对分别在S1和S2中。
利用递归分别计算前两种情况，最后分析第三种情况。然后合并，分析比较出三者中的最小值。

设计：

#include <iostream>
#include <cmath>
#include <algorithm>
#include <fstream>
#define MAX 100000
using namespace std;

//点结构
struct point {
	double x, y;
}p[MAX];

int a[MAX];			//保存排序索引

//距离
double distance(point a, point b) {
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

bool cmpx(const point a, const point b) {
	return a.x < b.x;
}

bool cmpy(int &a, int &b) {
	return p[a].y < p[b].y;
}

//最近点对
double closeset(point p[], int low, int high) {
	if (low == high)
		return MAX;
	if (low + 1 == high)
		return distance(p[low], p[high]);
	int mid = (low + high) / 2;
	double ans = min(closeset(p, low, mid), closeset(p, mid + 1, high));
	int i, j, s = 0;
	for (i = low; i <= high; i++) {
		if (p[mid].x - ans <= p[i].x && p[i].x <= p[mid].x + ans)
			a[s++] = i;
	}
	sort(a, a + s, cmpy);

	for (i = 0; i < s; i++) {
		int k = i + 7 > s ? s : i + 7;
		for (j = i + 1; j < k; j++) {
			if (p[a[j]].y - p[a[i]].y > ans)
				break;
			ans = min(distance(p[a[i]], p[a[j]]), ans);
		}
	}
	return ans;
}

int main() {
	int n, i = 0;
	cout << "输入点的个数：";
	cin >> n;
	point *p = new point[n];
	ifstream infile;
	infile.open("点.txt", ios::in);	//相对路径
	if (!infile.is_open()) {
		cout << "Error: Open file failure." << endl;
		system("pause");
		return 0;
	}

	while (!infile.eof() && i < n) {
		infile >> p[i].x >> p[i].y;
		i++;
	}
	infile.close();
	sort(p, p + n, cmpx);
	cout << "最近距离为：" << closeset(p, 0, n - 1) / 2 << endl;

	system("pause");
	return 0;
}

分析：
时间复杂度为O(nlogn)