Codeforces Gym-101116-F (Flight Plan)

题目链接：

$\quad$Gym-101116

题目大意：

$\quad$假设地球是一个球体，给出起点的经纬度$A(x_1, y_1),$ 终点的经纬度$B(x_2, y_2)$。
$\quad x :$ 纬度，北纬为正，南纬为负；
$\quad y :$ 经度，东经为正，西经为负。
$\quad$令：
$\qquad D_1$为两点的距离；
$\qquad D_2$为起点$A$先沿着纬线，再沿着经线 走到终点$B$的距离。
$\quad 求D_1, D_2$。

数据范围：

$\quad 1 \leq T \leq 10000$
$\quad \lvert x_1 \rvert, \lvert x_2 \rvert \leq 90, \quad \lvert y_1 \rvert , \lvert y_2 \rvert \leq 180$

解题思路：

$\quad$这道题不重要，重要的是怎么求地球上两点的距离。

$\quad$图的话，将就看这个吧！懒得画了，，，只看图就好，，，
$\quad$

$\quad$球面坐标和直角坐标存在一一对应的关系：球面有一点$M(r, \varphi, \theta)$，
$\quad$
$\quad$则，点$M$的直角坐标和球面坐标的关系为：

$$\begin{cases} x = r \cos\varphi \cos\theta \\ y = r \cos\varphi \sin\theta \\ z = r \sin\varphi \end{cases}$$
$\quad 设 A(r, \alpha, \beta), B(r, \varphi, \theta)$，
$\quad 则有\vec{OA} = (r \cos\alpha \cos\beta, r \cos\alpha \sin\beta, r \sin\alpha), \vec{OB} = (r \cos\varphi \cos\theta, r \cos\varphi \sin\theta, r \sin\varphi)$
$\quad$由此可得，夹角$ang = \arccos(\frac{\vec{OA} \cdot \vec{OB}}{\lvert \vec{OA} \rvert \cdot \lvert \vec{OB} \rvert} )$

$\qquad$化简得：$\color{navy}{ang = \arccos(\cos\alpha \cos\varphi \cos(\beta - \theta) + \sin\alpha \sin\varphi)}$
$\qquad$夹角出来了，距离就好办了。

$\quad$代码：

//化为弧度制
double rad(double d) {
    return d * PI / 180.0;
}

//求两点之间的距离
double dis(Point A, Point B) {
    A.x = rad(A.x); A.y = rad(A.y);
    B.x = rad(B.x); B.y = rad(B.y);
    double ang = acos(cos(A.x) * cos(B.x) * cos(A.y - B.y) + sin(A.x) * sin(B.x));
    double s = ang * R;    //R为地球半径
    return s;
}

$\quad$我还看到一个公式，但我不知道怎么推出来的，若有好心大佬，求解答。

double dis(Point A, Point B) {
    double Theta = A.x - B.x;
    double Beta = A.y - B.y;
    double ang = 2 * asin(sqrt(sin(Theta / 2) * sin(Theta / 2) + cos(A.x) * cos(B.x) * sin(Beta / 2) * sin(Beta / 2)));
    double s = ang * R;
    return s;
}

$\quad$距离的问题解决了，其他都不重要了。不对，这道题还没完呢！，，，
$\quad$有一个坑点：当一个点在东半球，另一个在西半球时，需要特别注意一下角度，稍微想想。

$\quad$AC代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
#include <set>
#include <iostream>
using namespace std;
typedef long long LL;
const int inf = 1 << 30;
const LL INF = 1LL << 60;
const double PI = acos(-1.0);
const double R = 6371.0;

int T;
struct Point {
    double x, y;
};
Point a, b, c;

//化为弧度制
double rad(double d) {
    return d * PI / 180.0;
}

//求弧长
double arc(double Theta, double r) {
    return Theta * r;
}

//求两点之间的距离
double dis(Point A, Point B) {
    A.x = rad(A.x); A.y = rad(A.y);
    B.x = rad(B.x); B.y = rad(B.y);
    double ang = acos(cos(A.x) * cos(B.x) * cos(A.y - B.y) + sin(A.x) * sin(B.x));

    /*这个公式也可以,但我不知道是怎么推的
    double Theta = A.x - B.x;
    double Beta = A.y - B.y;
    double ang = 2 * asin(sqrt(sin(Theta / 2) * sin(Theta / 2) + cos(A.x) * cos(B.x) * sin(Beta / 2) * sin(Beta / 2)));
    */
    double s = ang * R;
    return s;
}

int main()
{
    scanf("%d", &T);
    while(T--) {
        scanf("%lf %lf %lf %lf", &a.x, &a.y, &b.x, &b.y);
        //这个稍微想想就明白了,一个小坑点
        if(abs(a.y - b.y) > 180.0) {
            if(a.y < 0.0) a.y += 360.0;
            else if(b.y < 0.0) b.y += 360.0;
        }
        c.x = a.x; c.y = b.y;
        double r = R * cos(rad(a.x));
        printf("%.10lf %.10lf\n", dis(a, b), dis(b, c) + arc(rad(abs(a.y - b.y)), r));
    }
    return 0;
}