HDU1007(最近点对问题)

题意不难理解,就是找到最近的两个点,计算其距离,除以2就是所求的圆的半径...

思路很简单,运用分治的思想,先划分区间,分别找到左右区间中的最近点对,再合并区间,找到区间间的最近点对...

注意如果用qsort()进行排序可能会超时...

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <algorithm>

using namespace std;

int N,top;
struct node
{
	double x,y;
}point[100005],stack[100005];


int cmpx(struct node a,struct node b)
{
	return (a.x<b.x);
}

int cmpy(struct node a,struct node b)
{
	return (a.y<b.y);
}

double min(double left,double right)
{
	return left>right?right:left;
}

double dis(double x1,double y1,double x2,double y2)
{
	return sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));
}

double findpoint(int left,int right)
{
	double dl,dr,dm,d;
	int mid,i,j;
	mid=(left+right)/2;
	if (right-left==1)//处理两个点的情况,这种情况下可以直接得出距离
		return dis(point[left].x,point[left].y,point[right].x,point[right].y);
	else if (right-left==2)//处理三个点的情况,也可以计算最短距离
	{
		dl=dis(point[left].x,point[left].y,point[mid].x,point[mid].y);
		dr=dis(point[mid].x,point[mid].y,point[right].x,point[right].y);
		dm=dis(point[left].x,point[left].y,point[right].x,point[right].y);
		return min(min(dl,dr),dm);
	}
	else//按照分治的思想,二分处理其余情况,划分区间使问题规模缩小到两点以及三点情况,
	{
		dl=findpoint(left,mid);//左区间内的最短距离
		dr=findpoint(mid+1,right);//右区间内的最短距离
		d=min(dl,dr);//合并区间时要处理[mid-d,mid+d]范围内的点,找出这些点的最短距离,若找到的dm小于d,则刷新d的值
		top=0;
		for(i=left;i<=right;i++)
		{
			if (point[i].x>=point[mid].x-d&&point[i].x<=point[mid].x+d)
			{
				stack[top++]=point[i];//将[mid-d,mid+d]范围内的点转存到stack数组中,便于后续的排序
			}
		}
		sort(stack,stack+top,cmpy);//根据y值升序排序
		for (i=0;i<top;i++)
		{
			for(j=i+1;j<top&&j<=top+11;j++)//只需计算i往后的前11个点
			{
				if(fabs(stack[j].y-stack[i].y)>d) break;//若stack[j].y与stack[i].y的差已经大于d,那么接下来不可能找到与i点距离小于d的,break掉
				dm=dis(stack[i].x,stack[i].y,stack[j].x,stack[j].y);//计算dm
				if (dm<d) d=dm;//刷新d值
			}
		}
		return d;
	}
}
int main(void)
{
	int i;
	while(scanf("%d",&N)==1&&N!=0)
	{
		for (i=0;i<N;i++)
		{
			scanf("%lf %lf",&point[i].x,&point[i].y);
		}
		sort(point,point+N,cmpx);//根据x值升序排序
		printf("%.2lf\n",findpoint(0,N-1)/2);
	}
	return 0;
}