//一条队列对于第一个人的情况

//概率p1:队列保持不变

//p2:第一个人到队尾

//p3:第一个人出队

//p4:系统崩溃

//队列中有N个人,Tomato 在第M位置,求系统崩溃。Tomato 前面的人数为小于k个人的情况的概率

//dp[i][j] 表示队列中有i个人,Tomato 在第j个位置时出现目标状态的概率

//能够非常easy得到递推公式

//dp[i][1] = p2/(1-p1)*dp[i][i] + p4/(1-p1);

//dp[i][j] = p2/(1-p1)*dp[i][j-1] + p3/(1 -p1)*dp[i-1][j-1] + p4/(1-p1);

//dp[i][j] = p2/(1-p1)*dp[i][j-1] + p3/(1-p1)*dp[i-1][j-1] ;

//仅仅有三种情况发生在Tomato 身上,成功出队,系统崩溃时大于k系统崩溃时小于等于k

//所以初始条件dp[1][1] = p4/(p3+p4)

//然后每次先通过迭代求出dp[i][i] ,得到答案

//须要注意的是二维数组会超空间,所以要用滚动数组

#include<cstdio>

#include<cstring>

#include<iostream>

#include<cmath>

using namespace std ;

const int maxn = 2010 ;

double dp[2][maxn] ;

int n, m , k ;

int main()

{

   // freopen("in.txt" ,"r" ,stdin) ;

    double p1 , p2 , p3 , p4 ,tmp , temp;

    int i , j ;

    while(~scanf("%d%d%d%lf%lf%lf%lf" , &n ,&m , &k , &p1 , &p2 , &p3 , &p4))

    {

        memset(dp , 0  , sizeof(dp)) ;

        if(p4 == 0)

        {

            puts("0.00000") ;

            continue ;

        }

        dp[1][1] = p4/(p3+p4);

        for( i = 2;i <= n;i++)

        {

            tmp = 0 ;

            temp = 1 ;

            for( j = i;j > 0 ;j--)

            {

                if(j == 1)tmp+= p4/(1-p1)*temp ;

                else if(j > k)tmp+=p3/(1-p1)*dp[(i-1)%2][j-1]*temp ;

                else tmp+=(p3/(1 -p1)*dp[(i-1)%2][j-1] + p4/(1-p1))*temp ;

                temp*=p2/(1-p1);

            }

            dp[i%2][i] = tmp/(1- temp) ;

            for( j = 1;j < i ;j++)

            {

                if(j == 1)dp[i%2][1] = p2/(1-p1)*dp[i%2][i] + p4/(1-p1);

                else if(j <= k)dp[i%2][j] = p2/(1-p1)*dp[i%2][j-1] + p3/(1 -p1)*dp[(i-1)%2][j-1] + p4/(1-p1);

                else dp[i%2][j] = p2/(1-p1)*dp[i%2][j-1] + p3/(1-p1)*dp[(i-1)%2][j-1] ;

            }

        }

        printf("%.5lf\n" , dp[n%2][m]) ;

    }

}

hdu4089Activation 概率dp的更多相关文章

  1. Codeforces 28C [概率DP]

    /* 大连热身D题 题意: 有n个人,m个浴室每个浴室有ai个喷头,每个人等概率得选择一个浴室. 每个浴室的人都在喷头前边排队,而且每个浴室内保证大家都尽可能均匀得在喷头后边排队. 求所有浴室中最长队 ...

  2. HDU 4405 Aeroplane chess (概率DP)

    题意:你从0开始,要跳到 n 这个位置,如果当前位置是一个飞行点,那么可以跳过去,要不然就只能掷骰子,问你要掷的次数数学期望,到达或者超过n. 析:概率DP,dp[i] 表示从 i  这个位置到达 n ...

  3. POJ 2096 Collecting Bugs (概率DP)

    题意:给定 n 类bug,和 s 个子系统,每天可以找出一个bug,求找出 n 类型的bug,并且 s 个都至少有一个的期望是多少. 析:应该是一个很简单的概率DP,dp[i][j] 表示已经从 j ...

  4. POJ 2151 Check the difficulty of problems (概率DP)

    题意:ACM比赛中,共M道题,T个队,pij表示第i队解出第j题的概率 ,求每队至少解出一题且冠军队至少解出N道题的概率. 析:概率DP,dp[i][j][k] 表示第 i 个队伍,前 j 个题,解出 ...

  5. 概率DP light oj 1030

    t组数据 n块黄金 到这里就捡起来 出发点1 到n结束  点+位置>n 重掷一次 dp[i] 代表到这里的概率 dp[i]=(dp[i-1]+dp[i-2]... )/6  如果满6个的话 否则 ...

  6. hdu 4050 2011北京赛区网络赛K 概率dp ***

    题目:给出1-n连续的方格,从0开始,每一个格子有4个状态,左右脚交替,向右跳,而且每一步的步长必须在给定的区间之内.当跳出n个格子或者没有格子可以跳的时候就结束了,求出游戏的期望步数 0:表示不能到 ...

  7. [转]概率DP总结 by kuangbin

    概率类题目一直比较弱,准备把kuangbin大师傅总结的这篇题刷一下! 我把下面的代码换成了自己的代码! 原文地址:http://www.cnblogs.com/kuangbin/archive/20 ...

  8. SGU 422 Fast Typing(概率DP)

    题目大意 某人在打字机上打一个字符串,给出了他打每个字符出错的概率 q[i]. 打一个字符需要单位1的时间,删除一个字符也需要单位1的时间.在任意时刻,他可以花 t 的时间检查整个打出来的字符串,并且 ...

  9. HDU 4050 wolf5x(动态规划-概率DP)

    wolf5x Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Sub ...

随机推荐

  1. [POJ3233]Matrix Power Series 分治+矩阵

    本文为博主原创文章,欢迎转载,请注明出处 www.cnblogs.com/yangyaojia [POJ3233]Matrix Power Series 分治+矩阵 题目大意 A为n×n(n<= ...

  2. HDU——T 1054 Strategic Game

    http://acm.hdu.edu.cn/showproblem.php?pid=1054 Time Limit: 20000/10000 MS (Java/Others)    Memory Li ...

  3. cogs 2056. 无平方因子数

    2056. 无平方因子数 ★☆   输入文件:non.in   输出文件:non.out   简单对比时间限制:1 s   内存限制:256 MB [题目描述] 给出正整数n,m,区间[n,m]内的无 ...

  4. 数据库联表统计查询 Group by & INNER JOIN

    原数据表 视频信息表  tab_video_info 播放记录表  tab_play_record 需求 统计播放量(已经开始播放)最多的前20个视频: SELECT a.video_id, SUM( ...

  5. HMM(隐马尔可夫),这篇讲的不错

    http://www.cnblogs.com/skyme/p/4651331.html

  6. HDU 2879

    利用x<n的信息,可以证得当n为素数时,he[n]=2;同时,若n 为素数,则有HE[N^K]=2;因为若等式成立则有n|x(x-1).抓住这个证即可. 至于符合积性函数,想了很久也没想出来,看 ...

  7. bzoj 1600 &amp; Usaco 月赛 2008 建造栅栏 题解

    [原题] 1600: [Usaco2008 Oct]建造栅栏 Time Limit: 5 Sec  Memory Limit: 64 MB Submit: 785  Solved: 443 [Subm ...

  8. ios+openflow 问题

    环境:xcode5.1+ios7.1 需求:A试图 的scroll加入 B视图:[A addSubview:B.view] 问题: 1.B视图载入到A视图上了,但Openflow的图片未载入.后经调试 ...

  9. 使用Xcode-debug模式和release模式

    在开发过程中,我们常常须要用到NSLog输出一些信息,甚至有的开发过程.必须在控制台查看输出,有经验的程序猿通过控制台输出就能知道整个数据交互的一个流程. 可是一个公布的程序,里面带有太多的NSLog ...

  10. Linux下安装ipython与jupyter

    IPython从Python发展而来,更倾向于科学计算.互联网数据分析更喜欢用. 首先切换root用户: su - root pip3自动安装ipython [root@hear ~]# ipytho ...