hdu1007最近点对

最新推荐文章于 2019-12-18 20:08:00 发布

原创最新推荐文章于 2019-12-18 20:08:00 发布 · 347 阅读

0 ·

CC 4.0 BY-SA版权

hdu 专栏收录该内容

21 篇文章

订阅专栏

本文介绍了一种解决最近邻点对问题的有效算法。该算法通过递归地将平面划分成子区域，并利用分治策略寻找最短距离。特别地，文章详细展示了如何在每个步骤中使用排序和比较来确定任意两点间的最小距离。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
using namespace std;
struct D
{
double x,y;
} a[100100];
int tmp[100100];

double cmp(D a,D b)
{
return a.x<b.x;
}
double cmpy(int p,int q)
{
return a[p].y<a[q].y;
}
double dis(D a,D b)
{
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double dfs(int left,int right)
{
double d;
if(left==right-1)
return dis(a[left],a[right]);
if(left==right-2)
{
d = min(dis(a[left],a[right-1]),dis(a[left+1],a[right]));
d = min(d,dis(a[left],a[right]));
return d;
}
int mid = (left+right)/2;
//printf("%d %d %lf\n",left,right,d);
d = min(dfs(left,mid),dfs(mid+1,right));
int k = 0;
for(int i = left; i <= right; i++)
{
if(a[i].x>=a[mid].x-d&&a[i].x<=a[mid].x+d)
tmp[k++] = i;
}
sort(tmp,tmp+k,cmpy);
for(int i = 0; i < k; i++)
{
for(int j = i+1; j < k; j++)
{
if(fabs(a[tmp[j]].y-a[tmp[i]].y)<d)
{
double ds = dis(a[tmp[i]],a[tmp[j]]);
if(ds<d)
d = ds;
}
else break;
}
}
return d;
}
int main()
{
int n;
while(scanf("%d",&n)&&n)
{
for(int i = 0; i < n; i++)
scanf("%lf%lf",&a[i].x,&a[i].y);
sort(a,a+n,cmp);
printf("%.2lf\n",dfs(0,n-1)/2);
memset(tmp,0,sizeof(tmp));
}
return 0;
}