ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

POJ2151Check the difficulty of problems

2019-05-29 21:02:00  阅读:282  来源: 互联网

标签:概率 POJ2151Check int 32 problems long difficulty include define


 题目: 概率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

 

标签:概率,POJ2151Check,int,32,problems,long,difficulty,include,define
来源: https://www.cnblogs.com/Willems/p/10946176.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有