POJ 2096 找bug 期望dp

题目大意:

一个人受雇于某公司要找出某个软件的bugs和subcomponents,这个软件一共有n个bugs和s个subcomponents,每次他都能同时随机发现1个bug和1个subcomponent,问他找到所有的bugs和subcomponents的期望次数。

这道题目要用期望dp来进行统计

假设已经找到i个bug和j个subcomponents,这个状态记为dp[i][j],那么下次查找会出现4种状态:dp[i][j],dp[i+1][j],dp[i][j+1],dp[i+1][j+1]

这四种状态的概率为 i*j*1.0/n/s , (n-i)*j*1.0/n/s , i*(s-j)*1.0/n/s , (n-i)*(s-j)*1.0/n/s

那么每到下一种 状态步数要加1,那么我们领dp[i][j]为到达全部找完bug和subcomponents的数学期望,亦可理解为,在这个状态还需多少步能完成任务

这样我们就可以得到公式:

dp[i][j] = 1 + i*j*1.0/n/s * dp[i][j] +  (n-i)*j*1.0/n/s * dp[i+1][j] + i*(s-j)*1.0/n/s * dp[i][j+1] + (n-i)*(s-j)*1.0/n/s * dp[i+1][j+1]

那么我们就可以通过后面的状态不断递推来找到前面的状态的值就行了

我们可以很容易知道dp[n][s] = 0;

为了方便计算,不用排除特殊情况,那么我们就将dp[n+1][1]~dp[n+1][s]   dp[1][1+s]~dp[n][s+1]设为0作为初始化

然后从后往前不断完善这个矩阵的dp[][]值

最后输出dp[0][0]就是我们要求的值

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <cstring>
 4 using namespace std;
 5 #define N 1005
 6 double dp[N][N];
 7 int main()
 8 {
 9     int n,m;
10     while(~scanf("%d%d",&n,&m)){
11         for(int i=1;i<=n+1;i++)
12             dp[i][m+1]=0;
13
14         for(int i=1;i<=m+1;i++)
15             dp[n+1][i]=0;
16
17         dp[n][m]=0;
18
19         for(int i=n;i>=0;i--){
20             for(int j=m;j>=0;j--){
21                 if(i==n&&j==m)
22                     continue;
23                 double tmp = 1 + i*(m-j)*1.0/n/m * dp[i][j+1] + (n-i)*j*1.0/m/n * dp[i+1][j] + (n-i)*(m-j)*1.0/m/n * dp[i+1][j+1];
24                 dp[i][j] = tmp / (1-i*j*1.0/m/n);
25             }
26         }
27
28         printf("%.4f\n",dp[0][0]);
29     }
30     return 0;
31 }
时间: 2024-10-10 14:05:28

POJ 2096 找bug 期望dp的相关文章

poj 2096 Collecting Bugs——期望DP

题目:http://poj.org/problem?id=2096 f[ i ][ j ] 表示收集了 i 个 n 的那个. j 个 s 的那个的期望步数. #include<cstdio> #include<cstring> #include<algorithm> #define db double using namespace std; const int N=1005; db n,s,f[N][N]; int main() { scanf("%lf%l

POJ - 2096 Collecting Bugs 期望dp

期望是一个 DAG 模型,逆着递推即可~ #include <cstdio> #define N 1005 #define setIO(s) freopen(s".in","r",stdin) using namespace std; double f[N][N]; int main() { int n,s,i,j; scanf("%d%d",&n,&s); for(int i=n;i>=0;--i) { for

poj 2096 Collecting Bugs (概率dp)

题意:有s个系统,n种bug,一个程序员每天可以发现一个bug,求发现存在s个系统,n种bug的天数的期望 思路:定义dp[i][j]是已经发现i种bug,j个系统的期望 dp[i+1][j+1] 表示在一个新的系统中发现新bug 它的概率为 (n-i)/n*(s-j)/s dp[i+1][j]     表示在一个已经发现过的系统中发现了一种新的bug 概率为 (n-i)/n*j/s dp[i][j+1]     表示在一个新系统中发现一种已经发现过的bug 概率为 i/n*(s-j)/s dp

POJ 2096 Collecting Bugs(dp 期望)

题目链接:http://poj.org/problem?id=2096 Description Ivan is fond of collecting. Unlike other people who collect post stamps, coins or other material stuff, he collects software bugs. When Ivan gets a new program, he classifies all possible bugs into n ca

poj 2096 Collecting Bugs (概率dp 天数期望)

题目链接 题意: 一个人受雇于某公司要找出某个软件的bugs和subcomponents,这个软件一共有n个bugs和s个subcomponents,每次他都能同时随机发现1个bug和1个subcomponent,问他找到所有的bugs和subcomponents的期望次数. 分析: 期望倒着推,概率正着推. dp[i][j]表示已经找到i种bug,并存在于j个子系统中,要达到目标状态的天数的期望.显然,dp[n][s]=0,因为已经达到目标了.而dp[0][0]就是我们要求的答案.dp[i][

Poj 2096 Collecting Bugs (概率DP求期望)

C - Collecting Bugs Time Limit:10000MS     Memory Limit:64000KB     64bit IO Format:%I64d & %I64u Submit Status Description Ivan is fond of collecting. Unlike other people who collect post stamps, coins or other material stuff, he collects software b

poj 2096 Collecting Bugs(期望)

http://poj.org/problem?id=2096 程序的bug有n个子集,s个种类.每一个bug属于每个子集的概率为1/n,每一个bug属于每个种类的概率为1/s,问每个子集且每个种类都有bug的期望. 求期望,设dp[i][j]表示已有bug属于i个子集,j个种类的期望,现已知终态为dp[n][s] = 0,dp[i][j]可由逆推而得: dp[i][j],新的bug属于已有的i个子集j个分类,概率为i/n * j/s; dp[i][j+1],新的bug属于已有的i个子集但不属于已

POJ 2096 Collecting Bugs(概率DP)

解题思路: 题目比较难理解,大致题意就是一共有N种bugs,分别属于S个子系统,求找到N种BUG并且每种子系统的bug都被找到所需要的天数的数学期望. dp[i][j]表示找到 i 种bug 属于 j 个子系统到 目标状态所需要的数学期望.dp[i][j]可以由四种状态转移而来. (i * j) / (n * s) * d[i][j]; (n - i) * j / (n * s) * dp[i+1][j]; i * (s - j)/(n * s) * dp[i][j + 1]; (n - i)

POJ 2096 Collecting Bugs 概率dp(水

题目链接:点击打开链接 题意: 点击打开链接 对于这里的dp做法是: 写一个状态x,然后把从x转移出去的方程写出来,即 x = y1+y2+··· 其中所有的yi都是已知的. 这样我们就会得到一个方程是从未知到已知. 但是dp是由已知到未知.所以我们再呵呵回来.. #include <cstdio> #include <iostream> #include <cstring> #include <queue> #include <algorithm&g