题意：给10^5个点的坐标，求出最小的两个点之间的距离

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <string>
#include <vector>
#include <queue>   
#include <map>             
#define LL long long    
#define inf 999999999
#define N 100100
using namespace std;   //hdu 1007 （分治法）
/*
题意：给10^5个点的坐标，求出最小的两个点之间的距离
思路:用分治法（要么最小距离的两个点都在一边，要么一个在一边，另一个在另一边）
*/
struct zuobiao
{
    double x, y;
}a[N], b[N];
bool cmp(zuobiao a1, zuobiao b1)
{
    return a1.x<b1.x;
}
bool cmp1(zuobiao a1, zuobiao b1)
{
    return a1.y>b1.y;
}
double dis(zuobiao a1, zuobiao b1)
{
    return sqrt((a1.x-b1.x)*(a1.x-b1.x)+(a1.y-b1.y)*(a1.y-b1.y));
}
double d;
void f(int l, int r)
{
    if(l==r||!d)return ;
    double k;
    if(r-l==1)
    {
        k=dis(a[l], a[r]);
        if(d>k)
            d=k;
        return ;
    }
    int t, j, g, mid=(l+r)>>1;
    f(l, mid);
    f(mid+1, r);
    for(t=mid+1, g=0; t<=r; ++t)   //处理点在两边的情况
    {
        if(a[t].x-a[mid].x>=d)
            break;
        else
        {
            b[g].x=a[t].x;
            b[g++].y=a[t].y;
        }
    }
    for(t=mid; t>=l; --t)
    {
        if(a[mid+1].x-a[t].x>=d)
            break;
        else 
        {
            b[g].x=a[t].x;
            b[g++].y=a[t].y;
        }
    }
    sort(b, b+g, cmp1);
    for(j=0; j<g; ++j)
    {
        for(t=j+1; t<g; ++t)
        {
            if(b[j].y-b[t].y>=d)
                break;
            k=dis(b[j], b[t]);
            if(d>k)
                d=k;
        }
    }
    return ;
}

int main()
{
    int n, t;
    while(scanf("%d", &n)&&n)
    {
        for(t=0; t<n; ++t)
        {
            scanf("%lf%lf", &a[t].x, &a[t].y);
        }
        sort(a, a+n, cmp);
        d=inf;
        f(0, n-1);
        printf("%.2lf\n", d/2);
    }
    return 0;
}