概率/期望dp专题

已于 2022-08-30 22:00:40 修改
阅读量670 收藏 4
点赞数
分类专栏： ACM学习笔记 文章标签： 算法 c++ 蓝桥杯
于 2022-08-29 14:59:16 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_34416123/article/details/126585094
版权

Happy Running（推公式）

画一个二维的坐标图，推公式即可。

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;
 
int main()
{
    int T; scanf("%d", &T);
    while(T--)
    {
        double k, x;
        scanf("%lf%lf", &k, &x);
        if(k == x) puts("0.50");
        else if(k < x) printf("%.2f\n", 0.5 * (x * x - k * k) / (x * x) + 0.5);
        else if(x >= k - x) printf("%.2f\n", 0.5 * pow(2 * x - k, 2) / (x * x));
        else puts("0.00");
    }
	return 0;
}

恶意竞争（期望dp）

期望dp入门题

概率dp一般是设已经算了多少的情况下，剩余的期望。这样的话结果就是0，然后递推到初始状态。方程可能会原地踏步，自己推自己，这个时候就需要移项。

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

const int N = 1e3 + 10;
double dp[N][N];
int n, s;
 
int main()
{
    scanf("%d%d", &n, &s);
    for(int i = n; i >= 0; i--)         //逆推到0
        for(int j = s; j >= 0; j--)
        {
            if(i == n && j == s) continue;   //注意这里，不加的话会除以0
            double a = 1.0 * i / n, b = 1.0 * j / s;
            dp[i][j] = 1 / (1 - a * b) * (a * b + (1 - a) * b * (1 + dp[i + 1][j])
                        + a * (1 - b) * (1 + dp[i][j + 1]) + 
                        (1 - a) * (1 - b) * (1 + dp[i + 1][j + 1]));
        }
    printf("%.10f\n", dp[0][0]);
    
	return 0;
}

带富翁（期望dp）

期望dp一般倒推比较好想，正推也行，但会比较复杂。倒推指的是从当前状态到结束状态的期望，而不是像其他dp那样从初始状态到当前状态的期望

dp[i]表示从i位置走到n位置的期望分数

那么dp[n] = a[n]

dp[i] = a[i] + dp[j] / len

j是所有能达到的点，len是能到达的点的个数

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

const int N = 1e2 + 10;
double dp[N], a[N];
int n;
 
int main()
{
    scanf("%d", &n);
    _for(i, 1, n) scanf("%lf", &a[i]);

    dp[n] = a[n];
    for(int i = n - 1; i >= 1; i--)
    {
        dp[i] = a[i];
        int len = min(i + 6, n) - i;
        _for(j, i + 1, i + len)
            dp[i] += dp[j] / len;
    }
    printf("%.10f\n", dp[1]);
    
	return 0;
}

筛子游戏（期望dp+处理f[0]）

这题的难点在于dp方程里面有一个f[0]，这个f[0]是我们所求的，也是一个常数

处理的技巧就是设f[i] = A[i]f[0] + B[i],代入dp方程，得到A和B的递推式子

注意题目求的是大于n的,也就是f[n+1]=0,所以要从A[n],B[n]开始求。

递推求得A[0],B[0],由f[0] = A[0]f[0] + B[0]解出f[0]即可。

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

const int N = 600;
int n, a, b, c, k1, k2, k3, K;
double p[20], A[N], B[N];
 
int main()
{
    scanf("%d%d%d%d%d%d%d", &n, &k1, &k2, &k3, &a, &b, &c);
    _for(i, 1, k1)
        _for(j, 1, k2)
            _for(k, 1, k3)
            {
                if(i == a && j == b && k == c) p[0]++;
                else p[i + j + k]++;
            }
    K = k1 + k2 + k3;
    _for(i, 0, K) p[i] /= k1 * k2 * k3;

    for(int i = n; i >= 0; i--)
    {
        A[i] = p[0];
        B[i] = 1;
        _for(k, 1, K) 
        {
            A[i] += A[i + k] * p[k];
            B[i] += B[i + k] * p[k];
        }
    }
    printf("%.10f\n", B[0] / (1 - A[0]));
    
	return 0;
}

Balls（概率dp）

概率dp正推，期望dp逆推。这题先算概率，其实手推就可以发现是1+n/2了，因为每次期望增加是0.5。用概率dp的话，可以dp[i][j]表示i个黑球，j个白球的概率，然后由dp[i][j]更新dp[i+1][j]和dp[i][j + 1]即可。算出概率后枚举白球个数算期望，也可以得出答案就是1+n/2

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;
 
int main()
{
    int n; scanf("%d", &n);
    printf("%.7f\n", 1 + n * 0.5);
    
	return 0;
}

珂学送分（期望dp）

同样倒推，用dp[i]表示以i为起点所能得到的期望，这样dp[n] = 1

然后枚举当前的分段即可,可以预处理加速。

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

const int N = 1e5 + 10;
int a[N], n, m, x, p[N], mx;
double dp[N], s[N];
 
int main()
{
    scanf("%d%d", &n, &m);
    _for(i, 1, n) scanf("%d", &a[i]), mx = max(mx, a[i]);
    while(m--)
    {
        scanf("%d", &x);
        if(x < mx)
        {
            puts("YNOI is good OI!");
            continue;
        }

        int sum = 0, r = 0;
        _for(l, 1, n)
        {
            while(r + 1 <= n && sum + a[r + 1] <= x) r++, sum += a[r];
            p[l] = r;
            sum -= a[l];
        }

        s[n + 1] = s[n + 2] = 0;
        s[n] = dp[n] = 1;
        for(int i = n - 1; i >= 1; i--)
        {
            int l = i, r = p[i];
            dp[i] = 1 + 1.0 / (r - l + 1) * (s[l + 1] - s[r + 2]);
            s[i] = s[i + 1] + dp[i];
        }
        printf("%.2f\n", dp[1]);
    }
    
	return 0;
}

P5104 红包发红包（期望）

手推一下，发现是一直除以2

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

typedef long long ll;
const int mod = 1e9 + 7;
ll w, n, k;

ll binpow(ll a, ll b)
{
    ll res = 1;
    for(; b; b >>= 1)
    {
        if(b & 1) res = res * a % mod;
        a = a * a % mod;
    }
    return res;
}

ll inv(ll x) { return binpow(x, mod - 2); }
 
int main()
{
    scanf("%lld%lld%lld", &w, &n, &k);
    printf("%lld\n", w * inv(binpow(2, k)) % mod);
    
	return 0;
}

P6154 游走（期望）

这题不算是期望dp，用期望的定义可知答案就是所有路径的长度之和除以所有路径条数

dp一下即可

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

typedef long long ll;
const int N = 1e5 + 10;
const int mod = 998244353;
ll dp[N], cnt[N];
vector<int> g[N];
int d[N], n, m;

ll binpow(ll a, ll b)
{
    ll res = 1;
    for(; b; b >>= 1)
    {
        if(b & 1) res = res * a % mod;
        a = a * a % mod;
    }
    return res;
}

ll inv(ll x) { return binpow(x, mod - 2); }
 
int main()
{
    scanf("%d%d", &n, &m);
    while(m--)
    {
        int u, v;
        scanf("%d%d", &u, &v);
        g[u].push_back(v);
        d[v]++;
    }

    queue<int> q;
    _for(i, 1, n) 
    {
        cnt[i] = 1;
        if(!d[i]) q.push(i);
    }
    while(!q.empty())
    {
        int u = q.front(); q.pop();
        for(int v: g[u])
        {
            cnt[v] = (cnt[v] + cnt[u]) % mod;
            dp[v] = (dp[v] + dp[u] + cnt[u]) % mod;
            if(--d[v] == 0) q.push(v);
        }
    }

    ll a = 0, b = 0;
    _for(i, 1, n)
    {
        a = (a + dp[i]) % mod;
        b = (b + cnt[i]) % mod;
    }
    printf("%lld\n", a * inv(b) % mod);
    
	return 0;
}

P1297 [国家集训队]单选错位（期望）

计算每道题对答案的贡献即可

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

const int N = 1e7 + 10;
int a[N], n, A, B, C;
 
int main()
{
    scanf("%d%d%d%d%d", &n, &A, &B, &C, a + 1);
    for (int i = 2; i <= n; i++)
        a[i] = ((long long) a[i - 1] * A + B) % 100000001;
    for (int i = 1; i <= n; i++)
        a[i] = a[i] % C + 1;

    double ans = 0;
    _for(i, 1, n)
    {
        int j = (i == n) ? 1 : i + 1;
        ans += 1.0 * min(a[i], a[j]) / (1LL * a[i] * a[j]);
    }
    printf("%.3f\n", ans);
    
	return 0;
}

SP1026 FAVDICE - Favorite Dice（期望dp）

dp[i]表示从已经出现i种面时，还需要扔色子的次数，dp[n] = 0, dp[0]即是答案

注意写方程时注意当前这次操作需要1次，不要漏写1

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

const int N = 1e3 + 10;
double dp[N];
int n;
 
int main()
{
    int T; scanf("%d", &T);
    while(T--)
    {
        scanf("%d", &n);
        dp[n] = 0;
        for(int i = n - 1; i >= 0; i--)
        {
            double a = 1.0 * i / n;
            dp[i] = a / (1 - a) + 1 + dp[i + 1];
        }
        printf("%.2f\n", dp[0]);
    }
	return 0;
}

P4316 绿豆蛙的归宿（期望dp）

这题直接暴力搜出所有的路径竟然不会T，能过

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

typedef long long ll;
const int N = 1e5 + 10;
vector<pair<int, int>> g[N];
double dp[N], ans;
int n, m;

void dfs(int u, int fa, ll sum, double p)
{
    if(u == n) ans += sum * p;
    for(auto x: g[u])
    {
        int v = x.first, w = x.second;
        if(v == fa) continue;
        dfs(v, u, sum + w, p / g[u].size());
    }
}
 
int main()
{
    scanf("%d%d", &n, &m);
    while(m--)
    {
        int u, v, w;
        scanf("%d%d%d", &u, &v, &w);
        g[u].push_back({v, w});
    }

    dfs(1, 0, 0, 1);
    printf("%.2f\n", ans);
   
	return 0;
}

当然，正解同样是逆推，用dp[u]表示从u到n的期望

#include<bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

const int N = 1e5 + 10;
vector<pair<int, int>> g[N];
double dp[N], ans;
int d[N], n, m;
 
int main()
{
    scanf("%d%d", &n, &m);
    while(m--)
    {
        int u, v, w;
        scanf("%d%d%d", &u, &v, &w);
        g[u].push_back({v, w});
        d[v]++;
    }

    queue<int> q;
    vector<int> topo;
    _for(i, 1, n)
        if(!d[i])
            q.push(i);
    while(!q.empty())
    {
        int u = q.front(); q.pop();
        topo.push_back(u);
        for(auto x: g[u])
            if(--d[x.first] == 0)
                q.push(x.first);
    }

    for(int i = topo.size() - 1; i >= 0; i--)
    {
        int u = topo[i];
        for(auto x: g[u])
        {
            int v = x.first, w = x.second;
            dp[u] += 1.0 / g[u].size() * (w + dp[v]);
        }
    }
    printf("%.2f\n", dp[1]);
   
	return 0;
}

Alex Su (*^▽^*)
关注
专栏目录
浅谈 概率与期望 DP
Monica
09-11 1542
概率与期望DP，一直都不会啊，感觉好难完全没法思考…期望题一般是逆推，当然也有一些是顺推，然而现在我只做过一些水得很的期望题，好菜啊…UVa 11021 麻球繁衍 蓝书p140题解：现在有k只麻球，每只麻球都只能活一天，但是可能会在死前爆出0到n-1值麻球，几率分别为p0,p1,p2...pn−1p_0,p_1,p_2...p_{n-1}，求m天后所有麻球都死了的概率，因为每只麻球是独立的，这只麻球
概率dp 求期望
++ ‘s blog
07-23 1031
先引用大牛文章解释原理
概率及期望DP小结
weixin_30745553的博客
07-26 222
概念 \(PS\)：不需要知道太多概念，能拿来用就行了。 定义 样本（\(\omega\)）：一次随机试验产生的一个结果。 样本空间（\(\Omega\)）：一个随机试验的所有可能的结果的全体，即\(\Omega=\{\omega\}\)。 事件（\(A\)）：某一类结果，即\(A\subset\Omega\)。 基本事件（\(s\)）：各个互斥的事件即为基本事件。 我们借助样本空间...
期望DP
Morning_Glory_JR的博客
07-22 1342
转移方程 几种常见设转移方程数组的方法 设f[i]f[i]f[i]表示的是变成由iii状态变成 最终 状态的期望 按照题意直接设 把选择的东西加入数组，如f[i][j]f[i][j]f[i][j]表示第iii个物品选jjj个的期望或f[i][j]f[i][j]f[i][j]表示有iii个AAA物品，jjj个BBB物品的期望 求转移方程 应优先考虑进行操作后当前状态会怎么样，而不是如何变成当前状...
概率与期望DP
nameofcsdn的博客
01-24 1865
目录 一，概率DP CSU 1123: PK武林盟主 CSU 1342：Double CSU 1725 加尔鲁什·地狱咆哮对阵虚灵大盗拉法姆 二，期望DP LIghtOJ 1038 Race to 1 Again HDU 4405 Aeroplane chess（飞行棋） POJ 2096 Collecting Bugs SGU 495 Kids and Prizes 一，概率DP CSU 1123: PK武林盟主 题目： Description 枫之羽认为自己很强，想当武林盟主
kuangbin专题 概率/期望DP总结
mzip的博客
05-26 440
题目链接: https://vjudge.net/contest/76505#overview A- A Dangerous Maze LightOJ - 1027 题意: 现在在你面前有N个门。每个门要么把你传送出迷宫，要么把你传送来的位置且你的记忆也会回到初始的时候。现在给出每个门传送的时间t，若t为正数，说明该门花费时间t可以将你传送出迷宫，若 t 为负数，说明该门花费时间 t 将你传送...
期望与概率DP专题/高斯消元
Qianyri的博客
08-19 1434
高斯消元 总结：概率和期望问题的求解都要依靠递推式，概率问题按照递推式进行dp，期望问题根据递推式化简的难易程度通过化简再dp或高斯消元求解，该类问题大部分仍属于各个类型的dp求解问题 2018.9.3-概率与期望dp学习 Discovering Gold 一排1到n的格子，每个格子上有黄金 ai ,你最开始在 1 号,每一次投骰子决定到哪一个格子，超出1~n范围则重新投掷，你到了哪个格子...
期望dp和概率dp的区别
最新发布
07-29
- *1* [概率/期望dp专题](https://blog.csdn.net/qq_34416123/article/details/126585094)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_...
【专题总结】概率&期望DP
NExPlain的专栏
07-26 1137
下面是对一系列求概率&期望问题的总结：
kuangbin专题二十一概率&期望总结
qq_34921856的博客
05-30 754
做的时候觉得很难，很多题都是看题解做的，所以放了一段时间。这几天把kuangbin聚聚博客上的3篇论文仔细看了看，再看这些题就感觉思路很显然了。但仍然有几题不能想出来，故这里补一下总结，这也算kuangbin专题的最后一篇博客了。 首先将这几篇论文贴上来，论文真的真的很难懂。。认真看看但不要太深入即可。(高中生太强了。。 《信息学竞赛中概率问题求解初探》 《浅析竞赛中一类数学期望问题的解决方...
期望与概率dp
Jackzhan的博客
08-09 783
概率 概率加法 互不相交的事件，P(A或B发生)=P(A)+P(B) 如果A和B所涵盖的结果有交集，那么P(A或B发生)=P(A)+P(B)-P(A与B同时发生) 概率乘法 在两个 互不干扰的事中，事件A在其中一件事中，事件B在另外一件事中 那么P(A与B同时发生)=P(A)P(B) e.g. 掷两个骰子，P(第一个点数小于3且第二个点数为偶数)=(2/6)(3/6)=1/6 期望 事件A...
概率与期望DP小结。
hyesuixin的博客
07-19 1179
  好久没更新博客，来水一蛤，近三天做了二十道概率题，也算是摸到了一点门道。 首先关于一个很多博主的说到的点，也是我个人认为很重要的一个点，就是概率要顺着推，期望倒着推，我比较赞同以下这篇博客的观点：https://blog.csdn.net/nameofcsdn/article/details/52082746 其次就是关于概率DP状态的设计，概率DP的题往往设计状态比较直观，直接按照题意来...
【POJ 2096】Collecting Bugs（概率与期望+dp）
reverie_mjp的博客
11-15 372
空伤悲，独愁怅，怎解此刻凄凉？
期望dp和概率dp
solemntee的博客
03-09 470
概率正着推，期望倒着推。 期望dp 一、无后效性的期望dp 例1.1 Favorite Dice 甩一个n面的骰子，问每一面都被甩到的次数期望是多少。 令dp[i]dp[i]dp[i]表示目前已经甩到iii个数字，从当前状态到每一面都甩到的期望次数。 首先我们知道dp[n]=0dp[n]=0dp[n]=0，因为此时已经每一面都甩到了。 对于dp[i]dp[i]dp[i]，当前已经掷出iii个字母，掷出重复字母的概率是i/ni/ni/n，掷出新字母的概率是(n−i)/n(n-i)/n(n−i)/n 所以dp
概率期望dp
weixin_30446197的博客
08-13 145
1.codevs1060 搞笑世界杯 题目描述Description 搞笑世界杯的球票出售方式也很特别,它们只准备了两种球票.A 类 票 免费球票 B 类票 双倍价钱球票.购买时由工作人员通过掷硬币决定,投到正面 的买A类票, 反面的买B类票.并且主办方不可能倒贴钱,所以他们总是准备了同样多的A类票和B类票 你和你的朋友想计算一下排在队尾的两个人同时拿到一种票的概率是多少 (包...
期望DP初步
dk810510的博客
01-14 164
感觉期望DP这种东西像是玄学… 主要总结说一点基础性的东西, 或许对于理解题目的做法会有一点帮助. 首先是关于独立事件, 互斥事件的概念. 通俗地说, 就是对于两个事件A, B, 假如满足发生了其中一个就不会发生另一个, 则称A, B为互斥事件; 假如A与B的发生没有任何关系, 可能都发生, 可能只有一个发生, 也可能都不发生, 则称A, B为独立事件. 下文的...
期望&概率dp总结
十分残念的博客
01-25 4954
总算刷完kuangbin期望&amp;概率专题了，下面总结一下心得和题解！1.期望dp期望dp通常逆推，即从结果推向初始状态，也可以用记忆化搜索进行dp；E=Σp1*(E1+X1)+Σp2*(E+X2)其中E为当前状态的期望，E1为下一个状态的期望，p1和X1分别为将当前状态转移到下一个状态的概率和花费，p2和X2分别为保持当前状态的概率和花费。最后化简为E=(Σp1*(E1+X1)+Σp2*X2...
概率DP/期望DP总结
热门推荐
myjs999的博客
07-12 1万+
概述 一般来说，概率DP找到正确的状态定义后，转移是比较容易想到的。但状态一定是“可数”的，把有范围的整数作为数组下标。事实上，将问题直接作为状态是最好的。如问“n人做XX事的期望次数”，则设计状态为f[i]表示i个人做完事的期望。转移一般是递推，即从上一个状态转移得（填表）或转移向下一个状态（刷表）。 练习题 涂格子1 n个格子，每次随机涂一个，求涂满m个格子的期望次数。 如概述所...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值