计算机几何基础，计算多边形的重心，HDU-1115

我们已经能算多边形的面积，之后我们在往下引申一下，就要求多边形的重心。

给你2分钟思考一下，应该如何计算多边形的重心呢？

所以说这个是能推广的

上面的图片你是否看得懂呢？

下面我用更简单的思路说一下

给你一个多边形，然后呢如果多边形的质量仅仅只在定点，并且均匀分布，那么就可以用求三角形重心的方法求多边形的重心

但是如果多边形的质量是均匀分布在其内部的区域上的，也就是和面积有关，，我们就需要换另一种方式

我们首先把多边形分割成N-2个小三角形，之后求出每个三角形的重心，与面积。

每个三角形自己的重心乘自己的面积，然后加在一起，除以这个多边形的总面积，就可以求出重心的位置了

下面我们来开一道例题     HDU-1115

代码如下

 

#include<iostream>
#include<stdio.h>
using namespace std;
struct centre{
	double x,y;
};
int cas,n;
double Area(centre p0,centre p1,centre 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;	
}
int main(){
	int n;
	centre p0,p1,p2;
	scanf("%d",&n);
	while (n--)
	{
		double sum_x=0,sum_y=0,sum_area=0;
		double area;
		int m;
		scanf("%d",&m);
		scanf("%lf%lf",&p0.x,&p0.y);
		scanf("%lf%lf",&p1.x,&p1.y);
		for(int i=2;i<m;i++)
		{
			scanf("%lf%lf",&p2.x,&p2.y);
			area=Area(p0,p1,p2);
			sum_area+=area;
			sum_x+=(p0.x+p1.x+p2.x)*area;
			sum_y+=(p0.y+p1.y+p2.y)*area;//不再这里/3是因为会损失精度
			p1=p2;
		}
		printf("%.2lf %.2lf\n",sum_x/sum_area/3,sum_y/sum_area/3);
	}
	return 0;
}