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

#include<stdio.h>
#include<string.h>

double dp[1010][31][31],p[1010][31],S[1010][31];

int main()
{
    int n,m,T;
    while(~scanf("%d%d%d",&m,&n,&T)){
        if(m==0&&n==0&&T==0) break;
    memset(dp,0,sizeof(dp));
    memset(S,0,sizeof(S));
    for(int i=1 ; i<=n ; i++)
    for(int j=1 ; j<=m ; j++)
    scanf("%lf",&p[i][j]);

    //dp[i][j][k]:第i个队伍前j道题目完成了k个

    for(int i=1 ; i<=n ; i++)
    {dp[i][0][0]=1;
        for(int j=1 ; j<=m ; j++)
        {
            dp[i][j][0]=dp[i][j-1][0]*(1-p[i][j]);
        }
    }
    for(int i=1 ; i<=n ; i++)
    {
        for(int j=1 ; j<=m ; j++)
        {
            for(int k=1 ; k<=j ; k++)
            {
                dp[i][j][k]=dp[i][j-1][k-1]*p[i][j] + dp[i][j-1][k]*(1-p[i][j]);
            }
        }
    }
   //  S[i][M]:第i队伍完成M道题的概率
     for(int i=1 ; i<=n ; i++)
     {
         S[i][0]=dp[i][m][0];
         for(int j=1 ; j<=m ; j++)
         {
             S[i][j]=S[i][j-1]+dp[i][m][j];
         }
     }
     double ans1=1;
     for(int i=1 ; i<=n ; i++)
     {
             ans1*=(1-S[i][0]);
     }
     double ans2=1;
     for(int i=1 ; i<=n ; i++)
     {
         ans2*=(S[i][T-1]-S[i][0]);
     }
     printf("%.3f\n",ans1-ans2);
    }
}