hdoj2036(矢量求任意多边形的面积)

hdoj2036(矢量求任意多边形的面积)

改革春风吹满地

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

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
  
  

思路：

把这个图形按照顺时针求解面积：

S=S(p2,p1,0)+S(p3,p2,0)+S(p4,p3,0)+S(p5,p4,0)

此时s=s1+s2;,s1为待求面积

S是顺时针求得的，在这个过程中s2也是按照顺时针求得的并计算在了S内.

即：S2=S(p7,p1,0)+S(p6,p7,0)+S(p5,p6,0)

在求k=S(p6,p5,0)+S(P7,p6,0)+S(p1,p7,0)(相当于是按逆时针求得的k=-S2)

最后可推的：任意多边形面积公式：

S1=S(p2,p1,0)+S(p3,p2,0)+S(p4,p3,0)+S(p5,p4,0)+S(p6,p5,0)+S(P7,p6,0)+S(p1,p7,0)

 for(i=1;i<n;i++)
k+= (map[i].a*map[i-1].b)-map[i-1].a*map[i].b;
k+= (map[0].a*map[n-1].b)-map[n-1].a*map[0].b;//最后一个点和第一个点以及原点构成的三角形的面积 
这是这题的易错点，算矢量积时一定只能选定一个方向进行。
(x1,y1)*(x0,y0)和(x0,y0)*(x1,y1)的结果是不一样的！！！！ 
如：
计算(x0,y0),(x1,y1),(0,0)的面积时 
k=(x1*y0)-(x0*y1)   （后一个点的坐标，在前一个坐标的前面，即:(x1,y1)*(x0,y0)不是(x0,y0)*(x1,y1)） 
下面都要保证这个顺序求解：
k=(x1*y0)-(x0*y1)   +   (x2*y1)-(x1*y2)  +   (x3*y2)-(x2*y3) +......
  +  (x0*x(n-1))-(x(n-1)*y0)  (最后一个点和第一个点的坐标进行向量运算时,第一个点一定在最后一个前面)

公式很简单：k=(x1*y0)-(x0*y1)   +   (x2*y1)-(x1*y2)  +   (x3*y2)-(x2*y3) +...... +  (x0*x(n-1))-(x(n-1)*y0)

                      s=k/2;

My   solution:/2016.1.15/

#include<stdio.h>
struct stu
{
	int a,b;
}map[110];
int n;
double area()
{
	int i,j;
	double k=0;
	for(i=1;i<n;i++)
	k+= (map[i].a*map[i-1].b)-map[i-1].a*map[i].b;
	k+= (map[0].a*map[n-1].b)-map[n-1].a*map[0].b;//这个容易忽略,最后一个点和第一个点以及原点构成的三角形的面积 
	k/=2;
	if(k<0)
	k=-k;
	return k;
}
int main()
 {
 	int i,j,k;
 	double s;
    while(scanf("%d",&n)&&n)
    {
    	for(i=0;i<n;i++)
    	scanf("%d%d",&map[i].a,&map[i].b);
    	s=area();
    	printf("%.1lf\n",s);
	}
 }