多边心重心[nyoj 3 多边形重心问题, hdu 1115 Lifting the Stone]

Nyoj 3 多边形重心：点击打开链接

Hud 1115 Lifting the Stone：点击打开链接

题目意思都很简单，就是求一个多变型的重心。以前没有接触过多边形重心的求解问题。还好学习了，一下。

计算几何中：
三角形的重心： x = (xa+xb+xc)/3, y = (ya+yb+yc)/3;
四边形的重心：作一对角线，将它分成两个三角形分别求出重心与面积 (x1,y1) ,s1 ; (x2, y2), s2 则该四边形的重心为： x = (x1*s1+x2*s2)/(s1+s2), y = (y1*s1+y2*s2)/(s1+s2);
五边形则分为一个三角形与一个四边形……

任意多边形中直接取任一点（一般为原点）把多边形分为n-2个三角形 分别求重心
x=∑si*xi/∑si
y=∑si*yi/∑si
si为每块三角形的有向面积，这里说的是“有向面积”！无论是否为凸多边形，这个公式都是成立的。

以上就是求解多变型重心的问题，有了以上的知识+叉积求解三角形的面积，这个问题就手到擒来了。

code：

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;

const int N = 1e6 + 5;
const double eps = 1e-6;
struct Point {
    double x, y;
}p[N];

struct Vector{
    double x, y;
};

double Cross(Point a, Point b, Point c){
    Vector A, B;
    A.x = b.x - a.x; A.y = b.y - a.y;
    B.x = c.x - a.x; B.y = c.y - a.y;
    return A.x * B.y - B.x * A.y;
}

int main(){
    int T;
    scanf("%d", &T);
    while(T --){
        int n;
        scanf("%d", &n);
        for(int i = 0; i < n; i ++){
            scanf("%lf %lf", &p[i].x, &p[i].y);
        }
        double ansx = 0 , ansy = 0, Tarea = 0, allarea = 0;
        for(int i = 1; i < n - 1; i ++){
            Tarea = Cross(p[0], p[i], p[i + 1]);
            allarea += Tarea;
            ansx += Tarea * (p[0].x + p[i].x + p[i + 1].x);
            ansy += Tarea * (p[0].y + p[i].y + p[i + 1].y);
        }
        ansx /= (allarea * 3.0);ansy /= (allarea * 3.0);
        printf("%.2lf %.2lf\n",ansx + eps, ansy + eps);
    }
    return 0;
}

在nyoj上没有卡精度，但是在poj和hdu上都卡精度问题了。

这充分的体现了计算几何的素养。不是你的算法不到位，而是你是计算几何素养不到位。

在出现-0.000上是怎么处理的。等等，很多的问题，有待以后的遇见和解决。