题意:
输入数据包含多个测试实例,每个测试实例占一行,每行的开始是一个整数n(3<=n<=100),它表示多边形的边数(当然也是顶点数),然后是按照逆时针顺序给出的n个顶点的坐标(x1, y1, x2, y2... xn, yn),为了简化问题,这里的所有坐标都用整数表示。输入数据中所有的整数都在32位整数范围内,n=0表示数据的结束,不做处理。对于每个测试实例,请输出对应的多边形面积,结果精确到小数点后一位小数。每个实例的输出占一行。
示例输入:
3 0 0 1 0 0 1
4 1 0 0 1 -1 0 0 -1
0
示例输出:
0.5
2.0
错误解决方案:
按逆时针给出若干个顶点的坐标,设有5个顶点,得出以下多边形:
从顶点0出发,有两条线将其分成了(0,1,2)(0,2,3)(0,3,4)三个三角形,只要求出这个三个三角形的面积和就是这个多边形的面积。已知三个顶点求这个这角形的面积用海伦公式:s=sqrt(p*(p-a)*(p-b)*(p-c)),p=(a+b+c)/2。
那个问题来了,要是下面这种点是凹进去的,怎么办?
那么看下正确的解决方案。
#include<stdio.h>
#include<math.h>
typedef struct point
{
double x;
double y;
}Point;
double Area(int a,int b,int c,Point p[100])
{
double aa,bb,cc,pp;
aa=sqrt((p[a].x-p[b].x)*(p[a].x-p[b].x)+(p[a].y-p[b].y)*(p[a].y-p[b].y));
bb=sqrt((p[b].x-p[c].x)*(p[b].x-p[c].x)+(p[b].y-p[c].y)*(p[b].y-p[c].y));
cc=sqrt((p[c].x-p[a].x)*(p[c].x-p[a].x)+(p[c].y-p[a].y)*(p[c].y-p[a].y));
pp=(aa+bb+cc)/2.0;
return sqrt(pp*(pp-aa)*(pp-bb)*(pp-cc));
}
int main()
{
Point p[100];
int n,i;
double areaSum;
while(scanf("%d",&n)&&n)
{
areaSum=0;
for(i=0;i<n;i++)
scanf("%lf%lf",&p[i].x,&p[i].y);
for(i=2;i<n;i++)
{
areaSum+=Area(0,i-1,i,p);
}
printf("%.1lf\n",areaSum);
}
}
正确解决方案:
#include<stdio.h>
#include<math.h>
typedef struct point
{
double x;
double y;
}Point;
int main()
{
Point p[100];
int n,i;
double areaSum;
while(scanf("%d",&n)&&n)
{
areaSum=0;
for(i=0;i<n;i++)
scanf("%lf%lf",&p[i].x,&p[i].y);
for(i=1;i<n;i++)
{
areaSum+=(p[i-1].x*p[i].y-p[i].x*p[i-1].y)/2.0;//交叉相乘
}
areaSum+=(p[n-1].x*p[0].y-p[0].x*p[n-1].y)/2.0;//首尾交叉相乘
printf("%.1lf\n",areaSum);
}
}