《CSP认证真题》2020-9-4 4.星际旅行

最新推荐文章于 2022-09-16 17:15:13 发布

Coder_小庞

最新推荐文章于 2022-09-16 17:15:13 发布

阅读量285

点赞数 1

分类专栏：算法 CSP 文章标签：程序设计

本文链接：https://blog.csdn.net/qq2297329956/article/details/115431531

版权

算法同时被 2 个专栏收录

22 篇文章 1 订阅

订阅专栏

CSP

4 篇文章 0 订阅

订阅专栏

该博客主要介绍了CSP认证中的一道题目，涉及到星际旅行中的两个星球之间的最短路径计算。文章详细阐述了如何利用距离公式和几何原理，计算两星球之间的直线距离以及考虑行星半径后的实际路径，并给出了解题的C++代码实现。

摘要由CSDN通过智能技术生成

《CSP认证真题》2020-9-4 4.星际旅行

文章目录

《CSP认证真题》2020-9-4 4.星际旅行

在这里插入图片描述

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>

using namespace std;

const int N = 110, M = 2010;

int n, m;
double R;
double o[N], p[M][N];
double d[M], rd[M];
double ans[M];

inline double sqr(double x)
{
    return x * x;
}

int main()
{
    scanf("%d%d%lf", &n, &m, &R);
    for (int i = 0; i < n; i ++ ) scanf("%lf", &o[i]);
    for (int i = 0; i < m; i ++ )
    {
        double s = 0;
        for (int j = 0; j < n; j ++ )
        {
            scanf("%lf", &p[i][j]);
            s += sqr(p[i][j] - o[j]);
        }
        d[i] = sqrt(s);
        rd[i] = sqrt(s - sqr(R));
    }

    for (int i = 0; i < m; i ++ )
        for (int j = 0; j < i; j ++ )
        {
            double s = 0;
            for (int k = 0; k < n; k ++ ) s += sqr(p[i][k] - p[j][k]);
            double c = sqrt(s), a = d[i], b = d[j];
            double p = (a + b + c) / 2;
            double area = sqrt(p * (p - a) * (p - b) * (p - c));
            double h = area * 2 / c;
            if (h >= R || sqr(b) >= sqr(a) + s || sqr(a) >= sqr(b) + s)
            {
                ans[i] += c, ans[j] += c;
                continue;
            }
            double angle1 = acos((sqr(a) + sqr(b) - s) / (2 * a * b));
            double angle2 = acos(R / a);
            double angle3 = acos(R / b);
            double t = (angle1 - angle2 - angle3) * R + rd[i] + rd[j];
            ans[i] += t, ans[j] += t;
        }
    for (int i = 0; i < m; i ++ )
        printf("%.12lf\n", ans[i]);
    return 0;
}
*/