結果

問題 No.2829 GCD Divination
コンテスト
ユーザー northward
提出日時 2024-08-02 23:01:05
言語 Rust
(1.94.0 + proconio + num + itertools)
コンパイル:
/usr/bin/rustc_custom
実行:
./target/release/main
結果
MLE  
実行時間 -
コード長 1,174 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 3,949 ms
コンパイル使用メモリ 195,080 KB
実行使用メモリ 807,424 KB
最終ジャッジ日時 2026-04-03 19:45:58
合計ジャッジ時間 71,027 ms
ジャッジサーバーID
（参考情報） 		judge1_0 / judge4_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other MLE * 35
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#![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 2..=10_000_000 {
        if sieve[p].is_empty() {
            for q in (p..10_000_010).step_by(p) {
                sieve[q].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 d == k {
                continue;
            }

            if sieve[k].len() == 1 {
                s += dp[d];
                cnt += 1;
            } else if d == sieve[k][0] {
                let m = sieve[k][1];

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

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

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

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