求纬度线上任意两点的球面距离 (10分)

题目描述

已知地球的平均半径为6371.393千米，并已知位于同一纬度上两点的纬度、经度值（单位为度），编写程序计算这两点之间的弦长和球面距离分别是多少千米。要求所有的实数使用double类型，圆周率的近似值取3.141592653589793。 提示：可以用正负来区分北纬与南纬、东经与西经。即北纬为正，南纬为负；东经为正，西经为负。

输入格式:

在第一行中输入一个位于区间[-90,90]内代表纬度的实数；在第二行中输入两个位于区间[-180,180]内代表经度的实数，以空格分隔。

输出格式:

输出弦长和球面距离的值，保留3位小数，以空格分隔，最后换行。

输入样例:

32.6
-56.1 163.8

输出样例:

10090.980 11645.760

源代码

#include <stdio.h>
#include <math.h>
#include <string.h>
#define PI 3.141592653589793
#define R 6371.393
int main()
{
    double latitude,longitude1,longitude2,l;
    scanf("%lf %lf %lf",&latitude,&longitude1,&longitude2);
    l=fabs(longitude2-longitude1);
    if(l>180)
        l=360-l;
    double arclength,length;
    latitude=latitude*2*PI/360.0;
    l=fabs(longitude2-longitude1)*2*PI/360.0;
    arclength=R*acos(cos(latitude)*cos(latitude)*cos(fabs(l))+sin(latitude)*sin(latitude));
    length=2*R*sin(arclength/(2*R));
    printf("%.3lf %.3lf\n",length,arclength);
    return 0;
}
//弧长为  arclength=R*arccos[cosb1*cosb2*cos(a1-a2)+sinb1*sinb2]   R是地球半径   a是经度  b是纬度
//弦长公式  弦长和弧长的关系式   C = arcsin(L/2R)×2r  弧度制公式   弦长为L，弧长为C，半径长为R