結果

問題 No.1411 Hundreds of Conditions Sequences
ユーザー akakimidori
提出日時 2021-02-23 18:54:50
言語 Rust
(1.83.0 + proconio)
結果
AC  
実行時間 115 ms / 2,000 ms
コード長 2,023 bytes
コンパイル時間 15,623 ms
コンパイル使用メモリ 382,460 KB
実行使用メモリ 26,188 KB
最終ジャッジ日時 2024-09-22 13:50:04
合計ジャッジ時間 23,977 ms
ジャッジサーバーID
(参考情報)
judge2 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 62
権限があれば一括ダウンロードができます

ソースコード

diff #

fn read() -> Vec<usize> {
    let mut s = String::new();
    use std::io::*;
    std::io::stdin().read_to_string(&mut s).unwrap();
    let mut it = s.trim().split_whitespace();
    it.next();
    it.flat_map(|s| s.parse()).collect()
}

const MOD: usize = 1_000_000_007;

fn pow(mut r: usize, mut n: usize) -> usize {
    let mut t = 1;
    while n > 0 {
        if n & 1 == 1 {
            t = t * r % MOD;
        }
        r = r * r % MOD;
        n >>= 1;
    }
    t
}

fn main() {
    let a = read();
    let w = *a.iter().max().unwrap();
    let mut p = (0..=w).collect::<Vec<_>>();
    for i in (2..).take_while(|p| p * p <= w) {
        if p[i] == i {
            let mut j = i * i;
            while let Some(p) = p.get_mut(j) {
                *p = i;
                j += i;
            }
        }
    }
    let factorize = |mut n: usize| -> Vec<(usize, usize)> {
        let mut f = vec![];
        while n > 1 {
            let p = p[n];
            let mut c = 0;
            while n % p == 0 {
                n /= p;
                c += 1;
            }
            f.push((p, c));
        }
        f
    };
    let mut dp = vec![(0, 0); w + 1];
    for a in a.iter() {
        for (p, c) in factorize(*a) {
            let po = &mut dp[p];
            if c > po.0 {
                *po = (c, po.0);
            } else if c > po.1 {
                po.1 = c;
            }
        }
    }
    let all = a.iter().fold(1, |s, a| s * *a % MOD);
    let lcm = dp.iter().enumerate().fold(1, |s, (k, p)| s * k.pow(p.0 as u32) % MOD);
    use std::io::Write;
    let out = std::io::stdout();
    let mut out = std::io::BufWriter::new(out.lock());
    for a in a.iter() {
        let all = all * pow(*a, MOD - 2);
        let mut lcm = lcm;
        for (p, c) in factorize(*a) {
            let (x, y) = dp[p];
            if c == x {
                lcm = lcm * pow(pow(p, MOD - 2), x - y) % MOD;
            }
        }
        let ans = (all - lcm + MOD) % MOD;
        writeln!(out, "{}", ans).ok();
    }
}
0