poj 2151 Check the difficulty of problems

　　题意：有m到题，t个队伍，给出第 i 支队作出 第 j 道题的概率 Pij。求每个队至少作出一道题，作出最多的那个队伍至少作出 n 道题的概率。

　　dp[i][j][k]表示第 i 支队前 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][k]表示第 i 支队作出题目个数 0~k 的概率。

　　s[i][k]=dp[i][m][0]+dp[i][m][1]+...+dp[i][m][k].

　　s[i][m]-s[i][0]表示第 i 支队至少作出一道题的概率，s[i][n-1]-s[i][0]表示第 i 支队作出题目个数 为1~n-1的概率。

　

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

int m,t,n;
double dp[1001][35][35],s[1001][35],p[1001][35];
int main()
{
    while(scanf("%d%d%d",&m,&t,&n) && (m+t+n))
    {
        memset(dp,0,sizeof(dp));
        memset(s,0,sizeof(s));
        for(int i=1;i<=t;i++)
            for(int j=1;j<=m;j++)
                scanf("%lf",&p[i][j]);
        for(int i=1;i<=t;i++)   dp[i][0][0]=1.0;
        for(int i=1;i<=t;i++)
        {
            for(int j=1;j<=m;j++)
                dp[i][j][0]=dp[i][j-1][0]*(1-p[i][j]);
            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][0]=dp[i][m][0];
            for(int k=1;k<=m;k++)
                s[i][k]=s[i][k-1]+dp[i][m][k];
        }
        double p1=1.0;
        for(int i=1;i<=t;i++)
            p1*=(s[i][m]-s[i][0]);
        double p2=1.0;
        for(int i=1;i<=t;i++)
            p2*=(s[i][n-1]-s[i][0]);
        printf("%.3lf\n",p1-p2);
    }
    return 0;
}