結果

問題 No.1731 Product of Subsequence
ユーザー phspls
提出日時 2022-10-17 23:35:19
言語 Rust
(1.83.0 + proconio)
結果
AC  
実行時間 816 ms / 2,000 ms
コード長 1,189 bytes
コンパイル時間 13,131 ms
コンパイル使用メモリ 397,896 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-28 12:23:48
合計ジャッジ時間 21,330 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

use std::collections::HashMap;

const MOD: usize = 1e9 as usize + 7;

fn gcd(a: usize, b: usize) -> usize {
    if b == 0 { return a; }
    gcd(b, a%b)
}

fn main() {
    let mut nk = String::new();
    std::io::stdin().read_line(&mut nk).ok();
    let nk: Vec<usize> = nk.trim().split_whitespace().map(|s| s.parse().unwrap()).collect();
    let n = nk[0];
    let k = nk[1];
    let mut a = String::new();
    std::io::stdin().read_line(&mut a).ok();
    let a: Vec<usize> = a.trim().split_whitespace().map(|s| s.parse().unwrap()).collect();

    let mut dp = HashMap::new();
    dp.insert(1usize, 1usize);
    for i in 0..n {
        let mut next_dp = HashMap::new();
        for (&val, &cnt) in dp.iter() {
            let gcdval = gcd(a[i], k);
            let gcdval = gcd(gcdval*val, k);
            *next_dp.entry(val).or_insert(0usize) += cnt;
            *next_dp.entry(val).or_insert(0usize) %= MOD;
            *next_dp.entry(gcdval).or_insert(0usize) += cnt;
            *next_dp.entry(gcdval).or_insert(0usize) %= MOD;
        }
        dp = next_dp;
    }
    let result = MOD + *dp.get(&k).unwrap_or(&0usize) - if k == 1 { 1 } else { 0 };
    println!("{}", result % MOD);
}
0