結果

問題 No.1151 チャレンジゲーム
コンテスト
ユーザー tokusakurai
提出日時 2020-08-08 09:26:19
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 122 ms / 2,000 ms
コード長 2,409 bytes
コンパイル時間 2,256 ms
コンパイル使用メモリ 196,192 KB
最終ジャッジ日時 2025-01-12 18:46:57
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 50
権限があれば一括ダウンロードができます
コンパイルメッセージ
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