結果

問題 No.162 8020運動
ユーザー krotonkroton
提出日時 2015-03-06 01:46:54
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 1,749 bytes
コンパイル時間 1,406 ms
コンパイル使用メモリ 161,912 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-24 09:52:27
合計ジャッジ時間 2,242 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 2 ms
6,944 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 2 ms
6,944 KB
testcase_08 AC 2 ms
6,944 KB
testcase_09 AC 2 ms
6,944 KB
testcase_10 AC 2 ms
6,940 KB
testcase_11 AC 2 ms
6,940 KB
testcase_12 AC 2 ms
6,940 KB
testcase_13 AC 2 ms
6,940 KB
testcase_14 AC 1 ms
6,940 KB
testcase_15 AC 1 ms
6,940 KB
testcase_16 AC 2 ms
6,940 KB
testcase_17 AC 2 ms
6,940 KB
testcase_18 AC 2 ms
6,944 KB
testcase_19 AC 2 ms
6,944 KB
testcase_20 AC 2 ms
6,944 KB
testcase_21 AC 2 ms
6,944 KB
testcase_22 AC 2 ms
6,944 KB
testcase_23 AC 2 ms
6,944 KB
testcase_24 AC 2 ms
6,944 KB
testcase_25 AC 1 ms
6,944 KB
testcase_26 AC 2 ms
6,940 KB
testcase_27 AC 2 ms
6,940 KB
testcase_28 AC 2 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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