poj 2151 Check the difficulty of problems 概率DP

简单的概率DP！！！

结果=所以的情况-所以满足且<=n-1的情况！

代码如下：

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<cstring>
#include<vector>
#define ll __int64
#define pi acos(-1.0)
#define MAX 1001
using namespace std;
double p[MAX][32],dp[50][32],s[MAX][32],p1,p2;
int main(){
    int n,m,t,i,j,k;
    while(cin>>m>>t>>n){
        if(m==0&&t==0&&n==0) break;
        for(i=1;i<=t;i++)
            for(j=1;j<=m;j++)
                cin>>p[i][j];
        for(i=1;i<=t;i++){
            dp[0][0]=1;
            for(j=1;j<=m;j++) dp[j][0]=dp[j-1][0]*(1-p[i][j]);
            for(j=1;j<=m;j++)
                for(k=1;k<=j;k++){
                    dp[j][k]=dp[j-1][k-1]*p[i][j]+dp[j-1][k]*(1-p[i][j]);
                }
            s[i][0]=dp[m][0];
            for(j=1;j<=m;j++) s[i][j]=s[i][j-1]+dp[m][j];
        }
        p1=p2=1;
        for(i=1;i<=t;i++){
            p1*=(1-s[i][0]);
            p2*=(s[i][n-1]-s[i][0]);
        }
        printf("%.3lf
",p1-p2);
    }
    return 0;
}