hdoj1007

3 篇文章 0 订阅
#include"iostream"
#include"algorithm"
#include"stdio.h"
using namespace std;
const int vexnum=110000;
struct point
{
	double x,y;
};
bool cmpx(point &a,point&b)
{
	return a.x<b.x;
}
bool cmpy(point &a,point &b)
{
	return a.y<b.y;
}
double space(point &a,point &b)
{
	return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}
double minp(double a,double b)
{
	return a<b?a:b;
}
double f(int l,int r,point*a,point *mid)
{
	if(r-l==1) return space(a[r],a[l]);
	if(r-l==2) return minp(minp(space(a[l],a[l+1]),space(a[l],a[r])),space(a[l+1],a[r]));
	int m=(l+r)/2;
	double d=minp(f(l,m,a,mid),f(m+1,r,a,mid));
	int i,cnd=0,j;
	for(i=l;i<=r;i++)
	{
		if(a[i].x<a[m].x+d&&a[i].x>a[m].x-d)
			mid[cnd++]=a[i];
	}
	sort(mid,mid+cnd,cmpy);
	for(i=0;i<cnd;i++)
		for(j=i+1;j<cnd;j++)
		{
			if(mid[j].y-mid[i].y>d)
				break;
			d=minp(d,space(mid[i],mid[j]));
		}
		return d;
	
}
int main()
{
	//freopen("a.txt","r",stdin);
	point a[vexnum],mid[vexnum];
	int n,i;
	while(cin>>n&&n)
	{
		for(i=0;i<n;i++)
			cin>>a[i].x>>a[i].y;
		sort(a,a+n,cmpx);
		double min=f(0,n-1,a,mid);
		 printf("%0.2f\n",sqrt(min)/2);
	}
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值