UVA 11817 Tunnelling the Earth

最新推荐文章于 2018-08-29 20:20:17 发布

zxy_snow

最新推荐文章于 2018-08-29 20:20:17 发布

阅读量1.1k

点赞数

CC 4.0 BY-SA版权

分类专栏： URAL 计算几何文章标签： system

本文链接：https://blog.csdn.net/zxy_snow/article/details/6704225

计算几何同时被 2 个专栏收录

134 篇文章

订阅专栏

URAL

41 篇文章

订阅专栏

本文探讨了在地球表面计算两点间距离的方法，包括球面距离和直线距离的计算公式，并通过实例展示了算法的应用。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

经度维度计算球面距离以及直线距离。模板模板。。

#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 pi = acos(-1.0);
const double R = 6371009;
double angle_3d(double lng1, double lat1, double lng2, double lat2)
{
    double dlng = fabs(lng1 - lng2)*pi/180;
    while( dlng >= 2*pi )
        dlng -= pi*2;
    if( dlng > pi )
        dlng = 2*pi - dlng;
    lat1 *= pi/180;
    lat2 *= pi/180;
    return acos(cos(lat1)*cos(lat2)*cos(dlng) + sin(lat1)*sin(lat2));
}

int main()
{
    int ncases;
    double la,lb,wa,wb;
    
    scanf("%d", &ncases);
    
    while( ncases-- )
    {
        scanf("%lf%lf%lf%lf",&la, &wa, &lb, &wb);
        double ang = angle_3d(wa, la, wb, lb);
        double l1 = ang*R;
        double l2 = R*sqrt(2 - 2*cos(ang));
        printf("%.0lf\n", l1 - l2);
    }

return 0;
}