Cod sursa(job #709473)

Utilizator feelshiftFeelshift feelshift Data 8 martie 2012 09:56:28
Problema Infasuratoare convexa Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 2.62 kb
// http://infoarena.ro/problema/infasuratoare
#include <fstream>
#include <iomanip>
#include <cmath>
#include <algorithm>
#include <vector>
using namespace std;

#define x first
#define y second

ifstream in("infasuratoare.in");
ofstream out("infasuratoare.out");

const int oo = 0x3f3f3f3f;

int nrPoints,selectedPoints;
pair<double,double> minPoint = make_pair(oo,oo);
vector< pair<double,double> > points,convexHull;

double dist(pair<double,double> one,pair<double,double> two);
bool leftTurn(pair<double,double> one,pair<double,double> two,pair<double,double> three);
double determinant(pair<double,double> one,pair<double,double> two,pair<double,double> three);

class compare {
public:
	bool operator() (pair<double,double> one,pair<double,double> two) {
		/*double oneAtan = atan2(one.y - minPoint.y,one.x - minPoint.x);
		double twoAtan = atan2(two.y - minPoint.y,two.x - minPoint.x);

		if(oneAtan != twoAtan)
			return (oneAtan < twoAtan);
		else
			return (dist(minPoint,one) < dist(minPoint,two));*/
		return (leftTurn(minPoint,one,two));
	}
};

int main() {
	in >> nrPoints;

	pair<double,double> point;
	for(int i=1;i<=nrPoints;i++) {
		in >> point.x >> point.y;

		if(point.x < minPoint.x || (point.x == minPoint.x && point.y < minPoint.y))
			minPoint = point;

		points.push_back(point);
	}

	//out << "Initial point: " << minPoint.x << " " << minPoint.y << "\n";

	points.erase(find(points.begin(),points.end(),minPoint));
	sort(points.begin(),points.end(),compare());

	selectedPoints = 2;
	points.insert(points.begin(),minPoint);	
	
	convexHull.push_back(minPoint);
	convexHull.push_back(*(points.begin()+1));

	for(int i=1;i<nrPoints;i++) {
		while(!leftTurn(convexHull[selectedPoints-2],convexHull[selectedPoints-1],points[i]))
			convexHull.pop_back(), selectedPoints--;
		convexHull.push_back(points[i]), selectedPoints++;
	}

	/*for(unsigned int i=0;i<points.size();i++)
		out << points[i].x << " " << points[i].y << "\n";
	out << "\n\n";*/

	//out << "Convex hull:\n";
	out << selectedPoints << fixed << setprecision(6) << "\n";
	for(unsigned int i=0;i<convexHull.size();i++)
		out << convexHull[i].x << " " << convexHull[i].y << "\n";
	//out << "\n";

	return (0);
}

double dist(pair<double,double> one,pair<double,double> two) {
	return (sqrt(pow(one.x - two.x,2) + pow(one.y - two.y,2)));
}

bool leftTurn(pair<double,double> one,pair<double,double> two,pair<double,double> three) {
	return (determinant(one,two,three) > 0);
}

double determinant(pair<double,double> one,pair<double,double> two,pair<double,double> three) {
	return (one.x * two.y + three.x * one.y + two.x * three.y
		- three.x * two.y - one.x * three.y - two.x * one.y);
}