poj 2151

简单的概率dp，用到补集的概念

1.首先要明确p1="每个队至少过一题 "与p2="至少有一个队过N题"不是两个相互独立的事件，求他们同时发生的概率不能直接p1*p2

2.第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])

3.令ans[i][j]=dp[i][m][0]+...+dp[i][m][j]

问题可以转化为：所有至少做出一道的概率(p1) - 所有选手做出的题数 j>=1 && j<=n 的概率(p2)

5.p1=(ans[1][m]-ans[1][0])*(ans[2][m]-ans[2][0])*...*(ans[t][m]-ans[t][0]); p2=(ans[1][n-1]-ans[1][0])*(ans[2][n-1]-ans[2][0])*...*(ans[t][n-1]-ans[t][0]);

#include <iostream>
using namespace std;

int m,t,n,s[1050][35];
double p[1050][35],dp[1050][35][35];
int main()
{
	int 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])
	while(scanf("%d%d%d",&m,&t,&n)!=EOF)
	{
		if(!m && !t && !n)
			break;
		for(i=1;i<=t;i++)
			for(j=1;j<=m;j++)
				scanf("%lf",&p[i][j]);
		for(i=1;i<=t;i++)
		{
			dp[i][0][0]=1;
			for(j=1;j<=m;j++)
				for(k=0;k<=j;k++)
				{
					dp[i][j][k]=dp[i][j-1][k]*(1-p[i][j]);
					if(k!=0)
						dp[i][j][k]+=dp[i][j-1][k-1]*p[i][j];
				}
		}
		double tmp1,tmp2,tmp,ans;
		tmp=1;
		for(i=1;i<=t;i++)
		{
			tmp1=0;
			for(j=1;j<=m;j++)
				tmp1+=dp[i][m][j];
			tmp*=tmp1;
		}
		ans=tmp;
		tmp=1;
		for(i=1;i<=t;i++)
		{
			tmp1=0;
			for(j=1;j<n;j++)
				tmp1+=dp[i][m][j];
			tmp*=tmp1;
		}
		ans-=tmp;
		printf("%.3lf\n",ans);	
	}
	return 0;
}