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;
}