UVa 11021 Tribles / 概率DP

k只麻球 每活一天就会死 死之前有pi的概率生出i只麻球 求m天死光的的概率

f[i] = p0 + p1*dp(i-1) + p2*dp(i-1)^2+....+pn-1*dp(i-1)^n-1

好像是倒着推的

pj*dp(i-1)^j 表示1个麻球生了j个麻球在i-1天后死光 1个麻球在i-1天死的概率为dp(i-1) j个麻球就是j个dp(i-1)相乘

#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
const int maxn = 1010;
double dp[maxn];
double p[maxn];

int n, m, k;
int main()
{
	int T;
	int cas = 1;
	scanf("%d", &T);
	while(T--)
	{
		scanf("%d %d %d", &n, &k, &m);
		for(int i = 0; i < n; i++)
			scanf("%lf", &p[i]);
		dp[0] = 0;
		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);
		}
		printf("Case #%d: %.7lf\n", cas++, pow(dp[m], k));
	}
	return 0;
}