【CF113D】Museum【概率期望】【高斯消元】

题意:一张 nnn 个点的无向连通图,两个人开始时分别在 a,ba,ba,b。每次在 uuu 时会以 ppp 的概率原地不动,1−p1-p1p 的概率等概率随机选择到一个相邻的点,当两人在同一点时停止。分别求在每个点相遇的概率。

n≤22n\leq 22n22

网上一堆 “从起点走到 (i,j)(i,j)(i,j) 的概率”,看得我一脸懵逼……

很容易分析出转移矩阵 MMM,然后相当于求这个东西:

lim⁡t→+∞MtV\lim_{t\to +\infin}M^{t}Vt+limMtV

但这个并没有通用的求法,因为矩阵的特征向量有无数多个。

算法一

比较直观的解法。

f(i,j)f(i,j)f(i,j) 表示 (i,j)(i,j)(i,j) 这个状态到达次数的期望,即这个状态在所有世界线中出现次数的平均值。

由于最终点只有可能出现 000 次或 111 次,所以它的期望次数就是概率。

而这个期望随便消一下就可以算出来。

算法二

比较本质的解法。

考虑枚举一个最终状态 (s,s)(s,s)(s,s),在此条件下求 f(i,j)f(i,j)f(i,j) 表示最终到这个状态的概率,令 f(i,i)=[i=s]f(i,i)=[i=s]f(i,i)=[i=s],就可以消元了。但这样是 O(n7)O(n^7)O(n7) 的,无法通过。

考虑一次性把每个点作为终点的 nnn 个答案算出来,即构建出 (n2−n)×(n2)(n^2-n)\times (n^2)(n2n)×(n2) 的矩阵。这样有 n2n^2n2 个未知数但只有 n2−nn^2-nn2n 个方程,无法解出,但可以求出 f(a,b)f(a,b)f(a,b) 关于 f(1,1),f(2,2),…,f(n,n)f(1,1),f(2,2),\dots,f(n,n)f(1,1),f(2,2),,f(n,n) 的线性表达,就可以求出答案了。

复杂度 O(n6)O(n^6)O(n6)

所以算法一算法二以及上面那个假算法写出来都一样的……

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <cmath>
#include <vector>
using namespace std;
vector<int> e[25];
double p[25],a[505][505];
int n,m,sa,sb;
inline int id(int x,int y)
{if (x==y) return n*n-n+x;return (x-1)*(n-1)+y-(y>x);
}
void gauss(int n,int m)
{for (int i=1;i<=n;i++){int pos=i;for (int j=i+1;j<=n;j++) if (fabs(a[j][i])>fabs(a[pos][i])) pos=j;if (pos>i) swap(a[i],a[pos]);for (int j=1;j<=n;j++)if (j!=i){double t=a[j][i]/a[i][i];for (int k=i;k<=m;k++)a[j][k]-=t*a[i][k];}}
}
int main()
{scanf("%d%d%d%d",&n,&m,&sa,&sb);if (sa==sb){for (int i=1;i<=n;i++) if (i==sa) printf("%.10f ",1.0);else printf("%.10f ",0.0);return 0;}for (int i=1;i<=m;i++){int u,v;scanf("%d%d",&u,&v);e[u].push_back(v),e[v].push_back(u);}for (int i=1;i<=n;i++) scanf("%lf",&p[i]);for (int u=1;u<=n;u++)for (int v=1;v<=n;v++)if (u!=v){int s=id(u,v);double pu=1.0/e[u].size(),pv=1.0/e[v].size();for (int i=0;i<=(int)e[u].size();i++)for (int j=0;j<=(int)e[v].size();j++){double t=(i<(int)e[u].size()? (1-p[u])*pu:p[u])*(j<(e[v].size())? (1-p[v])*pv:p[v]);a[s][id(i<(int)e[u].size()? e[u][i]:u,j<(int)e[v].size()? e[v][j]:v)]=t;}a[s][s]-=1;}gauss(n*n-n,n*n);int s=id(sa,sb);for (int i=n*n-n+1;i<=n*n;i++) printf("%.10f ",-a[s][i]/a[s][s]);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314847.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

人生如戏,别太入戏

这里是Z哥的个人公众号每周五早8点 按时送达当然了&#xff0c;也会时不时加个餐&#xff5e;我的第「81」篇原创敬上大家好&#xff0c;我是Z哥。对&#xff0c;就是脑袋很大的那个。大到什么程度呢&#xff1f;我给新来的小伙伴们说说。我昨天还打算配副新眼镜来着&#xff0…

Codeforces Round #588 (Div. 2) F. Konrad and Company Evaluation 图论 + 建反图 好题

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一张nnn个点mmm条边的图&#xff0c;其中每个点iii初始编号为iii&#xff0c;边是有向的&#xff0c;方向为从编号大的指向编号小的。定义一个贡献为存在某三个点a,b,ca,b,ca,b,c有两条边为a−>b,b−…

【SDOI2017】硬币游戏【KMP】【概率期望】【高斯消元】

题意&#xff1a;给 nnn 个长度为 mmm 的 01 串&#xff0c;一个 01 串初始为空&#xff0c;不断随机一个字符加在后面&#xff0c;当出现给定的 nnn 个串中的一个时停止。分别求在 nnn 个串处停止的概率。 考场思路历程&#xff1a; 显然建出 AC 自动机&#xff08;flag&…

C#只用属性来解决兔子,不用方法和循环

属性在调用者看来就像一个普通的变量&#xff0c;但作为类的设计者&#xff0c;你可以利用属性来隐藏你类中的一些字段&#xff0c;使外界只能通过属性来访问你的字段&#xff0c;你可以通过属性来限制外界对你的字段的存取&#xff0c;就利用get、set。如果想让用户随意存取你…

Codeforces Round #730 (Div. 2) D2. RPD and Rap Sheet (Hard Version) 交互 + k进制的转换

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 定义a⊕kba\oplus_k ba⊕k​b为a,ba,ba,b在kkk进制下的不进位加法。系统会随机生成一个数xxx&#xff0c;你猜这个数&#xff0c;假设当前猜的数为yyy&#xff0c;如果猜对了就返回111&#xff0c;否则要猜的…

【ZJOI2019】线段树【线段树上dp】【大讨论】

题意&#xff1a;有一个 [1,n][1,n][1,n] 的线段树和 mmm 个区间赋值操作&#xff0c;求任取一个操作的子集并按顺序在线段树上跑后线段树上有 lazy 标记的点的个数之和 模 998244353998244353998244353。 n,m≤105n,m\leq 10^5n,m≤105 真线段树上 dp 考虑线段树的情况很复杂…

.netcore项目docker化,以及docker之间通信

简言&#xff1a;最近刚完成公司的新系统&#xff0c;系统使用的是微服务架构&#xff0c;由于领导说要将服务docker化。下面将我的研究结果分享出来&#xff0c;如若有错误的地方&#xff0c;还请各位大佬多多指点。目录&#xff1a;  什么是docker&#xff1f;使用docker有…

Codeforces Round #592 (Div. 2) F. Chips 构造 + 细节

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 恶心的构造题&#xff0c;思路很简单但是代码细节很多&#xff0c;搞了半天。 根据题目的性质不难发现&#xff0c;如果有两个相同颜色的球相邻&#xff0c;那么他们的颜色永远不会改变。 …

.NET CORE 多语言实现方案

根据市场需求&#xff0c;基于.NET CORE平台开发的RoadFlow工作流平台需要实现多语言版本。经过方案讨论和比对&#xff0c;决定采用.NET自带的本地化功能来实现多语言。话不多说&#xff0c;直接上实现方式。首先修改Startup.cs在public void ConfigureServices(IServiceColle…

【CTSC2018】暴力写挂【边分治】【虚树】

题意&#xff1a;给两棵基于同一点集的带边权树&#xff0c;记 lca⁡(x,y),depth⁡(x)\operatorname{lca}(x,y),\operatorname{depth}(x)lca(x,y),depth(x) 为第一棵树上的 lca、到根的边长度之和&#xff0c;lca⁡′(x,y),depth⁡′(x)\operatorname{lca}(x,y),\operatorname{…

hdu 6153 A Secret kmp + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个串a,ba,ba,b&#xff0c;让你求对于bbb的每个后缀&#xff0c;设其长度为lenlenlen&#xff0c;其在aaa中出现的次数为cntcntcnt&#xff0c;那么他的贡献为len∗cntlen*cntlen∗cnt&#xff0c;求每…

中国程序员,请挺起你的腰杆!

这两天发生一个事&#xff0c;登月50周年之际&#xff0c;阿波罗11号制导计算机&#xff08;AGC&#xff09;指令和登月模块的源代码在Github上发布公开了&#xff0c;大量中国人前往围观&#xff0c;把issues区当成了论坛版块灌水留言。猎奇起哄本为消遣作乐&#xff0c;有不妥…

【BZOJ3252】攻略【闵可夫斯基和】【堆启发式合并】

题意&#xff1a;给一棵 nnn 个点的有点权的树&#xff0c;你需要找 kkk 条根到叶子的路径&#xff0c;使得路径并集的权值和最大。 n≤2105n\leq 2\times 10^5n≤2105 其实就是个贪心&#xff0c;只是从这个角度更自然一点( 先有个显然的 dp&#xff0c;设 f(u,k)f(u,k)f(u,…

Codeforces Round #592 (Div. 2) G. Running in Pairs 构造(水)

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 史上最水GGG题&#xff0c;没有之一。 考虑最小的情况如何构造&#xff0c;显然就是让a,ba,ba,b都1−n1-n1−n依次排列即可&#xff0c;这样的最小值为n∗(n1)2\frac{n*(n1)}{2}2n∗(n1)​…

Async,Await和ConfigureAwait的关系

在.NET Framework 4.5中&#xff0c;async / await关键字已添加到该版本中&#xff0c;简化多线程操作&#xff0c;以使异步编程更易于使用。为了最大化利用资源而不挂起UI&#xff0c;你应该尽可能地尝试使用异步编程。虽然async / await让异步编程更简单&#xff0c;但是有一…

【洛谷P5114】八月脸【边分治】【闵可夫斯基和】

题意&#xff1a;一棵 nnn 个点的树&#xff0c;每个点有两个权值 ai,bia_i,b_iai​,bi​&#xff0c;有黑白两种颜色。mmm 次询问&#xff0c;每次给定一个 kkk,求一条端点异色的路径&#xff0c;使得 k∑ai∑bik\sum a_i\sum b_ik∑ai​∑bi​ 最大化。 n≤2105n\leq 2\times…

Educational Codeforces Round 75 (Rated for Div. 2) E2. Voting (Hard Version) 贪心

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; n≤2e5,m≤n,p≤1e9n\le2e5,m\le n,p\le 1e9n≤2e5,m≤n,p≤1e9 思路&#xff1a; 首先需要发现一些性质&#xff0c;假设preipre_iprei​代表所有mj<im_j< imj​<i的pjp_jpj​和。可以发现&#…

【招聘(深圳)】敢为软件技术有限公司 .Net 工程师

敢为软件介绍&#xff1a;深圳市敢为软件技术有限公司致力于打造面向未来的最酷最具影响力的高科技公司。敢为软件在万物互联、行业应用、大屏可视化及自然交互等方面形成完整的闭环&#xff0c;是物联网领域连接能力最强、应用领域最广、用户体验最好的高科技公司。在万物互联…

【WC2018】通道【边分治】【虚树】【树的直径】

题意&#xff1a;给三棵基于同一点集的带边权的树&#xff0c;边权非负&#xff0c;求两点间三棵树上距离之和的最大值。 n≤105n\leq 10^5n≤105 一句话题解&#xff1a;在第一棵树上做边分治&#xff0c;丢到第二棵树上建虚树&#xff0c;在虚树上根据第三棵树的直径dp。 首…

DevOps案例研究:知人善任——Google敏捷核心文化

内容来源&#xff1a;DevOps案例深度研究-Google敏捷实践战队&#xff0c;本文只展示部分PPT及研究成果&#xff0c;更多细节请关注案例分享会。本文内容贡献者&#xff1a;陈霁、刘翀、谈佳婧、张霖。阅读干货前先感受一下热烈的氛围~一、Google如何快速交付原型1.1 Savioke公…