一道高中就应该掌握的概率题,却想了很久而不得其解。不是因为求概率难,而是不知道怎么用算法实现。没想到dp可以如此轻松的解决刚开始思想上的种种难关,dp真的很强大啊。(ps:这道题很好,值得认真的做一下!)
#include<iostream>
#include<cstdio>
using namespace std;
double team[1001][31];
double dp[1001][31][31];//dp[i][j][k]表示第i个队前j道题做出k道的概率
int m,t,n;
int main()
{
int i,j,k;
double ans1,ans2,tmp;
while(scanf("%d%d%d",&m,&t,&n),m+t+n)
{
for(i=1;i<=t;i++)
{
for(j=1;j<=m;j++)
scanf("%lf",&team[i][j]);
}
for(i=1;i<=t;i++)//初始化做0道的概率
{
dp[i][0][0]=1.0;
for(j=1;j<=m;j++)
dp[i][j][0]=(1-team[i][j])*dp[i][j-1][0];
}
for(i=1;i<=t;i++)
{
for(j=1;j<=m;j++)
for(k=1;k<=j;k++)
dp[i][j][k]=dp[i][j-1][k-1]*team[i][j]+dp[i][j-1][k]*(1-team[i][j]);
}
ans1=1.0;//每队至少做一道的概率
for(i=1;i<=t;i++)
ans1*=(1-dp[i][m][0]);
ans2=1.0;//每队做1~n-1道的概率
for(i=1;i<=t;i++)
{
tmp=0.0;
for(j=1;j<n;j++)
tmp+=dp[i][m][j];
ans2*=tmp;
}
printf("%.3lf\n",ans1-ans2);
}
return 0;
}