【BZOJ4008】【HNOI2015】亚瑟王（DP，期望概率）

最新推荐文章于 2018-09-02 20:29:39 发布

原创最新推荐文章于 2018-09-02 20:29:39 发布 · 286 阅读

1 ·

CC 4.0 BY-SA版权

BZOJ 同时被 3 个专栏收录

215 篇文章

订阅专栏

省选

90 篇文章

订阅专栏

83 篇文章

订阅专栏

本文介绍了一种计算特定游戏模型中伤害期望值的方法，采用动态规划算法解决了一个简化版亚瑟王游戏中的问题，旨在帮助玩家最大化游戏中的伤害输出。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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 求出这一套卡牌在一局游戏中能造成的伤害的期望值。

Solution

考虑DP。
我们设 $f_{i,j}$ 表示前 $i$ 张牌有 $j$ 张发动的概率， $g_i$ 表示第 $i$ 张牌发动的概率，那么答案就是 $\sum_{i=1}^n g_i d_i$ 。

我们现在来考虑如何转移。
1. 从 $f_{i-1,j}$ 转移，那么第 $i$ 张牌有 $r-j$ 次机会被发动但是最终没有发动，所以 $f_{i,j}+=f_{i-1,j}\times (1-p_i)^{r-j}$
2. 从 $f_{i-1,j_1}$ 转移，那么第 $i$ 张牌有 $r-j+1$ 次机会被发动，其中有一次成功了，所以 $f_{i,j}+=f_{i-1,j-1}\times(1-(1-p_i)^{r-j+1})$
3. $g_i=\sum_{j=0}^r f_{i-1,j}\times (1-(1-p_i)^{r-j})$

然而zjp有更简短的转移，真的太爷了！

Code

/************************************************
 * Au: Hany01
 * Date: Apr 9th, 2018
 * Prob: [BZOJ4008][HNOI2015] 亚瑟王
 * Email: hany01@foxmail.com
************************************************/

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;
typedef pair<int, int> PII;
#define File(a) freopen(a".in", "r", stdin), freopen(a".out", "w", stdout)
#define rep(i, j) for (register int i = 0, i##_end_ = (j); i < i##_end_; ++ i)
#define For(i, j, k) for (register int i = (j), i##_end_ = (k); i <= i##_end_; ++ i)
#define Fordown(i, j, k) for (register int i = (j), i##_end_ = (k); i >= i##_end_; -- i)
#define Set(a, b) memset(a, b, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define x first
#define y second
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define ALL(a) (a).begin(), (a).end()
#define SZ(a) ((int)(a).size())
#define INF (0x3f3f3f3f)
#define INF1 (2139062143)
#define Mod (1000000007)
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define y1 wozenmezhemecaia

template <typename T> inline bool chkmax(T &a, T b) { return a < b ? a = b, 1 : 0; }
template <typename T> inline bool chkmin(T &a, T b) { return b < a ? a = b, 1 : 0; }

inline int read()
{
    register int _, __; register char c_;
    for (_ = 0, __ = 1, c_ = getchar(); c_ < '0' || c_ > '9'; c_ = getchar()) if (c_ == '-') __ = -1;
    for ( ; c_ >= '0' && c_ <= '9'; c_ = getchar()) _ = (_ << 1) + (_ << 3) + (c_ ^ 48);
    return _ * __;
}

const int maxn = 233, maxr = 133;

int main()
{
#ifdef hany01
    File("bzoj4008");
#endif

    static double p[maxn][maxr], g[maxn], f[maxn][maxr];
    static int d[maxn];

    for (static int T = read(); T --; )
    {
        register int n = read(), r = read();
        For(i, 1, n) {
            scanf("%lf", &p[i][1]), d[i] = read(), p[i][0] = 1, p[i][1] = 1 - p[i][1];
            For(j, 2, r) p[i][j] = p[i][j - 1] * p[i][1];
        }

        Set(f, 0), Set(g, 0);
        f[1][0] = p[1][r], f[1][1] = g[1] = 1 - p[1][r];
        register double Ans = g[1] * d[1];
        For(i, 2, n) {
            For(j, 0, r) {
                g[i] += f[i - 1][j] * (1 - p[i][r - j]);
                f[i][j] += f[i - 1][j] * p[i][r - j];
                if (j) f[i][j] += f[i - 1][j - 1] * (1 - p[i][r - j + 1]);
            }
            Ans += g[i] * d[i];
        }
        printf("%.10lf\n", Ans);
    }

    return 0;
}
//自在飞花轻似梦，无边丝雨细如愁。
//    -- 秦观《浣溪沙·漠漠轻寒上小楼》