結果

問題 No.162 8020運動
ユーザー nebukuro09
提出日時 2018-06-14 09:41:50
言語 D
(dmd 2.109.1)
結果
MLE  
実行時間 -
コード長 1,559 bytes
コンパイル時間 2,457 ms
コンパイル使用メモリ 159,888 KB
実行使用メモリ 609,564 KB
最終ジャッジ日時 2024-06-13 01:13:37
合計ジャッジ時間 14,274 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2 MLE * 1
other -- * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.stdio, std.array, std.string, std.conv, std.algorithm;
import std.typecons, std.range, std.random, std.math, std.container;
import std.numeric, std.bigint, core.bitop;

void main() {
    immutable int T = 14;
    auto N = readln.chomp.to!int;
    auto P = readln.split.map!(a => a.to!int / 100.0L).array;
    auto trans = new real[int][](1<<T);

    for (int m1 = 0; m1 < (1 << T); ++m1) {
        for(int m2 = m1; m2 >= 0; --m2) {
            m2 &= m1;
            real tmp = 1;
            for (int i = 0; i < T; ++i) {
                if ((1 << i) & m1) {
                    bool l = (i != 0) && ((1 << (i-1)) & m1);
                    bool r = (i != T-1) && ((1 << (i+1)) & m1);
                    real safe;
                    if (l && r)
                        safe = 1 - P[2];
                    else if (l || r)
                        safe = 1 - P[1];
                    else
                        safe = 1 - P[0];
                    tmp *=  ((1 << i) & m2) ? safe : 1 - safe;
                }
            }
            trans[m1][m2] = tmp;
        }
    }

    auto dp = new real[][](80-N+1, 1<<T);
    foreach (i; 0..80-N+1) dp[i].fill(0);
    dp[0][(1<<T)-1] = 1;

    foreach (i; 0..80-N) {
        foreach (m1; 0..(1<<T)) {
            if (dp[i][m1] == 0)
                continue;
            foreach (m2; trans[m1].keys)
                dp[i+1][m2] += dp[i][m1] * trans[m1][m2];
        }
    }

    real ans = 0;
    foreach (mask; 0..(1<<T))
        ans += dp[80-N][mask] * mask.popcnt;
    writefln("%.09f", ans*2);
}
0