UVA 10897 Travelling Distance

还是经度维度求距离。。。水题。。纯粹熟练模板。。

#include <map>
#include <set>
#include <queue>
#include <stack>
#include <math.h>
#include <time.h>
#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <limits.h>
#include <string.h>
#include <string>
#include <algorithm>
#define MID(x,y) ( ( x + y ) >> 1 )
#define L(x) ( x << 1 )
#define R(x) ( x << 1 | 1 )
#define BUG puts("here!!!")
#define STOP system("pause")

using namespace std;

const double r = 6371.01;
const double pi = acos(-1.0);
void get(double &lat, double  &lng)
{
	double x, y, z;
	char s[10];
	scanf("%lf%lf%lf%s",&x, &y, &z, s);
	lat = x + y/60 + z/3600;
	if( s[0] == 'S' ) lat *= -1;
	lat *= pi/180;
	
	scanf("%lf%lf%lf%s",&x, &y, &z, s);
	lng = x + y/60 + z/3600;
	if( s[0] == 'W' ) lng *= -1;
	lng *= pi/180;
}
double angle_3d(double lng1, double lat1, double lng2, double lat2)
{           //经度，纬度，经度，纬度 
    return acos(cos(lat1)*cos(lat2)*cos(lng1 - lng2) + sin(lat1)*sin(lat2));
}

int main()
{
	int ncases;
	double lat1, lng1, lat2, lng2;
	
	scanf("%d", &ncases);
	
	while( ncases-- )
	{
		get(lat1, lng1);
		get(lat2, lng2);
		
		double ang = angle_3d(lng1, lat1, lng2, lat2);
		
		double len = ang*r;
		
		printf("%.2lf\n", len);
	}

return 0;
}