【HDU1115】多边形重心的计算

1.如何计算一个多边形的重心？如果是三角形我们知道他就是对应点坐标的算术平均值。推广成多边形，其实也是类似的，只不过是点的加权平均值。

2.计算多边形的重心，我们可以把多边形分解成n-2个三角形，然后计算这n-2个三角形的重心，之后再计算这些重心的加权平均值。对于每个三角形的重心对于整个多边形的重心的贡献，就是他的权。而这个权就是这个三角形的面积。

                     

OK，知道了这些，那么就好办了。直接写代码了。但是这里面有一个坑点，不知道是我自己写的有问题，还是这里面的确存在这个坑点。就是在计算三角形的面积的时候，我们可以选择海伦公式（毕竟，编程时候最方便），也可以用叉积计算。但是我使用海伦公式的时候前几次都WA了，实在是找不到哪里错了，就换成了叉积。然后A了。

#include<bits/stdc++.h>
using namespace std;
#pragma warning(disable:4996)
struct P
{
	P() {};
	double x, y;
	
};
//double dis(P a, P b)
//{
//	return sqrt(pow(a.x - b.x, 2) + pow(a.y - b.y,2));
//}
//double Area(P a, P b, P c)
//{
//	double l1 = dis(a, b);
//	double l2 = dis(a, c);
//	double l3 = dis(b, c);
//	double p = (l1 + l2 + l3) / 2;
//	double s = sqrt(p*(p - l1)*(p - l2)*(p - l3));
//	return s;
//}
double Area(P p1, P  p2, P p3) 
{
	return ((p2.x - p1.x)*(p3.y - p1.y) - (p3.x - p1.x)*(p2.y - p1.y)) / 2;
}
int main()
{
	int T;
	scanf("%d", &T);
	while (T--)
	{
		int n;
		scanf("%d", &n);
		P a ,b, c;
		double sumX = 0;
		double sumY = 0;
		double sumArea = 0;
		scanf("%lf%lf%lf%lf", &a.x, &a.y, &b.x, &b.y);
		for (int i = 2; i < n; i++) 
		{
			scanf("%lf%lf", &c.x, &c.y);
			double area = Area(a, b, c);
			sumX += (a.x + b.x + c.x)*area;
			sumY += (a.y + b.y + c.y)*area;
			sumArea += area;//计算所有权值
			b = c;
		}
		double ansx = sumX / sumArea / 3;//求的多边形重心
		double ansy = sumY / sumArea / 3;
		printf("%.2lf %.2lf\n", ansx, ansy);
	}
}