[HDU 2036]改革春风吹满地

Posted 2020-09-30 NaVi_Awson

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

题解

计算几何裸题。叉积裸题。

将多边形拆成三角形，叉积计算平行四边形面积，答案累计，最后/2就好了。

注意题目给出的点是按逆时针顺序，注意有向面积的正负。

矢量叉积：

　　计算矢量叉积是与直线和线段相关算法的核心部分。设矢量P = ( x1, y1 )，Q = ( x2, y2 )，则矢量叉积定义为由(0,0)、p1、p2和p1+p2所组成的平行四边形的带符号的面积，即：P × Q = x1*y2 - x2*y1，其结果是一个标量。显然有性质 P × Q = - ( Q × P ) 和 P × ( - Q ) = - ( P × Q )。一般在不加说明的情况下，本文下述算法中所有的点都看作矢量，两点的加减法就是矢量相加减，而点的乘法则看作矢量叉积。

　　叉积的一个非常重要性质是可以通过它的符号判断两矢量相互之间的顺逆时针关系：

• 若 P × Q > 0 , 则P在Q的顺时针方向。
• 若 P × Q < 0 , 则P在Q的逆时针方向。
• 若 P × Q = 0 , 则P与Q共线，但可能同向也可能反向。

 

 1 #include<set>
 2 #include<map>
 3 #include<ctime>
 4 #include<cmath>
 5 #include<queue>
 6 #include<stack>
 7 #include<cstdio>
 8 #include<string>
 9 #include<vector>
10 #include<cstring>
11 #include<cstdlib>
12 #include<iostream>
13 #include<algorithm>
14 #define LL long long
15 #define RE register
16 #define IL inline
17 using namespace std;
18 
19 struct Point
20 {
21     double x,y;
22     Point (){}
23     Point (double _x,double _y){x=_x;y=_y;}
24     Point operator - (const Point &B)
25     const{
26         return Point(x-B.x,y-B.y);
27     }
28     int operator * (const Point &B)
29     const{
30         return x*B.y-B.x*y;
31     }
32 }a,b,c;
33 int n;
34 
35 int main()
36 {
37     scanf("%d",&n);
38     while (n)
39     {
40         double ans;
41         scanf("%lf%lf%lf%lf%lf%lf",&a.x,&a.y,&b.x,&b.y,&c.x,&c.y);
42         ans=(a-c)*(b-c);
43         for (RE int i=4;i<=n;i++)
44         {
45             b=c;
46             scanf("%lf%lf",&c.x,&c.y);
47             ans+=(a-c)*(b-c);
48         }
49         ans=ans/2;
50         printf("%.1lf\n",ans);
51         scanf("%d",&n);
52     }
53     return 0;
54 }