結果

問題 No.1358 [Zelkova 2nd Tune *] 語るなら枚数を...
ユーザー rokahikou1rokahikou1
提出日時 2021-01-23 20:57:03
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 385 ms / 2,000 ms
コード長 3,521 bytes
コンパイル時間 2,525 ms
コンパイル使用メモリ 172,216 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-08-30 05:20:12
合計ジャッジ時間 5,027 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 1 ms
4,384 KB
testcase_04 AC 1 ms
4,384 KB
testcase_05 AC 2 ms
4,384 KB
testcase_06 AC 6 ms
4,384 KB
testcase_07 AC 5 ms
4,380 KB
testcase_08 AC 5 ms
4,384 KB
testcase_09 AC 5 ms
4,380 KB
testcase_10 AC 4 ms
4,384 KB
testcase_11 AC 345 ms
4,384 KB
testcase_12 AC 303 ms
4,380 KB
testcase_13 AC 385 ms
4,384 KB
testcase_14 AC 276 ms
4,384 KB
testcase_15 AC 301 ms
4,380 KB
testcase_16 AC 361 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define rep(i, n) for(int(i) = 0; (i) < (n); (i)++)
#define FOR(i, m, n) for(int(i) = (m); (i) < (n); (i)++)
#define ALL(v) (v).begin(), (v).end()
#define LLA(v) (v).rbegin(), (v).rend()
#define SZ(v) (v).size()
#define INT(...)                                                               \
    int __VA_ARGS__;                                                           \
    read(__VA_ARGS__)
#define LL(...)                                                                \
    ll __VA_ARGS__;                                                            \
    read(__VA_ARGS__)
#define DOUBLE(...)                                                            \
    double __VA_ARGS__;                                                        \
    read(__VA_ARGS__)
#define CHAR(...)                                                              \
    char __VA_ARGS__;                                                          \
    read(__VA_ARGS__)
#define STRING(...)                                                            \
    string __VA_ARGS__;                                                        \
    read(__VA_ARGS__)
#define VEC(type, name, size)                                                  \
    vector<type> name(size);                                                   \
    read(name)
using namespace std;
using ll = long long;
using pii = pair<int, int>;
using pll = pair<ll, ll>;
using Graph = vector<vector<int>>;
template <typename T> struct edge {
    int to;
    T cost;
    edge(int t, T c) : to(t), cost(c) {}
};
template <typename T> using WGraph = vector<vector<edge<T>>>;
const int INF = 1 << 30;
const ll LINF = 1LL << 60;
const int MOD = 1e9 + 7;
template <class T> inline vector<T> make_vec(size_t a, T val) {
    return vector<T>(a, val);
}
template <class... Ts> inline auto make_vec(size_t a, Ts... ts) {
    return vector<decltype(make_vec(ts...))>(a, make_vec(ts...));
}

void read() {}
template <class T> inline void read(T &a) { cin >> a; }
template <class T, class S> inline void read(pair<T, S> &p) {
    read(p.first), read(p.second);
}
template <class T> inline void read(vector<T> &v) {
    for(auto &a : v)
        read(a);
}
template <class Head, class... Tail>
inline void read(Head &head, Tail &...tail) {
    read(head), read(tail...);
}
template <class T> inline void chmax(T &a, T &b) { (a < b ? a = b : a); }
template <class T> inline void chmin(T &a, T &b) { (a > b ? a = b : a); }

ll extgcd(ll a, ll b, ll &x, ll &y) {
    if(b == 0) {
        x = 1;
        y = 0;
        return a;
    }
    ll d = extgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}

ll floor_div(ll a, ll b) {
    assert(b != 0);
    if(b < 0)
        a = -a, b = -b;
    return a >= 0 ? a / b : (a - (b - 1)) / b;
}

ll ceil_div(ll a, ll b) {
    assert(b != 0);
    if(b < 0)
        a = -a, b = -b;
    return a >= 0 ? (a + b - 1) / b : a / b;
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    INT(t);
    while(t--) {
        VEC(ll, a, 3);
        LL(y);
        sort(LLA(a));
        ll n = a[0], k = a[1], h = a[2];
        ll res = 0;
        ll g = __gcd(k, h);
        k /= g, h /= g;
        ll p, q;
        extgcd(k, h, p, q);
        for(int i = 0; i <= y / n; i++) {
            ll s = y - n * i;
            if(s % g != 0)
                continue;
            s /= g;
            res += max(0LL, floor_div(s * q, k) - ceil_div(-s * p, h) + 1);
            res %= MOD;
        }
        cout << res << endl;
    }
}
0