球面距离计算公式：d(x1,y1,x2,y2)=r*arccos(sin(x1)*sin(x2)+cos(x1)*cos(x2)*cos(y1-y2))

x1,y1是纬度\经度的弧度单位，r为地球半径

 

把有用的信息从恶心的字符串中提取出来，然后套公式，再处理下恶心的精度，得出的距离需要四舍五入。

1 #include 
       
         2 #include
        
          3 #include
         
           4 #include
          
            5 #include
           
             6 #include
            
              7 #include
             
               8 #include
              
                9 #include
               
                10 using namespace std;11 #define N 10000012 #define LL long long13 #define INF 0xfffffff14 const double eps = 1e-8;15 const double pi = acos(-1.0);16 const double inf = ~0u>>2;17 const double r = 6875/2.0;18 char s[110];19 int dcmp(double x)20 {21 if(fabs(x)