結果

問題 No.1358 [Zelkova 2nd Tune *] 語るなら枚数を...
ユーザー tonyu0tonyu0
提出日時 2021-01-22 23:45:19
言語 Rust
(1.77.0 + proconio)
結果
WA  
実行時間 -
コード長 2,217 bytes
コンパイル時間 22,061 ms
コンパイル使用メモリ 376,668 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-08 19:02:47
合計ジャッジ時間 33,035 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 17 ms
5,248 KB
testcase_01 WA -
testcase_02 AC 0 ms
5,248 KB
testcase_03 AC 0 ms
5,248 KB
testcase_04 AC 1 ms
5,248 KB
testcase_05 AC 1 ms
5,376 KB
testcase_06 AC 22 ms
5,376 KB
testcase_07 AC 24 ms
5,376 KB
testcase_08 AC 54 ms
5,376 KB
testcase_09 AC 15 ms
5,376 KB
testcase_10 AC 53 ms
5,376 KB
testcase_11 TLE -
testcase_12 WA -
testcase_13 AC 1,990 ms
5,376 KB
testcase_14 AC 1,750 ms
5,376 KB
testcase_15 AC 1,758 ms
5,376 KB
testcase_16 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

use std::io::*;
const MOD: i64 = 1_000_000_07;
fn ext_gcd(a: i64, b: i64, x: &mut i64, y: &mut i64) -> i64 {
    if b == 0 {
        *x = 1;
        *y = 0;
        a
    } else {
        let d = ext_gcd(b, a % b, y, x);
        *y -= (a / b) * *x;
        d
    }
}

fn calc(a: i64, b: i64, c: i64, y: i64) -> i64 {
    let mut ans = 0;
    for i in 0..=y / a {
        let d = y - i * a;

        // d = b * x + c * y
        let mut x0 = 0;
        let mut y0 = 0;
        let g = ext_gcd(b, c, &mut x0, &mut y0);
        if d % g != 0 {
            continue;
        }
        x0 *= d / g;
        y0 *= d / g;
        let bs = b / g;
        let cs = c / g;
        if x0 < 0 {
            // while x0 < 0 {
            //     x0 += cs;
            //     y0 -= bs;
            // }
            let u = (-x0 + cs - 1) / cs;
            y0 -= u * bs;
            // while y0 >= 0 {
            //     ans += 1;
            //     y0 -= bs;
            // }
            if y0 >= 0 {
                ans = (ans + (y0 / bs + 1) % MOD) % MOD;
            }
        } else {
            // while y0 < 0 {
            //     y0 += bs;
            //     x0 -= cs;
            // }
            let u = (-y0 + bs - 1) / bs;
            x0 -= u * cs;
            // while x0 >= 0 {
            //     ans += 1;
            //     x0 -= cs;
            // }
            if x0 >= 0 {
                ans = (ans + (x0 / cs + 1) % MOD) % MOD
            }
        }
    }
    ans
}

fn main() {
    let mut s: String = String::new();
    std::io::stdin().read_to_string(&mut s).ok();
    let mut itr = s.trim().split_whitespace();
    let t: usize = itr.next().unwrap().parse().unwrap();
    for _ in 0..t {
        let n: i64 = itr.next().unwrap().parse().unwrap();
        let k: i64 = itr.next().unwrap().parse().unwrap();
        let h: i64 = itr.next().unwrap().parse().unwrap();
        let y: i64 = itr.next().unwrap().parse().unwrap();
        let nc = y / n;
        let kc = y / k;
        if nc <= 1000000 {
            println!("{}", calc(n, k, h, y));
        } else if kc <= 1000000 {
            println!("{}", calc(k, h, n, y));
        } else {
            println!("{}", calc(h, n, k, y));
        }
    }
}
0