結果

問題 No.1151 チャレンジゲーム
ユーザー kimiyukikimiyuki
提出日時 2020-08-07 23:20:34
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 21 ms / 2,000 ms
コード長 2,577 bytes
コンパイル時間 2,113 ms
コンパイル使用メモリ 206,764 KB
実行使用メモリ 4,400 KB
最終ジャッジ日時 2023-10-25 04:48:08
合計ジャッジ時間 4,203 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,348 KB
testcase_01 AC 2 ms
4,348 KB
testcase_02 AC 2 ms
4,348 KB
testcase_03 AC 2 ms
4,348 KB
testcase_04 AC 2 ms
4,348 KB
testcase_05 AC 2 ms
4,348 KB
testcase_06 AC 2 ms
4,348 KB
testcase_07 AC 2 ms
4,348 KB
testcase_08 AC 2 ms
4,348 KB
testcase_09 AC 2 ms
4,348 KB
testcase_10 AC 2 ms
4,348 KB
testcase_11 AC 2 ms
4,348 KB
testcase_12 AC 2 ms
4,348 KB
testcase_13 AC 8 ms
4,348 KB
testcase_14 AC 7 ms
4,348 KB
testcase_15 AC 7 ms
4,348 KB
testcase_16 AC 7 ms
4,348 KB
testcase_17 AC 8 ms
4,348 KB
testcase_18 AC 7 ms
4,348 KB
testcase_19 AC 7 ms
4,348 KB
testcase_20 AC 7 ms
4,348 KB
testcase_21 AC 7 ms
4,348 KB
testcase_22 AC 8 ms
4,348 KB
testcase_23 AC 8 ms
4,348 KB
testcase_24 AC 7 ms
4,348 KB
testcase_25 AC 7 ms
4,348 KB
testcase_26 AC 8 ms
4,348 KB
testcase_27 AC 7 ms
4,348 KB
testcase_28 AC 20 ms
4,400 KB
testcase_29 AC 21 ms
4,400 KB
testcase_30 AC 20 ms
4,400 KB
testcase_31 AC 21 ms
4,400 KB
testcase_32 AC 21 ms
4,400 KB
testcase_33 AC 21 ms
4,400 KB
testcase_34 AC 20 ms
4,400 KB
testcase_35 AC 20 ms
4,400 KB
testcase_36 AC 21 ms
4,400 KB
testcase_37 AC 21 ms
4,400 KB
testcase_38 AC 21 ms
4,400 KB
testcase_39 AC 21 ms
4,400 KB
testcase_40 AC 21 ms
4,400 KB
testcase_41 AC 21 ms
4,400 KB
testcase_42 AC 21 ms
4,400 KB
testcase_43 AC 21 ms
4,400 KB
testcase_44 AC 21 ms
4,400 KB
testcase_45 AC 20 ms
4,400 KB
testcase_46 AC 21 ms
4,400 KB
testcase_47 AC 20 ms
4,400 KB
testcase_48 AC 15 ms
4,400 KB
testcase_49 AC 20 ms
4,400 KB
testcase_50 AC 19 ms
4,400 KB
testcase_51 AC 21 ms
4,400 KB
testcase_52 AC 20 ms
4,400 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define REP(i, n) for (int i = 0; (i) < (int)(n); ++ (i))
#define REP3(i, m, n) for (int i = (m); (i) < (int)(n); ++ (i))
#define REP_R(i, n) for (int i = (int)(n) - 1; (i) >= 0; -- (i))
#define REP3R(i, m, n) for (int i = (int)(n) - 1; (i) >= (int)(m); -- (i))
#define ALL(x) ::std::begin(x), ::std::end(x)
using namespace std;

auto solve(int n, const vector<int> & a) {
    vector<int> pow_3(n + 1);
    pow_3[0] = 1;
    REP (i, n) {
        pow_3[i + 1] = pow_3[i] * 3;
    }
    vector<array<double, 2> > dp(pow_3[n]);
    REP_R (state, pow_3[n]) {
        int available = 0;
        int score = 0;
        {  // unpack
            int s = state;
            REP (i, n) {
                if (s % 3 == 0) {
                    available |= 1 << i;
                } else if (s % 3 == 1) {
                    score += a[i];
                } else {
                    score -= a[i];
                }
                s /= 3;
            }
        }
        if (not available) {
            dp[state][false] = score > 0;
            dp[state][true] = score <= 0;
        } else {
            REP (turn, 2) {
                REP (i, n) if (available & (1 << i)) {
                    double prob = 1;
                    REP (j, n) if (available & (1 << j)) {
                        double p = 1.0 / a[i];
                        double q = 1.0 / a[j];
                        if (a[i] == 1) {
                            double x = 1 - dp[state + (turn + 1) * pow_3[i]][not turn];
                            prob = min(prob, x);
                        } else if (a[j] == 1) {
                            double x = p * (1 - dp[state + (turn + 1) * pow_3[i]][not turn]);
                            double y = (1 - p) * dp[state + (not turn + 1) * pow_3[j]][turn];
                            prob = min(prob, x + y);
                        } else {
                            double x = p * (1 - dp[state + (turn + 1) * pow_3[i]][not turn]);
                            double y = (1 - p) * q * dp[state + (not turn + 1) * pow_3[j]][turn];
                            double r = (1 - p) * (1 - q);
                            prob = min(prob, (x + y) / (1 - r));
                        }
                    }
                    dp[state][turn] = max(dp[state][turn], prob);
                }
            }
        }
    }
    return dp[0][false];
}

int main() {
    int n; cin >> n;
    vector<int> a(n);
    REP (i, n) {
        cin >> a[i];
    }
    auto ans = solve(n, a);
    cout << setprecision(18) << ans << endl;
    return 0;
}
0