Uva 11021 Tribles(概率dp)

思路：重要的是看清楚题意，给（0–n-1） 个概率分别是细菌死亡后分裂（0–n-1）个细菌的概率。根据题意设dp[i]为一个细菌i天后死绝的概率。那么dp[i]就等于这个细菌不分裂的概率，加上分裂一次的概率*dp[i-1] (就是在i-

1 )天死亡的概率， 加上分裂两次的概率*dp[i-1] *dp[i-1] (两个分裂出来的细菌同时在i-1天内死亡的概率)。

整理一下得到递推式子为：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;


double p[20000];

double dp[20000];

double power(double a, int n)    //快速幂 a^n % p
{
    double ans = 1;
    while(n)
    {
        if(n & 1) ans = ans * a;
        a = a * a ;
        n >>= 1;
    }
    return ans;
}
int main()
{


    int t;
    scanf("%d", &t);
    for(int _ = 1;_<=t;_++)
    {

        int n, k, m;
        scanf("%d%d%d", &n, &k, &m);
        for(int i = 0;i<n;i++) scanf("%lf", &p[i]);
        dp[1] = p[0];

        for(int i = 2;i<=m;i++)
        {
            dp[i] = 0;

            for(int j = 0;j<n;j++)
            {
                dp[i] += p[j] * pow(dp[i-1], j);
                //cout<<"?";
            }
        }



        printf("Case #%d: %.7f\n", _, power(dp[m],k));

    }
    return 0;
}