-概率DP

已于 2022-05-27 18:31:01 修改
阅读量229 收藏
点赞数 2
分类专栏： # 动态规划 文章标签： 算法 c++
于 2022-05-08 22:15:16 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/PURSUE__LIFE/article/details/124655713
版权

目录

1、走路1

2、走路2

3、走路3

4、瓜子

5、走路4

6、守卫者的挑战

1、走路1

// problem： 走路1

#include <bits/stdc++.h>
using namespace std;
#define ll long long
typedef pair<int, int> PII;
double f[105];
vector<int> c[105];
int n, m;
int main(){
	scanf("%d %d",&n, &m);
	for(int i = 1; i <= m; ++i){
		int x, y;
		scanf("%d %d", &x, &y);
		c[x].push_back(y);
	}
	memset(f, 0, sizeof(f));
	f[1] = 1;
	for(int i = 1; i < n; ++i){
		int l = c[i].size();
		for(int j = 0; j < l; ++j)
			f[c[i][j]] += f[i] / l;
	}

	printf("%.10lf\n", f[n]);

	return 0;
}
/*
相信某些人可能存在某些疑问，为什么i号城市是从 1 - n 开始枚举的，可能1号城市
一开始只走向了 5 号城市呀。 为什么不用 dfs搜索？  用dfs搜索明显更好，但是如果你认为 循环
枚举有问题的话， 其实你这还是没有理解好题意。
就算一开始只走向了 5 号城市，那说明 1 -> 2 是走不通的，概率自然也为 0
*/

2、走路2

// problem :  走路2
/*
知识点：
(a / b) mod p  ==>  a * b^(p-2)  mod p         费马小定理

b^(p-2)  mod p  可以用快速幂求
*/
#include <bits/stdc++.h>
using namespace std;
#define ll long long
typedef pair<int, int> PII;
ll f[105];
vector<int> c[105];
int n, m;
const int p = 1e9 + 7;

ll rp(ll now, int k){
	ll res = 1;
	for(; k; k >>= 1, now *= now, now %= p)
		if(k & 1)
			res *= now, res %= p;
	return res;
}

int main(){
	scanf("%d %d",&n, &m);
	for(int i = 1; i <= m; ++i){
		int x, y;
		scanf("%d %d", &x, &y);
		c[x].push_back(y);
	}
	memset(f, 0, sizeof(f));
	f[1] = 1;
	for(int i = 1; i < n; ++i){
		int l = c[i].size();
		ll k = rp(l, p - 2);
		for(int j = 0; j < l; ++j)
			f[c[i][j]] += f[i] * k,
			f[c[i][j]] %= p;
	}

	printf("%lld\n", f[n]);

	return 0;
}

3、走路3

// problem :  走路3

#include <bits/stdc++.h>
using namespace std;
#define ll long long
typedef pair<int, int> PII;
ll f[105][105];  //位于  i 号城市，走了 j 条路的概率
vector<int> c[105];
int n, m;
const int p = 1e9 + 7;

ll rp(ll now, int k){
	ll res = 1;
	for(; k; k >>= 1, now *= now, now %= p)
		if(k & 1)
			res *= now, res %= p;
	return res;
}

int main(){
	scanf("%d %d",&n, &m);
	for(int i = 1; i <= m; ++i){
		int x, y;
		scanf("%d %d", &x, &y);
		c[x].push_back(y);
	}
	memset(f, 0, sizeof(f));
	f[1][0] = 1; // f 存的是概率 还不是期望
	for(int i = 1; i < n; ++i){
		int l = c[i].size();
		for(int j = 0; j < n; ++j)
			if(f[i][j])
				for(int k = 0; k < l; ++k)
					f[c[i][k]][j + 1] += f[i][j] * rp(l, p - 2),
					f[c[i][k]][j + 1] %= p;
	}

	ll ans = 0;
	for(int i  = 1; i < n; ++i)
		ans += f[n][i] * i, ans %= p;
	printf("%lld\n", ans);
	return 0;
}

可以进行优化， 从后往前求

 用f[x] 表示从x号城市出发  期望  经过多少条高速公路能够到达n号城市，f[n] = 0;
  x - > y  (1步)    y - > n    (f[y])
 求概率，起点到终点    求期望，终点到起点

// problem :  走路3  快速版
// 用f[x] 表示从x号城市出发  期望  经过多少条高速公路能够到达n号城市，f[n] = 0;
//   x - > y  (1步)    y - > n    (f[y])
// 求概率，起点到终点    求期望，终点到起点

#include <bits/stdc++.h>
using namespace std;
#define ll long long
typedef pair<int, int> PII;
ll f[105];  //位于  i 号城市，走了 j 条路的期望
vector<int> c[105];
int n, m;
const int p = 1e9 + 7;

ll rp(ll now, int k){
	ll res = 1;
	for(; k; k >>= 1, now *= now, now %= p)
		if(k & 1)
			res *= now, res %= p;
	return res;
}

int main(){
	scanf("%d %d",&n, &m);
	for(int i = 1; i <= m; ++i){
		int x, y;
		scanf("%d %d", &x, &y);
		c[x].push_back(y);
	}
	memset(f, 0, sizeof(f));
	for(int i = n - 1; i ; --i){  // 数据保证所有城市都可以到 n 号城市。
		int l = c[i].size();
		f[i] = 1;
		for(int j = 0; j < l; ++j)
			f[i] += f[c[i][j]] * rp(l, p - 2),
			f[i] %= p;
	}

	printf("%lld\n", f[1]);
	return 0;
}

4、瓜子

// problem : 瓜子 

#include <bits/stdc++.h>
using namespace std;
#define ll long long
typedef pair<int, int> PII;
ll f[2001][4005], a[4005];
const int p = 998244353;
int n;
ll rp(ll now, int k){
	ll res = 1;
	for(; k; k >>= 1, now *= now, now %= p)
		if(k & 1)
			res *= now, res %= p;
	return res;
}

int main(){
	scanf("%d",&n);
	for(int i = 1; i <= 2 * n; ++i)
		a[i] = rp(i, p - 2);
	// f[0][j] = 0;
	for(int i = 1; i <= n; ++i)
		for(int j = 0; j <= 2 * (n - i); ++j){
			f[i][j] = 1;
			f[i][j] += f[i - 1][j + 2] * i % p * a[i + j] % p;
			f[i][j] %= p;
			if(j)
				f[i][j] += f[i][j - 1] * j % p * a[i + j] % p,f[i][j] %= p;

		}
	printf("%lld\n", f[n][0]);
	return 0;
}

5、走路4

 这题的x和y没有大小之间的关系，所以可能存在环。如果存在环，就不满足最优子结构，就不用使用DP做。

该题算是模板题，采用 高斯消去

// problem : 走路4     有环  不满足最优子结构  不能使用动态规划
// 高斯消元 帮你忙！

#include <bits/stdc++.h>
using namespace std;
#define ll long long
typedef pair<int, int> PII;
int n, m;
const int p = 1e9 + 7;
vector<int> c[105];
ll f[105][105], v[105], a[105];
ll rp(ll now, int k){
	ll res = 1;
	for(; k; k >>= 1, now *= now, now %= p)
		if(k & 1)
			res *= now, res %= p;
	return res;
}
inline void gsxy(){
	for(int i = 1; i <= n; ++i){
		for(int j = i; j <= n; ++j)
			if(f[j][i]){
				for(int k = i; k <= n; ++k)
					swap(f[i][k], f[j][k]);
				swap(v[j], v[i]);
			}
		for(int j = i + 1; j <= n; ++j)
			if(f[j][i]){
				ll delta = f[j][i] * rp(f[i][i], p - 2) % p;
				for(int k = i; k <= n; ++k){
					f[j][k] -= f[i][k] * delta % p;
					if(f[j][k] < 0)
						f[j][k] += p;
				}
				v[j] -= v[i] * delta % p;
				if(v[j] < 0)
					v[j] += p;
			}
	}
	for(int i = n; i; --i){
		for(int j = i + 1; j <= n; ++j){
			v[i] -= f[i][j] * a[j] % p;
			if(v[i] < 0)
				v[i] += p;
		}
		a[i] = v[i] * rp(f[i][i], p - 2) % p;
	}
}
int main(){
	scanf("%d %d", &n, &m);
	for(int i = 1; i <= m; ++i){
		int x, y;
		scanf("%d %d", &x, &y);
		c[x].push_back(y);
	}
	for(int i = 1; i < n; ++i){
		f[i][i] = 1; v[i] = 1;
		int l = c[i].size();
		for(int j = 0; j < l; ++j)
			f[i][c[i][j]] = p - 1 * rp(l, p - 2);
	}
	f[n][n] = 1;
	gsxy();
	printf("%lld\n", a[1]);

	return 0;
}

6、守卫者的挑战

232. 守卫者的挑战 - AcWing题库https://www.acwing.com/problem/content/234/

分析：第一次做概率DP，题目懂，但是上手就不会写。而且这题的背包可能一开始会变成负数，最后再变成正数，最终也是满足条件的。考虑到这种情况，我们就要将背包的容量进行偏移。

如果有n个任务，那么最极端的情况下，背包会成 -n 的状态，所以我们的偏移量设成 -n 就够了；再考虑到一种情况，如果我们的k很大，其实相当于是浪费背包空间的一种情况，因为我们有n个任务，所以我们真实的背包也只需要n的容量即可。并且有一个最重要的前提就是找到dp的状态方程表示：

f[i][j][k] 表示经过 i 关，通过了 j 关，还剩下 k 的背包容量。

#include<bits/stdc++.h>
using namespace std;
int n,l,k;
const int N = 205;
double f[N][N][N<<1],p[N];
int main(){
    scanf("%d %d %d",&n,&l,&k);
    if(k>n)k = n;
    f[0][0][n+k] = 1;// 偏移n个单位,防止背包先出现负数，再出现正数的情况
    
    for(int i = 1;i<=n;++i){
        scanf("%lf",&p[i]);
        p[i]/=100;
    }
    for(int i = 1;i<=n;++i){
        int a;scanf("%d",&a);// 第i关通过后获得的背包
        for(int j = 0;j<i;++j){// 通关的个数 
            for(int v = 0;v<=n+n;++v){
                if(v+a>=0)f[i][j+1][min(n+n,v+a)] += p[i] * f[i-1][j][v];
                f[i][j][v] += (1-p[i]) * f[i-1][j][v];
            }
        }
    }
    double ans = 0;
    for(int i = l;i<=n;++i)
    for(int j = n;j<=n+n;++j)ans += f[n][i][j];
    printf("%.6lf\n",ans);
    
    return 0;
}

xingxg.
关注
  • 2
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
tyvj P1864 [Poetize I]守卫者的挑战
我要吃熊猫
07-05 617
数学期望DP
【动态规划】守卫者的挑战
wu_yihao的专栏
08-20 1561
守卫者的挑战 (guard.pas/c/cpp) 题目描述 打开了黑魔法师Vani的大门，队员们在迷宫般的路上漫无目的地搜寻着关押applepi的监狱的所在地。突然，眼前一道亮光闪过。“我，Nizem，是黑魔法圣殿的守卫者。如果你能通过我的挑战，那么你可以带走黑魔法圣殿的地图……”瞬间，队员们被传送到了一个擂台上，最初身边有一个容量为的包包。 擂台赛一共有项挑战，各项挑战依次进行。第项挑战
AcWing 232 守卫者的挑战(概率dp)
weixin_48024524的博客
07-28 98
但是再看一下，只有n次比赛，每次只加1，所以不用开那么大，开到400就可以了。这时候考虑数据，因为k很大，所以要是开200*200*2000肯定爆。重要的是double初始化不要用memeset。然后考虑dp的转移方程三个状态，i，j,k。进行i比赛，j次成功，还剩k个空位的意思。然后再仔细想一想，转移方程就是。...
概率dp
中转站
03-17 131
为了求解问题的期望值，我们需要用到概率公式求解； 但当情况复杂时，无法直接计算得到； 所有引入dp的概念来计算概率： 即通过前一步概率得到下一步概率，把计算结果的总概率分解成计算部分的概率。 用 dp 来dfs时，如果是从初始递推到末状态，会经过许多无效状态； 我们可以从末状态 反推 到初始状态。 dp[i][j]的定义是从[i][j]到末状态的期望。 dp的反推公式是： d...
Acwing 232. 守卫者的挑战
qq_35975367的博客
07-28 244
Acwing 232. 守卫者的挑战 题意： 有n个挑战，一开始背包容量为k，每次挑战有p[i]的概率成功，成功的话会得到一个大小为1的地图碎片或者是提升背包容量X，所有的地图碎片必须装在包里，问最后带地图离开的概率 题解： 设f(i,j,k)表示前i个挑战赢了j次，剩余k次的概率 这么开数组的话会爆炸(200 * 200 * 2000),但是我们知道n个挑战也就最多只有x个地图碎片可以拿，那么只要容量大于n就必然成功。所以记f(i,j,n)表示前i个挑战赢了j次，剩余容量等于n的概率，既然是剩余空间，值域
算法-动态规划- 概率 DP 与期望 DP（包含源程序）.rar
09-16
算法-动态规划- 概率 DP 与期望 DP（包含源程序）.rar
概率DP
09-04
2014清华梦之队NOIP考前提高班概率DP课件
数学期望和概率DP--2021.09.06(E).pdf
09-06
数学期望和概率DP--2021.09.06(E).pdf
区间DP概率DP树形DP插头DP
09-05
区间DP概率DP树形DP插头DP，每种DP一道典型例题，有助于初学者
概率dp 树形dp经典题目加解析
04-28
概率dp 树形dp经典题目加解析........高斯求期望...............
守卫者的挑战
weixin_34198881的博客
05-30 111
守卫者的挑战 您初始有一个容量为k的包，现在有n项挑战依次进行，第i项挑战成功的概率为\(p_i\%\),每次挑战成功会得到一个\(a_i\),\(a_i=-1\)表示得到一个地图碎片，否则表示得到一个容量为\(a_i\)的包，胜利的条件是至少成功l次，且得到的背包恰好能容纳所有的地图碎片（一个地图碎片占据容量1），问胜利的概率，\(0<=k<=2000,0<=n<=200...
0x3B 守卫者的挑战 （概率dp）
Hacheylight
02-04 238
蛤很明显题目都是叫你去求概率了 /cy 这个问题特别像背包 于是我们能够列出一个 dpdpdp 状态： dp[i][j][k]dp[i][j][k]dp[i][j][k] 表示前 iii 个回合赢了 jjj 次背包空间还剩下 kkk 的概率 显然分成赢了和输了两种情况： 赢了： dp[i+1][j+1][k+a[i+1]]+=dp[i][j][k]∗p[i+1]dp[i+1][j+1][k+a[i...
bzoj3029 守卫者的挑战（概率dp）
zmy_wky的博客
03-20 268
Description 　　打开了黑魔法师Vani的大门，队员们在迷宫般的路上漫无目的地搜寻着关押applepi的监狱的所在地。突然，眼前一道亮光闪过。“我，Nizem，是黑魔法圣殿的守卫者。如果你能通过我的挑战，那么你可以带走黑魔法圣殿的地图……”瞬间，队员们被传送到了一个擂台上，最初身边有一个容量为K的包包。 　　擂台赛一共有N项挑战，各项挑战依次进行。第i项挑战有一个属性aiaia_i，...
Codevs P1997 守卫者的挑战
Starlight
10-22 493
Codevs P1997 守卫者的挑战题目描述 Description　　打开了黑魔法师Vani的大门，队员们在迷宫般的路上漫无目的地搜寻着关押applepi的监狱的所在地。突然，眼前一道亮光闪过。“我，Nizem，是黑魔法圣殿的守卫者。如果你能通过我的挑战，那么你可以带走黑魔法圣殿的地图……”瞬间，队员们被传送到了一个擂台上，最初身边有一个容量为K的包包。 　　擂台赛一共有项挑战，各项挑战依次进
bzoj3029 守卫者的挑战
YHaX的博客
07-02 1500
数学期望dp
bzoj 3029 守卫者的挑战
heheda_is_an_OIer的博客
11-16 688
水题，唯一值得提的是我们不用维护背包剩余空间为O（K）的情况，只要维护最大到n的情况就行了。因为只要剩余空间到了n，无论如何也不会装满。#include#include#include#include#include#define md#define ll long long#define inf (int) 1e9#define eps 1e-8#define N 250using namesp
BZOJ2435 守卫者的挑战 - 期望dp
Mys_C_K的博客
08-02 399
传送门 题解： 考虑dp。dp[i][j][k]表示当前第i场，已经赢了j场，背包容量为k（可能为负）的概率。 转移显然。可以用n+1表示背包容量无穷大（因为显然碎片最多n个，容量太大就变成无穷大了）。 代码： #include //BZOJ 3029 #include #include #define MAXN 210 using namespace std; double dp[
守卫者的挑战【概率DP】【处理内存、偏移量】
谁是凶手1703
10-27 233
守卫者的挑战 题意 n项挑战，每项挑战，有获胜的概率，每项挑战的收益有正数、-1.问这n项挑战中，获胜l场以上，并且收益为>=0的概率。 思路 n个物品中选择>=k个，01背包的思路，然后有物品的个数限制，多加一维度表示拿了几个物品，有收益的限制，再多加一维。因此状态表示为：f[i] [j] [k] :前i项挑战，获胜j场，收益为k的概率。 题意要求，转移过程中可以有负数，意味着收益可以为 负数。题目要求我们求最后收益>=0的概率，因此我们把一定收益>=0的收益概率压缩到收益为n（为
期望dp和概率dp的区别
最新发布
07-29
期望dp和概率dp是两种不同的动态规划方法。 期望dp是指通过计算每个状态的期望值来求解最终的期望。在期望dp中，我们通常定义dp\[i\]表示在第i个状态时的期望值，然后通过状态转移方程来更新dp数组，最终得到最终...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

xingxg.

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值