UVA 10002 Center of Masses【求凸包重心】_根据凸包算法求出的点计算中心点-CSDN博客

本文链接：https://blog.csdn.net/J_Dark/article/details/8892278

该博客介绍了如何解决UVA 10002题，即在给定无序的凸多边形点后，计算其质心。解题方法包括构建凸包，然后计算凸包的质心，当顶点数少于3时结束。参考链接提供了一种二维几何平面中求凸多边形重心的策略。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目大意：给出凸多边形上点（注意！不是按照一定顺序给出，可能是乱序），求该凸多边形的质心，另外当凸多边形顶点数小于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 <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
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<point> p;
vector<int> 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<nodeNum; i++){
		cin >> 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<nodeNum; i++){
		while( top && turnRight(p[CH[top-1]], p[CH[top]], p[i]) )
		{
		   top--;
		}
		CH[++top] = i;
	}

	int len = top;
	//从排序最高点到到起点0fab反向作凸包右链  过程2
	CH[++top] = nodeNum-2;
	for(int i=nodeNum-3; 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<top ; i++){
       area = Area(p[CH[sPos]],p[CH[sPos+1]],p[CH[i]]) ;
       sum_area += area ;
       sum_x += (p[CH[sPos]].x + p[CH[sPos+1]].x + p[CH[i]].x) * area ;
       sum_y += (p[CH[sPos]].y + p[CH[sPos+1]].y + p[CH[i]].y) * area ;
       p[CH[sPos+1]] = p[CH[i]] ;
    }
	printf ("%.3lf %.3lf\n", sum_x/sum_area/3, sum_y/sum_area/3);
}

///
int main(){
	while(cin >> nodeNum && nodeNum >= 3)
	{
		Input();
		ComputeCH(); //计算凸包
		ComputeGC();  //计算凸多边形重心并输出
	}
	return 0;
}