結果

問題 No.162 8020運動
コンテスト
ユーザー GOTKAKO
提出日時 2025-12-19 23:13:31
言語 C++17
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 793 ms / 5,000 ms
コード長 1,432 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,055 ms
コンパイル使用メモリ 203,360 KB
実行使用メモリ 80,512 KB
最終ジャッジ日時 2025-12-19 23:13:55
合計ジャッジ時間 22,960 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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

int main(){
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);

    int A; cin >> A;
    double p0,p1,p2; cin >> p0 >> p1 >> p2;
    p0 /= 100,p1 /= 100,p2 /= 100;

    int n = 14,n2 = 1<<n;
    vector<double> dp(n2);
    vector<vector<pair<int,double>>> trans(n2);
    for(int i=0; i<n2; i++){
        for(int k=i; ; k=(k-1)&i){
            double mul = 1;
            for(int p=0; p<n; p++) if(i&(1<<p)){
                int ok = (i>>(p+1))&1;
                if(p) ok += (i>>(p-1))&1;
                if(k&(1<<p)){
                    if(ok == 0) mul *= (1-p0);
                    else if(ok == 1) mul *= (1-p1);
                    else mul *= (1-p2);
                }
                else{
                    if(ok == 0) mul *= p0;
                    else if(ok == 1) mul *= p1;
                    else mul *= p2;
                }
            }
            trans.at(i).push_back({k,mul});
            if(k == 0) break;
        }
    }
    dp.at(n2-1) = 2;
    while(A < 80){
        vector<double> next(n2);
        for(int i=0; i<n2; i++){
            if(dp.at(i) == 0) continue;
            double d = dp.at(i);
            for(auto [k,mul] : trans.at(i)) next.at(k) += d*mul;
        }
        swap(dp,next),A++;
    }
    double answer = 0;
    for(int i=0; i<n2; i++) answer += __popcount(i)*dp.at(i);
    cout << fixed << setprecision(20) << answer << endl;
}
0