UVA11021 Tribbles

题目大意：n个麻球，第一天有k个，麻球生命期为一天，临近死亡前会有i的几率生出Pi个麻球。问m天后麻球全部死亡概率

设f[i]表示i天后一个麻球全部死亡的概率

有f[1] = P0

f[i] = P0 + P1 * f[1] + P2 * f[2]^2 + ... + Pi * f[i]^i + ... +Pn * f[n] ^ n

即：一个麻球在第一天结束生i个麻球，i个麻球会从第二天开始执行第一天的麻球的行为，因此i天后全部死亡对第二天的麻球来说就是i-1天后全部死亡。由于麻球死亡相互独立，只需i次方即可

最终答案f[m]^k

蛇皮自己写的pow被卡精度

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <cstdlib>
 #include <algorithm>
 #include <queue>
 #include <vector>
 #include <cmath>
 #define min(a, b) ((a) < (b) ? (a) : (b))
 #define max(a, b) ((a) > (b) ? (a) : (b))
 #define abs(a) ((a) < 0 ? (-1 * (a)) : (a))
 inline void swap(int &a, int &b)
 {
     long long tmp = a;a = b;b = tmp;
 }
 inline void read(int &x)
 {
     x = ;char ch = getchar(), c = ch;
     while(ch < '' || ch > '') c = ch, ch = getchar();
     while(ch <= '' && ch >= '') x = x *  + ch - '', ch = getchar();
     if(c == '-') x = -x;
 }
 const int INF = 0x3f3f3f3f;
 const int MAXN =  + ;
 int t,n,k,m;
 double p[MAXN], dp[MAXN];
 /*double pow(double a, int b)
 {
     double r = 1, base = a;
     for(;b;b >>= 1)
     {
         if(b & 1) r *= base;
         base *= base;
     }
     return r;
 }*/
 int main()
 {
     read(t);
     int ca = ;
     for(;t;-- t)
     {
         ++ ca;
         memset(dp, , sizeof(dp));
         read(n), read(k), read(m);
         -- n;
         for(register int i = ;i <= n;++ i) scanf("%lf", &p[i]);
         dp[] = p[];
         for(register int i = ;i <= m;++ i)
             for(register int j = ;j <= n;++ j)
                 dp[i] += p[j] * pow(dp[i - ], j);
         printf("Case #%d: %.7lf\n", ca, pow(dp[m], k));
     }
     return ;
 } 

UVA11021