UVA 11021繁衍模型+概率计算
1 /*UVA 11021 2 繁衍问题: 3 总结:不同个体都可以看成相互独立的事件,即每个个体后代单独存货 4 输入: n (1<= n<=1000) , k (0<= k<=1000) and m (0<= m<=1000) 5 n:可能生产数 k:初始生命体数 m:m天都全部死亡 6 p[0] ,p[1]...p[n-1]:p[i]产生i个后代的概率 7 输出:m天后全部死亡的概率 8 递推公式:f[x]:1个生命体x天后死亡的概率 9 kf[k,x]:k个生命体x天后死亡的概率 10 f[x]=p[0]+p[1]*kf[1,x-1]+p[2]*kf[2,x-1].....p[n-1]*kf[n-1,x-1]; 11 kf[k,x]=f[x]^k;相互独立 12 边界条件:f[1]=p[0]; 13 */ 14 #include<iostream> 15 #include<stdio.h> 16 #include<string.h> 17 #include<algorithm> 18 #include<stdlib.h> 19 #include<math.h> 20 #include<queue> 21 #include<vector> 22 #include<map> 23 24 using namespace std; 25 26 27 int k,n,m; 28 double p[1005],f[1005]; 29 /*double f(double d) 30 { 31 if (d==1) return p[0]; 32 double ans=p[0]; 33 for(int i=1;i<=n-1;i++) ans+=p[i]*pow(f(d-1),i); 34 return ans; 35 } 迭代次数多,RE*/ //注意 36 int main() 37 { 38 int t; 39 cin>>t; 40 for(int cas=1;cas<=t;cas++) 41 { 42 cin>>n>>k>>m; 43 for(int i=0;i<n;i++) cin>>p[i]; 44 f[1]=p[0]; 45 for(int i=2;i<=m;i++) 46 { 47 f[i]=p[0]; 48 for(int j=1;j<=n-1;j++) f[i]+=p[j]*pow(f[i-1],j); 49 } 50 cout<<"Case #"<<cas<<": "; 51 52 printf("%.7lf\n",pow(f[m],k)); 53 } 54 }