ACM计算几何知识点详解与模板代码1

1. 求解多边形面积

• 题目描述

给出一个简单多边形(没有缺口)，它的边要么是垂直的，要么是水平的。要求计算多边形的面积。
多边形被放置在一个X-Y的卡笛尔平面上，它所有的边都平行于两条坐标轴之一。然后按逆时针方向给出各顶点的坐标值。所有的坐标值都是整数(因此多边形的面积也为整数)

输入

第 一行给出多边形的顶点数n(n≤100)。接下来的n行每行给出多边形一个顶点的坐标值X和Y(都为整数并且用空格隔开)。顶点按逆时针方向逐个给出。并 且多边形的每一个顶点的坐标值-200≤x，y≤200。多边形最后是靠从最后一个顶点到第一个顶点画一条边来封闭的。

输出

仅有一行，包含一个整数，表示多边形的面积。

样例输入

10
0 0
4 0
4 1
3 1
3 3
2 3
2 2
1 2
1 3
0 3

样例输出

9

解决思路：

将点的问题转化为向量的问题是计算几何中的常见思路，首先我们会想到高数中如何求解平行四边形的面积？那就是用相邻的两边构成两个向量，然后进行叉乘运算即可，接着我们会想到如何求三角形面积，那就是平行四边形的一半。所以现在你应该想到如何去做多边形的面积了，那就是将其分为很多三角形。现在，问题来了，我们该如何分割呢？

其实这也非常简单，我们任取一点作为定点，分别做与这个点相连的向量，接着让相邻点的向量分别相乘即可。

具体证明过程可以直接观看上图，下面直接给出结论

#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
const double eps=1e-5;
struct vec//定义一个向量结构体，便于计算
{
    double x,y;
    vec(double xx=0,double yy=0)
    {
        x=xx;
        y=yy;
    }
};
double operator*(vec a,vec b)//按C++代码风格重载运算符
{
    return a.x*b.y-a.y*b.x;
}
void solveArea(vec a[],int n)
{
    int i,j;
    double area=0;
    for(i=0;i<n;i++)
    {
        area+=a[i]*a[(i+1)%n];
    }
    double k=0.f;
    k=abs(0.5*area);
    cout<<k<<endl;
}
int main()
{
    int n,i,j;
    double c,b;
    vec *a=new vec[100];
    while(cin>>n)
    {
        for(i=0;i<n;i++)
        {
            cin>>a[i].x>>a[i].y;
        }
         solveArea(a,n);
    }
}