HDU 1007 Quoit Design（分治）

题目链接

题意

给你n个点，求平面最近点对，答案除二。

思路

暴力骗、分治皆可，仅供自己回忆向题解，甚看。
网上题解挺多，这里存个分治代码。
下面随便讲讲，没看别的题解估计看不懂我再讲啥，看了也看不懂。
先按x升序排序分治，合并操作时，按y升序排序，每个点最多只查几个点（具体懒得算）。
因为合并时满足，俩区间各自点最小距离为d，那么相邻区间合并查找，对于左边区间一个点，右边区间可能更新答案的点也不会超三个，但是合并排序了，查找略多一小点。

类似题

POJ 3714 要保证最近点对不在同一个点集，那么只需求dis时，判断是否在相同点集，是返回无穷大即可。

代码

#include <bits/stdc++.h>
using namespace std;

struct Node
{
    double x, y;
    bool operator < (const Node &_a) const
    {
        return x < _a.x;
    }
}e[100005];
int tmp[100005];

bool cmp(int a, int b)
{
    return e[a].y < e[b].y;
}

double dis(int a, int b)
{
    return sqrt((e[a].x-e[b].x)*(e[a].x-e[b].x)+(e[a].y-e[b].y)*(e[a].y-e[b].y));
}

double f(int l, int r)
{
    if(l+2 == r) return min(dis(l,r), min(dis(l+1,r), dis(l,l+1)));
    if(l+1 == r) return dis(l,r);
    int mid = (l+r) >> 1;
    double ans = min(f(l,mid), f(mid+1,r));
    int k = 0;
    for(int i = l; i <= r; ++i) if(abs(e[i].x - e[mid].x) < ans) tmp[k++] = i;
    sort(tmp, tmp+k, cmp);
    for(int i = 0; i < k; ++i)
    {
        for(int j = i+1; j < k && e[tmp[j]].y-e[tmp[i]].y < ans; ++j)
        {
            ans = min(ans, dis(tmp[i], tmp[j]));
        }
    }
    return ans;
}

int main()
{
#ifndef ONLINE_JUDGE
	freopen("1.in","r",stdin);
	freopen("1.out","w",stdout);
#endif
    int n;
    while(scanf("%d",&n), n)
    {
        for(int i = 1; i <= n; ++i) scanf("%lf%lf",&e[i].x,&e[i].y);
        sort(e+1,e+1+n);
        printf("%.2f\n",f(1,n)/2.0);
    }
	return 0;
}