~~這題是在幹嘛啊？怕不是又是a*b-a-b~~

然而萬萬沒想到，這是到圖論題

設\(dis[i]\)為在\(%x\)意義下，能到達的樓層為i的最小值

也就是說只用\(y, z\)能到達的樓層在\(%x\)意義下的最小值

不難推出方程

\[dis[(i + y) \% x] = min(dis[(i + y) \% x], dis[i] + y)\]

\[dis[(i + z) \% x] = min(dis[(i + z) \% x], dis[i] + z)\]

看到這兩個柿子，不難想到圖論中的最短路，所以我們可以用最短路來求出\(0~x\)的\(dis\)值

求除了dis以後有什麽用呢？

跟據dis定義，我們可以通過y, z到達的最小樓層，以該樓層為起點（設該點為s），我們可以跳到\(s + x\)

總共我們可以跳到\((H - dis[i]) / x + 1\)層

不難證明，每個樓層是不會被重復統計的

於是我們就可以以優秀的復雜度完成此題了

#include<bits/stdc++.h>
using namespace std;
#define il inline
#define re register
#define int long long
il int read() {
    re int x = 0, f = 1; re char c = getchar();
    while(c < '0' || c > '9') { if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - 48, c = getchar();
    return x * f;
}
#define rep(i, s, t) for(re int i = s; i <= t; ++ i)
#define mem(k, p) memset(k, p, sizeof(k))
#define maxn 100005
int n, m, x, y, z, dis[maxn], vis[maxn], ans;
queue<int>q;
il void SPFA() {
    mem(dis, 63), q.push(1 % x), dis[1 % x] = 1;
    while(!q.empty()) {
        int u = q.front(); q.pop(), vis[u] = 0;
        int v = (u + y) % x;
        if(dis[v] > dis[u] + y) {
            dis[v] = dis[u] + y;
            if(!vis[v]) vis[v] = 1, q.push(v);
        }
        v = (u + z) % x;
        if(dis[v] > dis[u] + z) {
            dis[v] = dis[u] + z;
            if(!vis[v]) vis[v] = 1, q.push(v);
        }
    }
}
signed main() {
    n = read(), x = read(), y = read(), z = read();
    SPFA();
    rep(i, 0, x - 1) if(dis[i] <= n) ans += (n - dis[i]) / x + 1;
    printf("%lld", ans);
    return 0;
}
 

跳樓機（同余類最短路）