正题

题目链接:https://www.luogu.com.cn/problem/CF605E

题目大意

给出nnn个点的一张完全有向图，每一天iii到jjj的路径有pi,jp_{i,j}pi,j的概率出现。
询问从111出发走到nnn在最优策略下的期望天数。

1≤n≤103,pi,i=11\leq n\leq 10^3,p_{i,i}=11≤n≤103,pi,i=1

解题思路

设gig_igi表示从iii到nnn的答案，那么ggg满足式子
gi=∑j=1ngj∏gk<gj(1−pj,k)×pi,j×∏gk<gj(1−pi,k)g_i=\sum_{j=1}^{n}\frac{g_{j}}{\prod_{g_k<g_j}(1-p_{j,k})}\times p_{i,j}\times \prod_{g_k<g_j}(1-p_{i,k})gi=j=1∑n∏gk<gj(1−pj,k)gj×pi,j×gk<gj∏(1−pi,k)
注意到里面有条件gk<gjg_k<g_jgk<gj，但是我们并不知道ggg，好像就死循环了。

但是我们可以知道gng_ngn肯定是最小的而且是000，然后剩下第二小的只会受到gng_ngn的影响，也就是这个里面有一个的ggg是完整的。

可以猜测这个第二小的肯定是剩下的ggg里面最小的那个，因为如果不是，那么显然不满足我们最小化的条件，从一个更大的走到了本应该可以更小的。

所以我们每次剩下的点中找到一个最小的gj∏gk<gj(1−pj,k)\frac{g_{j}}{\prod_{g_k<g_j}(1-p_{j,k})}∏gk<gj(1−pj,k)gj来更新它对其他点的影响即可。

时间复杂度O(n2)O(n^2)O(n2)

code

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1100;
int n;bool v[N];
double p[N][N],prod[N],f[N];
int main()
{scanf("%d",&n);for(int i=1;i<=n;i++){for(int j=1;j<=n;j++){scanf("%lf",&p[i][j]);p[i][j]/=100.0;}}if(n==1)return puts("0")&0;for(int i=1;i<n;i++)prod[i]=1-p[i][n],f[i]=1;v[n]=1;while(1){double low=1e9;int x;for(int i=1;i<=n;i++)if(f[i]/(1-prod[i])<low&&!v[i])low=f[i]/(1-prod[i]),x=i;if(x==1)return printf("%.10lf",low)&0;v[x]=1;for(int i=1;i<=n;i++){if(v[i])continue;f[i]+=low*p[i][x]*prod[i];prod[i]*=(1-p[i][x]);}}return 0;
}

CF605E-Intergalaxy Trips【期望dp】相关推荐

luogu P4745 [CERC2017]Gambling Guide（期望DP + 最短路实现）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 P4745 [CERC2017]Gambling Guide(期望DP + 最短路实现) Weblin ...
[SCOI2008]奖励关（期望dp）
你正在玩你最喜欢的电子游戏,并且刚刚进入一个奖励关.在这个奖励关里,系统将依次随机抛出k次宝物,每次你都可以选择吃或者不吃(必须在抛出下一个宝物之前做出选择,且现在决定不吃的宝物以后也不能再吃). 宝 ...
【bzoj4318】OSU! 期望dp
题目描述 osu 是一款群众喜闻乐见的休闲软件. 我们可以把osu的规则简化与改编成以下的样子: 一共有n次操作,每次操作只有成功与失败之分,成功对应1,失败对应0,n次操作对应为1个长度为n的0 ...
【loj6342】跳一跳期望dp
题目描述一个人从 $1$ 开始向 $n$ 跳,在 $i$ 时会等概率跳到 $i,i+1,...,n$ 之一.求从 $1$ 跳到 $n$ 的期望步数. $n\le 10^7$ . 题解期望dp傻逼题 ...
【高斯消元】兼【期望dp】例题
[总览] 高斯消元基本思想是将方程式的系数和常数化为矩阵,通过将矩阵通过行变换成为阶梯状(三角形),然后从小往上逐一求解. 如:$3X_1 + 2X_2 + 1X_3 = 3$ $ ...
Luogu P3251 [JLOI2012]时间流逝期望dp
题面题面题解期望$dp$好题! 今年$ZJOI$有讲过这题... 首先因为$T$只有$50$,大力$dfs$后发现,可能的状态数最多只有$20w$左右,所以我们就可以大力 ...
【bzoj5197】[CERC2017]Gambling Guide 期望dp+堆优化Dijkstra
题目描述给定一张n个点,m条双向边的无向图. 你要从1号点走到n号点.当你位于x点时,你需要花1元钱,等概率随机地买到与x相邻的一个点的票,只有通过票才能走到其它点. 每当完成一次交易时,你可以选择 ...
【BZOJ】1076 [SCOI2008]奖励关期望DP+状压DP
[题意]n种宝物,k关游戏,每关游戏给出一种宝物,可捡可不捡.每种宝物有一个价值(有负数).每个宝物有前提宝物列表,必须在前面的关卡取得列表宝物才能捡起这个宝物,求期望收益.k<=100,n&l ...
string （KMP+期望DP）
Time Limit: 1000 ms Memory Limit: 256 MB Description 给定一个由且仅由字符 'H' , 'T' 构成的字符串$S$. 给定一个最初为空的字符串 ...
概率期望题（期望 DP）做题记录
概率期望题(期望 DP)做题记录 P3830 [SHOI2012]随机树难点在于第二问:生成树的期望深度. 不 wei zhuo 捏,设 $dp_{i,j}$ 表示已经有了 $i$ 个叶子结 ...

CF605E-Intergalaxy Trips【期望dp】

正题

题目大意

解题思路

code

CF605E-Intergalaxy Trips【期望dp】相关推荐

最新文章

热门文章