uva 11817 - Tunnelling the Earth

题意：从地球上的一个点到另一个点，求两点的球面距离和直线距离之差。假定地球是正球体，半径为6371009米。

 

#include<iostream>
#include<cmath>
#define r 6371009
#define pi 2.0*asin(1.0)

using namespace std;

double ang(double lt1,double lt2,double lg1,double lg2)
{
    return acos(cos(lt1)*cos(lt2)*cos(lg1-lg2)+sin(lt1)*sin(lt2));
}
double dis_line(double lt1,double lt2,double lg1,double lg2)
{
    return r*sqrt(2-2*(cos(lt1)*cos(lt2)*cos(lg1-lg2)+sin(lt1)*sin(lt2)));
}
double dis_sphere(double lt1,double lt2,double lg1,double lg2)
{
    return r*ang(lt1,lt2,lg1,lg2);
}
int main()
{
    int t;
    double lt1,lt2,lg1,lg2;
    cin>>t;
    while(t--)
    {
        cin>>lt1>>lg1>>lt2>>lg2;
        lt1*=pi/180.0;
        lt2*=pi/180.0;
        lg1*=pi/180.0;
        lg2*=pi/180.0;
        long long l=0.5+dis_sphere(lt1,lt2,lg1,lg2)-dis_line(lt1,lt2,lg1,lg2);
        cout<<l<<endl;
    }
    return 0;
}