POJ 1654 Area（任意多边形面积）

转载请注明出处，谢谢 http://blog.csdn.net/ACM_cxlove?viewmode=contents           by---cxlove

题目：从一个点出发，8个方向，给出每一步的方向，求出走过的路径形成的多边形的面积。不过题目说了会首尾相接，而且不会有交叉，不过显然不一定是凸多边形。

http://poj.org/problem?id=1654

其实利用向量叉积求解多边形面积，也不只是适用于凸边形。出现凹的部分，刚好一正一负，抵消了。

不过这题讨厌的是精度。

double是不能过的，必须要用整数，又不能溢出，用64位整数，输出的时候面积只需要除以2，所以只需要判断奇偶。。。。。啊啊啊WA好多次不只因为精度，八个方向对应的坐标变化错了好几次，东西不分啊 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#include<string>
#include<algorithm>
#include<queue>
#define LL __int64
#define eps 1e-7
#define N 2000000
#define MOD 1000000007
#define inf 1<<30
#define zero(a) (fabs((double)(a))<eps)
using namespace std;
struct Point{
    int x,y;
}p[1000005];
int n;
int way[9][2]={1,-1,1,0,1,1,0,-1,0,0,0,1,-1,-1,-1,0,-1,1};
LL xmul(Point p0,Point p1,Point p2){
    return (LL)(p1.x-p0.x)*(p2.y-p0.y)-(LL)(p2.x-p0.x)*(p1.y-p0.y);
}
LL area(){
    LL ans=0;
    for(int i=1;i<=n;i++)
        ans+=xmul(p[0],p[i],p[i+1]);
    return ans<0?-ans:ans;
}
char str[1000005];
int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%s",str);
        p[0].x=0;p[0].y=0;
        for(n=0;str[n+1];n++){
            int dir=str[n]-'1';
            p[n+1].x=p[n].x+way[dir][0];
            p[n+1].y=p[n].y+way[dir][1];
        }
        LL ans=area();
        if(ans%2==0) printf("%I64d\n",ans/2);
        else printf("%I64d.5\n",ans/2);
    }
    return 0;
}