結果

問題 No.1358 [Zelkova 2nd Tune *] 語るなら枚数を...
ユーザー firiexpfiriexp
提出日時 2021-01-23 17:51:03
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,307 bytes
コンパイル時間 953 ms
コンパイル使用メモリ 102,140 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-06-10 01:29:20
合計ジャッジ時間 10,115 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 8 ms
6,812 KB
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <algorithm>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <numeric>
#include <bitset>
#include <cmath>

static const int MOD = 1000000007;
using ll = long long;
using u32 = unsigned;
using u64 = unsigned long long;
using namespace std;

template<class T> constexpr T INF = ::numeric_limits<T>::max() / 32 * 15 + 208;

template<typename T>
T extgcd(T a, T b, T &x ,T &y){
    for (T u = y = 1, v = x = 0; a; ) {
        ll q = b/a;
        swap(x -= q*u, u);
        swap(y -= q*v, v);
        swap(b -= q*a, a);
    }
    return b;
}

template<class T>T divfloor(T x, T y){
    if((x < 0) != (y < 0)) return (x+y-1)/y;
    else return x/y;
}

template<class T>T divceil(T x, T y){
    if((x < 0) != (y < 0)) return x/y;
    else return (x+y-1)/y;
}

void solve(){
    ll a, b, c, y;
    cin >> a >> b >> c >> y;
    if(a > c) swap(a, c);
    if(b > c) swap(b, c);
    ll ans = 0, g = gcd(a, b);
    a /= g; b /= g;
    for (ll x2 = y; x2 >= 0; x2 -= c) {
        if(x2 % g) continue;
        ll x = x2/g, p, q;
        extgcd(a, b, p, q);
        p *= x; q *= x;
        ans += max(0LL, divfloor(q, a)-divceil(-p, b)+1);
    }
    cout << ans%MOD << "\n";
}

int main() {
    int t;
    cin >> t;
    while(t--){
        solve();
    }
    return 0;
}
0