BZOJ3534:[SDOI2014]重建(矩阵树定理)
Description
T国有N个城市,用若干双向道路连接。一对城市之间至多存在一条道路。
在一次洪水之后,一些道路受损无法通行。虽然已经有人开始调查道路的损毁情况,但直到现在几乎没有消息传回。
幸运的是,此前T国政府调查过每条道路的强度,现在他们希望只利用这些信息估计灾情。具体地,给定每条道路在洪水后仍能通行的概率,请计算仍能通行的道路恰有N-1条,且能联通所有城市的概率。
Input
输入的第一行包含整数N。
接下来N行,每行N个实数,第i+l行,列的数G[i][j]表示城市i与j之
间仍有道路联通的概率。
输入保证G[i][j]=G[j][i],且G[i][j]=0;G[i][j]至多包含两位小数。
Output
输出一个任意位数的实数表示答案。
你的答案与标准答案相对误差不超过10^(-4)即视为正确。
Sample Input
0 0.5 0.5
0.5 0 0.5
0.5 0.5 0
Sample Output
HINT
1 < N < =50
数据保证答案非零时,答案不小于10^-4
Solution
题目即让求:($T$是生成树,$e$是边)
$\sum_{T}\prod_{e\in T} p_e \prod_{e\notin T}(1-p_e)$
把第二个$\prod$变一下
$\sum_{T}\prod_{e\in T} p_e \frac{\prod_{e}(1-p_e)}{\prod_{e\in T}(1-p_e)}$
也就是
$\prod_{e}(1-p_e)\sum_{T}\prod_{e\in T} \frac{p_e}{1-p_e}$。
其中行和矩阵$-$边权矩阵的行列式的值$=\sum_{T}\prod_{e\in T} w_e$,其中$w$是边权。
然后就可以高斯消元求解了。
Code
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<cmath> 5 #define N (59) 6 #define eps (1e-10) 7 using namespace std; 8 9 int n; 10 double a[N][N],f[N][N],ans=1; 11 12 void Gauss() 13 { 14 int w=1; 15 for (int i=1; i<=n-1; ++i) 16 { 17 int num=i; 18 for (int j=i+1; j<=n-1; ++j) 19 if (fabs(f[j][i])>fabs(f[num][i])) num=j; 20 if (num!=i) swap(f[num],f[i]), w=-w; 21 for (int j=i+1; j<=n-1; ++j) 22 { 23 double t=f[j][i]/f[i][i]; 24 for (int k=i; k<=n-1; ++k) 25 f[j][k]-=t*f[i][k]; 26 } 27 } 28 for (int i=1; i<=n-1; ++i) ans*=f[i][i]; 29 for (int i=1; i<=n; ++i) 30 for (int j=i+1; j<=n; ++j) 31 ans*=1-a[i][j]; 32 printf("%.10lf\n",ans*w); 33 } 34 35 int main() 36 { 37 scanf("%d",&n); 38 for (int i=1; i<=n; ++i) 39 for (int j=1; j<=n; ++j) 40 { 41 scanf("%lf",&a[i][j]); 42 if (a[i][j]==1) a[i][j]-=eps; 43 if (i==j) continue; 44 f[i][j]=-a[i][j]/(1-a[i][j]); 45 f[i][i]-=f[i][j]; 46 } 47 Gauss(); 48 }
转载于:https://www.cnblogs.com/refun/p/10375388.html
BZOJ3534:[SDOI2014]重建(矩阵树定理)相关推荐
- P3317-[SDOI2014]重建【矩阵树定理,数学期望】
正题 题目链接:https://www.luogu.com.cn/problem/P3317 题目大意 nnn个点若干条边.告诉你每条边出现的概率,求刚好出现一颗生成树的概率是多少. 解题思路 矩阵树 ...
- 【学习笔记】矩阵树定理(Matrix-Tree)
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 目录 一.矩阵树定理 二.常用定理 三.例题 1. Luogu P6178 [模板]Matrix-Tr ...
- 最小生成树、矩阵树定理、Prufer序列总结
Kruskal算法 按边权排序,从小到大合并不在同一集合两点即可 Prim算法 每次加入一个到当前已选点集最近的点 P2619 [国家集训队]Tree I 考虑二分,每次给白边加上一个mid,通过这种 ...
- 矩阵树定理2020HDU多校第6场j-Expectation[位运算+期望]
矩阵树定理 用于求解图上面生成树的个数,生成树的个数等于基尔霍夫矩阵的任何一个N-1阶主子式的行列式的绝对值 矩阵树模板 struct Matrix_Tree {ll a[N][N];Matrix_T ...
- Luogu P4336 [SHOI2016]黑暗前的幻想乡(容斥,矩阵树定理,子集反演)
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 Luogu P4336 [SHOI2016]黑暗前的幻想乡(容斥,矩阵树定理) Problem n≤1 ...
- 图论数学:矩阵树定理
运用矩阵树定理进行生成树计数 给定一个n个点m条边的无向图,问生成树有多少种可能 直接套用矩阵树定理计算即可 矩阵树定理的描述如下: 首先读入无向图的邻接矩阵,u-v G[u][v]++ G[v][u ...
- [CF917D]Stranger Trees[矩阵树定理+解线性方程组]
题意 给你 \(n\) 个点的无向完全图,指定一棵树 \(S\),问有多少棵生成树和这棵树的公共边数量为 \(k\in[0,n-1]\) \(n\leq 100\) 分析 考虑矩阵树定理,把对应的树边 ...
- HDU多校6 - 6836 Expectation(矩阵树定理+高斯消元求行列式)
题目链接:点击查看 题目大意:给出一张由 n 个点和 m 条边组成的无向图,对于任意一个生成树,其权值为 n - 1 条边的边权进行二进制的 and 运算,现在需要在图中任意选择一个生成树,问期望权值 ...
- Wannafly挑战赛23F-计数【原根,矩阵树定理,拉格朗日插值】
正题 题目链接:https://ac.nowcoder.com/acm/contest/161/F 题目大意 给出nnn个点的一张图,求它的所有生成树中权值和为kkk的倍数的个数.输出答案对ppp取模 ...
最新文章
- ICRA2021| 自动驾驶相关论文汇总
- 【Kaldi安装问题】Nvidia库cub安下载失败
- OBJECT_ID 与objectproperty
- 3.IDA-数据显示窗口(导出窗口、导入窗口、String窗口、...窗口)
- linux安装grub界面,系统安装: GRUB安装,配置及使用汇总
- 比较OpenBLAS,Intel MKL和Eigen的矩阵相乘性能
- 接口的抽象与实现(概述)
- IFIX 需要权限打开某个画面
- python有道字典_Python基于有道实现英汉字典功能
- 数据结构:B树和B+树的插入、删除图文详解
- 被单位开除,已经交了14年的养老保险,该怎么办?
- linux进程名称最大长度,linux – 进程名称长度的最大允许限制是多少?
- [Bzoj1009][HNOI2008]GT考试(动态规划)
- 大学生游戏静态HTML网页作业--美丽中国
- NOJ——[1480] 懒惰的风纪委Elaine
- ESP8266 AP服务器设置
- 每日一个小技巧:1分钟告诉你截图翻译软件哪个好用
- 计算机如果没有什么 就无法启动,电脑开机没有任何反应
- Qt编写4K/8K大分辨率播放器(8K占用1%CPU)
- 静态路由和默认路由的原理与基础配置命令