結果

問題 No.162 8020運動
ユーザー nebukuro09nebukuro09
提出日時 2018-06-14 10:10:29
言語 C++17(clang)
(17.0.6 + boost 1.83.0)
結果
AC  
実行時間 2,053 ms / 5,000 ms
コード長 1,696 bytes
コンパイル時間 2,014 ms
コンパイル使用メモリ 165,220 KB
実行使用メモリ 200,320 KB
最終ジャッジ日時 2024-05-07 17:46:26
合計ジャッジ時間 44,092 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 950 ms
200,316 KB
testcase_01 AC 955 ms
200,312 KB
testcase_02 AC 1,764 ms
200,192 KB
testcase_03 AC 1,636 ms
200,188 KB
testcase_04 AC 2,048 ms
200,192 KB
testcase_05 AC 2,031 ms
200,192 KB
testcase_06 AC 2,041 ms
200,320 KB
testcase_07 AC 2,041 ms
200,192 KB
testcase_08 AC 2,053 ms
200,064 KB
testcase_09 AC 938 ms
200,316 KB
testcase_10 AC 949 ms
200,192 KB
testcase_11 AC 951 ms
200,184 KB
testcase_12 AC 956 ms
200,188 KB
testcase_13 AC 946 ms
200,188 KB
testcase_14 AC 958 ms
200,192 KB
testcase_15 AC 946 ms
200,316 KB
testcase_16 AC 1,110 ms
200,192 KB
testcase_17 AC 957 ms
200,188 KB
testcase_18 AC 943 ms
200,192 KB
testcase_19 AC 1,748 ms
200,320 KB
testcase_20 AC 947 ms
200,192 KB
testcase_21 AC 1,814 ms
200,320 KB
testcase_22 AC 950 ms
200,192 KB
testcase_23 AC 1,819 ms
200,064 KB
testcase_24 AC 1,754 ms
200,192 KB
testcase_25 AC 1,850 ms
200,184 KB
testcase_26 AC 933 ms
200,064 KB
testcase_27 AC 1,523 ms
200,312 KB
testcase_28 AC 1,554 ms
200,316 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define REP(i,n) for (int i=0;i<(n);i++)
#define REP2(i,m,n) for (int i=m;i<(n);i++)
typedef long long ll;
typedef long double ld;

unordered_map<int, double> trans[1<<14];
double dp[2][1<<14];

int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);

    const int T = 14;
    int N;
    double P[3];
    cin >> N;
    REP(i, 3) cin >> P[i], P[i] /= 100.0;

    for (int m1 = 0; m1 < (1 << T); ++m1) {
        for(int m2 = m1; m2 >= 0; --m2) {
            m2 &= m1;
            double 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);
                    double 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;
        }
    }

    REP(i, 2) REP(j, 1<<T) dp[i][j] = 0;
    dp[0][(1<<T)-1] = 1;
    int cur = 0, tar = 1;

    REP(i, 80-N) {
        REP(j, 1<<T)
            dp[tar][j] = 0;
        REP(m1, 1<<T) {
            if (dp[cur][m1] == 0)
                continue;
            for (auto itr: trans[m1])
                dp[tar][itr.first] += dp[cur][m1] * itr.second;
        }
        swap(cur, tar);
    }

    double ans = 0;
    REP(mask, 1<<T)
        ans += dp[cur][mask] * __builtin_popcount(mask);
    cout.precision(20);
    cout << fixed  << ans*2 << endl;
}
0