多边形重心问题 NYOJ 3 思路+代码

多边形重心问题

时间限制: 3000 ms  |  内存限制: 65535 KB
难度: 5
描述
在某个多边形上,取n个点,这n个点顺序给出,按照给出顺序将相邻的点用直线连接, (第一个和最后一个连接),所有线段不和其他线段相交,但是可以重合,可得到一个多边形或一条线段或一个多边形和一个线段的连接后的图形; 
如果是一条线段,我们定义面积为0,重心坐标为(0,0).现在求给出的点集组成的图形的面积和重心横纵坐标的和;
输入
第一行有一个整数0<n<11,表示有n组数据;
每组数据第一行有一个整数m<10000,表示有这个多边形有m个顶点;
输出
输出每个多边形的面积、重心横纵坐标的和,小数点后保留三位;
样例输入
3
3
0 1
0 2
0 3
3
1 1
0 0
0 1
4
1 1
0 0
0 0.5
0 1
样例输出
0.000 0.000
0.500 1.000
0.500 1.000
上传者

iphxer


思路:

首先得知道两个定理。

已知三角形各顶点A(Xi,Yi)( i=1,2,3)
S=((x2-x1)*(y3-y1)-(x3-x1)*(y2-y1));
重心坐标(Xg,Yg)为
Xg=(x1+x2+x3)/3;
Yg=(y1+y2+y3)/3;
有两种方法,我用的第二种。

算法一:一般适合凸多边形


      n边多边形可以分成n-2个三角形,将这些三角形看做质点(质点的位置是三角形的重心x1,x2,..,质量是面积s1,s2,..),
      那么多边形就由这些质点组成,质点坐标以其质量为权的加权算术平均数即是多边形重心坐标x。


x=(x1*s1+x2*s2+...)/(s1+s2+...)
y=(y1*s1+y2*s2+...)/(s1+s2+...)

s=s1+s2+...






算法二:任意多边形


       将算法一改进,n边多边形中每两个点(有顺序)加上原点可构成n个三角形,
       将这些三角形看做质点(质点的位置是三角形的重心x1,x2...,质量是面积(有正负)s1,s2,...),
       那么多边形就由这些质点组成,质点坐标以其质量为权的加权算术平均数即是多边形重心坐标x。
       多边形的面积s是这n个三角形面积(有正负)的代数和的绝对值。


x=(x1*s1+x2*s2+...)/(s1+s2+...)
x=(x1*s1+x2*s2+...)/(s1+s2+...)

s=|s1+s2+...|

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

double cross(double a[2],double b[2]){   //求向量a,b的叉积大小
	return a[0]*b[1]-a[1]*b[0];
}

int main(void)
{
	int ncase;
	cin>>ncase;
	while(ncase--)
        {
		int i,k;
		double S,tS,a[10002][2],sx,sy;
		cin>>k;
		for(i=1;i<=k;i++)
			cin>>a[i][0]>>a[i][1];
		S=0.;sx=0.;sy=0.;          //S面积,xy横纵坐标和
		for(i=1;i<=k;i++)
		{
			tS=cross(a[i],a[i%k+1])/2.;//最后一个与第一个相连 取余
			S+=tS;
			sx+=tS*(a[i][0]+a[i%k+1][0])/3;//s*x1
			sy+=tS*(a[i][1]+a[i%k+1][1])/3;
		}
		if(fabs(S)<1e-7)//精度判断
			puts("0.000 0.000");
		else
			printf("%.3lf %.3lf\n",fabs(S),(sx+sy)/S);
	}
	return 0;
}





  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值