Cod sursa(job #765571)

Utilizator BlaugranasEnal Gemaledin Blaugranas Data 8 iulie 2012 10:40:04
Problema Cele mai apropiate puncte din plan Scor 90
Compilator cpp Status done
Runda Arhiva educationala Marime 0.89 kb
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
struct P
{long long x,y;};
P a[100001],b[100001];
int i,n;

int A(P a,P b)
{return a.x<b.x||(a.x==b.x&&a.y<b.y);}

long long B(P p,P q)
{return (p.x-q.x)*(p.x-q.x)+(p.y-q.y)*(p.y-q.y);}

long long D(int p,int q)
{int m=(p+q)/2,i,j,l,k;
long long r=1LL<<62;
if(q-p<=1)
       return r;
if(q-p==2)
       return B(a[p],a[p+1]);
if(q-p==3)
       return min(B(a[p],a[p+2]),min(B(a[p],a[p+1]),B(a[p+1],a[p+2])));
r=min(D(p,m),D(m,q));
for(j=0,i=p;i<=q&&j<8;i++)
if(abs(a[i].x-a[m].x)<=r)
       b[++j]=a[i];
for(l=1;l<j;l++)
for(k=l+1;k<=j;k++)
       r=min(r,B(b[l],b[k]));
return r;}

int main()
{freopen("cmap.in","r",stdin);
freopen("cmap.out","w",stdout);
scanf("%d",&n);
for(i=1;i<=n;i++)
       scanf("%lld%lld",&a[i].x,&a[i].y);
sort(a+1,a+n+1,A);
printf("%.6lf\n",sqrt(D(1,n)));
return 0;}