叉积求多边形面积
题意:以二维坐标的原点为起点,每次只走一步,每步只走根号2的长度或1的长度,且每步走的方向八个:上/下/左/右/左上/左下/右上/右下,而且每个方向用一个数字代替,1 2 3 4 6 7 8 9,5代表终止,即最后一步回到了原点,构成了多边形。现在给定一个以数字组成的字符串,字符串的每一个数字字符代表了走的方向,问走的路线构成的多边形的面积是多少?
算法:这题的字符串长度可以达到10^6个,而且多边形不分凹凸,这样的话,就很难求它的面积。但可以用“叉乘”的特性来解些题。如|a×b|=|a|*|b|*sinθ,若a和b为三解形两边的向量,那么这个三角形的面积就是S=|a×b|/2;利用向量叉乘的性质我们可以求出n个边构成的多边形的面积,但这些点的顺序必须是确定的。
在计算多边形面积时,只需要在平面中任意取一点,再取多边形中相邻两点,把所有的计
算结果相加,得出答案的绝对值就是多边形的面积,下面举几个例子:
例a : S = Soab + Sobc + Socd + Soda ,4个三角形都是顺时针,面积为他们和的绝对值。
例b : S =Soab+Sobc+Socd+Soda ,4个三角形中后3个是顺时针,后3个的绝对值的面积Sobcda,但是Soab为逆时针,符号与后三个相反,4个相加刚好把多余的抵消。例c,d:同理。
所以中心点可以任意选,对于些题,我们设中心点为(0,0),而且题中的输入的字符串的每个数字代表一个点(5除外),而且是有序的,这就符合要求了。
注意:1.字符串用数组,而且这些点也用数组存的话,此题必然会超内存。所以,最好的方法是过程中取点,过程中输入字符,这样只需要两个变量就可以解决问题了。
2.还有一点得注意,对于此题中多边形的面积,必定是0.5的整数倍,因为它由多个1*1的三角形组成,而在求sum时,sum=2*s,所以sum必定是一个整数,所以,我们可以用__int64类型的变量来存结果,如果用double可能存在精度丢失的情况,一般会wa。最后我们只需要判断sum%2的值就能确定输出是整数还是要加0.5了。
代码:
#include<stdio.h>
int a[10][2]={{0,0},{-1,-1},{0,-1},{1,-1},{-1,0},{0,0},{1,0},{-1,1},{0,1},{1,1}};//方向向量,想起了广搜的方向向量。。。
int main(){
__int64 area;
int i,n;
char ch;
int x1,y1,x2,y2;
scanf("%d",&n);
getchar();
while(n--){
area=0;x1=0;y1=0;
while((ch=getchar())!='\0'){
if(ch=='5') break;
x2=x1+a[ch-'0'][0];
y2=y1+a[ch-'0'][1];
area+=x1*y2-x2*y1;
x1=x2;
y1=y2;
}
if(area<0)
area=-area;
printf(area%2==0?"%I64d\n":"%I64d.5\n",area/2);
}
}