【CF113D】Museum
Solution
额题意的话大概就是给一个无向图然后两个人给两个出发点,每个点每分钟有\(p[i]\)的概率停留,问这两个人在每个点相遇的概率是多少
如果说我们知道最后在哪里相遇,处理起来会比较方便。注意到\(n\)比较小,所以考虑枚举最后相遇的房间,然后可以看成求两个人同时从最后相遇的房间出发,走到\(X\)和\(Y\)的概率,这样初始状态什么的比较好表示,就会好求很多了(跟。。hnoi2013游走有点像?)
考虑\(dp\) ,假设我们现在枚举到终点是\(t\),\(f[i][j]\)表示第一个人在\(i\)这个点,第二个人在\(j\)这个点的概率是多少
那么有初始状态:
1.\(i==j==t\)的情况,\(f[i][j]=1\)
2.\(i==j!=t\)的情况,\(f[i][j]=0\)
考虑从第\(i\)个房间走向一个相邻房间的概率\(k[i]\),记第\(i\)个点的度数为\(du[i]\),那么有
\]
(有\(1-p[i]\)的概率离开这个房间,并且走向每个相邻房间的概率是一样的)
所以我们可以得到转移:
f[i][j]&=p[i]*p[j]*f[i][j]&(都停留)\\
&+k[i]*p[j]*\sum f[u][j]&(第一个人走一步)\\
&+k[j]*p[i]*\sum f[i][v]&(第二个人走一步)\\
&+k[i]*k[j]*\sum f[u][v]&(两个人都走一步)
\end{aligned}
\]
其中\(u\)满足\(i\)和\(u\)在原图中有一条边直接相连,\(v\)满足\(v\)和\(j\)在原图中有一条边直接相连
然后移下项,对于每一对\((i,j)\)我们都可以得到一条长得像这样的式子:
\]
将\(f[i][j]\)看成未知数,就可以直接用高斯消元来做了
注意,移项的时候,因为在初始化中\(f[i][i]\)的值已经定下来了,所以应该当成常数来看
代码大概长这个样子
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#define db double
using namespace std;
const int N=25;
int mp[N][N],du[N];
db p[N],a[N*N][N*N],k[N],f[N],num[N][N];
int n,m,X,Y,mx;
db ans[N];
void fillmatrix(int t);
db gauss(int n);
void solve(int t);
void Fill(int x,int y);
int Id(int x,int y){return num[x][y];}
int main(){
#ifndef ONLINE_JUDGE
freopen("a.in","r",stdin);
#endif
int x,y;
scanf("%d%d%d%d",&n,&m,&X,&Y);
for (int i=1;i<=m;++i){
scanf("%d%d",&x,&y);
mp[x][y]=mp[y][x]=1;
++du[x]; ++du[y];
}
if (X==Y){
for (int i=1;i<=n;++i)
if (i==X) printf("%.10lf\n",1.0);
else printf("%.10lf\n",0.0);
return 0;
}
for (int i=1;i<=n;++i) scanf("%lf",p+i);
for (int i=1;i<=n;++i) k[i]=(1.0-p[i])/(1.0*du[i]);
mx=0;
for (int i=1;i<=n;++i)
for (int j=1;j<=n;++j)
if (i!=j) num[i][j]=++mx;
for (int i=1;i<=n;++i)
solve(i);
for (int i=1;i<=n;++i)
printf("%.10lf ",ans[i]);
printf("\n");
}
void solve(int t){
fillmatrix(t);
ans[t]=gauss(mx);
}
void fillmatrix(int t){
int id;
memset(a,0,sizeof(a));
memset(f,0,sizeof(f));
f[t]=1;
for (int i=1;i<=n;++i){
for (int j=1;j<=n;++j){
if (i!=j) Fill(i,j);
}
}
}
void Fill(int x,int y){
int id=Id(x,y);
a[id][id]=p[x]*p[y]-1.0;
for (int i=1;i<=n;++i){
if (mp[i][x]){
if (i!=y)
a[id][Id(i,y)]+=k[x]*p[y];
else
a[id][mx+1]-=k[x]*p[y]*f[y];
}
if (mp[i][y]){
if (i!=x)
a[id][Id(x,i)]+=k[y]*p[x];
else
a[id][mx+1]-=k[y]*p[x]*f[x];
}
for (int j=1;j<=n;++j){
if (mp[i][x]&&mp[j][y]){
if (i!=j)
a[id][Id(i,j)]+=k[x]*k[y];
else
a[id][mx+1]-=k[x]*k[y]*f[j];
}
}
}
}
db gauss(int n){
int id;
db tmp;
for (int i=1;i<=n;++i){
id=i;
for (int j=i+1;j<=n;++j)
if (fabs(a[j][i])>fabs(a[id][i])) id=j;
if (id!=i)
for (int j=1;j<=n+1;++j) swap(a[i][j],a[id][j]);
for (int j=i+1;j<=n;++j){
tmp=a[j][i]/a[i][i];
for (int k=i;k<=n+1;++k)
a[j][k]-=tmp*a[i][k];
}
}
for (int i=n;i>=1;--i){
for (int j=n;j>i;--j)
a[i][n+1]-=a[i][j]*a[j][n+1];
a[i][n+1]/=a[i][i];
}
return a[Id(X,Y)][mx+1];
}
【CF113D】Museum的更多相关文章
- 【254】◀▶IEW-Unit19
Unit 19 Technology Communication I.名词性从句在雅思写作中的运用 英语中哪些位置可以放名词? 1)主语 2)宾语 3)表语 4)同位语 名词的位置放一个句子=名词性从 ...
- Python高手之路【六】python基础之字符串格式化
Python的字符串格式化有两种方式: 百分号方式.format方式 百分号的方式相对来说比较老,而format方式则是比较先进的方式,企图替换古老的方式,目前两者并存.[PEP-3101] This ...
- 【原】谈谈对Objective-C中代理模式的误解
[原]谈谈对Objective-C中代理模式的误解 本文转载请注明出处 —— polobymulberry-博客园 1. 前言 这篇文章主要是对代理模式和委托模式进行了对比,个人认为Objective ...
- 【原】FMDB源码阅读(三)
[原]FMDB源码阅读(三) 本文转载请注明出处 —— polobymulberry-博客园 1. 前言 FMDB比较优秀的地方就在于对多线程的处理.所以这一篇主要是研究FMDB的多线程处理的实现.而 ...
- 【原】Android热更新开源项目Tinker源码解析系列之一:Dex热更新
[原]Android热更新开源项目Tinker源码解析系列之一:Dex热更新 Tinker是微信的第一个开源项目,主要用于安卓应用bug的热修复和功能的迭代. Tinker github地址:http ...
- 【调侃】IOC前世今生
前些天,参与了公司内部小组的一次技术交流,主要是针对<IOC与AOP>,本着学而时习之的态度及积极分享的精神,我就结合一个小故事来初浅地剖析一下我眼中的“IOC前世今生”,以方便初学者能更 ...
- Python高手之路【三】python基础之函数
基本数据类型补充: set 是一个无序且不重复的元素集合 class set(object): """ set() -> new empty set object ...
- Python高手之路【一】初识python
Python简介 1:Python的创始人 Python (英国发音:/ˈpaɪθən/ 美国发音:/ˈpaɪθɑːn/), 是一种解释型.面向对象.动态数据类型的高级程序设计语言,由荷兰人Guido ...
- 【开源】简单4步搞定QQ登录,无需什么代码功底【无语言界限】
说17号发超简单的教程就17号,qq核审通过后就封装了这个,现在放出来~~ 这个是我封装的一个开源项目:https://github.com/dunitian/LoTQQLogin ————————— ...
随机推荐
- 【转】自动化测试 - Appium + Python史上最全最简环境搭建步骤
一,为什么是Appium借一张图: 1.1 Appium优点 l 开源 l 跨架构:NativeApp.Hybird App.Web App l 跨设备:Android.iOS.Firefox ...
- Linux 安装FastDFS<单机版>(使用Mac远程访问)
阅读本文需要先阅读安装FastDFS<准备> 一 编译环境 yum install gcc-c++ yum -y install libevent yum install -y pcre ...
- NO.08--VUE之自定义组件添加原生事件
前几篇给大家分享了我的业余的“薅羊毛”的经历,回归正题,讲回vue吧: 许多vue新手在工作开发中会遇到一个问题,直接使用 button 添加原生事件是没有问题的,但是使用自定义组件添加原生事件时,就 ...
- Update类型_JDBC的方法_JAVA方法_Loadrunner脚本
java vuser JDBC 参数化的方法 如果不进行参数化 直接把32 33行去掉 ,sql 值写到valuers 中就行了 下面这是 insert,delete,update 三种方法 ...
- python数据可视化——matplotlib 用户手册入门:使用指南
参考matplotlib官方指南: https://matplotlib.org/tutorials/introductory/usage.html#sphx-glr-tutorials-introd ...
- 树形dp入门两题
题目描述 小明对数学饱有兴趣,并且是个勤奋好学的学生,总是在课后留在教室向老师请教一些问题.一天他早晨骑车去上课,路上见到一个老伯正在修剪花花草草,顿时想到了一个有关修剪花卉的问题.于是当日课后,小明 ...
- hbase实战——(1.1 nosql介绍)
什么是nosql NoSQL(NoSQL = Not Only SQL),意思是不仅仅是SQL的扩展,一般指的是非关系型的数据库. 随着互联网web2.0网站的兴起,传统的关系数据库在应付web2.0 ...
- USACO 2.3.3 Zero Sum 和为零(深搜枚举)
Description 请考虑一个由1到N(N=3, 4, 5 ... 9)的数字组成的递增数列:1 2 3 ... N. 现在请在数列中插入“+”表示加,或者“-”表示减,抑或是“ ”表示空白,来将 ...
- c# webBrowser清除缓存问题
1.webBrowser的浏览器为ie. 2.通过清除ie缓存即可. 3.代码调用如下: public enum ShowCommands : int { SW_HIDE = , SW_SHOWNOR ...
- c# 加载图片 正在被占用问题
问题情境:图片文件加载到pdf中,程序没有退出,再次加载该图片文件,提示被占用. 解决办法: 1.加载文件会锁定该文件,fromfile方法会导致占用内存较大,不使用该方法. FileStream f ...