poj 3407 Brookebond s'en va en guerre... && poj 2354 Titanic

最新推荐文章于 2020-03-06 15:47:33 发布

原创最新推荐文章于 2020-03-06 15:47:33 发布 · 562 阅读

0 ·

CC 4.0 BY-SA版权

学生时代刷题归档专栏收录该内容

146 篇文章

订阅专栏

本文探讨了一个看似法语的英语题目，实际涉及求解两点在球面上的距离，利用经纬度坐标进行计算。文章详细介绍了球面距离公式，并通过代码实现了解题过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意：这奇葩的题题目看着是不是很不像英语。。。没错。。这尼玛好像是法语。。。。orz。。，大致意思就是求两个点的球面距离，不过给的是经纬度坐标

思路:公式题。。套公式就好了

/*
球面距离公式
dis = r*acos(sin(a1)*sin(a2)+cos(a1)*cos(a2)*cos(b1-b2)) ;
 
其中a1,a2表示维度（北正南负），b1,b2表示经度（东正西负）。R为球的半径
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<algorithm>
using namespace std;
char s[100],s1[100];
const double pi=acos(-1.0);
double solve(double a,double c,char s[])
{
	a=a+c/60;
	if(s[0]=='S' || s[0]=='W')
		a=-a;
	a=a*pi/180;
	return a;
}
int main()
{
	double a,b,c,d;
	double p1x,p2x,p1y,p2y;
	while(~scanf("%lf%lf%s%lf%lf%s",&a,&b,s1,&c,&d,s))
	{
		p1x=solve(a,b,s1);
		p1y=solve(c,d,s);
		scanf("%lf%lf%s%lf%lf%s",&a,&b,s1,&c,&d,s);
		p2x=solve(a,b,s1);
		p2y=solve(c,d,s);
		double xx= acos(sin(p1x) * sin(p2x) + cos(p1x) * cos(p2x) * cos(p2y - p1y)) * 6370;
		printf("%.3lf\n",xx);
	}
	return 0;
}

poj 2354

尼玛的。。这题输入格式为什么这么坑。。

/*
球面距离公式
dis = r*acos(sin(a1)*sin(a2)+cos(a1)*cos(a2)*cos(b1-b2)) ;

其中a1,a2表示维度（北正南负），b1,b2表示经度（东正西负）。R为球的半径
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<algorithm>
using namespace std;
char s[1000],s1[1000];
const double pi=acos(-1.0);
const double eps=1e-8;
double solve(double a,double b,double c,char s[])
{
	a=a+c/3600+b/60;
	if(s[2]=='S' || s[2]=='W')
		a=-a;
	a=a*pi/180;
	return a;
}
int dblcmp(double k)
{
	if(fabs(k)<eps)
		return 0;
	return k>0?1:-1;
}
int main()
{
	double a,b,c;
	double a1,b1,a2,b2;
	while(gets(s))
	{
		gets(s);
		gets(s);
		scanf("%lf^%lf'%lf",&a,&b,&c);
		gets(s);
		a1=solve(a,b,c,s);
		scanf("and %lf^%lf'%lf",&a,&b,&c);
		gets(s);
		b1=solve(a,b,c,s);



		gets(s);
		scanf("%lf^%lf'%lf",&a,&b,&c);
		gets(s);
		a2=solve(a,b,c,s);
		scanf("and %lf^%lf'%lf",&a,&b,&c);
		gets(s);
		b2=solve(a,b,c,s);
		gets(s);
		double xx=6875.0/2;
		xx=xx*acos(sin(a1)*sin(a2)+cos(a1)*cos(a2)*cos(b1-b2));
		printf("The distance to the iceberg: %.2lf miles.\n",xx);
		if(dblcmp(xx-100+0.005)<0)
			printf("DANGER!\n");
	}
	return 0;
}