任意多边形面积

最新推荐文章于 2021-04-08 17:59:49 发布

classmonster

最新推荐文章于 2021-04-08 17:59:49 发布

阅读量449

点赞数

分类专栏： c&&c++

c&&c++ 专栏收录该内容

50 篇文章 0 订阅

订阅专栏

给定多边形的顶点坐标（有序），让你来求这个多边形的面积，你会怎么做？
我们知道，任意多边形都可以分割为N个三角形，所以，如果以这为突破点，那么我们第一步就是把给定的多边形，分割为数个三角形，分别求面积，最后累加就可以了，把多边形分割为三角形的方式多种多样，在这里，我们按照如下图的方法分割：

图1

S点作为起始点（点1），a->e依次作为点2,3……。
一个三角形的面积是怎样的呢？
根据线性代数的知识，我们有如下的三角形面积公式，称之为有向面积(signed area)：

将这个行列式以第三列展开可以得到：

这就是以点1、2、3构成的三角形的有向面积（点如果是顺时针给出，有向面积为负，逆时针给出，有向面积为正），那么继续我们的工作，通过三角形的面积公式，来得到多边形的面积公式：
对于图1而言，多边形的面积就是：
S(1->6)=S(1,2,3)+S(1,3,4)+S(1,4,5)+S(1,5,6)
这里我们不免有些疑问，第一，图1所给出的是凸多边形，那这种算法对于非凸多边形是否同样适用呢？比如下面这个最简单的凸多边形的图形：

图2

用刚才的划分方法的话，就会出现一个诡异的问题，那就是有一个三角形出现在了图形的外面，而另外一个又超出了多边形的范围（划分为了Sab，Sbc两个图形），那么这样再用刚才的公式求面积，结果还是正确的么？
S(1->4)=S(1,2,3)+S(1,3,4)
先公布结论，这个式子是正确的，等等，为什么？还记得刚才我提到了那个“有向面积”的概念么？忘了的话，请回头看看加重了的字。
请注意从图中看，Sab点为顺时针排列，Sbc点为逆时针排列，面积从数值上就是从Sab这个超过范围的大三角形中去掉Sbc这个小三角形，最后的结果神奇的就是多边形Sabc的面积，那么这个结论能否推广到任意多边形呢？

图3

在这里不做证明，下面给出的公式，就是任意多边形的面积公式：

有了以上公式我们就可以计算任意多边形的面积了；

下面看一道题。

Description

“ 改革春风吹满地,
不会AC没关系;
实在不行回老家，
还有一亩三分地。
谢谢!（乐队奏乐）”

话说部分学生心态极好，每天就知道游戏，这次考试如此简单的题目，也是云里雾里，而且，还竟然来这么几句打油诗。
好呀，老师的责任就是帮你解决问题，既然想种田，那就分你一块。
这块田位于浙江省温州市苍南县灵溪镇林家铺子村，多边形形状的一块地，原本是linle 的，现在就准备送给你了。不过，任何事情都没有那么简单，你必须首先告诉我这块地到底有多少面积，如果回答正确才能真正得到这块地。
发愁了吧？就是要让你知道，种地也是需要AC知识的！以后还是好好练吧...

Input

输入数据包含多个测试实例，每个测试实例占一行，每行的开始是一个整数n(3<=n<=100)，它表示多边形的边数（当然也是顶点数），然后是按照逆时针顺序给出的n个顶点的坐标（x1, y1, x2, y2... xn, yn）,为了简化问题，这里的所有坐标都用整数表示。
输入数据中所有的整数都在32位整数范围内，n=0表示数据的结束，不做处理。

Output

对于每个测试实例，请输出对应的多边形面积，结果精确到小数点后一位小数。
每个实例的输出占一行。

Sample Input

    
    
     
     3 0 0 1 0 0 1
4 1 0 0 1 -1 0 0 -1
0

Sample Output

代码如下：

#include<iostream>
#include<algorithm>
#include<iomanip>
using namespace std;
int x[105],y[105];
int s(double x1,double y1,double x2,double y2)
{
	int s=x1*y2-x2*y1;
	return s;
}
int main()
{
	int a;
	while(cin>>a)
	{
		double z=0;
		if(a==0) break;
		
		for(int i=0;i<a;i++)
		{
			cin>>x[i]>>y[i];
		}
		x[a]=x[0];
		y[a]=y[0];
		for(int i=0;i<a;i++)
		{
			z=z+s(x[i],y[i],x[i+1],y[i+1]);
		
		}
		cout<<fixed<<setprecision(1)<<z/2<<endl;
	}
	
	return 0;
}