题目:泰坦尼克(⊙v⊙);通过无线电测量冰山的位置,给出船和冰山的精度纬度,求距离。
分析:计算几何、大地坐标系。利用公式计算即可;圆心角:acos(cos(la1)*cos(la2)*cos(lo1-lo2)+sin(la1)*sin(la2)。
注意:1.输入格式好恶心;2.用c++救AC了g++无限WA。。。
#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cmath>
using namespace std;
double dist( double l1, double d1, double l2, double d2 )
{
double r = 6875.0/2;
double p = acos(-1.0);
l1 *= p/180;d1 *= p/180;
l2 *= p/180;d2 *= p/180;
return r*acos(cos(l1)*cos(l2)*cos(d1-d2)+sin(l1)*sin(l2));
}
int main()
{
char temp[100];
double d,m,s,l1,l2,d1,d2,dis;
for ( int i = 0 ; i < 9 ; ++ i ) scanf("%s",temp);
scanf("%lf^%lf'%lf\" %s",&d,&m,&s,temp);
l1 = d+m/60+s/3600;
if ( temp[0] == 'S' ) l1 *= -1;
scanf("%s",temp);
scanf("%lf^%lf'%lf\" %s.",&d,&m,&s,temp);
d1 = d+m/60+s/3600;
if ( temp[0] == 'W' ) d1 *= -1;
for ( int i = 0 ; i < 5 ; ++ i ) scanf("%s",temp);
scanf("%lf^%lf'%lf\" %s",&d,&m,&s,temp);
l2 = d+m/60+s/3600;
if ( temp[0] == 'S' ) l2 *= -1;
scanf("%s",temp);
scanf("%lf^%lf'%lf\" %s.",&d,&m,&s,temp);
d2 = d+m/60+s/3600;
if ( temp[0] == 'W' ) d2 *= -1;
scanf("%s",temp);
dis = dist(l1,d1,l2,d2);
printf("The distance to the iceberg: %.2lf miles.\n",dis);
if ( floor(dis+0.005) < 100 ) printf("DANGER!\n");
return 0;
}