POJ 2354 大地坐标系

#include <cstdio>
#include <cstdlib>
#include <cmath>
const int maxn = 100;
char str[maxn];
double d, m, s, L1, L2, D1, D2;
double Dis(double l1, double d1, double l2, double d2)
{
    double r = 6875.0 / 2;
    double p = acos(-1.0);
    l1 *= p / 180.0; l2 *= p / 180.0;
    d1 *= p / 180.0; d2 *= p / 180.0;
    double d = r * sqrt(2 - 2 * (cos(l1) * cos(l2) * cos(d1 - d2) + sin(l1) * sin(l2)));
    return 2 * asin(d / (2 * r)) * r;
}
int main()
{
    while (~scanf("%*s%*s%*s%*s%*s%*s%*s%*s%*s"))
    {
        scanf("%lf^%lf'%lf\"%s", &d, &m, &s, &str);
        L1 = d + m / 60 + s / 3600; if (str[0] == 'S') L1 *= -1;
        scanf("%*s");
        scanf("%lf^%lf'%lf\"%s", &d, &m, &s, &str);
        D1 = d + m / 60 + s / 3600; if (str[0] == 'W') D1 *= -1;
        scanf("%*s%*s%*s%*s%*s");
        scanf("%lf^%lf'%lf\"%s", &d, &m, &s, &str);
        L2 = d + m / 60 + s / 3600; if (str[0] == 'S') L2 *= -1;
        scanf("%*s");
        scanf("%lf^%lf'%lf\"%s", &d, &m, &s, &str);
        D2 = d + m / 60 + s / 3600; if (str[0] == 'W') D2 *= -1;
        scanf("%*s");
        double ans = Dis(L1, D1, L2, D2);
        printf("The distance to the iceberg: %.2lf miles.\n", ans);
        if (floor(ans * 100 + 0.5) < 100 * 100) printf("DANGER!\n")                                              ;
    }
    return 0;
}

和uva 11817 基本相同。