HDU4089/Uva1498 Activation 概率DP(好题)
题意:Tomato要在服务器上激活一个游戏,一开始服务器序列中有N个人,他排在第M位,每次服务器会对序列中第一位的玩家进行激活,有四种结果:
1.有p1的概率会激活失败,这时候序列的状态是不变的。
2.有p2的概率第一位的玩家会连接错误,这时候序列中第一位的玩家会成为最后一位,其他玩家相对位置不变。
3.有p3的概率第一位的玩家激活成功,这时候第一位的玩家会离开序列。
4.有p4的概率服务器崩溃,发生这件事之后所有玩家都不能激活了。
求Tomato遇到服务器崩溃并且在服务器崩溃时处于前K位的概率。
解法:这道题是一道非常好的概率DP题。必须一步一步来全部理解之后自己写一遍。
首先是设计状态dp[i]j[]代表队伍有i个人主角现在拍第j的概率。
写出状态转移方程(搭配题意写):
j==1 ,dp[i][j]=p1dp[i][j]+p2dp[i][i]
1<j<=k ,dp[i][j]=p1dp[i][j]+p2dp[i][j-1]+p3dp[i-1][j-1]+p4
j>k ,dp[i][j]=p1dp[i][j]+p2dp[i][j-1]+p3dp[i-1][j-1]
初步化简上诉式子:
其中 pp=1/(1-p1)
j==1 ,dp[i][j]=pp*p2dp[i][i]
1<j<=k ,dp[i][j]=pp*p2dp[i][j-1]+pp*p3dp[i-1][j-1]+pp*p4
j>k ,dp[i][j]=pp*p2dp[i][j-1]+pp*p3dp[i-1][j-1]
深度抽象上诉式子:
j==1 ,dp[i][j]=k[j] * dp[i][i] + b[j]
1<j<=k ,dp[i][j]=k[j] * dp[i][j-1] + b[j]
j>k ,dp[i][j]=k[j] * dp[i][j-1] + b[j]
(为什么是这样:我们仔细观察初步化简的式子:每一个方程右边除了dp[i][i]/dp[i][j-1]是与i同阶的,后面的都是i-1阶的,因为我们是一层一层做dp的,也就是说后面的可以看成是常数,所以这里的kj和bj具体是什么我们先根据初步化简的式子求出来放到数组里当成常数然后就不用管了)
我们仔细观察上诉深度化简的式子,发现其实第i阶就是有n个变量dp[i][j],然后有n个不同的方程,那么就一定可以用高斯消元解出方程!!!
但是再次仔细观察发现:欸!这n个方程结构很简单恰好形成了个环形!!那么我们就可以利用这点直接先解出一个未知数,然后递推出全部的未知数!
怎么解?
我们用笨方法模拟一下:假设只有3个未知数
dp1=k1dp3+c1
dp2=k2( k1dp3 + c1 ) +c2
dp3=k3( k2( k1dp3 + c1 ) +c2 ) + c3
我们尝试把dp3化简: dp3=k3k2k1dp3 + k3k2c1 + k3c2 + c3
那么规律不就出来了吗!
于是我们解出dp[i][i]再解出dp[i][1],然后就可以愉快地递推后面的值了。
#include<bits/stdc++.h>
using namespace std;
const int N=2e3+;
const double eps=1e-;
int n,m,K;
double p1,p2,p3,p4,pp,p14;
double k[N],b[N],dp[N][N]; void init() {
for (int i=;i<=n;i++) k[i]=b[i]=;
for (int i=;i<=n;i++) for (int j=;j<=n;j++) dp[i][j]=;
} int main()
{
while (cin>>n && n) {
cin>>m>>K;
cin>>p1>>p2>>p3>>p4;
init();
if (fabs(p4)<eps) {
printf("%.5lf\n",); continue;
}
pp=1.0/(1.0-p1); p14=1.0/(1.0-p1-p4); for (int i=;i<=n;i++) k[i]=pp*p2; dp[][]=p4/(-p1-p2);
for (int i=;i<=n;i++) {
b[]=pp*p4;
for (int j=;j<=K;j++) b[j]=pp*p3*dp[i-][j-]+pp*p4;
for (int j=K+;j<=i;j++) b[j]=pp*p3*dp[i-][j-]; double kk=1.0,bb=;
for (int j=i;j;j--) {
bb+=kk*b[j];
kk=k[j]*kk;
} dp[i][i]=bb/(1.0-kk);
dp[i][]=k[]*dp[i][i]+b[];
for (int j=;j<i;j++) dp[i][j]=k[j]*dp[i][j-]+b[j];
}
printf("%.5lf\n",dp[n][m]);
}
return ;
}
HDU4089/Uva1498 Activation 概率DP(好题)的更多相关文章
- 13年山东省赛 The number of steps(概率dp水题)
转载请注明出处: http://www.cnblogs.com/fraud/ ——by fraud The number of steps Time Limit: 1 Sec Me ...
- Codeforces 148D 一袋老鼠 Bag of mice | 概率DP 水题
除非特别忙,我接下来会尽可能翻译我做的每道CF题的题面! Codeforces 148D 一袋老鼠 Bag of mice | 概率DP 水题 题面 胡小兔和司公子都认为对方是垃圾. 为了决出谁才是垃 ...
- [HDU 4089]Activation[概率DP]
题意: 有n个人排队等着在官网上激活游戏.Tomato排在第m个. 对于队列中的第一个人.有以下情况: 1.激活失败,留在队列中等待下一次激活(概率为p1) 2.失去连接,出队列,然后排在队列的最后( ...
- [HDU4089]Activation(概率DP)
HDU4089 题意:有n个人排队等着在官网上激活游戏.Tomato排在第m个. 对于队列中的第一个人.有一下情况: 1.激活失败,留在队列中等待下一次激活(概率为p1) 2.失去连接,出队列,然后排 ...
- HDU 4089 Activation 概率DP 难度:3
http://acm.hdu.edu.cn/showproblem.php?pid=4089 这道题中一共有两个循环: 1.事件1 如果一直落在Activation failed事件上,那么就会重新继 ...
- poj 3071 Football (概率DP水题)
G - Football Time Limit:1000MS Memory Limit:65536KB 64bit IO Format:%I64d & %I64u Submit ...
- hdu4405Aeroplane chess 概率dp水题
//从0到n有n+1个格子 //对于格子i,掷一次骰子的数为x.那么能够从位置i到位置i+x //格子之间有连线,假设格子a和b有连线,那么从a到b不用掷骰子 //求从0到n的骰子掷的次数的期望 // ...
- 【BZOJ1426】收集邮票 概率DP 论文题 推公式题
链接: #include <stdio.h> int main() { puts("转载请注明出处[辗转山河弋流歌 by 空灰冰魂]谢谢"); puts("网 ...
- [转]概率DP总结 by kuangbin
概率类题目一直比较弱,准备把kuangbin大师傅总结的这篇题刷一下! 我把下面的代码换成了自己的代码! 原文地址:http://www.cnblogs.com/kuangbin/archive/20 ...
随机推荐
- NIO编程之多客户端聊天系统
1. 服务端 import java.io.IOException; import java.net.InetSocketAddress; import java.nio.ByteBuffer; im ...
- CSS基础知识总结之css样式引用的三种方式
在html中增加css样式有三种: 1.在标签中增加style属性: <!DOCTYPE html> <html lang="en" xmlns="ht ...
- PHP获取时间排除周六、周日的两个方法
//方法一: <?php $now = time(); //指定日期用法 $now = strtotime('2014-01-08') ; $day = 3600*24; $total = 12 ...
- CF1182 D Complete Mirror——思路
题目:http://codeforces.com/contest/1182/problem/D 很好的思路是从度数为1的点和直径来入手. 找一条直径.看看直径的两个端点是否合法. 如果都不合法,那么根 ...
- LRU management
LRU management 字典树用来查找值,实现map<string,int>操作 tips:tot必须从一开始QAQ #include<bits/stdc++.h> us ...
- php二维数组对某一列进行数据组装
demo示例: $data = array( '0' => array('id' => 10, 'username' => '陈一'), '1' => array('id' = ...
- 通过 homebrew下载速度过慢的解决方案
请移步 https://www.cnblogs.com/jingxiaoniu/p/11123377.html 进行操作
- Linux(RHEL7)下安装vsftp服务
1.安装vsftp(没有配置yum源的先配置yum源) yum install -y vsftpd 2.启动ftp服务 systemctl start vsftpd.service 3.打开防火墙 f ...
- Repository模式--采用EF Fluent API使用EntityTypeConfiguration分文件配置Model映射关系
EF中类EntityTypeConfiguration是一个很有用的类,在nopCommerence中就使用这个类来分文件分文件配置Model映射关系.今天我就来谈谈Repository模式在Enti ...
- MySQL 导入导出 CSV 文件
导入 导出 清空表中的所有数据 注意事项 常见问题 ERROR 1290 (HY000): The MySQL server is running with the --secure-file-pri ...