hdu1007 Quoit Design 【计算几何+分治】

题目大意：

给出平面上n个点，求最近两点距离的一半。

解题思路：

很经典的平面分治题。

假设我们把所有点按 $x$ 坐标平均分成了左右两个部分，设平分直线为$x=x_0$，那么最近点对$（p,q）$的距离就是下面二者的最小值：
(1)$p,q$同属于左半边或右半边时的最近点对距离。
(2)$p,q$属于不同区域时的最近点对距离。

对于情况(1)，是与原题完全相同的子问题，可以递归处理。设左右两部分求出的最小距离为$d$。

对于情况(2)，我们就只用考虑距离小于$d$的点对了.
不妨令$p$在$x=x_0$左边，$q$在$x=x_0$右边。
那么有$x_0-d<x_p<x_0<x_q<x_0+d$。

考虑满足上述条件的一个点$p$，那么能与它更新答案的点$q$需要满足
$x_0<x_q<x_0+d，y_p-d<y_q<y_p$
对于满足条件的点$q$同理。

可以证明对于每个点，另一侧满足条件的点最多有三个，如图：

右侧正方形中最多能有三个点，不然右侧点对最小距离就会小于d了。

这样一来，我们只需在分治同时把左右不分的点按y值归并排序，处理中间点时每个点就可以O(1)转移了。

所以总复杂度为$O(logn)$。

细节详见代码

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

const int N=100005;
struct point
{
    double x,y;
    point(){}
    point(double _x,double _y):
        x(_x),y(_y){}
    inline friend point operator - (const point &a,const point &b)
    {return point(a.x-b.x,a.y-b.y);}
    inline double dis(){return sqrt(x*x+y*y);}
}p[N];
int n;

inline bool cmp(const point &a,const point &b)
{
    if(a.x==b.x)return a.y<b.y;
    return a.x<b.x;
}

double solve(int l,int r)
{
    if(l==r)return 1e20;
    int mid=l+r>>1;
    double x0=(p[mid].x+p[mid+1].x)/2;
    double d=min(solve(l,mid),solve(mid+1,r));

    static point a[N],b[N],c[N];
    int bn=0,cn=0,l1=l,l2=mid+1;
    for(int i=l;i<=r;i++)
    {
        if(l1<=mid&&(l2==r+1||p[l1].y<=p[l2].y))
        {
            a[i]=p[l1++];
            if(x0-d<a[i].x)b[++bn]=a[i];
        }
        else
        {
            a[i]=p[l2++];
            if(x0+d>a[i].x)c[++cn]=a[i];
        }
    }
    for(int i=l;i<=r;i++)p[i]=a[i];

    for(int i=1,j=1;i<=bn||j<=cn;)
    {
        if(i<=bn&&(j==cn+1||b[i].y<c[j].y))
        {
            for(int k=j-1;k;k--)
            {
                if(b[i].y-d>=c[k].y)break;
                d=min(d,(b[i]-c[k]).dis());
            }
            i++;
        }
        else 
        {
            for(int k=i-1;k;k--)
            {
                if(c[j].y-d>=b[k].y)break;
                d=min(d,(b[k]-c[j]).dis());
            }
            j++;
        }
    }
    return d;
}

int main()
{
    //freopen("lx.in","r",stdin);
    while(scanf("%d",&n)&&n)
    {
        for(int i=1;i<=n;i++)
            scanf("%lf%lf",&p[i].x,&p[i].y);
        sort(p+1,p+n+1,cmp);
        printf("%.2f\n",solve(1,n)/2);
    }
    return 0;
}