hdu4089Activation 概率dp
//一条队列对于第一个人的情况
//概率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的更多相关文章
- Codeforces 28C [概率DP]
/* 大连热身D题 题意: 有n个人,m个浴室每个浴室有ai个喷头,每个人等概率得选择一个浴室. 每个浴室的人都在喷头前边排队,而且每个浴室内保证大家都尽可能均匀得在喷头后边排队. 求所有浴室中最长队 ...
- HDU 4405 Aeroplane chess (概率DP)
题意:你从0开始,要跳到 n 这个位置,如果当前位置是一个飞行点,那么可以跳过去,要不然就只能掷骰子,问你要掷的次数数学期望,到达或者超过n. 析:概率DP,dp[i] 表示从 i 这个位置到达 n ...
- POJ 2096 Collecting Bugs (概率DP)
题意:给定 n 类bug,和 s 个子系统,每天可以找出一个bug,求找出 n 类型的bug,并且 s 个都至少有一个的期望是多少. 析:应该是一个很简单的概率DP,dp[i][j] 表示已经从 j ...
- POJ 2151 Check the difficulty of problems (概率DP)
题意:ACM比赛中,共M道题,T个队,pij表示第i队解出第j题的概率 ,求每队至少解出一题且冠军队至少解出N道题的概率. 析:概率DP,dp[i][j][k] 表示第 i 个队伍,前 j 个题,解出 ...
- 概率DP light oj 1030
t组数据 n块黄金 到这里就捡起来 出发点1 到n结束 点+位置>n 重掷一次 dp[i] 代表到这里的概率 dp[i]=(dp[i-1]+dp[i-2]... )/6 如果满6个的话 否则 ...
- hdu 4050 2011北京赛区网络赛K 概率dp ***
题目:给出1-n连续的方格,从0开始,每一个格子有4个状态,左右脚交替,向右跳,而且每一步的步长必须在给定的区间之内.当跳出n个格子或者没有格子可以跳的时候就结束了,求出游戏的期望步数 0:表示不能到 ...
- [转]概率DP总结 by kuangbin
概率类题目一直比较弱,准备把kuangbin大师傅总结的这篇题刷一下! 我把下面的代码换成了自己的代码! 原文地址:http://www.cnblogs.com/kuangbin/archive/20 ...
- SGU 422 Fast Typing(概率DP)
题目大意 某人在打字机上打一个字符串,给出了他打每个字符出错的概率 q[i]. 打一个字符需要单位1的时间,删除一个字符也需要单位1的时间.在任意时刻,他可以花 t 的时间检查整个打出来的字符串,并且 ...
- HDU 4050 wolf5x(动态规划-概率DP)
wolf5x Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Sub ...
随机推荐
- [POJ3233]Matrix Power Series 分治+矩阵
本文为博主原创文章,欢迎转载,请注明出处 www.cnblogs.com/yangyaojia [POJ3233]Matrix Power Series 分治+矩阵 题目大意 A为n×n(n<= ...
- HDU——T 1054 Strategic Game
http://acm.hdu.edu.cn/showproblem.php?pid=1054 Time Limit: 20000/10000 MS (Java/Others) Memory Li ...
- cogs 2056. 无平方因子数
2056. 无平方因子数 ★☆ 输入文件:non.in 输出文件:non.out 简单对比时间限制:1 s 内存限制:256 MB [题目描述] 给出正整数n,m,区间[n,m]内的无 ...
- 数据库联表统计查询 Group by & INNER JOIN
原数据表 视频信息表 tab_video_info 播放记录表 tab_play_record 需求 统计播放量(已经开始播放)最多的前20个视频: SELECT a.video_id, SUM( ...
- HMM(隐马尔可夫),这篇讲的不错
http://www.cnblogs.com/skyme/p/4651331.html
- HDU 2879
利用x<n的信息,可以证得当n为素数时,he[n]=2;同时,若n 为素数,则有HE[N^K]=2;因为若等式成立则有n|x(x-1).抓住这个证即可. 至于符合积性函数,想了很久也没想出来,看 ...
- bzoj 1600 & Usaco 月赛 2008 建造栅栏 题解
[原题] 1600: [Usaco2008 Oct]建造栅栏 Time Limit: 5 Sec Memory Limit: 64 MB Submit: 785 Solved: 443 [Subm ...
- ios+openflow 问题
环境:xcode5.1+ios7.1 需求:A试图 的scroll加入 B视图:[A addSubview:B.view] 问题: 1.B视图载入到A视图上了,但Openflow的图片未载入.后经调试 ...
- 使用Xcode-debug模式和release模式
在开发过程中,我们常常须要用到NSLog输出一些信息,甚至有的开发过程.必须在控制台查看输出,有经验的程序猿通过控制台输出就能知道整个数据交互的一个流程. 可是一个公布的程序,里面带有太多的NSLog ...
- Linux下安装ipython与jupyter
IPython从Python发展而来,更倾向于科学计算.互联网数据分析更喜欢用. 首先切换root用户: su - root pip3自动安装ipython [root@hear ~]# ipytho ...