結果
問題 | No.1358 [Zelkova 2nd Tune *] 語るなら枚数を... |
ユーザー |
|
提出日時 | 2021-02-16 17:23:59 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 600 ms / 2,000 ms |
コード長 | 1,207 bytes |
コンパイル時間 | 2,084 ms |
コンパイル使用メモリ | 197,908 KB |
最終ジャッジ日時 | 2025-01-18 21:30:41 |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 17 |
ソースコード
#include <bits/stdc++.h> using namespace std; const long long MOD = 1e9 + 7; long long extgcd(long long a, long long b, long long &x, long long &y) { long long d = a; if (b != 0) { d = extgcd(b, a%b, y, x); y -= a/b * x; } else { x = 1; y = 0; } return d; } long long solve(long long a, long long b, long long c, long long n) { long long ret = 0; long long x, y; long long g = extgcd(a, b, x, y); for (int i = 0; i * c <= n; i++) { long long r = n - i * c; if (r == 0) ret = (ret + 1) % MOD; else if (r % g == 0) { long long xt = x * (r / g); long long yt = y * (r / g); long long l = xt / (b / g); xt -= l * (b / g); yt += l * (a / g); if (xt < 0) { xt += b / g; yt -= a / g; } if (yt >= 0) { ret += yt / (a / g) + 1; ret %= MOD; } } } return ret; } int main() { ios_base::sync_with_stdio(0); cin.tie(0); int T; cin >> T; for (int i = 0; i < T; i++) { long long n, k, h, y; cin >> n >> k >> h >> y; vector<long long> v{n, k, h}; sort(v.begin(), v.end()); cout << solve(v[0], v[1], v[2], y) << endl; } return 0; }