#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;

struct point
{
 double x;
 double y;
};

//求多邊形的重心算法
//說明:
//求多邊形重心并不是簡(jiǎn)單的把求三角形的重心公式推廣就行了
//我的算法是在平面上取一點(diǎn)(一般取原點(diǎn), 這樣可以減少很多計(jì)算, 而且使思路更清晰^_^)
//這樣就得到了N個(gè)三角形OP[i]P[i+1](其中點(diǎn)的順序要為逆時(shí)針的),
//分別求出這N個(gè)三角形的重心Ci和面積Ai(注意此處面積是又向面積, 就是用叉乘求面積時(shí)保留其正負(fù)號(hào))
//在求出A = A1+A2+...+AN(同樣保留正負(fù)號(hào)的代數(shù)相加)
//最終重心C = sigma(Ai+Ci)/A;
point gravity(point *p, int n)
{
 double area = 0;
 point center;
 center.x = 0;
 center.y = 0;

 for (int i = 0; i < n-1; i++)
 {
  area += (p[i].x*p[i+1].y - p[i+1].x*p[i].y)/2;
  center.x += (p[i].x*p[i+1].y - p[i+1].x*p[i].y) * (p[i].x + p[i+1].x);
  center.y += (p[i].x*p[i+1].y - p[i+1].x*p[i].y) * (p[i].y + p[i+1].y);
 }

 area += (p[n-1].x*p[0].y - p[0].x*p[n-1].y)/2;
 center.x += (p[n-1].x*p[0].y - p[0].x*p[n-1].y) * (p[n-1].x + p[0].x);
 center.y += (p[n-1].x*p[0].y - p[0].x*p[n-1].y) * (p[n-1].y + p[0].y);

 center.x /= 6*area;
 center.y /= 6*area;

 return center;
}