题目概述
给定地球上两点的经纬度,求出两点间球面距离
设地球为直径6875的标准球体
时限
1000ms/3000ms
输入
输入共九行,第四行为点1的纬度(角度,分,秒)及南北纬,第五行为点1的经度及东西经,第七,八行为点2的,输入只有一组
限制
没有限制
输出
第一行为字符串
The distance to the iceberg: # miles.
其中#为保留两位小数的浮点数,为所求距离,若#<100,则在第二行输出
DANGER!
样例输入
Message #513.
Received at 22:30:11.
Current ship’s coordinates are
41^46’00” NL
and 50^14’00” WL.
An iceberg was noticed at
41^14’11” NL
and 51^09’00” WL.
===
样例输出
The distance to the iceberg: 52.04 miles.
DANGER!
讨论
计算几何,解析几何,给出经纬度算距离是有公式的,按照题解代码中的变量如下
dis=r⋅arccos(cos(y1)⋅cos(y2)⋅cos(x2−x1)+sin(y1)⋅sin(y2))
其中y是纬度,x是经度,r是地球半径
这个公式具体的推导过程可以参考
http://www.360doc.com/content/14/0117/10/325430_345890919.shtml
计算并不复杂,但是要写成一个公式确实需要不少化简
实现方面,这个题由于最后的DANGER!所以对精度有点奇怪的要求,先+0.005消除舍入误差,然后-100和0.0001比较,不这样处理的基本都会错
题解状态
212K,0MS,C++,1124B
题解代码
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define INF 0x3f3f3f3f
#define MAXN 103
#define memset0(a) memset(a,0,sizeof(a))
#define EPS 1e-8
char s[20];//string 存放没用的输入的字符串
const double r = 6875 / 2.0, pi = atan(1.0) * 4;//地球半径 题目给的是直径 pi 角度转弧度用
void fun()
{
int a, b, c;
scanf("%s%s%s%s%s%s%s%s%s%d^%d\'%d\"%s", s, s, s, s, s, s, s, s, s, &a, &b, &c, s);//input
double y1 = (a + b / 60.0 + c / 3600.0)*(s[0] == 'N' ? 1 : -1) / 180 * pi;//角度 分 秒的转换都是60进制
scanf("%s%d^%d'%d\"%s", s, &a, &b, &c, s);//input
double x1 = (a + b / 60.0 + c / 3600.0)*(s[0] == 'E' ? 1 : -1) / 180 * pi;
scanf("%s%s%s%s%s%d^%d'%d\"%s", s, s, s, s, s, &a, &b, &c, s);//input
double y2 = (a + b / 60.0 + c / 3600.0)*(s[0] == 'N' ? 1 : -1) / 180 * pi;
scanf("%s%d^%d'%d\"%s", s, &a, &b, &c, s);//input
double x2 = (a + b / 60.0 + c / 3600.0)*(s[0] == 'E' ? 1 : -1) / 180 * pi;
double dis = r*acos(cos(y1)*cos(y2)*cos(x2 - x1) + sin(y1)*sin(y2));
scanf("%s", s);//input
printf("The distance to the iceberg: %.2lf miles.\n%s", dis, dis + 0.005 - 100 < 0.0001 ? "DANGER!\n" : "");//output
}
int main(void)
{
//freopen("vs_cin.txt", "r", stdin);
//freopen("vs_cout.txt", "w", stdout);
fun();
}
EOF