题意:给你n个点,让你求出n个点里距离最近的两点距离。
思路:分治法先把点按照坐标x进行排序,划分为左右两部分,
最近点对有三种情况,都在左半部分,都在右半部分,两点分别在左边和右边
递归求出都在左边和都在右边的情况,选一个最小值m1。
然后求出两边都有的情况。分别在两边的情况可以对点进行进一步的筛选,将点的下标
保存在t[],对筛选出的点按照y坐标进行排序,方便求出最小值。
代码如下:
#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<math.h>
using namespace std;
struct point
{
double x,y;
};
point p[100005];
int t[100005];
double min(double a,double b)
{
return (a>b)?b:a;
}
double dis(point a,point b)
{
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
bool cmp(const point &a,const point &b)
{
return a.x<b.x;
}
bool cmp1(int a,int b)
{
return p[a].y<p[b].y;
}
double spit(int a,int b)
{
int i,j,c;
double m1;
if(a+1==b)
return dis(p[a],p[b]);
if(a+2==b)
return min(dis(p[a],p[a+1]),min(dis(p[a],p[b]),dis(p[b],p[a+1])));
c=(a+b)>>1;
m1=min(spit(a,c),spit(c+1,b));
int k=0;
for(i=a;i<=b;i++)
{
if(p[i].x>=p[c].x-m1&&p[i].x<=p[c].x+m1)
t[k++]=i;
}
sort(t,t+k,cmp1);
for(i=0;i<k;i++)
for(j=i+1;j<k;j++)
{
if(p[t[j]].y-p[t[i]].y>m1) break;
m1=min(m1,dis(p[t[i]],p[t[j]]));
}
return m1;
}
int main()
{
int n;
int i;
double r;
while(scanf("%d",&n)){
if(n==0) break;
for(i=0;i<n;i++) scanf("%lf %lf",&p[i].x,&p[i].y);
sort(p,p+n,cmp);
r=spit(0,n-1);
r=r/2;
printf("%.2lf\n",r);
}
return 0;
}