結果

問題 No.1151 チャレンジゲーム
ユーザー tokusakuraitokusakurai
提出日時 2020-08-08 09:26:19
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 96 ms / 2,000 ms
コード長 2,409 bytes
コンパイル時間 2,058 ms
コンパイル使用メモリ 204,300 KB
実行使用メモリ 10,236 KB
最終ジャッジ日時 2024-10-01 06:54:29
合計ジャッジ時間 5,303 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 1 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 1 ms
5,376 KB
testcase_04 AC 1 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 2 ms
5,376 KB
testcase_13 AC 26 ms
6,784 KB
testcase_14 AC 19 ms
6,656 KB
testcase_15 AC 27 ms
6,912 KB
testcase_16 AC 20 ms
6,864 KB
testcase_17 AC 21 ms
6,872 KB
testcase_18 AC 23 ms
6,868 KB
testcase_19 AC 20 ms
6,872 KB
testcase_20 AC 20 ms
6,944 KB
testcase_21 AC 25 ms
6,944 KB
testcase_22 AC 19 ms
6,816 KB
testcase_23 AC 24 ms
6,912 KB
testcase_24 AC 25 ms
6,820 KB
testcase_25 AC 21 ms
6,820 KB
testcase_26 AC 21 ms
6,816 KB
testcase_27 AC 24 ms
6,820 KB
testcase_28 AC 48 ms
9,976 KB
testcase_29 AC 56 ms
10,120 KB
testcase_30 AC 60 ms
10,064 KB
testcase_31 AC 52 ms
10,036 KB
testcase_32 AC 52 ms
9,992 KB
testcase_33 AC 62 ms
10,016 KB
testcase_34 AC 59 ms
9,988 KB
testcase_35 AC 67 ms
10,216 KB
testcase_36 AC 70 ms
10,136 KB
testcase_37 AC 65 ms
10,084 KB
testcase_38 AC 45 ms
10,060 KB
testcase_39 AC 62 ms
9,996 KB
testcase_40 AC 51 ms
10,008 KB
testcase_41 AC 74 ms
9,864 KB
testcase_42 AC 69 ms
10,044 KB
testcase_43 AC 62 ms
9,996 KB
testcase_44 AC 38 ms
10,076 KB
testcase_45 AC 61 ms
10,040 KB
testcase_46 AC 63 ms
9,996 KB
testcase_47 AC 50 ms
10,032 KB
testcase_48 AC 10 ms
10,160 KB
testcase_49 AC 96 ms
9,992 KB
testcase_50 AC 33 ms
10,140 KB
testcase_51 AC 57 ms
10,196 KB
testcase_52 AC 54 ms
10,236 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:36:9: warning: '*sum[<unknown>]' may be used uninitialized [-Wmaybe-uninitialized]
   36 |     int all = sum[(1<<N)-1];
      |         ^~~

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define rep(i, n) for(int i = 0; i < n; i++)
#define rep2(i, x, n) for(int i = x; i <= n; i++)
#define rep3(i, x, n) for(int i = x; i >= n; i--)
#define elif else if
#define sp(x) fixed << setprecision(x)
#define pb push_back
#define eb emplace_back
#define all(x) x.begin(), x.end()
#define sz(x) (int)x.size()
using ll = long long;
using ld = long double;
using pii = pair<int, int>;
using pil = pair<int, ll>;
using pli = pair<ll, int>;
using pll = pair<ll, ll>;
const ll MOD = 1e9+7;
//const ll MOD = 998244353;
const int inf = (1<<30)-1;
const ll INF = (1LL<<60)-1;
const ld EPS = 1e-10;
template<typename T> bool chmax(T &x, const T &y) {return (x < y)? (x = y, true) : false;};
template<typename T> bool chmin(T &x, const T &y) {return (x > y)? (x = y, true) : false;};

int main(){
    int N;
    cin >> N;
    int A[N];
    rep(i, N) cin >> A[i];
    int sum[1<<N];
    rep(i, 1<<N){
        sum[i] = 0;
        rep(j, N) if(i&(1<<j)) sum[i] += A[j];
    }
    int all = sum[(1<<N)-1];
    ld dp[1<<N][201][2];
    rep(i, 1<<N) rep(j, 201) rep(k, 2) dp[i][j][k] = 0.0;
    rep(j, 201){
        if(j*2 > all) rep(k, 2) dp[0][j][k] = 1.0;
    }
    rep(i, 1<<N){
        if(i == 0) continue;
        rep(j, all-sum[i]+1){
            ld memo1[N][N], memo2[N][N];
            fill(memo1[0], memo1[N], -1.0), fill(memo2[0], memo2[N], -1.0);
            rep(k, N) rep(l, N){
                if(!(i&(1<<k)) || !(i&(1<<l))) continue;
                ld s = 1.0/ld(A[k]), t = 1.0/ld(A[l]);
                ld p = dp[i^(1<<k)][j+A[k]][1], q = dp[i^(1<<l)][j][0];
                ld x = (s*p+(1-s)*t*q)/(s+t-s*t);
                ld y = t*q+(1-t)*x;
                memo1[k][l] = x, memo2[k][l] = y;
            }
            int res[N];
            rep(k, N){
                res[k] = -1;
                rep(l, N){
                    if(memo2[k][l] < -EPS) continue;
                    if(res[k] == -1) res[k] = l;
                    if(memo2[k][l] < memo2[k][res[k]]) res[k] = l;
                }
            }
            int a = -1;
            rep(k, N){
                if(res[k] == -1) continue;
                if(a == -1) a = k; 
                if(memo1[k][res[k]] > memo1[a][res[a]]) a = k;
            }
            dp[i][j][0] = memo1[a][res[a]], dp[i][j][1] = memo2[a][res[a]];
        }
    }
    cout << sp(10) << dp[(1<<N)-1][0][0] << endl;
}
0