Cod sursa(job #1142878)

Utilizator federerUAIC-Padurariu-Cristian federer Data 14 martie 2014 12:52:01
Problema Infasuratoare convexa Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.18 kb
#include<fstream>
#include<algorithm>
#include<vector>
#include<stack>
#include<iomanip>
#define Nmax 120010
#define x first
#define y second
using namespace std;
typedef pair<double, double> Point;
vector<Point>V;
ifstream fin("infasuratoare.in");
ofstream fout("infasuratoare.out");
Point S[Nmax];
int top;

double cross_product(Point A, Point B, Point C)
{
	return (B.x - A.x) * (C.y - A.y) - (B.y - A.y) * (C.x - A.x);
}

bool cmp(Point &a, Point &b)
{
	return cross_product(V[0], a, b)<0;
}

void Sort()
{
	int pos = 0, i;
	for (i=1; i < V.size();++i)
	if (V[i] < V[pos])
		pos = i;
	sort(V.begin() + 1, V.end(), cmp);
}

void ConvexHull()
{
	Sort();
	top = 0;
	S[++top] = V[0];
	S[++top] = V[1];
	for (int i = 2; i < V.size(); ++i)
	{
		while (cross_product(S[top-1], S[top], V[i])>0 && top >= 2)
			top--;
		S[++top] = V[i];
	}
}

int main()
{
	int N;
	double a, b;
	fin >> N;
	for (int i = 1; i <= N; ++i)
	{
		fin >> a >> b;
		Point p = make_pair(a, b);
		V.push_back(p);
	}
	ConvexHull();
	fout << top << '\n';
	fout << fixed;
	for (int i = 1; i <= top; ++i)
		fout << setprecision(12)<<S[i].x << ' ' << S[i].y << '\n';

	return 0;
}