Cod sursa(job #2554471)

Utilizator vladth11Vlad Haivas vladth11 Data 22 februarie 2020 22:42:53
Problema Cele mai apropiate puncte din plan Scor 95
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 2.25 kb
//Scuzati bulaneala
#include <bits/stdc++.h>
#define debug(x) cerr << #x << " " << x << "\n"

using namespace std;

pair <int,int> v[100001];
double distanta(pair <int,int> n, pair <int,int> m)
{
    pair <double,double> a = {n.first,n.second},b = {m.first,m.second};
    return sqrt((b.first - a.first) * (b.first - a.first) + (b.second - a.second) * (b.second - a.second));
}
class InParser {
private:
	FILE *fin;
	char *buff;
	int sp;

	char read_ch() {
		++sp;
		if (sp == 4096) {
			sp = 0;
			fread(buff, 1, 4096, fin);
		}
		return buff[sp];
	}

public:
	InParser(const char* nume) {
		fin = fopen(nume, "r");
		buff = new char[4096]();
		sp = 4095;
	}

	InParser& operator >> (int &n) {
		char c;
		while (!isdigit(c = read_ch()) && c != '-');
		int sgn = 1;
		if (c == '-') {
			n = 0;
			sgn = -1;
		} else {
			n = c - '0';
		}
		while (isdigit(c = read_ch())) {
			n = 10 * n + c - '0';
		}
		n *= sgn;
		return *this;
	}

	InParser& operator >> (long long &n) {
		char c;
		n = 0;
		while (!isdigit(c = read_ch()) && c != '-');
		long long sgn = 1;
		if (c == '-') {
			n = 0;
			sgn = -1;
		} else {
			n = c - '0';
		}
		while (isdigit(c = read_ch())) {
			n = 10 * n + c - '0';
		}
		n *= sgn;
		return *this;
	}
};
set <pair <int,int> > s;
int main()
{
    InParser cin("cmap.in");
    ofstream cout("cmap.out");
    ios_base::sync_with_stdio(false);
    int n,i;
    cin >> n;
    for(i = 1; i <= n; i++)
    {
        cin >> v[i].first >> v[i].second;
        //s.insert(v[i]);
    }
    double minim = (double)(1 << 30);
    sort(v + 1,v + n + 1);
    int j = 1;
    for(i = 1; i <= n; i++)
    {

        while(v[i].first - v[j].first >= minim * minim)
        {
            s.erase({v[j].second,v[j].first});
            j++;
        }
        auto it1 = s.lower_bound({v[i].second - minim, v[i].first});
        auto it2 = s.upper_bound({v[i].second + minim, v[i].first});

        for (auto it = it1; it != it2; ++it)
        {
            int dx =  it->second;
            int dy = it->first;
            minim = min(minim, distanta({dx,dy},v[i]));
        }
        s.insert({v[i].second, v[i].first});
    }
    cout << fixed << setprecision(6) << minim;
    return 0;
}