poj 2151 概率DP

http://poj.org/problem?id=2151

概率还是不行。简单的一个式子还是能够推的

可是方程构造以及思路  全然不行啊....

參考小优的博客http://blog.csdn.net/lyy289065406/article/details/6648579

感觉思路比較好的

1、dp三维设定的意义

2、s[i][j]推出Ans  好精妙

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;
#define CL(a,b) memset(a,b,sizeof(a))
#define IN(s) freopen(s,"r",stdin)

const int MAXT = 1000+5;
const int MAXN = 35;

int n,m,t;
double p[1000+5][MAXN];
double dp[MAXT][MAXN][MAXN];
double s[MAXT][MAXN];

double solve()
{
    CL(dp,0);CL(s,0);
    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]=1.0;
            for(int k=1;k<=j;k++)dp[i][j][0]*=(1-p[i][k]);
            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 j=1;j<=m;j++)
            for(int k=0;k<=j;k++)
                s[i][j]+=dp[i][m][k];      //第i支队在前m题做出0~j题的概率
    }

    double ret1=1.0,ret2=1.0;
    for(int i=1;i<=t;i++)
    {
        ret1*=(s[i][m]-s[i][0]);
        ret2*=(s[i][n-1]-s[i][0]);
    }
    return ret1-ret2;
}

int main()
{
    //IN("poj2151.txt");
    while(~scanf("%d%d%d",&m,&t,&n) && m+t+n)
    {
        for(int i=1;i<=t;i++)
            for(int j=1;j<=m;j++)
                scanf("%lf",&p[i][j]);
        printf("%.3lf\n",solve());
    }
    return 0;
}