[hdu2036]改革春风吹满地

“ 改革春风吹满地,
不会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

一道计算几何裸题

每次选三个点叉一下基即可

因为三角形面积==叉基的二分之一的绝对值

证明：易证三角形面积==叉基的二分之一的绝对值

然后就做出来了

// It is made by XZZ
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
#define rep(a,b,c) for(rg int a=b;a<=c;a++)
#define drep(a,b,c) for(rg int a=b;a>=c;a--)
#define erep(a,b) for(rg int a=fir[b];a;a=nxt[a])
#define il inline
#define rg register
#define vd void
typedef long long ll;
typedef double db;
struct Point{
    db x,y;
    Point(){}
    Point(db _x,db _y){x=_x,y=_y;}
};
typedef Point P;
typedef Point Vector;
typedef Point V;
P operator + (P a,P b){return P(a.x+b.x,a.y+b.y);}
P operator - (P a,P b){return P(a.x-b.x,a.y-b.y);}
P operator * (P a,P b){return P(a.x*b.x,a.y*b.y);}
P operator * (P a,db b){return P(a.x*b,a.y*b);}
db operator ^ (P a,P b){return a.x*b.y-a.y*b.x;}
int main(){
    int n;db ans;
    P a,b,c;
    while(scanf("%d",&n),n){
        scanf("%lf%lf%lf%lf%lf%lf",&a.x,&a.y,&b.x,&b.y,&c.x,&c.y);
        ans=(a-b)^(a-c);
        rep(i,4,n)b=c,scanf("%lf%lf",&c.x,&c.y),ans+=(a-b)^(a-c);
        printf("%.1lf\n",fabs(ans/2.00));
    }
    return 0;
}

 更新：果然AC了