UVA 11021 麻球繁衍

Posted 2022-11-23 pangbi

题意：

     有K只麻球，每只生存一天就会死亡，每只麻球在死之前有可能生下一些麻球，生i个麻球的概率是pi，问m天后所有的麻球都死亡的概率是多少？

思路：
      涉及到全概率公式，因为麻球的各种活动都互不影响，所以现在只考虑一直麻球，我们假设f[i]是第i天全部都死亡的概率，那么f[i] = p0 + p1*f[i-1] + p2*f[i-1]^2 + ...pn-1*f[i - 1]^(n-1)，其中pjf(i-1)^j的含义是这个麻球生了j个后代，他们在i-1天后全部死亡，注意这j个后代的死亡是独立的，而每个死亡的概率都是f(i-1)^j也就是用前一天的全部死亡概率来代替今天的每一只死亡的概率，又因为今天的每只的生死概率什么的都是独立的，所以p2*f[i-1]^2可以理解成剩下2只，然后两只都死了，这样最后在第m天死光的概率就是f[m]，但是这个只是一只麻球的，所有麻球都死光是f[m]^k。

 1 #include <iostream>
 2 #include <cstring>
 3 #include <cstdio>
 4 #include <algorithm>
 5 #include <cmath>
 6 #define  maxn 1010
 7 using namespace std;
 8 int n,k,m;
 9 double p[maxn],f[maxn];
10 int main()
11     int T;
12     int i,j,k;
13     scanf("%d",&T);
14     for(int cas=1;cas<=T;cas++)
15         scanf("%d%d%d",&n,&k,&m);
16         for(i=0;i<n;i++)scanf("%lf",&p[i]);
17         f[0]=0;
18         f[1]=p[0];
19         for(i=2;i<=m;i++)
20             f[i]=0;
21             for(int j=0;j<n;j++)
22                 f[i]+=p[j]*pow(f[i-1],j);
23              
24         
25         printf("Case #%d: %.7lf\n",cas,pow(f[m],k));
26         
27         
28     
29     return 0;
30