牛客练习赛17 D 经纬度(数学题)

点击打开链接

 

解析:

...不知道的公式，推出来的结果...，完全懵逼

 

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
typedef long long int ll;
const double pi=acos(-1.0);
const double R =6371009;
double rad(double x)
{
    return x*pi/180.0;
}
double Dis(double lon1,double lat1,double lon2,double lat2)  //算直线距离
{
    lat1=rad(lat1);
    lat2=rad(lat2);
    double a=lat1-lat2;
    lon1=rad(lon1);
    lon2=rad(lon2);
    //double b=lon1-lon2;
    //double s=2*asin(sqrt(pow(sin(a/2),2)+cos(lat1)*cos(lat2)*pow(sin(b/2),2)));
    double s = R*sqrt(2-2*(cos(lon1)*cos(lon2)*cos(a)+sin(lon1)*sin(lon2)));
    //s=R*s;
    s=round(s*10000)/10000;
    return s;
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        double a,b,x,y;
        scanf("%lf%lf%lf%lf",&a,&b,&x,&y);
        double ans1=Dis(a,b,x,y);
        double ans2=2*asin((ans1*0.5)/R)*R;  //球面距离
        printf("%.0lf\n",ans2-ans1);
    }
}