hdoj2036(矢量求任意多边形的面积)
改革春风吹满地
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 26630 Accepted Submission(s): 13726
不会AC没关系;
实在不行回老家,
还有一亩三分地。
谢谢!(乐队奏乐)”
话说部分学生心态极好,每天就知道游戏,这次考试如此简单的题目,也是云里雾里,而且,还竟然来这么几句打油诗。
好呀,老师的责任就是帮你解决问题,既然想种田,那就分你一块。
这块田位于浙江省温州市苍南县灵溪镇林家铺子村,多边形形状的一块地,原本是linle 的,现在就准备送给你了。不过,任何事情都没有那么简单,你必须首先告诉我这块地到底有多少面积,如果回答正确才能真正得到这块地。
发愁了吧?就是要让你知道,种地也是需要AC知识的!以后还是好好练吧...
输入数据中所有的整数都在32位整数范围内,n=0表示数据的结束,不做处理。
每个实例的输出占一行。
3 0 0 1 0 0 1 4 1 0 0 1 -1 0 0 -1 0
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);
}
}