結果

問題 No.162 8020運動
ユーザー kmjp
提出日時 2015-02-09 02:36:03
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 118 ms / 5,000 ms
コード長 1,716 bytes
コンパイル時間 1,475 ms
コンパイル使用メモリ 161,632 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-23 16:01:21
合計ジャッジ時間 4,575 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
using namespace std;
double memo[70][20];
int P[3];
// N本の歯がY年後に何本残るか？
double func(int Y, int N) {
    
    // おめでとうございます。80歳までN本残りました。
    if(Y == 0)
        return N;
    
    // メモ参照
    if(memo[Y][N] >= 0)
        return memo[Y][N];
    
    if(N == 1) {
        // 歯が1本の時は確率が独特なので、別処理した方が楽です。
        // 1年残る確率×(Y-1)年残る確率
        memo[Y][N] = (100-P[0])/100.0 * func(Y-1, 1);
    }
    else {
        memo[Y][N] = 0;
        // 歯の壊れ方を2**N通り試す
        for(int mask=0; mask < 1<<N; mask++) {
            //この壊れ方をする確率は？
            double prob=1;
            for(int i = 0; i < N; i++) {
                double bad;
                if(i==0 || i==N-1) //両端だけ確率が異なる
                    bad = P[1]/100.0;
                else //両端以外
                    bad = P[2]/100.0;
                
                if(mask & (1<<i)) //歯が残る
                    prob *= 1-bad;
                else //残念虫歯でした
                    prob *= bad;
            }
            
            //連続する歯の残り方は？
            int left=0;
            for(int i = 0; i < N; i++) {
                if(mask & (1<<i)) {
                    //残っている歯が連続している
                    left++;
                }
                else {
                    //虫歯で残った歯が途切れた
                    if(left > 0)
                        memo[Y][N] += prob * func(Y-1,left);
                    left = 0;
                }
            }
            
            if(left > 0)
                memo[Y][N] += prob * func(Y-1,left);
        }
    }
    
    return memo[Y][N];
}
int main(int argc,char** argv){
    int x,y,T;
    
    cin >> T;
    cin >> P[0] >> P[1] >> P[2];
    
    // メモ初期化
    for(x = 0; x < 70; x++)
        for(y = 0; y < 20; y++)
            memo[x][y] = -1;
    
    // 上の歯と下の歯の期待値は同じなので、片方を求めて２倍すればよい
    printf("%.9lf\n" , 2 * func(80-T, 14));
    
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0