POJ 2151 Check the difficulty of problems(概率DP）

题目链接

不会做，概率论学的不好啊。看的这个题解，这里面好像有点错误，sum[i][j]数组，表示第i个队伍解出0,1,2...j的题目概率的和。这个问题主要是转化，和怎么计算，状态转移很常见，如果想不出怎么计算，肯定挂。初始化问题，让我错了几次，这是多久没写题了啊。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <map>
#include <cstring>
#include <cmath>
using namespace std;
#define eps 1e-7
double dp[1001][31][31];
double p[1001][31];
double sum[1001][31];
int main()
{
    int n,m,t,i,j,k;
    double ans,temp;
    while(scanf("%d%d%d",&m,&t,&n)!=EOF)
    {
        if(!m&&!t&&!n) break;
        ans = 1;
        temp = 1;
        memset(dp,0,sizeof(dp));
        memset(sum,0,sizeof(sum));
        memset(p,0,sizeof(p));
        for(i = 1; i <= t; i ++)
        {
            for(j = 1; j <= m; j ++)
            {
                scanf("%lf",&p[i][j]);
            }
        }
        for(i = 1; i <= t; i ++)
        {
            dp[i][1][1] = p[i][1];
            dp[i][1][0] = 1 - p[i][1];
            for(j = 2; j <= m; j ++)
            {
                dp[i][j][0] = dp[i][j-1][0]*(1 - p[i][j]);
                for(k = 1; k <= j; k ++)
                {
                    dp[i][j][k] = (1 - p[i][j])*dp[i][j-1][k] + p[i][j]*dp[i][j-1][k-1];
                }
            }
        }
        for(i = 1; i <= t; i ++)
        {
            sum[i][0] = dp[i][m][0];
            for(j = 1; j <= m; j ++)
            {
                sum[i][j] = dp[i][m][j] + sum[i][j-1];
            }
        }
        for(i = 1; i <= t; i ++)
        {
            ans = ans*(sum[i][m] - sum[i][0]);
        }
        for(i = 1; i <= t; i ++)
        {
            temp = temp*(sum[i][n-1] - sum[i][0]);
        }
        printf("%.3lf\n",ans-temp+eps);
    }
    return 0;
}