UVA 10002 Center of Masses【求凸包重心】

题目大意：给出凸多边形上点（注意！不是按照一定顺序给出，可能是乱序），求该凸多边形的质心，另外当凸多边形顶点数小于3，结束输入。

解题策略：凸多边形顶点无序给出—>求出凸包（使凸多边形以顶点逆序保存）—>求凸包质心（由于二维几何平面质量均匀）—>求凸包重心。

                     求凸多边形重心，参看链接：http://www.cnblogs.com/jbelial/archive/2011/08/08/2131165.html

/*
   UVA 10002 Center of Masses
   AC by J_Dark
   ON 2013/5/6 20:31
   Time 0.305s     跑得好快 = =
*/
#include 
#include 
#include 
#include 
#include 
using namespace std;
/////////////////////////////////////
struct point{
	double x, y;
	point(double a, double b){
		x = a;
		y = b;
	}
	double Distance(point t){
		return sqrt((x-t.x)*(x-t.x) + (y-t.y)*(y-t.y));
	}
};
vector p;
vector CH; //存放凸包顶点序号   模拟栈
int testCase, top, cc=0, nodeNum;
/////////////////////////////////////
void Input(){
	p.clear();
	CH.clear();
	CH.resize(nodeNum+5);
	double xx, yy;
	for(int i=0; i> xx >> yy;
		p.push_back(point(xx, yy));
	}
}

bool cmp(point a, point b){
	if(a.y == b.y)  return a.x < b.x;
	return a.y < b.y;
}

bool turnRight(point px1, point px2, point pp){
    const double eps = 1e-20;
	if((px2.x - px1.x)*(pp.y - px2.y) - (pp.x - px2.x)*(px2.y - px1.y) <= eps) return true;
	return false;
}
void ComputeCH(){
	sort(p.begin(), p.end(), cmp);
	CH[0] = 0;
	CH[1] = 1;

	top = 1;
	//从起点0到到排序最后点作凸包右链  过程1
	for(int i=2; i=0; i--){
		//top!=len, 不考虑已在过程1生成凸包上的点
		while( top!=len && turnRight(p[CH[top-1]], p[CH[top]], p[i]) )
		{
		   top--;
		}
		CH[++top] = i;
	}
}

/*计算凸多边形重心模板  该模板另外在求解的过程中，不
  需要考虑点的输入顺序是顺时针还是逆时针，相除后就抵消了。
*/
double Area(point p0, point p1, point p2){
    double area = 0;
    area =  p0.x * p1.y + p1.x * p2.y + p2.x * p0.y - p1.x * p0.y - p2.x * p1.y - p0.x * p2.y;
    return area / 2 ;
}
void ComputeGC(){
	int sPos = 0;
    double sum_x , sum_y , sum_area , area;
	sum_x = sum_y = sum_area = area = 0;
	//一定得按照序进行套用，我从最左下点开始,逆时针扫整个凸包
    for(int i=2; i> nodeNum && nodeNum >= 3)
	{
		Input();
		ComputeCH(); //计算凸包
		ComputeGC();  //计算凸多边形重心并输出
	}
	return 0;
}