Cod sursa(job #2353555)

Utilizator oanceadavidOancea David oanceadavid Data 24 februarie 2019 12:59:57
Problema Cele mai apropiate puncte din plan Scor 75
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 2.46 kb
#include <vector>
#include <iostream>
#include <algorithm>
#include <cassert>
#include <fstream>
#include <iomanip>
#include <cmath>
using namespace std;

unsigned long long distanta(const pair<int, int>& point_a, const pair<int, int>& point_b)
{
	return  ((unsigned long long)(point_a.first - point_b.first) * (unsigned long long)(point_a.first - point_b.first)
		+ (unsigned long long)(point_a.second - point_b.second) * (unsigned long long)(point_a.second - point_b.second));
}

int number_of_points(int left, int right)
{
	return right - left + 1;
}

unsigned long long cmap(const vector< pair<int, int> >& points, int left, int right)
{
	if (number_of_points(left, right) <= 3 || (points[left].first == points[right].first))
	{
		unsigned long long mini = numeric_limits<unsigned long long>::max();
		for (int i = left; i < right; i++)
		{
			for (int j = i + 1; j <= right; j++)
			{
				mini = min(mini, distanta(points[i], points[j]));
			}
		}
		return mini;
	}
	int midIdx_start = (right + left) / 2;
	int midIdx = midIdx_start;
	while (midIdx < right && points[midIdx_start].first == points[midIdx].first) midIdx++;
	midIdx--;
	// cout<<midIdx<<endl;
	unsigned long long mini = min(cmap(points, left, midIdx), cmap(points, midIdx + 1, right));
	int i, j;
	i = midIdx;
	j = midIdx + 1;

	assert(i >= left);
	assert(j <= right);
	int prev_x = points[i].first;
	int succ_x = points[j].first;
	long long distance = sqrt(mini);
	while (i >= left && (succ_x - points[i].first) < distance) i--;
	i++;
	while (j <= right && (points[j].first - prev_x) < distance) j++;
	j--;

	if (i > j) {
		//cout << i << " " << j << endl;
		//cout << " === " << points[i].first << " " << points[j].first << endl;
	}
	if (i < j) {
		//		cout << i << " " << j << endl;
				//cout << " --- " << points[i].first << " " << points[j].first << endl;
		for (int k = i; k <= midIdx; ++k) {
			for (int l = midIdx + 1; l <= j; ++l) {
				mini = min(mini, distanta(points[l], points[k]));
				//cout<<"k "<<k<<' '<<"l "<<l<<endl;
			}
		}
		//mini=min(cmap(points,i,j),mini);
	}
	return mini;
}

int main()
{
	int n;
	ifstream in("cmap.in");
	ofstream out("cmap.out");
	in.tie(NULL);
	ios::sync_with_stdio(false);
	in >> n;
	vector< pair<int, int> > points(n);
	for (int i = 0; i < n; i++)
	{
		in >> points[i].first >> points[i].second;
	}
	sort(begin(points), end(points));
	out << fixed << setprecision(6) << sqrt(cmap(points, 0, n - 1));
	return 0;
}