3270: 博物馆

Time Limit: 30 Sec  Memory Limit: 128 MB
Submit: 269  Solved: 147
[Submit][Status][Discuss]

Description

有一天Petya和他的朋友Vasya在进行他们众多旅行中的一次旅行,他们决定去参观一座城堡博物馆。这座博物馆有着特别的样式。它包含由m条走廊连接的n间房间,并且满足可以从任何一间房间到任何一间别的房间。
两个人在博物馆里逛了一会儿后两人决定分头行动,去看各自感兴趣的艺术品。他们约定在下午六点到一间房间会合。然而他们忘记了一件重要的事:他们并没有选好在哪儿碰面。等时间到六点,他们开始在博物馆里到处乱跑来找到对方(他们没法给对方打电话因为电话漫游费是很贵的)
不过,尽管他们到处乱跑,但他们还没有看完足够的艺术品,因此他们每个人采取如下的行动方法:每一分钟做决定往哪里走,有Pi 的概率在这分钟内不去其他地方(即呆在房间不动),有1-Pi 的概率他会在相邻的房间中等可能的选择一间并沿着走廊过去。这里的i指的是当期所在房间的序号。在古代建造是一件花费非常大的事,因此每条走廊会连接两个不同的房间,并且任意两个房间至多被一条走廊连接。
两个男孩同时行动。由于走廊很暗,两人不可能在走廊碰面,不过他们可以从走廊的两个方向通行。(此外,两个男孩可以同时地穿过同一条走廊却不会相遇)两个男孩按照上述方法行动直到他们碰面为止。更进一步地说,当两个人在某个时刻选择前往同一间房间,那么他们就会在那个房间相遇。
两个男孩现在分别处在a,b两个房间,求两人在每间房间相遇的概率。

Input

第一行包含四个整数,n表示房间的个数;m表示走廊的数目;a,b (1 ≤ a, b ≤ n),表示两个男孩的初始位置。
之后m行每行包含两个整数,表示走廊所连接的两个房间。
之后n行每行一个至多精确到小数点后四位的实数 表示待在每间房间的概率。
题目保证每个房间都可以由其他任何房间通过走廊走到。

Output

输出一行包含n个由空格分隔的数字,注意最后一个数字后也有空格,第i个数字代表两个人在第i间房间碰面的概率(输出保留6位小数)
注意最后一个数字后面也有一个空格

Sample Input

2 1 1 2
1 2
0.5
0.5

Sample Output

0.500000 0.500000

HINT

对于100%的数据有 n <= 20,n-1 <= m <= n(n-1)/2

Source

题解待更!

代码:

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<map>
#include<vector>
#define eps 1e-9
using namespace std;
int n,m,tot,aa,bb;
int d[405];
double a[405][405],p[405],f[405];
vector<int> e[405];
int id(int x,int y){return (x-1)*n+y;
}
void build(int x,int y)
{a[id(x,y)][id(x,y)]--;for (unsigned int i=0; i<e[x].size(); i++)for (unsigned int j=0; j<e[y].size(); j++){int xx=e[x][i],yy=e[y][j];if (xx!=yy){if (xx==x && yy==y) a[id(x,y)][id(xx,yy)]+=1.0*p[x]*p[y];else if (xx==x) a[id(x,y)][id(xx,yy)]+=1.0*p[xx]*(1-p[yy])/d[yy];else if (yy==y) a[id(x,y)][id(xx,yy)]+=1.0*(1-p[xx])*p[yy]/d[xx];else a[id(x,y)][id(xx,yy)]+=1.0*(1-p[xx])*(1-p[yy])/d[xx]/d[yy];}}
}
void gauss()
{double t; int to;for (int i=1; i<=tot; i++){t=0,to=0;for (int j=i; j<=tot; j++) if (fabs(a[j][i])>t) t=fabs(a[j][i]),to=j;if (to!=i) for (int j=1; j<=tot+1; j++) swap(a[to][j],a[i][j]);for (int j=i+1; j<=tot; j++) {if (fabs(a[j][i])<eps) continue;t=a[j][i]/a[i][i];for (int k=i; k<=tot+1; k++) a[j][k]-=t*a[i][k];}}for (int i=tot; i; i--){t=1.0*a[i][tot+1];for (int j=i+1; j<=tot; j++) t-=1.0*f[j]*a[i][j];f[i]=1.0*t/a[i][i];}
}
void init()
{scanf("%d%d%d%d",&n,&m,&aa,&bb);tot=n*n;a[id(aa,bb)][tot+1]=-1;for (int i=1; i<=n; i++) e[i].push_back(i);for (int i=1; i<=m; i++){int u,v; scanf("%d%d",&u,&v); d[u]++; d[v]++;e[u].push_back(v); e[v].push_back(u);}for (int i=1; i<=n; i++) scanf("%lf",&p[i]);for (int i=1; i<=n; i++) for (int j=1; j<=n; j++) build(i,j);
}
int main()
{init();gauss();for(int i=1;i<=n;i++){int t=id(i,i);printf("%.6lf",f[t]);if(i!=n)printf(" ");}return 0;
}

View Code

转载于:https://www.cnblogs.com/HQHQ/p/5793819.html

BZOJ3270: 博物馆相关推荐

  1. bzoj3270 博物馆

    Description 有一天Petya和他的朋友Vasya在进行他们众多旅行中的一次旅行,他们决定去参观一座城堡博物馆.这座博物馆有着特别的样式.它包含由m条走廊连接的n间房间,并且满足可以从任何一 ...

  2. 【BZOJ3270】【高斯消元】博物馆

    Description 有一天Petya和他的朋友Vasya在进行他们众多旅行中的一次旅行,他们决定去参观一座城堡博物馆.这座博物馆有着特别的样式.它包含由m条走廊连接的n间房间,并且满足可以从任何一 ...

  3. BZOJ3270 博物館 概率DP 高斯消元

    BZOJ3270 博物館 概率DP 高斯消元 @(XSY)[概率DP, 高斯消元] Description 有一天Petya和他的朋友Vasya在进行他们众多旅行中的一次旅行,他们决定去参观一座城堡博 ...

  4. 【权游相关】龙的简史 | 混乱博物馆

    世界上有讲不完的龙的故事,在欧洲,巨龙总是与骑士和公主联系起来.他们遥远地根植于古希腊和北欧的神话中,然而巨龙成为今天这幅双翼避日.口吐烈火的冷血形象,却是很晚近的事情. 本期混乱博物馆将带给你欧洲巨 ...

  5. 华盛顿多家博物馆因美政府停摆“闭门谢客”

    当地时间1月17日,位于华盛顿国家广场区域的多家博物馆,因美国政府停摆而持续处于关门状态. 游客在关门的博物馆门口停留. 博物馆门口摆放着关门的公告. (中新社记者 陈孟统 摄)

  6. 初中参观机器人博物馆的作文_2021考研英语二大作文模板整理

    陪伴你 一直到 考研结束 公众号又又又改版了,大家的订阅列表(也就是公众号的推文)不再按照时间顺序,而是按互动程度排序.希望大家可以多多互动(点个好看吧). 写在前面: ▲本文推荐的模板仅供参考,切勿 ...

  7. 谷歌等揭露「AI任务疑难」:存在局限的ImageNet等基准,就像无法代表「整个世界」的博物馆...

    来源:AI科技评论 作者:杏花 编辑:青暮 在日常生活中,我们需要一些「标准」来衡量个人的行为. 而在科研工作中,研究人员也需要一些「基准」来评估模型的性能. 因此,不管是普遍的「标准」还是特定的「基 ...

  8. 都说做ToB商业模式,VR AR博物馆内容的矿或许可以挖下

    技术和人文的结合,或许是一个可以挖掘的商业发展模式. [编者按]作者: 李亚倢 有一种习惯在国内可能不明显,如果去国外旅游,博物馆和标志性景观建筑总是要去一下的.因为在这些地方汇聚着整个城市的文化积淀 ...

  9. 现代浏览器博物馆_云旅游!Tableau 为你揭秘纽约现代艺术博物馆的珍贵馆藏

    互动话题见文末 "双节"假期刚刚结束,回归岗位的小伙伴们是否还对假期生活有点意犹未尽呢?近年来,大家出行已不仅是欣赏好风景以及逛吃逛吃了,许多爱好艺术的小伙伴也会踏进艺术馆,利用假 ...

最新文章

  1. 什么是Sprint?
  2. js html css淘宝足迹日历效果,js实现简单的日历显示效果函数示例
  3. springboot(七) 配置嵌入式Servlet容器
  4. win10只有edge浏览器能上网的解决方法
  5. java不进入for_为什么阿里巴巴Java开发手册中强制要求不要在foreach循环里进行元素的remove和add操作?...
  6. 动图处理_100+动图带你看懂百大加工工艺(成型、表面处理、链接、切割)
  7. SVN删除某版本解决办法
  8. 初始MySQL增删改查数据
  9. WPS快捷键之 EXCEL高级
  10. 最小二乘求解与正则化
  11. matlab中的uigetfile函数和setappdata函数
  12. 对企业形象识别系统(CIS)的了解
  13. 树莓派linux led字符设备驱动(互斥)
  14. 百度LBS开放平台个性化地图 制作一款独一无二的地图
  15. Excel如何批量添加超链接
  16. Eclipse WTP
  17. 初级调色转档CameraRaw
  18. Java小白 前端后端开发 环境搭建【jdk+idea+webstorm+maven+nodejs+vue+mysql】
  19. 织梦dede文章列表调用标签的用法和规则
  20. vue-3d-model:一个展示三维模型的 Vue 组件

热门文章

  1. [春节话题]春节联欢晚会的感受
  2. TypeError: the JSON object must be str, bytes or bytearray, not NoneType
  3. linux下面获取当前bing-国内版的壁纸
  4. 虚拟机没有域服务器,虚拟化现存域控制器
  5. 关联查询---Mybatis学习笔记(九)
  6. java基础知识陷阱
  7. 解决cc1plus.exe: out of memory allocating bytes
  8. python之从文件中按行读取数据
  9. python 决策树 math库 c45算法
  10. redis的常用命令及实例讲解