Pagini recente » Cod sursa (job #985969) | Cod sursa (job #2218946) | Cod sursa (job #1500535) | Cod sursa (job #1505244) | Cod sursa (job #1035969)
#include<stdio.h>
#include<algorithm>
#include<math.h>
using namespace std;
#define INF 2000000000
#define NMAX 100003
struct punct
{
int x,y;
};
int n,na;
punct v[NMAX],a[NMAX];
inline bool cmpx(const punct &a,const punct &b)
{return a.x<b.x;}
inline bool cmpy(const punct &a,const punct &b)
{return a.y<b.y;}
inline double dist(punct i,punct j)
{
return sqrt((i.x-j.x)*(i.x-j.x) + (i.y-j.y)*(i.y-j.y));
}
inline double min(double a,double b)
{
if(a-b>0.000001)
return b;
return a;
}
inline double caz3(int l,int r)
{
double sol=INF;
if(r-l+1>=2)
sol=dist(v[l],v[l+1]);
if(r-l+1==3)
{
sol=min(sol,dist(v[l],v[l+2]));
sol=min(sol,dist(v[l+1],v[l+2]));
}
return sol;
}
double rez(int l,int r)
{
double sol=INF;
/*if(r-l+1<=3)
{
sol=caz3(l,r);
return sol;
}*/
if(l>=r)
return INF;
if(l==r-1)
return dist(v[l],v[r]);
int med=(l+r)>>1,i,j,p=1;
sol=min(rez(l,med) , rez(med+1,r));
na=0;
for(i=l;i<=r;++i)
if(abs(v[med].x-v[i].x)<=sol)
a[++na]=v[i];
sort(a+1,a+1+na,cmpy); //
for(i=1;i<=na;++i)
{
while(a[i].y-a[p].y>sol)
++p;
for(j=p;j<i;++j)
sol=min(sol,dist(a[i],a[j]));
}
return sol;
}
int main()
{
freopen("cmap.in","r",stdin);
freopen("cmap.out","w",stdout);
int i;
double sol;
scanf("%d",&n);
for(i=1;i<=n;++i)
{
scanf("%d%d",&v[i].x,&v[i].y);
}
sort(v+1,v+1+n,cmpx); //
sol=rez(1,n);
printf("%.6lf",sol);
return 0;
}