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 0x7fffffffusing 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;}