hdu1007 Quoit Design 平面最近点对（分治）

题意：

就是平面最近点对←_←

思路：

先把点按x坐标排好序，利用分治的思想，比如要求区间[l, r]的最近点对，先求区间[l, mid]的最近点对，再求[mid + 1, r]的最近点对，设此时的最近距离为best。
但是best没有考虑一个点在左边，一个点在右边的情况。此时我们再枚举区间[mid - best, mid + best]中的所有点更新答案。
只不过这样裸枚举会超时，我们先把区间里的点放进一个temp数组，然后按y坐标排序，当两个点的y坐标之差超过best就没必要继续枚举了。
对了题目要求最后答案除以2。

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <climits>
#include <cmath>

#define MAXN 100005
#define INF 0x7fffffff

using namespace std;

const double eps = 1e-8;

struct Point{
    double x, y;
    Point() {}
    Point(double a, double b) : x(a), y(b) {}
    bool operator < (const Point &p) const {
        return x < p.x;
    }
};

int n;
Point p[MAXN];
Point temp[MAXN];

bool cmp(Point a, Point b) { return a.y < b.y; }

double dist(int a, int b) {
    double t = (p[a].x - p[b].x) * (p[a].x - p[b].x) + (p[a].y - p[b].y) * (p[a].y - p[b].y);
    return t;
}

double dist(Point a, Point b) {
    double t = (a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y);
    return t;
}

double dfs(int l, int r) {
    int mid = l + r >> 1;
    if (l + 1 == r) return dist(l, r);
    if (l + 2 == r) return min(dist(l, r), min(dist(l, l + 1), dist(l + 1, r)));
    double best = min(dfs(l, mid), dfs(mid + 1, r));
    int tot = 0;
    for (int i = l; i <= r; i++)
        if (p[i].x >= p[mid].x - best && p[i].x <= p[mid].x + best)
            temp[tot++] = p[i];
    sort(temp, temp + tot, cmp);
    for (int i = 0; i < tot; i++)
        for (int j = i + 1; j < tot; j++)
            if (temp[j].y - temp[i].y >= best) break;
            else best = min(best, dist(temp[i], temp[j]));
    return best;
}

int main() {
    while (~scanf("%d", &n) && n) {
        for (int i = 0; i < n; i++)
            scanf("%lf %lf", &p[i].x, &p[i].y);
        sort(p, p + n);
        double ans = sqrt(dfs(0, n - 1)) / 2;
        printf("%.2lf\n", ans);
    }
    return 0;
}