UVA 10245 The Closest Pair Problem 最近点问题 分治算法

题意，给出n个点的坐标，找出两点间最近的距离，如果小于10000就输出INFINITY。

纯暴力是会超时的，所以得另辟蹊径，用分治算法。

递归思路将点按坐标排序后，分成两块处理，最近的距离不是在两块中的一块中，就会存在于跨越中线的点对中。

查找跨越中线的点比较麻烦，之前已经求出两块中的最小距离，只要在x范围在[m-d,m+d]的点中找对，更新最小距离，最后返回最小距离即可。

代码：

 

 /*

 *   Author:        illuz <iilluzen[at]gmail.com>

 *   Blog:          http://blog.csdn.net/hcbbt

 *   File:          uva10245.cpp

 *   Lauguage:      C/C++

 *   Create Date:   2013-09-04 19:57:26

 *   Descripton:    UVA 10245 The Closest Pair Problem, partitation

 */

#include <cstdio>

#include <cstdlib>

#include <cmath>

#include <algorithm>

using namespace std;

#define rep(i, n) for (int i = 0; i < (n); i++)

#define sqr(x) ((x) * (x))

#define dis(i, j) (sqrt(sqr(p[i].x - p[j].x) + sqr(p[i].y - p[j].y)))



const int MAXN = 10010;

struct Point {

	double x;

	double y;

} p[MAXN];

int n, t[MAXN];

double Min;



bool cmp1(Point a, Point b) {

	return a.x < b.x;

}



bool cmp2(int a, int b) {

	return p[a].y < p[b].y;

}



double solve(int l, int r) {

	if (l == r) return 0xfffffff;

	if (l + 1 == r) return dis(l, r);

	int mid = (l + r) / 2;

	double d = min(solve(l, mid), solve(mid + 1, r));

	int c = 0;

	for (int i = l; i <= r; i++)

		if (fabs(p[mid].x - p[i].x) <= d)

			t[c++] = i;

	sort(t, t + c, cmp2);

	for (int i = 0; i < c; i++)

		for (int j = i + 1; j < c && p[t[j]].y - p[t[i]].y < d; j++) {

			double td = dis(t[i], t[j]);

			if (d - td > 1e-9) d = td;

		}

	return d;

}



int main() {

	while (scanf("%d", &n) && n) {

		rep(i, n) scanf("%lf%lf", &p[i].x, &p[i].y);

		sort(p, p + n, cmp1);

		Min = solve(0, n - 1);

		if (Min - 10000 <= 1e-9) printf("%.4lf\n", Min);

		else printf("INFINITY\n");

	}

	return 0;

}