多边形面积及重心求解算法 -- 计算几何

最新推荐文章于 2021-05-06 18:51:08 发布

Spikeeee-

最新推荐文章于 2021-05-06 18:51:08 发布

阅读量211

点赞数 1

分类专栏： # 计算几何 ACM

本文链接：https://blog.csdn.net/weixin_45919985/article/details/114092525

版权

ACM 同时被 2 个专栏收录

124 篇文章 2 订阅

订阅专栏

计算几何

4 篇文章 0 订阅

订阅专栏

叉乘应用：计算三角形面积
向量p=(a,b), q=(c,d)
p × q = ad － bc = - q × p
叉乘的大小等于于2倍三角形面积

测试样例：
在这里插入图片描述

#include <bits/stdc++.h>
#define ll long long
#define _for(i, a, b) for (int i = (a); i < (b); ++i)
#define sc scanf
#define pr printf
#define TLE                      \
    ios::sync_with_stdio(false); \
    cin.tie(0);
const int maxn = 1e5 + 5;
const int INF = 0x3f3f3f3f;
using namespace std;

int n;
double gx = 0.0, gy = 0.0;
double area = 0.0, _area =0.0;
struct Node {
	double x,y;
}c[maxn];

int main() {
	TLE;
	cin >> n;
	_for(i,0,n) cin >> c[i].x >> c[i].y;
	c[n].x=c[0].x;
	c[n].y=c[0].y;
	_for(i,0,n) {
		_area = (c[i].x*c[i+1].y - c[i].y*c[i+1].x)/2.0;
		area += _area;
		gx += _area * (c[i].x + c[i+1].x)/3.0;
		gy += _area * (c[i].y + c[i+1].y)/3.0;
	}
	if(area == 0.0) gx = gy = 0.0;
	else {
		gx = fabs(gx/area);
		gy = fabs(gy/area);
	}
	cout << "gx = " << gx << endl;
	cout << "gy = " << gy << endl;
	return 0;
}