#include "OrdenGraham.h"
#include <iostream>
 
using namespace std;
 
angulo::angulo()
{
 
};
angulo::~angulo()
{
 
};
 
float distsqr(const point &a, const point &b)
{
	return (a.x - b.x)*(a.x - b.x) + (a.y - b.y)*(a.y - b.y);
};
 
float dist(const point &a, const point &b)
{
	return sqrt(distsqr(a, b));
};
 
float cross(const point &a, const point &b, const point &c){
	float d = (b.x - a.x) * (c.y - a.y) - (b.y - a.y) * (c.x - a.x);
	if(fabs(d) < 1e-9) return 0.0;
	return d;
};
 
bool angleCmp(const point &self, const point &that)
{
    //angulo a;
	float t = cross(an.pivote, that, self);
	if(t < 0.0) return true;
	if(t == 0)
	{
		return (distsqr(an.pivote, self) < distsqr(an.pivote, that));
	}
	return false;
};
 
 
Vect graham(Vect p)
{
     angulo an;
     for(int i = 1; i < p.size(); ++i)
	 {
		if(p[i].y < p[0].y || (p[i].y == p[0].y && p[i].x < p[0].x ))
			swap(p[0], p[i]);
	 }
 
     an.pivote = p[0];
 
	 cout <<"orden"<<endl;
	 for(int i = 0; i < p.size(); ++i)
	 {
             cout << p[i].x << "  " << p[i].y << endl;
     }
     return p;
};
 

#ifndef ORDEN_GRAHAM
#define ORDEN_GRAHAM
 
#include <cmath>
#include <iostream>
#include <vector>
 
struct point
{
	float x, y;
	point() {}
	point(float X, float Y): x(X), y(Y) {}
};
 
typedef std::vector<point> Vect;
float distsqr(const point &a, const point &b);
float dist(const point &a, const point &b);
float cross(const point &a, const point &b, const point &c); 
bool angleCmp(const point &self, const point &that);
Vect graham(Vect p);
 
 
class angulo
{
      public:
             angulo();
             ~angulo();
             point pivote;
};
 
#endif