結果

問題 No.2829 GCD Divination
ユーザー naut3naut3
提出日時 2024-08-02 23:08:42
言語 Rust
(1.77.0 + proconio)
結果
MLE  
実行時間 -
コード長 1,179 bytes
コンパイル時間 13,401 ms
コンパイル使用メモリ 401,880 KB
実行使用メモリ 589,332 KB
最終ジャッジ日時 2024-08-02 23:09:03
合計ジャッジ時間 20,205 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 MLE -
testcase_01 -- -
testcase_02 -- -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#![allow(non_snake_case)]
#[allow(unused_imports)]
use proconio::{fastout, input, marker::*};

#[fastout]
fn main() {
    input! {
        N: usize,
    }

    let mut sieve = vec![vec![]; 10_000_010];

    for p in 2u32..=10_000_000 {
        if sieve[p as usize].is_empty() {
            for q in (p * 2..10_000_010).step_by(p as usize) {
                sieve[q as usize].push(p);
            }
        }
    }

    let mut dp = vec![0.; 10_000_010];

    for k in 2..=10_000_000 {
        let mut s = 0.;
        let mut cnt = 0;

        for &d in sieve[k].iter() {
            if sieve[k].len() == 1 {
                s += dp[d as usize];
                cnt += 1;
            } else if d == sieve[k][0] {
                let m = sieve[k][1];

                s += dp[d as usize] * (m - 1) as f32;
                cnt += m - 1;
            } else {
                let m = sieve[k][0];

                s += dp[d as usize] * (m - 1) as f32;
                cnt += m - 1;
            }
        }

        dp[k] = s + (1 * (k - 1 - cnt as usize)) as f32;
        dp[k] /= (k - 1) as f32;
        dp[k] += (k) as f32 / (k - 1) as f32;
    }

    println!("{}", dp[N] - 1.);
}
0