设为首页 加入收藏

TOP

poj 2096 Collecting Bugs(期望)
2015-07-20 17:44:55 来源: 作者: 【 】 浏览:3
Tags:poj 2096 Collecting Bugs 期望

?

?

程序的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个子集但不属于已有的j个分类,概率为i/n * (1 - j/s);

dp[i+1][j],新的bug不属于已有的i个子集但属于已有的j个分类,概率为(1-i/n)*j/s;

dp[i+1][j+1],新的bug不属于已有的i个子集也不属于已有的j个分类,概率为(1 - i/n) * (1 - j/s);

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

?

?

#include 
  
   
#include 
   
     #include
     #include 
     
       #include 
      
        #include 
       
         #include 
        
          #include 
         
           #include 
          
            #include 
           
             #include 
            
              #include 
             
               #include 
              
                //#define LL __int64 #define LL long long #define eps 1e-12 #define PI acos(-1.0) using namespace std; const int INF = 0x3f3f3f3f; const int maxn = 4010; double dp[1010][1010]; int main() { int n,s; while(~scanf(%d %d,&n,&s)) { dp[n][s] = 0; for(int i = n; i >= 0; i--) { for(int j = s; j >= 0; j--) { if(i == n && j == s) continue; dp[i][j] = (dp[i+1][j]*(n-i)*j + dp[i][j+1]*i*(s-j) + dp[i+1][j+1]*(n-i)*(s-j)+n*s)/(n*s - i*j)*1.0; } } printf(%.4f ,dp[0][0]); } return 0; } 
              
             
            
           
          
         
        
       
      
     
   
  


?

】【打印繁体】【投稿】【收藏】 【推荐】【举报】【评论】 【关闭】 【返回顶部
分享到: 
上一篇hdu 4791 dp预处理+二分 下一篇zoj 3103 Cliff Climbing(spfa )

评论

帐  号: 密码: (新用户注册)
验 证 码:
表  情:
内  容:

·MySQL 安装及连接-腾 (2025-12-25 06:20:28)
·MySQL的下载、安装、 (2025-12-25 06:20:26)
·MySQL 中文网:探索 (2025-12-25 06:20:23)
·Shell脚本:Linux Sh (2025-12-25 05:50:11)
·VMware虚拟机安装Lin (2025-12-25 05:50:08)