杭电2036

改革春风吹满地

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 48710    Accepted Submission(s): 24897

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

解题思路：

无论三角形的顶点位置如何，△PMN总可以用一个直角梯形（或矩形）和两个直角三角形面积的和差来表示
而在直角坐标系中，已知直角梯形和直角三角形的顶点的坐标，其面积是比较好求的。
下面以一种情形来说明这个方法，其它情形方法一样，表达式也一样（表达式最好加上绝对值，确保是正值）
如图情形（P在上方，M在左下，N在右下），过P作X轴的平行线L，作MA⊥L，NB⊥L（设P在A、B之间）
则A、B的坐标是A(c,b)，B(e,b)
所以PA＝a－c,PB＝e－a,AM＝b－d,BN＝b－f,AB＝e－c
所以S△PMN＝S梯形AMNB－S△PAM－S△PBN
＝(b－d＋b－f)(e－c)/2－(b－d)(a－c)/2－(b－f)(e－a)/2
＝(ad＋be＋cf－af－bc－de)/2

即：

三角形三顶点坐标分别为A(a,b),B(c,d),C(e,f)，那么这个三角形的面积为
S=1/2*三阶行列式，

三阶行列式为：

a  b  1

c  d  1

e  f  1

 

S=S123+S134+S145+S156+S167+S178

>在坐标系下，三角形面积公式如下。可按向量思路。

S=(double)(x2*y3+x1*y2+x3*y1-x3*y2-x2*y1-x1*y3)/2.0;

代码如下：

#include<iostream>
#include<cstdio>
using namespace std;
double do_area(int a,int b,int c,int d,int e,int f)
{
	return (a*d+c*f+b*e-e*d-c*b-a*f)/2.0;
}
 
int main()
{
	int i,n;
	struct xy
	{
		int x;
		int y;
	}a[105];
	double area;
	while(cin>>n)
	{
		area=0;
		if(n==0)  break;
		for(i=0;i<n;i++)
			cin>>a[i].x>>a[i].y;
		for(i=1;i<n-1;i++)
		{
			area+=do_area(a[0].x,a[0].y,a[i].x,a[i].y,a[i+1].x,a[i+1].y);
		}
		printf("%.1lf\n",area);
	}
	return 0;
}