#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <string>
#include <vector>
#include <queue>
#include <map>
#define LL long long
#define inf 999999999
#define N 100100
using namespace std; //hdu 1007 (分治法)
/*
题意:给10^5个点的坐标,求出最小的两个点之间的距离
思路:用分治法(要么最小距离的两个点都在一边,要么一个在一边,另一个在另一边)
*/
struct zuobiao
{
double x, y;
}a[N], b[N];
bool cmp(zuobiao a1, zuobiao b1)
{
return a1.x<b1.x;
}
bool cmp1(zuobiao a1, zuobiao b1)
{
return a1.y>b1.y;
}
double dis(zuobiao a1, zuobiao b1)
{
return sqrt((a1.x-b1.x)*(a1.x-b1.x)+(a1.y-b1.y)*(a1.y-b1.y));
}
double d;
void f(int l, int r)
{
if(l==r||!d)return ;
double k;
if(r-l==1)
{
k=dis(a[l], a[r]);
if(d>k)
d=k;
return ;
}
int t, j, g, mid=(l+r)>>1;
f(l, mid);
f(mid+1, r);
for(t=mid+1, g=0; t<=r; ++t) //处理点在两边的情况
{
if(a[t].x-a[mid].x>=d)
break;
else
{
b[g].x=a[t].x;
b[g++].y=a[t].y;
}
}
for(t=mid; t>=l; --t)
{
if(a[mid+1].x-a[t].x>=d)
break;
else
{
b[g].x=a[t].x;
b[g++].y=a[t].y;
}
}
sort(b, b+g, cmp1);
for(j=0; j<g; ++j)
{
for(t=j+1; t<g; ++t)
{
if(b[j].y-b[t].y>=d)
break;
k=dis(b[j], b[t]);
if(d>k)
d=k;
}
}
return ;
}
int main()
{
int n, t;
while(scanf("%d", &n)&&n)
{
for(t=0; t<n; ++t)
{
scanf("%lf%lf", &a[t].x, &a[t].y);
}
sort(a, a+n, cmp);
d=inf;
f(0, n-1);
printf("%.2lf\n", d/2);
}
return 0;
}
题意:给10^5个点的坐标,求出最小的两个点之间的距离
最新推荐文章于 2023-03-21 14:07:47 发布