UVA 11021繁衍模型+概率计算

/*UVA 11021
繁衍问题：
总结：不同个体都可以看成相互独立的事件，即每个个体后代单独存货
输入： n (1<= n<=1000) , k (0<= k<=1000) and m (0<= m<=1000)
       n：可能生产数     k：初始生命体数     m：m天都全部死亡
       p[0] ,p[1]...p[n-1]:p[i]产生i个后代的概率
输出：m天后全部死亡的概率
递推公式：f[x]:1个生命体x天后死亡的概率
          kf[k,x]:k个生命体x天后死亡的概率
          f[x]=p[0]+p[1]*kf[1,x-1]+p[2]*kf[2,x-1].....p[n-1]*kf[n-1,x-1];
          kf[k,x]=f[x]^k;相互独立
边界条件：f[1]=p[0];
*/
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<stdlib.h>
#include<math.h>
#include<queue>
#include<vector>
#include<map>

using namespace std;


int k,n,m;
double p[1005],f[1005];
/*double f(double d)
{
    if (d==1) return p[0];
    double ans=p[0];
    for(int i=1;i<=n-1;i++) ans+=p[i]*pow(f(d-1),i);
    return ans;
} 迭代次数多，RE*/    //注意
int main()
{
    int t;
    cin>>t;
    for(int cas=1;cas<=t;cas++)
    {
        cin>>n>>k>>m;
        for(int i=0;i<n;i++) cin>>p[i];
        f[1]=p[0];
        for(int i=2;i<=m;i++)
        {
            f[i]=p[0];
            for(int j=1;j<=n-1;j++) f[i]+=p[j]*pow(f[i-1],j);
        }
        cout<<"Case #"<<cas<<": ";

        printf("%.7lf
",pow(f[m],k));
    }
}