UVA

一道平面分治的板子题（做之前刚学会平面分治，太菜了）

把所有点按照x排序后，分成左半边的最近距离和右半边的最近距离，再以中点为圆心，目前的最近距离为半径，做这个圆内的枚举，求最短距离。

注意：sqrt()里面的数如果是int类型会比double类型慢，如果次数多了会T

代码如下：

#include<bits/stdc++.h>
using namespace std;
struct node{
    double x,y;
}a[10005];
bool cmp(node& a, node& b){
    return a.x < b.x;
}
double dis(int i,int j){
    return sqrt((a[i].x-a[j].x)*(a[i].x-a[j].x)+(a[i].y-a[j].y)*(a[i].y-a[j].y));
}
double erfen(int l,int r){
    if (l>r) return 40004;
    int mid=(l+r)>>1;
    int st=mid,ed=mid;
    double d=std::min(erfen(l,mid-1),erfen(mid+1,r));
    while (st>=l && a[st].x+d>a[mid].x) st--;
    while (ed<=r && a[ed].x-d<a[mid].x) ed++;
    for (int i=st+1; i<ed; i++){
        for (int j=i+1; j<ed; j++){
            d=std::min(dis(i,j),d);
        }
    }
    return d;
}
int main(){
    int n;
    while (scanf("%d",&n)){
        if (n==0) break;
        for (int i=1; i<=n; i++){
            scanf("%lf%lf",&a[i].x,&a[i].y);
        }
        sort(a+1,a+1+n,cmp);
        double ret=erfen(1,n);
        if (ret>10000) printf("INFINITY
");
        else printf("%.4f
",ret);
    }
}