結果

問題 No.162 8020運動
ユーザー mayoko_mayoko_
提出日時 2015-03-08 20:18:33
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 172 ms / 5,000 ms
コード長 1,925 bytes
コンパイル時間 843 ms
コンパイル使用メモリ 76,776 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-24 16:26:16
合計ジャッジ時間 4,307 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 4 ms
6,816 KB
testcase_01 AC 56 ms
6,944 KB
testcase_02 AC 104 ms
6,944 KB
testcase_03 AC 151 ms
6,944 KB
testcase_04 AC 150 ms
6,940 KB
testcase_05 AC 148 ms
6,944 KB
testcase_06 AC 151 ms
6,944 KB
testcase_07 AC 155 ms
6,940 KB
testcase_08 AC 156 ms
6,944 KB
testcase_09 AC 4 ms
6,940 KB
testcase_10 AC 103 ms
6,944 KB
testcase_11 AC 56 ms
6,944 KB
testcase_12 AC 24 ms
6,944 KB
testcase_13 AC 172 ms
6,944 KB
testcase_14 AC 10 ms
6,940 KB
testcase_15 AC 8 ms
6,940 KB
testcase_16 AC 29 ms
6,940 KB
testcase_17 AC 13 ms
6,940 KB
testcase_18 AC 136 ms
6,940 KB
testcase_19 AC 95 ms
6,944 KB
testcase_20 AC 109 ms
6,944 KB
testcase_21 AC 106 ms
6,940 KB
testcase_22 AC 97 ms
6,944 KB
testcase_23 AC 109 ms
6,940 KB
testcase_24 AC 100 ms
6,940 KB
testcase_25 AC 118 ms
6,944 KB
testcase_26 AC 4 ms
6,940 KB
testcase_27 AC 64 ms
6,940 KB
testcase_28 AC 138 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <sstream>
#include <string>
#include <vector>
#include <map>
#include <algorithm>
#include <iostream>
#include <utility>
#include <set>
#include <cctype>
#include <queue>
#include <stack>
#include <cstdio>
#include <cstdlib>
#include <cmath>

using namespace std;
typedef long long ll;

double memo[25][20]; // memo[Y][R]はあとY年残ってる時R本の歯が連続してたらそのうち何本残るかを返す
double P[3];

double dfs(int year, int rest) {
    if (year == 0) return (double)rest;
    if (rest <= 0) return 0;
    if (memo[year][rest] >= 0) return memo[year][rest];
    if (rest == 1) {
        return (1-P[0]) * dfs(year-1, rest);
    } else {
        double ret = 0;
        for (int s = 0; s < (1<<rest); s++) {
            // まずそのような状態になる確率を求める
            double p = 1;
            for (int i = 0; i < rest; i++) {
                if ((s>>i)&1) {
                    if (i == 0) p *= (1-P[1]);
                    else if (i == rest-1) p *= (1-P[1]);
                    else p *= (1-P[2]);
                } else {
                    if (i == 0) p *= P[1];
                    else if (i == rest-1) p *= P[1];
                    else p *= P[2];
                }
            }
            // 次に分割分だけdfsして期待値を出す
            for (int cur = 0; cur < rest; cur++) {
                if ((s>>cur)&1) {
                    int i = cur;
                    while ((s>>i)&1) i++;
                    ret += p*dfs(year-1, i-cur);
                    cur = i;
                }
            }
        }
        return memo[year][rest] = ret;
    }
}

int main(void) {
    int A;
    cin >> A;
    for (int i = 0; i < 25; i++) for (int j = 0; j < 20; j++) {
        memo[i][j] = -1;
    }
    for (int i = 0; i < 3; i++) {
        cin >> P[i];
        P[i] /= 100.0;
    }
    printf("%.10lf\n", dfs(80-A, 14) * 2);
    return 0;
}
0