結果

問題 No.58 イカサマなサイコロ
ユーザー tonyu0tonyu0
提出日時 2020-03-26 07:58:21
言語 Rust
(1.83.0 + proconio)
結果
AC  
実行時間 1 ms / 5,000 ms
コード長 1,173 bytes
コンパイル時間 11,270 ms
コンパイル使用メモリ 400,792 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2025-01-02 00:44:55
合計ジャッジ時間 12,234 ms
ジャッジサーバーID
(参考情報)
judge5 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

diff #

use std::io::*;

fn main() {
    let mut s: String = String::new();
    std::io::stdin().read_to_string(&mut s).ok();
    let mut itr = s.trim().split_whitespace();
    let n: usize = itr.next().unwrap().parse().unwrap();
    let k: usize = itr.next().unwrap().parse().unwrap();

    let mut dp1: Vec<f64> = vec![0.0; 101];
    let mut dp2: Vec<f64> = vec![0.0; 101];
    dp1[0] = 1.0;
    dp2[0] = 1.0;
    for _ in 0..n - k {
        let mut tmp = vec![0.0; 101];
        for i in 0..80 {
            for j in 1..7 {
                tmp[i + j] += dp1[i] / 6.0;
            }
        }
        dp1 = tmp;
    }
    for _ in 0..k {
        let mut tmp = vec![0.0; 101];
        for i in 0..80 {
            for j in 4..7 {
                tmp[i + j] += dp1[i] / 3.0;
            }
        }
        dp1 = tmp;
    }

    for _ in 0..n {
        let mut tmp = vec![0.0; 101];
        for i in 0..80 {
            for j in 1..7 {
                tmp[i + j] += dp2[i] / 6.0;
            }
        }
        dp2 = tmp;
    }

    let mut ans = 0.0;
    for i in 0..80 {
        for j in 0..i {
            ans += dp1[i] * dp2[j];
        }
    }
    println!("{:.10}", ans);
}
0