結果

問題 No.2902 ZERO!!
ユーザー Yukino DX.Yukino DX.
提出日時 2024-10-26 12:12:05
言語 Rust
(1.77.0 + proconio)
結果
AC  
実行時間 38 ms / 2,000 ms
コード長 3,509 bytes
コンパイル時間 13,367 ms
コンパイル使用メモリ 403,368 KB
実行使用メモリ 10,876 KB
最終ジャッジ日時 2024-10-26 12:12:55
合計ジャッジ時間 15,803 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,816 KB
testcase_01 AC 1 ms
6,816 KB
testcase_02 AC 38 ms
10,876 KB
testcase_03 AC 1 ms
6,816 KB
testcase_04 AC 9 ms
6,820 KB
testcase_05 AC 19 ms
6,816 KB
testcase_06 AC 37 ms
10,472 KB
testcase_07 AC 12 ms
6,820 KB
testcase_08 AC 34 ms
9,792 KB
testcase_09 AC 34 ms
9,556 KB
testcase_10 AC 25 ms
7,644 KB
testcase_11 AC 20 ms
6,820 KB
testcase_12 AC 3 ms
6,820 KB
testcase_13 AC 34 ms
9,736 KB
testcase_14 AC 27 ms
7,864 KB
testcase_15 AC 29 ms
8,696 KB
testcase_16 AC 18 ms
6,816 KB
testcase_17 AC 33 ms
9,588 KB
testcase_18 AC 32 ms
9,112 KB
testcase_19 AC 28 ms
8,388 KB
testcase_20 AC 25 ms
7,368 KB
testcase_21 AC 26 ms
7,792 KB
testcase_22 AC 14 ms
6,816 KB
testcase_23 AC 18 ms
6,816 KB
testcase_24 AC 1 ms
6,820 KB
testcase_25 AC 1 ms
6,820 KB
testcase_26 AC 1 ms
6,816 KB
testcase_27 AC 1 ms
6,824 KB
testcase_28 AC 1 ms
6,820 KB
testcase_29 AC 1 ms
6,816 KB
testcase_30 AC 1 ms
6,816 KB
testcase_31 AC 1 ms
6,820 KB
testcase_32 AC 1 ms
6,820 KB
testcase_33 AC 1 ms
6,816 KB
testcase_34 AC 1 ms
6,816 KB
testcase_35 AC 1 ms
6,816 KB
testcase_36 AC 1 ms
6,816 KB
testcase_37 AC 1 ms
6,820 KB
testcase_38 AC 1 ms
6,816 KB
testcase_39 AC 1 ms
6,820 KB
testcase_40 AC 1 ms
6,816 KB
testcase_41 AC 1 ms
6,816 KB
testcase_42 AC 1 ms
6,816 KB
testcase_43 AC 1 ms
6,820 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

use proconio::input;

fn main() {
    input! {
        n:usize,
    }

    let ps = eratosthenes_sieve(n);
    let mut f_ps = vec![0; ps.len()];
    for (i, &p) in ps.iter().enumerate() {
        let mut q = p;
        while n >= q {
            f_ps[i] += n / q;
            q *= p;
        }
    }

    let e_2 = f_ps[0];
    let mut muls = vec![1; e_2 + 1];
    for &f_p in f_ps.iter() {
        for j in 1..=f_p {
            muls[j] *= 1 + f_p / j;
            muls[j] %= MOD;
        }
    }

    let mut ans = 0;
    for i in 1..=e_2 {
        ans += (muls[i] + MOD - 1) % MOD;
        ans %= MOD;
    }

    println!("{}", ans);
}

const MOD: usize = 998244353;

use math::*;
mod math {
    #[allow(dead_code)]
    pub fn pow(mut n: usize, mut e: usize, m: usize) -> usize {
        let mut res = 1;
        n %= m;
        while 0 < e {
            if e & 1 != 0 {
                res *= n;
                res %= m;
            }

            n *= n;
            n %= m;
            e >>= 1;
        }

        res
    }

    #[allow(dead_code)]
    pub fn gcd(a: usize, b: usize) -> usize {
        if b == 0 {
            a
        } else {
            gcd(b, a % b)
        }
    }

    #[allow(dead_code)]
    pub fn lcm(a: usize, b: usize) -> usize {
        a / gcd(a, b) * b
    }

    #[allow(dead_code)]
    pub fn prime_factor(mut n: usize) -> std::collections::HashMap<usize, usize> {
        let mut pfacs = std::collections::HashMap::new();
        for d in 2.. {
            if d * d > n {
                break;
            }

            while n % d == 0 {
                *pfacs.entry(d).or_insert(0) += 1;
                n /= d;
            }
        }
        if n != 1 {
            *pfacs.entry(n).or_insert(0) += 1;
        }

        pfacs
    }

    #[allow(dead_code)]
    pub fn divisors(n: usize) -> Vec<usize> {
        let mut divs = vec![];
        for d in (1..).take_while(|&d| d * d <= n) {
            if n % d == 0 {
                divs.push(d);
                if n / d != d {
                    divs.push(n / d);
                }
            }
        }

        divs
    }

    #[allow(dead_code)]
    pub fn eratosthenes_sieve(n: usize) -> Vec<usize> {
        let mut is_primes = vec![true; n + 1];
        is_primes[0] = false;
        is_primes[1] = false;
        let mut primes = vec![];

        for i in 2..=n {
            if !is_primes[i] {
                continue;
            }

            for j in (2 * i..=n).step_by(i) {
                is_primes[j] = false;
            }
        }

        for (val, is_prime) in is_primes.into_iter().enumerate() {
            if is_prime {
                primes.push(val);
            }
        }

        primes
    }

    #[allow(dead_code)]
    pub fn base_n(mut base10: usize, base: usize, len: usize) -> Vec<usize> {
        let mut base_n = vec![];
        for _ in 0..len {
            base_n.push(base10 % base);
            base10 /= base;
        }
        base_n.reverse();

        base_n
    }

    #[allow(dead_code)]
    pub fn div_floor(mut a: i64, mut b: i64) -> i64 {
        if b < 0 {
            a = -a;
            b = -b;
        }

        if a >= 0 {
            a / b
        } else {
            (a + 1) / b - 1
        }
    }

    #[allow(dead_code)]
    pub fn div_ceil(mut a: i64, mut b: i64) -> i64 {
        if b < 0 {
            a = -a;
            b = -b;
        }
        if a > 0 {
            (a - 1) / b + 1
        } else {
            a / b
        }
    }
}
0