hdu 2036 改革春风吹满地 （叉乘法、海伦公式）

Problem 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
0.5
2.0

解析：理解题意，给出的多边形边数不确定，但是可以知道的一点是，n边形可以分解成n-2个两两相邻的三角形，题目要求多边形的面积，我们就只需要求出每个三角形面积，逐个面积相加。

那么问题来了，n给出的数据最大有100，这个时候是不是应该要有个方法有个公式什么的来重复每个三角形的面积计算就好了呢？答案是肯定的，并且还不止一个！！！

就通常的一种，海伦公式来说，在数学上有这么一句定义

三角形的面积S等于三角形周长C的一半分别与三边长a、b、c之差的乘积与周长C一半的乘积的开方

用公式来表示的话就是

p=(a+b+c)/2
S=sqrt(p*(p-a)*(p-b)*(p-c))

可是用海伦方法解决这个题目是行不通的，为什么？因为这个题目并没有说明多边形是凸多边形还是凹多边形，以上公式求的多边形面积只能用于凸多边形分解成若干个三角形

那么，接下来就要用到另一种方法来求了，它的名字叫叉乘法，详细说明如下：

#include<iostream>
#include<cstdio>
using namespace std;
int main()
{
    int n,i;
    int a[105],b[105];
    double s;
    while(cin>>n)
    {
        if(n==0) break;
        s=0;
        for(i=1;i<=n;i++)
          cin>>a[i]>>b[i];
        for(i=1;i<n;i++)
          s=s+(a[i]*b[i+1]-a[i+1]*b[i])*0.5;    //用叉乘法选择顺时针或者是逆时针先求前面n-3个三角形面积
        s=s+(a[i]*b[1]-a[1]*b[i])*0.5;   //最后一个n-2的三角形跟第一个三角形共一条边，单独求解
        printf("%.1lf\n",s);
    }
    return 0;
}

总结
叉乘法使用的时候简单方便，适用于任意多边形，凹凸都可
可是需要注意范围：

 多边形的顶点输入按顺时针或者逆时针排列

下面给出叉乘法求多边形模板

typedef struct
{
    double x,y;
} Point; 

double polygonarea(Point *polygon,int N)
{
    int i,j;
    double area = 0;
    for (i=0;i<N;i++)
    {
        j = (i + 1) % N;
        area += polygon[i].x * polygon[j].y;
        area -= polygon[i].y * polygon[j].x;
    }
    area /= 2;    
    return(area < 0 ? -area : area);
}