poj 2151 Check the difficulty of problems

   需要一定的概率基础，然后加上一定的dp基础就可以直接解决了，我概率也都忘得差不多了，竟然最后连pn是p1的子问题都分不清了，哎，糊涂啊。   i队做一题的概率非常好求（1减去一道题都没做出来的概率），然后所有队的概率相乘就是每队都至少做出一道题的概率，然后算出至少做出1道题，但不存在一个队做出的题数大于等于N的概率，然后2者相减（pn是p2的子问题），算pn的时候需要用到dp，主要的dp状态转移方程为：（dp[i][j]表示i道题做出j道的概率。）

for(int j=1;j<=M;++j)

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;

int M,T,N;
double p[1010][35],p1[1010],pn[1010];
double dp[35][35];

int main(void)
{
    while(scanf("%d %d %d",&M,&T,&N),M!=0||T!=0||N!=0)
    {
    for(int i=1;i<=T;++i)
     for(int j=1;j<=M;++j)
        scanf("%lf",p[i]+j);
     double ans = 1;
     for(int i=1;i<=T;++i)//算出每队至少做出来一题的概率
     {
         double tmp = 1;
         for(int j=1;j<=M;++j)
          tmp *= (1-p[i][j]);
         ans *= (1-tmp);
     }
     dp[0][0] = 1;
     for(int i=1;i<=T;++i)
     {
         for(int j=1;j<=M;++j)
         for(int k=0;k<=j;++k)
         {
             if(j==k)
             {
                dp[j][k] = dp[j-1][k-1]*p[i][j];
             }
             else if(k==0)
             {
                dp[j][0] = dp[j-1][0]*(1-p[i][j]);
             }
             else
             {
                 dp[j][k] = dp[j-1][k]*(1-p[i][j])+dp[j-1][k-1]*p[i][j];
             }
         }
         pn[i] = 0;
         for(int ii=1;ii<N;++ii)
         pn[i] += dp[M][ii];
     }
     double tmp = 1;
     for(int i=1;i<=T;++i)
     {
        tmp *= pn[i];
     }
     ans -= tmp;
     printf("%.3lf\n",ans);
    }

    return 0;
}

         for(int k=0;k<=j;++k)
         {
             if(j==k)
             {
                dp[j][k] = dp[j-1][k-1]*p[i][j];
             }
             else if(k==0)
             {
                dp[j][0] = dp[j-1][0]*(1-p[i][j]);
             }
             else
             {
                 dp[j][k] = dp[j-1][k]*(1-p[i][j])+dp[j-1][k-1]*p[i][j];
             }
         }