結果

問題 No.162 8020運動
コンテスト
ユーザー kroton
提出日時 2015-03-06 01:46:54
言語 C++11
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 1,749 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 965 ms
コンパイル使用メモリ 182,112 KB
実行使用メモリ 7,844 KB
最終ジャッジ日時 2026-03-10 12:13:06
合計ジャッジ時間 1,759 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge1_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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

double dp[100][100];
double P[10];

double solve(int A, int rem){
    if(A >= 80)return rem;
    if(rem == 0)return 0;

    double &res = dp[A][rem];
    if(res >= -0.5)return res;

    res = 0;

    if(rem == 1){
        res += (1 - P[0]) * solve(A + 1, 1);
        return res;
    }

    for(int s=0;s<rem;s++){
        double p = 1;
        if(s == 0 || s == rem - 1){
            p *= 1 - P[1];
        } else {
            p *= 1 - P[2];
        }

        if(s - 1 == 0){
            p *= P[1];
        }
        if(s - 1 > 0){
            p *= P[2];
        }

        if(s + 1 == rem - 1){
            p *= P[1];
        }
        if(s + 1 < rem - 1){
            p *= P[2];
        }

        res += p * solve(A + 1, 1);
    }

    for(int len=2;len<=rem;len++){
        for(int s=0;s+len<=rem;s++){
            int t = s + len - 1;

            double p = 1;

            if(s - 1 < 0){
                p *= 1 - P[1];
            }
            if(s - 1 == 0){
                p *= P[1] * (1 - P[2]);
            }
            if(s - 1 > 0){
                p *= P[2] * (1 - P[2]);
            }

            if(t + 1 > rem - 1){
                p *= 1 - P[1];
            }
            if(t + 1 == rem - 1){
                p *= (1 - P[2]) * P[1];
            }
            if(t + 1 < rem - 1){
                p *= (1 - P[2]) * P[2];
            }

            for(int k=0;k<len-2;k++)p *= 1 - P[2];

            res += p * solve(A + 1, len);
        }
    }

    return res;
}

int main(){
    memset(dp, -1, sizeof(dp));

    int A;
    cin >> A;

    cin >> P[0] >> P[1] >> P[2];
    P[0] /= 100;
    P[1] /= 100;
    P[2] /= 100;

    printf("%.10f\n", 2 * solve(A, 14));
    return 0;
}
0