Time Limit: 20 Sec  Memory Limit: 512 MBSec  Special Judge
Submit: 1952  Solved: 1159
[Submit][Status][Discuss]

Description

小 K 不慎被 LL 邪教洗脑了,洗脑程度深到他甚至想要从亚瑟王邪教中脱坑。

他决定,在脱坑之前,最后再来打一盘亚瑟王。既然是最后一战,就一定要打得漂
亮。众所周知,亚瑟王是一个看脸的游戏,技能的发动都是看概率的。作为一个非
洲人,同时作为一个前 OIer,小 K 自然是希望最大化造成伤害的期望值。但他已
经多年没写过代码,连 Spaly都敲不对了,因此,希望你能帮帮小 K,让他感受一
下当欧洲人是怎样的体验。 
本题中我们将考虑游戏的一个简化版模型。 
玩家有一套卡牌,共 n张。游戏时,玩家将 n 张卡牌排列成某种顺序,排列后
将卡牌按从前往后依次编号为 1 ~  n。本题中,顺序已经确定,即为输入的顺序。
每张卡牌都有一个技能。第 i 张卡牌的技能发动概率为 pi,如果成功发动,则会对
敌方造成di点伤害。也只有通过发动技能,卡牌才能对敌方造成伤害。基于现实因
素以及小K非洲血统的考虑,pi不会为 0,也不会为 1,即 0 < pi < 1。 
一局游戏一共有 r 轮。在每一轮中,系统将从第一张卡牌开始,按照顺序依次
考虑每张卡牌。在一轮中,对于依次考虑的每一张卡牌: 
1如果这张卡牌在这一局游戏中已经发动过技能,则 
1.1 如果这张卡牌不是最后一张,则跳过之(考虑下一张卡牌); 
否则(是最后一张),结束这一轮游戏。 
2否则(这张卡牌在这一局游戏中没有发动过技能),设这张卡牌为第 i 张 
2.1将其以 pi的概率发动技能。 
2.2如果技能发动,则对敌方造成 di点伤害,并结束这一轮。 
2.3如果这张卡牌已经是最后一张(即 i 等于n),则结束这一轮;否则,
考虑下一张卡牌。 
请帮助小 K 求出这一套卡牌在一局游戏中能造成的伤害的期望值。 

Input

输入文件的第一行包含一个整数 T,代表测试数据组数。

接下来一共 T 组数据。 
每组数据的第一行包含两个用空格分开的整数 n和r,分别代表卡牌的张数和
游戏的轮数。 
接下来 n行,每行包含一个实数和一个整数,由空格隔开,描述一张卡牌。第
i 行的两个数为 pi和 di,分别代表第 i 张卡牌技能发动的概率(实数)和技能发动
造成的伤害(整数)。保证 pi最多包含 4位小数,且为一个合法的概率。 

Output

对于每组数据,输出一行,包含一个实数,为这套卡牌在这一局游戏中造成的

伤害的期望值。对于每一行输出,只有当你的输出和标准答案的相对误差不超过
10^-8时——即|a-o|/a<=10-8时(其中a是标准答案,o是输出),你的输出才会被判为正确。
建议输出10 位小数。 

Sample Input

1
3 2
0.5000 2
0.3000 3
0.9000 1

Sample Output

3.2660250000

HINT

一共有 13 种可能的情况:

1.  第一轮中,第 1张卡牌发动技能;第二轮中,第 2张卡牌发动技能; 
概率为 0.15,伤害为5。 
2.  第一轮中,第 1张卡牌发动技能;第二轮中,第 3张卡牌发动技能; 
概率为 0.315,伤害为3。 
3.  第一轮中,第 1张卡牌发动技能;第二轮不发动技能; 
概率为 0.035,伤害为2。 
4.  第一轮中,第 2张卡牌发动技能;第二轮中,第 1张卡牌发动技能; 
概率为 0.075,伤害为5。 
5.  第一轮中,第 2张卡牌发动技能;第二轮中,第 3张卡牌发动技能; 
概率为 0.0675,伤害为4。 
6.  第一轮中,第 2张卡牌发动技能;第二轮不发动技能; 
概率为 0.0075,伤害为3。 
7.  第一轮中,第 3张卡牌发动技能;第二轮中,第 1张卡牌发动技能; 
概率为 0.1575,伤害为3。 
8.  第一轮中,第 3张卡牌发动技能;第二轮中,第 2张卡牌发动技能; 
概率为 0.04725,伤害为4。 
9.  第一轮中,第 3张卡牌发动技能;第二轮不发动技能; 
概率为 0.11025,伤害为1。 
10.  第一轮不发动技能;第二轮中,第 1张卡牌发动技能; 
概率为 0.0175,伤害为2。 
11.  第一轮不发动技能;第二轮中,第 2张卡牌发动技能; 
概率为 0.00525,伤害为3。 
12.  第一轮不发动技能;第二轮中,第 3张卡牌发动技能; 
概率为 0.011025,伤害为1。 
13.  第一轮不发动技能;第二轮亦不发动技能; 
概率为 0.001225,伤害为0。 
造成伤害的期望值为概率与对应伤害乘积之和,为 3.266025。 
对于所有测试数据, 1 <= T <= 444, 1 <= n <= 220, 0 <= r <= 132, 0 < pi < 1, 0 <= di <= 1000。  
除非备注中有特殊说明,数据中 pi与di均为随机生成。 
请注意可能存在的实数精度问题,并采取适当措施。 

Source

感觉期望的题都好绕啊qwq。。

首先根据期望的线性,每个位置对BOSS造成的伤害都是独立的。

根据期望的定义而且此题中每个位置只允许攻击一次,我们只需算出每个位置在$r$轮中出现的概率即可

这个概率为$1 - (1 - p[i])^r$,即总的概率减去一次都没出现的概率

但是这题有一个特殊限制“如果前面的打出了,那么本轮游戏结束”,也就是说如果$i - 1$位置打出了,那么$i$位置本次肯定是不能出现的

这样的话,如果$i$前面一共有$j$个位置打出了,也就是$i$有$j$次强制没被打出,那么$i$出现的概率为$1 - (1 - p[i])^{r - j}$

考虑这玩意儿如何计算,设$f[i][j]$表示前$i$个中强制打出了$j$个的概率

这显然是可以递推完成的,讨论一下这个位置是否打出就行了

那么$g[i] = \sum_i f[i][j] * (1 - (1 - p[i])^{r - j})$

最终答案为$ans = \sum g[i] * d[i]$

// luogu-judger-enable-o2
// luogu-judger-enable-o2
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#define LL long long
using namespace std;
const int MAXN = , INF = 1e9 + ;
inline int read() {
char c = getchar(); int x = , f = ;
while(c < '' || c > '') {if(c == '-') f = -; c = getchar();}
while(c >= '' && c <= '') x = x * + c - '', c = getchar();
return x * f;
}
int N, R;
double f[][],/*前i轮中,恰好有j张牌发动了攻击的概率*/ P[], D[], PowP[][], G[];
void GetPow() {
for(int i = ; i <= N; i++) {
PowP[i][] = ;
for(int j = ; j <= R; j++)
PowP[i][j] = PowP[i][j - ] * ( - P[i]);
}
}
main() {
int QWQ;
scanf("%d", &QWQ);
while(QWQ--) {
memset(f, , sizeof(f));
scanf("%d %d", &N, &R);
for(int i = ; i <= N; i++)
scanf("%lf %lf", &P[i], &D[i]);
GetPow();
f[][] = G[] = - PowP[][R]; f[][] = PowP[][R];
for(int i = ; i <= N; i++) {
G[i] = ;
for(int j = ; j <= min(i, R); j++) {
if(j) f[i][j] += f[i - ][j - ] * ( - PowP[i][R - j + ]);
//发动攻击,那么R-j+1轮中不会出现
if(i != j) f[i][j] += f[i - ][j] * PowP[i][R - j];
//第i张牌始终没有出现过
G[i] += f[i - ][j] * ( - PowP[i][R - j]);
}
}
double ans = ;
for(int i = ; i <= N; i++)
ans += G[i] * D[i];
printf("%.10lf\n", ans);
}
return ;
}

BZOJ4008: [HNOI2015]亚瑟王(期望dp)的更多相关文章

  1. P3239 [HNOI2015]亚瑟王 期望dp

    这个题一看就是期望dp,但是我有个问题,一个事件的期望等于他所有事件可能行乘权值的和吗...为什么我有天考试的时候就不对呢...求大佬解释一下. 至于这道题,f[i][j]代表前i个有j个发动技能,这 ...

  2. P3239 [HNOI2015]亚瑟王 期望 dp

    LINK:亚瑟王 Saber!Excalibur! 比较难的期望dp. 可以发现如果暴力枚举所有的局面复杂度很高 . 转换的思路则是 期望的线性性. 求出每张牌的期望累加即可. 考虑每张牌的期望=这张 ...

  3. 洛谷 P3239 [HNOI2015]亚瑟王(期望dp)

    题面 luogu 题解 一道复杂的期望\(dp\) 思路来源:__stdcall 容易想到,只要把每张牌打出的概率算出来就可以求出\(ans\) 设\(fp[i]\)表示把第\(i\)张牌打出来的概率 ...

  4. BZOJ4008. [HNOI2015]亚瑟王 期望概率dp

    看到这道题想什么? 一个好转移的状态由于T最多444所以把每个点控制在O(400000)以内,所以对于n和r最多乘一次因此猜f[n][r],f[r][n],首先一轮一轮的搞不好转移,那么先想一想f[n ...

  5. [HNOI2015]亚瑟王(期望+DP)

    题解 利用期望的线性性,可以把问题转化为求每一个卡牌造成期望的期望值. 然后我们就需要知道每一个卡牌发动技能的概率. 因为当某一张卡牌发动技能时这一轮会结束,这就很难直接计算了. 我们使用DP 设dp ...

  6. BZOJ [HNOI2015]亚瑟王 ——期望DP

    发现每张卡牌最后起到作用只和是否打出去了有关. 而且每张牌打出去的概率和之前的牌打出去的情况有关. 所以我们按照牌的顺序进行DP. 然后记录$i$张牌中打出$j$张的概率,然后顺便统计答案. 直接对系 ...

  7. [HNOI2015]亚瑟王[期望DP]

    也许更好的阅读体验 \(\mathcal{Description}\) 给出\(n\)个技能,每个技能按输入顺序有\(p[i]\)的概率释放并造成\(d[i]\)的伤害.每轮游戏从前往后顺序查看每个技 ...

  8. 概率DP——BZOJ4008 [HNOI2015]亚瑟王

    [HNOI2015]亚瑟王 Description 小 K 不慎被 LL 邪教洗脑了,洗脑程度深到他甚至想要从亚瑟王邪教中脱坑.他决定,在脱坑之前,最后再来打一盘亚瑟王.既然是最后一战,就一定要打得漂 ...

  9. 【BZOJ4008】[HNOI2015]亚瑟王 期望

    [BZOJ4008][HNOI2015]亚瑟王 Description 小 K 不慎被 LL 邪教洗脑了,洗脑程度深到他甚至想要从亚瑟王邪教中脱坑. 他决定,在脱坑之前,最后再来打一盘亚瑟王.既然是最 ...

随机推荐

  1. WSGI学习系列Pecan

    Pecan Introduce Pecan是一个轻量级的基于Python的Web框架, Pecan的目标并不是要成为一个“full stack”的框架, 因此Pecan本身不支持类似Session和D ...

  2. ecshop点击订购、加入按钮没反应的解决方法

    今天做ecshop站的时候,测试数据,发现点击订购.加入按钮都没反应,网上搜索,有些人说是修改了common.js,我将原始版本复原也没反映.后来重新安装ecshop,仔细研究发现,原来头部文件pag ...

  3. SUN巡检命令

    # hostname (主机名)# hostid# uname -X# uname -a # w (进程)# who# last# ps -eaf# /usr/ucb/ps -aux# prstat ...

  4. supermarket SSM

    1.数据库 2.整体空架构(jar包) 3.工程依赖 24节点 build节点 <properties> </project.build.sourceEncoding> < ...

  5. HTML知识点梳理

  6. jquery-validate插件

    jQuery Validation 插件 优点:1.表单验证非常简单方便,并且提供了许多配置项目2.国际化,可以自定义提示信息 命令行安装 //初始化bowerbower init //使用bower ...

  7. 一道Java集合框架题

    问题:某班30个学生的学号为20070301-20070330,全部选修了Java程序设计课程,给出所有同学的成绩(可用随机数产生,范围60-100),请编写程序将本班各位同学的成绩按照从低到高排序打 ...

  8. Python开发环境Wing IDE如何使用GTK和PyGObject

    Wing IDE是一个集成开发环境,可用于编辑.测试和调试使用PyGObject为GTK编写的Python代码.Wing IDE提供自动完成.调用提示.一个强大的调试器,以及许多其他功能,可帮助用户编 ...

  9. REP-0118:can not create temporary file(无法创建临时文件)

    解决办法: 查看一下注册表里面的reports_tmp 的路径 HKEY_LOCAL_MACHINE\SOFTWARE\Wow6432Node\ORACLE\KEY_DevSuiteHome1 .是不 ...

  10. react-native 视频播放器(很不错哦)

    第一步: npm i -S react-native-af-video-player(安装前:先安装: react-native-video.react-native-keep-awake.react ...