UVA1001 Say Chesse 解题报告

UVA1001 Say Chesse 解题报告

题目链接

https://vjudge.net/problem/UVA-1001

题目大意

无限大的三维立体奶酪里有n（0≤n≤100）个球形的洞。你的任务是帮助小老鼠A用最短的时间到达小老鼠O所在位置。奶酪里的移动速度为10秒一个单位，但是在洞里可以瞬间移动。洞和洞可以相交。输入n个球的位置和半径，以及A和O的坐标，求最短时间。

解题思路

考虑两个洞a和b之间的距离，假设他们的半径分别为r1和r2，因为洞内可以瞬移，所以a和b之间的距离为他们的欧几里得距离减去他们各自的半径，即sqrt((x[a] - x[b]) * (x[a] - x[b]) + (y[a] - y[b]) * (y[a] - y[b]) + (z[a] - z[b]) * (z[a] - z[b])) - r[a] - r[b]，这个值可能小于等于0，小于等于0则代表两个洞相交，则视作距离为0即可。至于起点和终点，我们可以将其视作半径为0的洞。我们将这n+2个点两两之间计算距离连边建图，最后跑一遍Floyd求最短路即可（跑Dijkstra也可以，但这题的数据量本身就非常小，Floyd代码更简短）。

代码

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using ull = unsigned long long;
using ld = long double;
#define endl '\n';
const int maxn = 1e2 + 10;
const ll INF = 0x3f3f3f3f3f3f3f;
const int mod = 1e9 + 7;
int n;
double d[maxn][maxn];
double x[maxn], y[maxn], z[maxn], r[maxn];

void init() {
    for (int i = 1; i <= n; i++)
            for (int j = 1; j <= n; j++)
                d[i][j] = INF;
}

double getDis(int a, int b) {
    double dis = sqrt((x[a] - x[b]) * (x[a] - x[b]) + (y[a] - y[b]) * (y[a] - y[b]) + (z[a] - z[b]) * (z[a] - z[b])) - r[a] - r[b];  // 洞内可以瞬移，所以减去半径ra和rb
    return dis > 0 ? dis : 0;   // 两个球体洞的圆心距离减去ra和rb如果小于0，则代表两个球体相交，因为在洞内是可以瞬移的，所以我们视作距离为0
}

void floyd() {
    for (int k = 1; k <= n; k++)
        for (int i = 1; i <= n; i++)
            for (int j = 1; j <= n; j++)
                if (d[i][k] < INF && d[k][j] < INF)
                    d[i][j] = min(d[i][j], d[i][k] + d[k][j]);
}

void solve() {
    int kase = 0;
    while (cin >> n, n != -1) {
        cout << "Cheese " << ++kase << ": Travel time = ";
        // 输入球形洞口的坐标半径
        for (int i = 1; i <= n; i++) {
            cin >> x[i] >> y[i] >> z[i] >> r[i];
        }
        n += 2;
        // 输入起点终点坐标半径，视为半径为0的洞
        for (int i = n - 1; i <= n; i++) {
            cin >> x[i] >> y[i] >> z[i];
            r[i] = 0;
        }
        // 初始化d[][]
        init();
        // 建图
        for (int i = 1; i <= n; i++) {
            for (int j = 1; j <= n; j++) {
                if (i != j)
                    d[i][j] = getDis(i, j);
            }
        }
        // 计算最短路
        floyd();
        // 输出起点到终点的距离
        d[n - 1][n] *= 10;
        cout << fixed << setprecision(0) << d[n - 1][n] << " sec\n";
    }
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    cout << fixed;
    cout.precision(18);

    solve();
    return 0;
}