POJ 2096(概率DP)

原创 2016年02月24日 16:48:31

题目链接:POJ 2096

解题思路:
总的来说状态方程还是蛮好想的,每次找bug的时候,有可能会出现四种情况,一是在已经找到bug的系统中找到了旧bug,二是在新的系统中找到旧bug,三是在已经找到bug的系统中找到了新bug,二是在新的系统中找到新bug。于是,状态dp[i][j]表示在 j 个系统中找到 i 个bug需要的天数期望,初始状态dp[0][0]=0,状态转移方程如下:

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

通过移向可以得到新的一条方程:

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

然而,这样做是会有问题的,问题在于最后一次递推,即 i = n 且 j = s 的时候,会出现除以0的情况,这样显然得不到正确结果,而且这也间接说明,这个状态方程其实是错的。因此我们需要对状态进行修改,即状态dp[i][j]表示在 j 个系统中找到 i 个bug后,距离目标状态的天数期望,初始状态dp[n][s]=0,状态转移方程如下:

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

最终转移方程如下:

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

AC代码:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

int n,s;
double dp[1005][1005];

int main()
{
    while(~scanf("%d %d",&n,&s))
    {
        double p1,p2,p3;

        for(int i=n;i>=0;i--)
        {
            for(int j=s;j>=0;j--)
            {
                if(i==n&&j==s)
                    continue;

                p1=double((n-i)*j)/n/s;             
                p2=double((i)*(s-j))/n/s;               
                p3=double((n-i)*(s-j))/n/s;

                dp[i][j]=(p1*dp[i+1][j] + p2*dp[i][j+1] + p3*dp[i+1][j+1] + 1) / (1.0-double(i*j)/n/s);
            }
        }
        printf("%.4f\n",dp[0][0]);
    }

    return 0;
}

总结:
1、概率DP,套基本的期望公式即可;
2、逆向递推,
3、此题精度有毒

相关文章推荐

POJ 2096 Collecting Bugs(概率dp)

Language: Default Collecting Bugs Time Limit: 10000MS   Memory Limit: 64000K Tot...

POJ 2096 Collecting Bugs 概率dp 求期望 入门

题意: 一个人很喜欢去发现bug。每一天能

[poj 2096]Collecting Bugs 概率DP

Collecting Bugs Time Limit: 10000MS Memory Limit: 64000K Total Submissions: 4767 Accepted:...

poj 2096 Collecting Bugs(概率dp)

题意:Ivan每天能找到一个bug,现在有一个工程,有n种bug,s个子系统。问在这个工程中找出n种bug并且每个子系统至少有1个bug的天数的期望是多少。 思路:用dp[i][j]表示当前找出了i...
  • qian99
  • qian99
  • 2014年03月17日 19:45
  • 435

POJ 2096 Collecting Bugs (概率dp)

POJ 2096 Collecting Bugs (概率dp)

POJ 2096 Collecting Bugs (概率DP)

Collecting Bugs Time Limit: 10000MS   Memory Limit: 64000K Total Submissions: 4038   Accepted:...

POJ 2096 Collecting Bugs 详细题解(概率DP求期望)

Collecting Bugs Time Limit: 10000MS   Memory Limit: 64000K Total Submissions: 5302 ...

poj2096(概率DP)

Collecting Bugs Time Limit: 10000MS   Memory Limit: 64000K Total Submissions: 1656   Accepted:...

poj2096——概率DP

Collecting Bugs Time Limit: 10000MS   Memory Limit: 64000K Total Submissions: 5371 ...

POJ 2096 Collecting Bugs 马尔可夫链 概率DP

AC自动机的专题刷完了...开始刷概率了... 题目大意: 就是现在要从一个软件中找出bug, 这个软件有s个子模块,一共有n种bug,现在每天你可以找到一个bug, 这个bug属于任何一...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:POJ 2096(概率DP)
举报原因:
原因补充:

(最多只允许输入30个字)