http://acm.nyist.net/JudgeOnline/problem.php?pid=3

一道计算几何求多边形重心问题,

题意:已知一多边形没有边相交,质量分布均匀。顺序给出多边形的顶点坐标,求其重心。

1,质量集中在顶点上。n个顶点坐标为(xi,yi),质量为mi,则重心
  X = ∑( xi×mi ) / ∑mi
  Y = ∑( yi×mi ) / ∑mi
  特殊地,若每个点的质量相同,则
  X = ∑xi / n
  Y = ∑yi / n

2,质量分布均匀。这个题就是这一类型,算法和上面的不同。
  特殊地,质量均匀的三角形重心:
  X = ( x0 + x1 + x2 ) / 3
  Y = ( y0 + y1 + y2 ) / 3

这一题这两种方法都用到了,首先是把该多边形划分成多个三角形,并求出每个三角形的面积和重心,然后以每个三角形的重心为顶点,构造出一个新的多边形,而该多边形和原多边形的重心相同,只是该多边形的质量都集中在顶点上了,因此可以用第一种方法来求,又因为质量和面积成正比,因此可以转化为求三角形面积,至于如何求三角形面积可以用叉积,注意当三角形在多边形外边时为负,这里判断三角形面积为正还是负的方法还可以用右手螺旋定则,如果是顺时针这位负,逆时针为正

AC代码:

#include<cstdio>
#include<cmath>
#include<string>
#define exp 1e-7
using namespace std;
typedef struct str
{
  double x;
  double y;
  str():x(0.0),y(0.0){}
}point;
point s[10005];
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
		int n;
		scanf("%d",&n);
		for(int i=0;i<n;++i)
		scanf("%lf%lf",&s[i].x,&s[i].y);
		double sum=0.0;
		point cur;
		for(int i=1;i<=n;++i )
		{
			double temp=(s[i%n].x*s[i-1].y-s[i-1].x*s[i%n].y)/2;//以原点为基准点
			sum+=temp;
			cur.x+=temp*(s[i%n].x+s[i-1].x)/3.0;
			cur.y+=temp*(s[i%n].y+s[i-1].y)/3.0;
		}
		if(fabs(sum-0.0)<exp) printf("0.000 0.000\n");
		else  printf("%.3lf %.3lf\n",sum,(cur.x+cur.y)/sum);
	}return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值