題目大意：
給你m,t,n，代表t個隊伍，做m道題，對於所有的隊伍，全部做對1道或者以上的題目，並且至少有1個隊伍做對n道題目的概率。
做法：
p[i][j]代表第i個隊伍做對第j道題目的概率。
dp[i][j][k]代表第i個隊伍，對於前j道題目做對k道題目的概率。
所有隊伍都做對1道及以上題目的概率bans*=(1-dp[i][m][0])//1-隊伍做對0道題目的概率的乘積。
對於所有隊伍都做對一道以及以上題目的概率的情況下：
至少有一個隊伍做對n到題目的概率=所有隊伍都做對1道題目的概率-所有隊伍都做對小於n道題目的概率。
代碼：
[html]
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<queue>
#include<stack>
#include<map>
#include<string>
#include<stdlib.h>
#define INF_MAX 0x7fffffff
#define INF 999999
#define max3(a,b,c) (max(a,b)>c?max(a,b):c)
#define min3(a,b,c) (min(a,b)<c?min(a,b):c)
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
struct node
{
int u;
int v;
int w;
bool friend operator < (node a, node b){
return a.w < b.w;
}
}edge[1001];
int gcd(int n,int m){if(n<m) swap(n,m);return n%m==0?m:gcd(m,n%m);}
int lcm(int n,int m){if(n<m) swap(n,m);return n/gcd(n,m)*m;}
double p[1001][51];
double dp[1001][31][31];
int main()
{
int m,t,n,i,j,k;
while(scanf("%d%d%d",&m,&t,&n)&&(m||n||t))
{
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++)
{
dp[i][0][j]=0;
}
}
for(i=1;i<=t;i++)
{
for(j=1;j<=m;j++)
{
dp[i][j][0]=dp[i][j-1][0]*(1-p[i][j]);
for(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];
}
}
}
double bans;//所有人都做對1道及以上題目的概率
bans=1;
for(i=1;i<=t;i++)
{
bans*=(1-dp[i][m][0]);
}
double ans;//至少有一個人做對n道題目的概率
ans=1;
for(i=1;i<=t;i++)
{
double ks;
ks=0;
for(j=1;j<n;j++)
{
ks+=dp[i][m][j];
}
ans*=ks;
}
ans=bans-ans;
printf("%.3f\n",ans);
}
return 0;
}
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<queue>
#include<stack>
#include<map>
#include<string>
#include<stdlib.h>
#define INF_MAX 0x7fffffff
#define INF 999999
#define max3(a,b,c) (max(a,b)>c?max(a,b):c)
#define min3(a,b,c) (min(a,b)<c?min(a,b):c)
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
struct node
{
int u;
int v;
int w;
bool friend operator < (node a, node b){
return a.w < b.w;
}
}edge[1001];
int gcd(int n,int m){if(n<m) swap(n,m);return n%m==0?m:gcd(m,n%m);}
int lcm(int n,int m){if(n<m) swap(n,m);return n/gcd(n,m)*m;}
double p[1001][51];
double dp[1001][31][31];
int main()
{
int m,t,n,i,j,k;
while(scanf("%d%d%d",&m,&t,&n)&&(m||n||t))
{
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++)
{
dp[i][0][j]=0;
}
}
for(i=1;i<=t;i++)
{
for(j=1;j<=m;j++)
{
dp[i][j][0]=dp[i][j-1][0]*(1-p[i][j]);
for(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];
}
}
}
double bans;//所有人都做對1道及以上題目的概率
bans=1;
for(i=1;i<=t;i++)
{
bans*=(1-dp[i][m][0]);
}
double ans;//至少有一個人做對n道題目的概率
ans=1;
for(i=1;i<=t;i++)
{
double ks;
ks=0;
for(j=1;j<n;j++)
{
ks+=dp[i][m][j];
}
ans*=ks;
}
ans=bans-ans;
printf("%.3f\n",ans);
}
return 0;
}
