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);
}
}