【HDU 2036】改革春风吹满地

问题描述
给出一n边形的顶点坐标(按逆时针方向),求此多边形面积。多组数据。
样例输入
3 0 0 1 0 0 1
4 1 0 0 1 -1 0 0 -1
0
样例输出
0.5
2.0
算法讨论
将多边形分解成多个三角形,利用叉积求出各个三角形面积,相加得答案。由于可能有凹多边形存在,所以不加绝对值。

const
  maxn=100;
var
  a:array[1..maxn,1..2] of longint;
  i,j,n:longint;
  s:real;

function cp(x1,y1,x2,y2,x,y:longint):real;
var
  ans:real;
begin
  ans:=0.5*((x1-x)*(y2-y)-(x2-x)*(y1-y));
  exit(ans)
end;

begin
  read(n);
  while n<>0 do
    begin
      fillchar(a,sizeof(a),0);
      s:=0;
      for i:=1 to n do
        read(a[i,1],a[i,2]);
      for i:=1 to n-2 do
        s:=s+cp(a[i+1,1],a[i+1,2],a[i+2,1],a[i+2,2],a[1,1],a[1,2]);
      writeln(s:0:1);
      read(n)
    end;
end.

【HDU 2036】改革春风吹满地_第1张图片
Pixiv ID:61813698

你可能感兴趣的:(Pascal,计算几何,叉积,计算几何)