結果

問題 No.1411 Hundreds of Conditions Sequences
ユーザー akakimidoriakakimidori
提出日時 2021-02-23 18:54:50
言語 Rust
(1.77.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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 1 ms
5,376 KB
testcase_03 AC 1 ms
5,376 KB
testcase_04 AC 1 ms
5,376 KB
testcase_05 AC 1 ms
5,376 KB
testcase_06 AC 1 ms
5,376 KB
testcase_07 AC 1 ms
5,376 KB
testcase_08 AC 1 ms
5,376 KB
testcase_09 AC 1 ms
5,376 KB
testcase_10 AC 1 ms
5,376 KB
testcase_11 AC 1 ms
5,376 KB
testcase_12 AC 103 ms
24,280 KB
testcase_13 AC 81 ms
23,308 KB
testcase_14 AC 54 ms
20,284 KB
testcase_15 AC 80 ms
23,204 KB
testcase_16 AC 23 ms
13,092 KB
testcase_17 AC 78 ms
23,052 KB
testcase_18 AC 101 ms
24,348 KB
testcase_19 AC 102 ms
24,656 KB
testcase_20 AC 54 ms
20,288 KB
testcase_21 AC 111 ms
24,772 KB
testcase_22 AC 56 ms
20,956 KB
testcase_23 AC 102 ms
24,276 KB
testcase_24 AC 15 ms
9,856 KB
testcase_25 AC 103 ms
24,388 KB
testcase_26 AC 76 ms
22,936 KB
testcase_27 AC 53 ms
20,492 KB
testcase_28 AC 85 ms
23,564 KB
testcase_29 AC 29 ms
15,204 KB
testcase_30 AC 70 ms
22,356 KB
testcase_31 AC 70 ms
21,812 KB
testcase_32 AC 109 ms
24,780 KB
testcase_33 AC 111 ms
24,912 KB
testcase_34 AC 111 ms
24,824 KB
testcase_35 AC 114 ms
24,796 KB
testcase_36 AC 112 ms
24,792 KB
testcase_37 AC 111 ms
24,784 KB
testcase_38 AC 109 ms
24,908 KB
testcase_39 AC 109 ms
24,916 KB
testcase_40 AC 113 ms
24,972 KB
testcase_41 AC 108 ms
24,944 KB
testcase_42 AC 107 ms
24,780 KB
testcase_43 AC 108 ms
24,740 KB
testcase_44 AC 110 ms
24,872 KB
testcase_45 AC 110 ms
24,908 KB
testcase_46 AC 111 ms
24,812 KB
testcase_47 AC 115 ms
24,904 KB
testcase_48 AC 113 ms
24,916 KB
testcase_49 AC 105 ms
24,876 KB
testcase_50 AC 114 ms
24,780 KB
testcase_51 AC 110 ms
24,912 KB
testcase_52 AC 91 ms
16,588 KB
testcase_53 AC 95 ms
16,588 KB
testcase_54 AC 95 ms
16,508 KB
testcase_55 AC 95 ms
16,592 KB
testcase_56 AC 93 ms
16,344 KB
testcase_57 AC 92 ms
16,584 KB
testcase_58 AC 98 ms
16,648 KB
testcase_59 AC 98 ms
16,592 KB
testcase_60 AC 98 ms
16,552 KB
testcase_61 AC 94 ms
16,484 KB
testcase_62 AC 59 ms
10,464 KB
testcase_63 AC 86 ms
26,188 KB
権限があれば一括ダウンロードができます

ソースコード

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