hdu 1007 Quoit Design 平面最近点对

难度：5

最近听说inline可以让非递归的函数更快，但是试了一下感觉不是很好。

平面最近点对，计算机和模板题

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;
const int maxn = 100010;
struct node {
    double x , y;
};
int n , a[maxn];
node p[maxn];
double cmpx(node a, node b) {
    return a.x < b.x;
}
double cmpy(int a , int b) {
    return p[a].y < p[b].y;
}
double dis(node a , node b) {
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double min3(double a,double b,double c) {
    return (a < b)?(a < c ? a : c ) : (b < c ? b : c);
}
double find(int l, int r) {
    if(r == l+1)
        return dis(p[l] , p[r]);
    if(r == l+2)
        return min3(dis(p[l],p[r]) , dis(p[l],p[l+1]) , dis(p[l+1],p[r]));
    int mid = (l + r) >> 1;
    double ans = min(find(l,mid) , find(mid+1,r));
    int cnt = 0;
    for(int i=l;i<=r;i++)
        if(p[i].x>=p[mid].x-ans && p[i].x<=p[mid].x+ans)
            a[cnt++] = i;
    sort(a , a+cnt,cmpy);
    for(int i=0;i<cnt;i++)
    for(int j=i+1;j<cnt;j++) {
        if(p[a[j]].y-p[a[i]].y >= ans) break;
        ans = min(ans , dis(p[a[i]] , p[a[j]]));
    }
    return ans;
}
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 , cmpx);
        printf("%.2lf\n" , find(0,n-1)/2.0);
    }
    return 0;
}