結果

問題 No.162 8020運動
ユーザー mayoko_mayoko_
提出日時 2015-03-08 20:18:33
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 134 ms / 5,000 ms
コード長 1,925 bytes
コンパイル時間 1,034 ms
コンパイル使用メモリ 73,456 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-09-06 22:15:44
合計ジャッジ時間 3,929 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 4 ms
4,380 KB
testcase_01 AC 51 ms
4,384 KB
testcase_02 AC 86 ms
4,380 KB
testcase_03 AC 133 ms
4,380 KB
testcase_04 AC 133 ms
4,376 KB
testcase_05 AC 134 ms
4,376 KB
testcase_06 AC 132 ms
4,376 KB
testcase_07 AC 132 ms
4,376 KB
testcase_08 AC 133 ms
4,380 KB
testcase_09 AC 4 ms
4,380 KB
testcase_10 AC 86 ms
4,380 KB
testcase_11 AC 51 ms
4,380 KB
testcase_12 AC 22 ms
4,384 KB
testcase_13 AC 133 ms
4,384 KB
testcase_14 AC 7 ms
4,380 KB
testcase_15 AC 8 ms
4,376 KB
testcase_16 AC 27 ms
4,380 KB
testcase_17 AC 12 ms
4,380 KB
testcase_18 AC 122 ms
4,380 KB
testcase_19 AC 79 ms
4,376 KB
testcase_20 AC 94 ms
4,376 KB
testcase_21 AC 94 ms
4,376 KB
testcase_22 AC 86 ms
4,376 KB
testcase_23 AC 94 ms
4,376 KB
testcase_24 AC 86 ms
4,380 KB
testcase_25 AC 103 ms
4,376 KB
testcase_26 AC 4 ms
4,380 KB
testcase_27 AC 58 ms
4,380 KB
testcase_28 AC 122 ms
4,380 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