假设地球是一个标准的球体，其半径是R，忽略地形对距离的影响。球面上的点使用经纬度表示，规定东经为正、西经为负、北纬为正、南纬为负。设A点的经度是α1、纬度是β1；B点的经度是α2、纬度是β2，则A和B之间的距离可按下式近似计算：

同时约定：东经为正,西经为负；南纬为90°+地理纬度值、北纬为90°-地理纬度值
D=Rarccos(C)
其中：C=sin(β1)sin(β2)cos(α1-α2)+cos(β1)cos(β2)；
特别注意：arccosC的单位是弧度.

假设R=6371km，编写函数计算地球表面任意两点之间的距离。

两个点的经纬度

两点之间的距离

1. #include <stdio.h>
2. #include <stdlib.h>
3. #include <math.h>
4. #define pi 3.14159
5. int main()
6. {
7.     double R=6371.0,a1,b1,a2,b2,d,c;
8.     scanf("%lf%lf%lf%lf",&a1,&b1,&a2,&b2);
9.     a1=a1/180*pi;
10.     b1=b1/180*pi;
11.     a2=a2/180*pi;
12.     b2=b2/180*pi;
13.     c=cos(b1)*cos(b2)*cos(a1-a2)+sin(b1)*sin(b2);
14.     d=acos(c)*R;
15.     printf("%.6lf",d);
16.     return 0;
17. }