zoj3598----球面三角形内角

最新推荐文章于 2018-02-27 09:59:52 发布

carter91620

最新推荐文章于 2018-02-27 09:59:52 发布

阅读量1.3k

点赞数 1

分类专栏：计算几何文章标签： struct c

本文链接：https://blog.csdn.net/dicpu/article/details/7442524

版权

计算几何专栏收录该内容

17 篇文章 0 订阅

订阅专栏

公式公式，公式的水题。。。。

#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define pi acos(-1.0)
#define eps 1e-8
struct TPoint
{
	double al,bk;
}p1,p2,p3;
double getangle(TPoint a,TPoint b)
{
	double al=fabs(a.al-b.al)*pi/180;
	while(al>=2*pi) al-=2*pi;
	if(al>=pi) al=2*pi-al;
	a.bk*=pi/180;
	b.bk*=pi/180;
	return cos(a.bk)*cos(b.bk)*cos(al)+sin(a.bk)*sin(b.bk);
}
int main()
{
	int re,ri;
	double ju,wu,r=10000;
	scanf("%d", &re);
	for(ri=1;ri<=re;ri++)
	{
		scanf("%lf%lf",&p1.al,&p1.bk);
		scanf("%lf%lf",&p2.al,&p2.bk);
		scanf("%lf%lf",&p3.al,&p3.bk);
		double ang1,ang2,ang3,l1,l2,l3,a,b,c;
		ang1=getangle(p1,p2);//cos
		ang2=getangle(p1,p3);
		ang3=getangle(p2,p3);
		a=sqrt(1-ang1*ang1); //sin
		b=sqrt(1-ang2*ang2);
		c=sqrt(1-ang3*ang3);

		double h1,h2,h3;
		h1=(ang1-ang2*ang3)/(b*c);
		h2=(ang2-ang1*ang3)/(a*c);
		h3=(ang3-ang1*ang2)/(a*b);

		double h;
		h=acos(h1)+acos(h2)+acos(h3);
		double w=h*180/pi;
		printf("%.2f\n",w);
	}
	return 0;
}