poj2151

这道题目的意思很简单，有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道题的概率）

下面是代码：

#include <stdio.h>
#include <stdlib.h>
#define Max 100
float dp[Max][Max];
int m,t,n;
int main(){
	while(scanf("%d%d%d",&m,&t,&n)){
		if(m==0 && t==0 && n==0)
			break;
		float pro1=1.0,pro2=1.0,temp;
		float result=1.0;
		while(t--){
			dp[0][0]=1.0;
			for(int i=1;i<=m;i++){
			    scanf("%f",&temp);
				pro1*=(1-temp);
				for(int j=0;j<=i;j++){
					float cov1,cov2;
					if(j-1<0)
						cov1=0.0;
					else
						cov1=dp[i-1][j-1]*temp;
					if(j>i-1)
						cov2=0.0;
					else
						cov2=dp[i-1][j]*(1-temp);
				    dp[i][j]=cov1+cov2;
				}
			}
			    result*=(1-pro1);
				pro1=1.0;
				float ac=0.0;
			    for(int i=1;i<n;i++)
				    ac+=dp[m][i];
				pro2*=ac;
		}
		result=result-pro2;
		printf("%.3f\n",result);
	}
	return 0;
}