HDU ACM 1007 Quoit Design 分而治之的方法，最近点

意甲冠军：给n坐标点。半一对点之间的距离所需的距离最近。

分析：分而治之的方法，最近点。

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

#define N 100005

double min(double a,double b)
{
	return a<b?
a:b;
}

struct POINT
{
	double x,y;
};
POINT point[N],*px[N],*py[N];

double dis(const POINT* p1,const POINT* p2)
{
	return sqrt(pow(p1->x-p2->x,2.0)+pow(p1->y-p2->y,2.0));
}

bool cmpx(const POINT* p1,const POINT* p2)
{
	return p1->x<p2->x;
}

bool cmpy(const POINT* p1,const POINT* p2)
{
	return p1->y<p2->y;
}

double core(int s,int e)  //分治法求最小点对核心代码
{
	int mid,i,j,cnt;
	double ans;

	if(s+1==e)        //仅仅有两个点的情况
		return dis(px[s],px[e]);
	if(s+2==e)        //仅仅有三个点的情况
		return min(dis(px[s],px[s+1]),min(dis(px[s+1],px[e]),dis(px[s],px[e])));
	mid=(s+e)>>1;
	ans=min(core(s,mid),core(mid+1,e));  //递归求解
	for(cnt=0,i=s;i<=e;i++)    //把x坐标在px[mid].x-ans~px[mid].x+ans范围内的点取出来
		if(px[i]->x>=px[mid]->x-ans && px[i]->x<=px[mid]->x+ans)
			py[cnt++]=px[i];
	sort(py,py+cnt,cmpy);        //按y值排序
	for(i=0;i<cnt;i++)
		for(j=i+1;j<cnt;j++)
		{
			if(py[j]->y-py[i]->y>=ans)
				break;
			ans=min(ans,dis(py[i],py[j]));
		}
	return ans;
}

int main()      
{
	int i,n;

	while(scanf("%d",&n)!=EOF && n)
	{
		for(i=0;i<n;i++)
		{
			scanf("%lf%lf",&point[i].x,&point[i].y);
			px[i]=&point[i];
		}
		sort(px,px+n,cmpx);       //按x坐标排序
		printf("%.2lf\n",core(0,n-1)/2.0);
	}
    return 0;      
}

版权声明：本文博主原创文章，博客，未经同意不得转载。