hdoj1007

由于这道题数据的特点，按照x排序会莫名其妙的TLE，按y排序就好了，也不需要什么辅助数组，一个p数组足矣，546ms AC

按照数组已经排序好的特点，如果从左边取一个点，它和右边取的一个点的y差值大于min(d)，就不需要继续从右边取点了，因为右边的点也肯定大于min(d)，这是只要将左边的点向后取一个就好了

AC喜洋洋~~~~被按x排序坑的吐血了~~

回过头来会找找更好的算法，摆脱这种不稳定性的~~~

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;

typedef struct {
	double x;
	double y;
} point;
point p[100000+10];
bool cmp(point a,point b){
	return a.y<b.y;
}
double dis(point a,point b){
	return pow(a.x-b.x,2)+pow(a.y-b.y,2);
}
double min(double a,double b){
	if(a<b) return a;
	else return b;
}
double divide(int s,int e){
	int m,i,j;
	double min1,mind,min2;
	if(s+1==e) return dis(p[s],p[e]);
	if(s+2==e) return min(dis(p[s],p[s+1]),min(dis(p[s],p[e]),dis(p[s+1],p[e])));
	m=(s+e)/2;
	min1=min(divide(s,m),divide(m+1,e));
	mind=sqrt(min1);
	for(i=s;i<=m;i++)
		for(j=m+1;j<=e;j++){                     
			if(p[j].y-p[i].y>mind) break;
			if(p[j].x-p[i].x>mind||p[i].x-p[j].x>mind) continue;
			min2=dis(p[i],p[j]);
			if(min2<min1) min1=min2;
		}
	return min1;
}

int main(){
	int n;
	int i;
	double radius;
	while(scanf("%d",&n)){
		if(n==0) break;
		for(i=1;i<=n;i++) scanf("%lf %lf",&p[i].x,&p[i].y);
		sort(p+1,p+n+1,cmp);
		radius=divide(1,n);
		radius=sqrt(radius)/2;
		printf("%.2lf\n",radius);
	}
	return 0;
}