Pagini recente » Cod sursa (job #2639724) | Cod sursa (job #2065131) | Cod sursa (job #1907784) | Cod sursa (job #2675868) | Cod sursa (job #2492231)
#include <bits/stdc++.h>
#include <limits>
using namespace std;
vector <pair<long,long> > v;
long long get_distance(pair <long,long> a,pair<long,long> b)
{
return (a.first-b.first)*(a.first-b.first)+(a.second-b.second)*(a.second-b.second);
}
long long divide(long long st,long long dr)
{
long long m, dist_min;
if(dr-st<=3)
{
long long dist = std::numeric_limits<long long>::infinity();
for(int i=st;i<=dr;i++)
{
for(int j=i+1;j<=dr;j++)
{
dist = min(dist,get_distance(v[i],v[j]));
}
}
return dist;
}
m = st + (dr - st) / 2;
dist_min = min(divide(st,m),divide(m,dr));
vector < pair<long,long> > band;
for(int i=st;i<=dr;i++)
{
if(abs(v[m].first-v[i].first)<dist_min)
{
band.push_back({v[i].second,v[i].first});
}
}
sort(band.begin(),band.end());
for(int i=0;i<band.size();i++)
{
for(int j=i+1;j<band.size() && j < i + 8;j++)
{
long long dist;
dist = get_distance(band[i],band[j]);
if(dist < dist_min)
{
dist_min = dist;
}
}
}
return dist_min;
}
int main()
{
ifstream fin("cmap.in");
ofstream fout("cmap.out");
int n;
long long x,y,dist_min;
fin>>n;
for(int i = 1; i <= n; i++) {
fin>>x>>y;
v.push_back({x,y});
}
sort(v.begin(),v.end());
dist_min = divide(0,n-1);
fout<<fixed<<setprecision(6)<<sqrt((double)dist_min);
fin.close();
fout.close();
return 0;
}