JZOJ4722 跳楼机巧妙地转换为最短路模型

最新推荐文章于 2022-05-06 14:38:28 发布

YxuanwKeith

最新推荐文章于 2022-05-06 14:38:28 发布

阅读量1.1k

点赞数 2

分类专栏：算法-最短路

本文链接：https://blog.csdn.net/YxuanwKeith/article/details/52267231

版权

算法-最短路专栏收录该内容

4 篇文章 0 订阅

订阅专栏

题目大意

有一个 $H$ 高的楼房，现在给你三个正整数 $x,y,z$ 表示，在每一层 $i$ ，你有4种选择，到 $i+x$ ， $i+y$ ， $i+z$ 层或第 $1$ 层。问能到多少个楼层。

$H \leq 10^{18}$
$x,y,z\leq 10^5$

解题思路

也就是说有多少 $p=ax+by+cx(p \leq H)$ 我们可一考虑一种一种操来做，假设我们一开始用前两个操作到了第 $i$ 层，然后后面全部用第三种操作，那么对答案的贡献就是 $\lfloor\frac{H-i}{z}\rfloor$ 。但是会用算重的情况，那么我们考虑用一种不会算重的方法来统计答案，设 $F_i(i \in[0,z))$ 表示在用前两个数能到达最小的位置 $F_i$ ，满足 $F_i \bmod{z} = i$ ，那么最后的答案就是 $\sum_{i=0}^{z-1}(\lfloor\frac{H-F_i}{z}\rfloor+1)$ 。

那么关键就是如何求出 $F$ 数组。那么我们考虑 $x$ 和 $y$ 在 $F$ 数组上的转移，不难发现：
$F_{(i+y) \bmod{z}} = F_i+y$
$F_{(i+z) \bmod(z)} = F_i+z$
那么我们把这两个点之间连边，跑一次SPFA就能得到 $F$ 数组了。

程序

我的程序与上面说的有点不一样， $Dis_i$ 表示达到 j <script type="math/tex" id="MathJax-Element-1803">j</script>位置最小的块是哪一块（每z分一块）

//YxuanwKeith
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>

using namespace std;
typedef long long LL;

const int MAXN = 5e5 + 5;

LL H, Dis[MAXN], Ans;
bool Flag[MAXN];
int S[3], D[MAXN];
int tot, Last[MAXN], Go[MAXN * 2], Len[MAXN * 2], Next[MAXN * 2];

void Link(int u, int v, int Dis) {
    Next[++ tot] = Last[u], Last[u] = tot, Go[tot] = v, Len[tot] = Dis; 
}

void MakeRoad(int x, int y) {
    int Dis = y / x, Lef = y % x;
    for (int i = 0; i < x; i ++) 
        Link(i, (i + Lef) % x, Dis + (i + Lef > x - 1));
}

void Spfa() {
    memset(Dis, 60, sizeof Dis);
    int l = 0, r = 1;
    Dis[0] = 0, D[1] = 0;
    while (l != r) {
        l = (l + 1) % MAXN;
        int Now = D[l];
        for (int p = Last[Now]; p; p = Next[p]) {
            int v = Go[p];
            if (Dis[v] > Dis[Now] + Len[p]) {
                Dis[v] = Dis[Now] + Len[p];
                if (!Flag[v]) {
                    Flag[v] = 1;
                    r = (r + 1) % MAXN;
                    D[r] = v;
                }
            }
        }
        Flag[Now] = 0;
    }
}

void Prepare() {
    sort(S, S + 3);
    int x = S[0], y = S[1], z = S[2];
    H --;
    MakeRoad(x, y);
    MakeRoad(x, z);
    Spfa();
}

void Solve() {
    int x = S[0];
    for (int i = 0; i < x; i ++) 
        Ans += max((H / x) + ((H % x) >= i) - Dis[i], 0ll);
    printf("%lld\n", Ans);
}

int main() {
    scanf("%lld", &H);
    scanf("%d%d%d", &S[0], &S[1], &S[2]);
    Prepare();
    Solve();
}