POJ2151Check the difficulty of problems

 题目: 概率dp||组合数学;

 

http://poj.org/problem?id=2151

 

一套题,有T个题,M个人应考,已知每个人做来某题的概率。问X的概率。X满足,每个考生至少做来一道题。至少有一人做的题不少于N道

思路:
p[i][j][k] 表示第i个考试前j个题会做k道的概率。再根据题意进行DP。
p[i][j][k] = p[i][j-1][k]*(1-a[i][j]) + p[i][j-1][k-1]*a[i][j];

然后用 每个人至少做出一道题的概率,减去 每个人至少做一道,至多做n-1道的概率; just answer。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#define LL long long
#define ULL unsigned long long
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define dep(i,j,k) for(int i=k;i>=j;i--)
#define INF 0x3f3f3f3f
#define mem(i,j) memset(i,j,sizeof(i))
#define make(i,j) make_pair(i,j)
#define pb push_back
using namespace std;
double p[1010][32][32],P1,P2,a[1010][32],tmp,res;
int main()
{
    int i,j,k,m,n,t;
    while(~scanf("%d%d%d",&m,&t,&n)){//t人 m题
        if(m==0&&t==0&&n==0) return 0;
        memset(p,0,sizeof(p));
        for(i=1;i<=t;i++)
         for(j=1;j<=m;j++)
           scanf("%lf",&a[i][j]);
        for(i=1;i<=t;i++)
         for(j=0;j<=m;j++)
          for(k=0;k<=j;k++){
               if(j==0&&k==0) p[i][j][k]=1;
               else if(j==k)  p[i][j][k]=p[i][j-1][k-1]*a[i][j];
               else if(k==0)  p[i][j][k]=p[i][j-1][k]*(1-a[i][j]);
               else p[i][j][k]=p[i][j-1][k-1]*a[i][j]+p[i][j-1][k]*(1-a[i][j]);
          }
        res=1;P1=1;
        for(i=1;i<=t;i++){
            P1*=(1-p[i][m][0]);
            tmp=0;
            for(j=1;j<n;j++) tmp+=p[i][m][j];
            res*=tmp;
        }
        res=P1-res;
        printf("%.3f\n",res); 
    }  return 0;
}
View Code

 

参考:https://www.cnblogs.com/hua-dong/p/8044969.html

 

转载于:https://www.cnblogs.com/Willems/p/10946176.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值