poj——1031——fence

最新推荐文章于 2019-09-21 00:52:51 发布

原创最新推荐文章于 2019-09-21 00:52:51 发布 · 1.3k 阅读

1 ·

CC 4.0 BY-SA版权

AC→poj 同时被 2 个专栏收录

11 篇文章

订阅专栏

计算几何

7 篇文章

订阅专栏

本文提供了一道来自POJ平台的问题1031的解答思路及代码实现，通过计算几何的方法解决了关于多边形面积的问题。

address：http://poj.org/problem?id=1031

My thought：说实话看见英语题的时候状态不咋好，后来勉强把题目看完。

不是太懂，想了好久，题目根本通不了，这英语太菜让我无语，我去度了一下，有的解法是用积分，有的说这是一道很水的题目。

我大概都看了看，有一些代码提示错误，下面这个代码是可以用的。题目中有一部分不是太懂，如果哪位朋友看懂了，可以留下言，谢谢。

#include<iostream>
#include<cmath>
#include<iomanip>
#include<algorithm>
using namespace std;
#define max(a,b) ((a)>(b)?(a):(b))   
#define min(a,b) ((a)<(b)?(a):(b))

const double PI=3.1415926535;
struct Point{
	double x,y;
};
Point point[205];
double beta(int k){
	//double result=(point[k].x)*(point[k+1].y)-(point[k+1].x)*(point[k].y);	
	double a=sqrt(point[k].x*point[k].x+point[k].y*point[k].y);	//点到原点的长度
	double deg1=acos(point[k].x/a);//与X轴的夹角
	if(point[k].y<0)
		deg1=2*PI-deg1;//求角的时候求出的是与X轴正半轴的夹角，如果在X轴下面，求出的角大小就是2*pi-这个角
	double b=sqrt(point[k+1].x*point[k+1].x+point[k+1].y*point[k+1].y);
	double deg2=acos(point[k+1].x/b);
	if(point[k+1].y<0)
		deg2=2*PI-deg2;
	double ans=deg2-deg1;
	if(ans>=PI)
		ans-=2*PI;
	else if(ans<=-PI)
		ans+=2*PI;
	return ans;
}

int main(){
	double k,h;
	int n,i;
	double angle,maximum,minimum,answer;
	while(cin>>k>>h>>n)
	{
		for(i=0;i<n;i++)
			cin>>point[i].x>>point[i].y;
		angle=0;
		maximum=0;
		minimum=0;
		point[n]=point[0];
		for(i=0;i<n;i++){
			angle+=beta(i);
			maximum=max(maximum,angle);
			minimum=min(minimum,angle);
		}
		
		angle=max(angle,-angle);
		if(fabs(angle-2*PI)<1e-4)
			answer=2*PI;
		else{
			answer=maximum-minimum;
			if(answer>2*PI)
				answer=2*PI;
		};
		answer*=k*h;	
		cout<<fixed<<setprecision(2)<<answer<<endl;
		//system("pause");
	}

	return 0;
}