hdu2036 改革春风吹满地

 

改革春风吹满地

 

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 30217    Accepted Submission(s): 15577

 

 

Problem 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

 

 

 

用海伦公式，需要考虑凸凹边形的情况，而且海伦公式计算量大，损失精度 

 

#include<cstdio>
#include<iostream>
#include<cmath>
using namespace std;

struct point  //点坐标 
{
	double x,y;
}p[101];

double dist(point p1,point p2) //求两点间的距离 
{
	return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
int judge(point p1,point p2,point p3)  //利用向量叉乘判断该分割的三角形面积是加还是减 
{
	double x1,y1,x2,y2;
	x1=p2.x-p1.x;
	y1=p2.y-p1.y;
	x2=p3.x-p2.x;
	y2=p3.y-p2.y;
	if(x1*y2-x2*y1>0) return 1;   //如果叉乘结果>0,即对应内角小于180度 
	else return -1;   //否则，对应内角大于180度 
}
int main()
{
	int n;
	while(cin>>n,n)
	{
		for(int i=1;i<=n;i++)
		  cin>>p[i].x>>p[i].y;
		double s=0;
		for(int i=2;i<=n-1;i++)
		{
			int temp=judge(p[1],p[i],p[i+1]);
			double a=dist(p[1],p[i]);
			double b=dist(p[1],p[i+1]);
			double c=dist(p[i],p[i+1]);
			double p=(a+b+c)/2;
			s+=sqrt(p*(p-a)*(p-b)*(p-c))*temp;//海伦公式求面积，根据判断结果决定是加还是减该面积 
		}
		printf("%.1lf\n",s);
	}
	return 0;
}

 

 

 

求多边形的面积最基本的还是把多边形分割成一个个三角形来计算。
计算三角形的面积可以用公式 S=(a*b*sin<a,b>)/2,其实就是三角形两条边向量的叉乘的 1/2
运用叉乘可以避免考虑凸凹边形的情况，也不必考虑顶点的输入时顺时针还是逆时针。 
叉乘的性质：设两向量P和Q 
1.P ×Q > 0 则 Q在 P的逆时针方向 
2.P ×Q < 0 则 Q在 P的顺时针方向
3.P ×Q = 0 则 Q和 P共线，方向可能相同也可能不相同 

 

<span style="font-size:14px;">#include<cstdio>
#include<iostream>
using namespace std;

struct point  //点坐标 
{
	double x,y;
}p[101];

double area(point p1,point p2,point p3) //求三角形面积
{
	return ((p2.x-p1.x)*(p3.y-p1.y)-(p3.x-p1.x)*(p2.y-p1.y))/2.0;
}
int main()
{
	int n;
	while(cin>>n,n)
	{
		for(int i=1;i<=n;i++)
		  cin>>p[i].x>>p[i].y;
		double s=0;
		for(int i=2;i<=n-1;i++)  //笔者这里以第一个顶点为分割的公共顶点
		   s+=area(p[1],p[i],p[i+1]);
		printf("%.1lf\n",s);
	}
	return 0;
} </span>