結果

問題 No.242 ビンゴゲーム
コンテスト
ユーザー t98slider
提出日時 2022-06-21 08:38:13
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
MLE  
実行時間 -
コード長 1,351 bytes
コンパイル時間 1,805 ms
コンパイル使用メモリ 171,596 KB
実行使用メモリ 532,816 KB
最終ジャッジ日時 2024-10-14 01:32:02
合計ジャッジ時間 8,513 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample MLE * 2 -- * 1
other -- * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

int main(){
    int n;
    cin >> n;
    vector<double> dp(1 << 25);
    dp[0] = 1;
    for(int i = 0; i < n; i++){
        int rem = 99 - i;
        vector<double> dp2(1 << 25);
        for(int j = 0; j < (1 << 25); j++){
            if(dp[j] == 0)continue;
            int cnt = 25 - __builtin_popcount(j);
            dp2[j] += dp[j] * (rem - cnt) / rem;
            if(cnt == 0)continue;
            for(int k = 0; k < 25; k++){
                if(j >> k & 1)continue;
                dp2[j | (1 << k)] += dp[j] * cnt / rem;
            }
        }
        swap(dp, dp2);
    }
    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[i];
    }
    printf("%.15lf\n", ans);
}
0