eoj1127 计算几何 任意多边形面积

题目：eoj1127

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

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

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

0.5
2.0

 

 

 

题目分析：

计算任意多边形面积。工具：向量叉乘。向量叉乘为三角形有向面积的两倍。因此可选取任意点，如原点与每条边的两个端点形成三角形，通过有向面积的正负向消可以得出任意多边形面积，可做过作图直观理解这一定理。

 

AC代码：

<span style="font-size:18px;">#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
 
using namespace std;
struct Point{			//点类
    double x,y;
    Point(int xx=0,int yy=0):x(xx),y(yy){}
};
struct Vec{			//向量类，可合并向量类和点类。但这样区分开来更好理解
    double x,y;
    Vec(){}
    Vec(const Point& a,const Point& b){		//将两点转化成向量
        x=a.x-b.x;
        y=a.y-b.y;
    }
    friend double operator*(Vec a,Vec b){		//重载叉乘运算
        return a.x*b.y-a.y*b.x;
    }
}v[105];
 
int main()
{
    int n;
    Point o={0.0,0.0};
    while(~scanf("%d",&n),n){
        for(int i=0;i<n;++i){
            double a,b;
            scanf("%lf%lf",&a,&b);
            v[i]={Point{a,b},o};
        }
        double ans=0.0;
        for(int i=0;i<n-1;++i){
            ans+=v[i]*v[i+1];
        }
        ans+=v[n-1]*v[0];
        printf("%.1f\n",fabs(0.5*ans));
    }
    return 0;
}
</span>