概率dp poj 2151

题意：

这道题目的意思很简单，有t个ACM队，m个题目，题目给出了每个队对每个题目做出的概率大小（0到1之间，包含0和1），要求每个队至少做出一道题（签到题），同时，要求获胜队必须至少能够做出n道题（获胜对不止一个），这到题目实际上就是一个dp+概率论知识，dp的推导也是概率论中的全概率公式推导出来的，之后就是如何利用概率论知识求解问题了。首先解释一下dp如何推导。我们令dp[i][j]为前i道题中做出j道题的概率，那么依据全概率公式，现在有两种情况：

1）第i道题做出来的，在此条件下dp[i][j]=dp[i-1][j-1]

2）第i道题没有做出来，在此条件下dp[i][j]=dp[i-1][j]

那么依据全概率公式dp[i][j]=dp[i-1][j-1]*p[i]+dp[i-1][j]*(1-p[i]),其中p[i]为第i题做出的概率。

状态转移方程确定后，那么就可以求解出前i个题目中做出j道题的概率了。

问题是所有的基本条件都已经确定了，那么题目要求解问题的概率该如何计算。这里，很容易认为每个队至少做出一道题和获胜队至少做出n道题是独立事件。其实不然。应该说他们之间存在包含于被包含的关系，这里分两种情况讨论：

我们将每个队都至少作出一道题记为事件A，获胜队至少作出n道题记为事件B，题目要求的就是p(AB)的结果

1）当n=1时，事件B包含事件A，那么p(AB)就转化为p(A),即为每个队都至少作出一道题的概率。

2）当n>1时，事件A包含事件AB，那么p(AB)就是p(A)-p(A-AB),这里补充说明一下，事件A-AB就是所有队都只作出少于n道题的概率（即作出1到n-1道题的概率）
————————————————
版权声明：本文为CSDN博主「阿杜dyh」的原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/i_want_to_be_a_god/article/details/25926243

 #include<cstdio>
 #include<algorithm>
 #include<math.h>
 #include<string.h>
 using namespace std;
 const int maxn=1e3+;
 double dp[maxn][][];  //表示第i只队伍在j道题的情况下赢k局的概率
 double p[maxn][];
 int main()
 {
     int n,t,m;
     while(scanf("%d%d%d",&n,&t,&m)!=EOF){
         if(!(n+t+m)) break;
         for(int i=;i<t;i++)
         for(int j=;j<=n;j++){
             scanf("%lf",&p[i][j]);
         }
         memset(dp,,sizeof(dp));
         for(int i=;i<t;i++){
             dp[i][][]=;
             for(int j=;j<=n;j++){
                 dp[i][j][]=dp[i][j-][]*(-p[i][j]);
                 for(int k=;k<=j;k++){
                     dp[i][j][k]=dp[i][j-][k-]*p[i][j]+dp[i][j-][k]*(-p[i][j]);
                 }
             }
         }
         double ans=;
         for(int i=;i<t;i++)
             ans*=-dp[i][n][];
         double temp=;
         for(int i=;i<t;i++){
             double sum=;
             for(int j=;j<m;j++)
                 sum+=dp[i][n][j];  //计算出做出的题在1~n-1的数量的概率；
             temp*=sum;             //将每一个做不出的概率相乘；
         }
         ans-=temp;         //减去不符合的情况
         printf("%.3f\n",ans);
     }
     return ;
 }