結果

問題 No.242 ビンゴゲーム
ユーザー t98slidert98slider
提出日時 2022-06-21 08:35:51
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
MLE  
実行時間 -
コード長 1,349 bytes
コンパイル時間 1,626 ms
コンパイル使用メモリ 174,484 KB
実行使用メモリ 814,796 KB
最終ジャッジ日時 2024-10-14 01:28:48
合計ジャッジ時間 4,339 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 MLE -
testcase_01 -- -
testcase_02 -- -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

int main(){
    int n;
    cin >> n;
    vector<vector<double>> dp(n + 1, vector<double>(1 << 25));
    dp[0][0] = 1;
    for(int i = 0; i < n; i++){
        int rem = 99 - i;
        for(int j = 0; j < (1 << 25); j++){
            if(dp[i][j] == 0)continue;
            int cnt = 25 - __builtin_popcount(j);
            dp[i + 1][j] += dp[i][j] * (rem - cnt) / rem;
            if(cnt == 0)continue;
            for(int k = 0; k < 25; k++){
                if(j >> k & 1)continue;
                dp[i + 1][j | (1 << k)] += dp[i][j] * cnt / rem;
            }
        }
    }
    vector<int> cand;
    int S = 0b11111;
    cerr << (bitset<25>(S));
    for(int i = 0; i < 5; i++){
        cand.push_back(S << 5 * i);
    }
    S = 0b0000100001000010000100001;
    for(int i = 0; i < 5; i++){
        cand.push_back(S << i);
    }
    cand.push_back(0);
    for(int i = 0; i < 5; i++){
        cand.back() |= 1 << i * 5 + i;
    }
    cand.push_back(0);
    for(int i = 0; i < 5; i++){
        cand.back() |= 1 << i * 5 + 5 - i;
    }
    double ans = 0;
    for(int i = 0; i < (1 << 25); i++){
        if(__builtin_popcount(i) < 5)continue;
        int cnt = 0;
        for(auto &&S: cand){
            cnt += ((S & i) == S);
        }
        ans += cnt * dp[n][i];
    }
    printf("%.15lf\n", ans);
}
0