bzoj3534 [Sdoi2014]重建

变形的$Martix-Tree$定理

发现我们要求的是$\prod_{i \in E}{p_{i}} * \prod_{i \notin E}{(1-p_{i})}$

然后呢? 矩阵树对重边也有效对吧。考虑带权图,发现建出来的矩阵的任何一个$n-1$阶主子式的行列式的值都是其所有生成树的边权之积的和

那么就可以搞了,考虑每一条边权为$\frac{p_{i}}{(1-p_{i})}$,最后再乘一个$\prod{(1-p_{i})}$即可

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<iostream>
 4 #include<algorithm>
 5 #include<cmath>
 6 #define eps 1e-10
 7 #define N 55
 8 using namespace std;
 9 int n;
10 double g[N][N],a[N][N],ans=1,tmp=1;
11 void gauss(){
12     n--;
13     for(int k=1;k<=n;k++){
14         int bst=k;
15         for(int i=k+1;i<=n;i++)
16             if(fabs(a[i][k])>fabs(a[bst][k]))bst=i;
17         if(bst!=k){
18             for(int i=k;i<=n;i++)
19                 swap(a[k][i],a[bst][i]);
20             ans*=-1;
21         }
22         for(int i=k+1;i<=n;i++){
23             double t=a[i][k]/a[k][k];
24             for(int j=k;j<=n;j++)
25                 a[i][j]-=t*a[k][j];
26         }
27         ans*=a[k][k];
28     }
29 }
30 int main(){
31     scanf("%d",&n);
32     for(int i=1;i<=n;i++)
33         for(int j=1;j<=n;j++)
34             scanf("%lf",&g[i][j]);
35     for(int i=1;i<=n;i++){
36         for(int j=i+1;j<=n;j++){
37             if(g[i][j]==1)g[i][j]-=eps;
38             tmp*=(1-g[i][j]);
39             g[i][j]=g[i][j]/(1-g[i][j]);
40             a[i][i]+=g[i][j];a[j][j]+=g[i][j];
41             a[i][j]=-g[i][j];a[j][i]=-g[i][j];
42         }
43     }
44     gauss();
45     printf("%0.10lf\n",ans*tmp);
46     return 0;
47 }

原文地址:https://www.cnblogs.com/Ren-Ivan/p/8185639.html

时间: 2024-11-05 21:46:42

bzoj3534 [Sdoi2014]重建的相关文章

[SDOI2014]重建

P3317 [SDOI2014]重建 题解 直接M-T肯定不对 推出的是对于所有树的生成概率的和,可以考虑行列式的期望,再交换求和号即可 同时乘上π(1-P)再变化初始的概率就有点厉害了 一种变化的技巧 代码: #include<bits/stdc++.h> #define il inline #define reg register int #define numb (ch^'0') #define ld long double using namespace std; typedef lo

【SDOI2014】【BZOJ3534】重建

Description T国有N个城市,用若干双向道路连接.一对城市之间至多存在一条道路. 在一次洪水之后,一些道路受损无法通行.虽然已经有人开始调查道路的损毁情况,但直到现在几乎没有消息传回. 辛运的是,此前T国政府调查过每条道路的强度,现在他们希望只利用这些信息估计灾情.具体地,给定每条道路在洪水后仍能通行的概率,请计算仍能通行的道路恰有N-1条,且能联通所有城市的概率. Input 输入的第一行包含整数N. 接下来N行,每行N个实数,第i+l行,列的数G[i][j]表示城市i与j之 间仍有

BZOJ 3534 [Sdoi2014]重建

题解:矩阵树定理 邻接矩阵-度数矩阵(期望下) 求出来的行列式为所有(生成树边权乘积)的和 每条边边权化为 c/(1-c),最后乘上π(1-c),对1边权特殊处理一下 问题:矩阵树定理不熟,不会证明 #include<iostream> #include<cstdio> #include<cstring> #include<cmath> using namespace std; const int maxn=60; const double eps=1e-9

luogu3317 [SDOI2014]重建

原来矩阵树定理对于边是概率的情况也是适用的qwqwq. ref #include <iostream> #include <cstdio> #include <cmath> using namespace std; int n; double w[55][55], uu; const double eps=1e-7; void gauss(){ for(int i=1; i<=n; i++){ int maxi=i; for(int j=i+1; j<=n;

luoguP3317 [SDOI2014]重建 变元矩阵树定理 + 概率

首先,我们需要求的是 $$\sum\limits_{Tree} \prod\limits_{E \in Tree} E(u, v) \prod\limits_{E \notin Tree} (1 - E(u, v))$$ 我们知道变元矩阵树定理 ---> 不知道请见此 我们自然希望要求和的事物只跟生成树的边有关 因此考虑把$\prod\limits_{E \notin Tree} (1 - E(u, v))$转化为$\prod\limits_{E} (1 - E(u, v)) * \frac{1

P3317 [SDOI2014]重建

\(\color{#0066ff}{ 题目描述 }\) T国有N个城市,用若干双向道路连接.一对城市之间至多存在一条道路. 在一次洪水之后,一些道路受损无法通行.虽然已经有人开始调查道路的损毁情况,但直到现在几乎没有消息传回. 幸运的是,此前T国政府调查过每条道路的强度,现在他们希望只利用这些信息估计灾情.具体地,给定每条道路在洪水后仍能通行的概率,请计算仍能通行的道路恰有N-1条,且能联通所有城市的概率. \(\color{#0066ff}{输入格式}\) 输入的第一行包含整数N. 接下来N行

P3317 [SDOI2014]重建 变元矩阵树定理 高斯消元

传送门:https://www.luogu.org/problemnew/show/P3317 这道题的推导公式还是比较好理解的,但是由于这个矩阵是小数的,要注意高斯消元方法的使用: #include <algorithm> #include <iterator> #include <iostream> #include <cstring> #include <cstdlib> #include <iomanip> #include

Matrix-Tree定理题表(已完成(ojbk))

矩阵树这个东西……并不懂什么基尔霍夫矩阵……背了一下结论……(顺便用这个东西加强了一下矩阵)(打板子的时候还是该取负取负,因为不取负才有可能是负数,最后答案一定是正数???(ryf说一定是这样))bzoj3534:[Sdoi2014]重建 矩阵树定理的一个小概念+板子(实数高斯消元)+处理精度(把0搞成一个很小的小数)bzoj4596:[Shoi2016]黑暗前的幻想乡 容斥一下+板子(求逆元)bzoj4031:[HEOI2015]小Z的房间 板子(辗转相除)bzoj1002:[FJOI2007

康复计划#5 Matrix-Tree定理(生成树计数)的另类证明和简单拓展

本篇口胡写给我自己这样的什么都乱证一通的口胡选手 以及那些刚学Matrix-Tree,大致理解了常见的证明但还想看看有什么简单拓展的人- 大概讲一下我自己对Matrix-Tree定理的一些理解.常见版本的证明.我自己的证明,以及简单的一些应用(比如推广到有向图.推广到生成树边权的乘积和什么的,非常基础). 应该看到这里的人都知道Matrix-Tree定理是干什么的吧-就是统计一个无向图的生成树个数,表示成一个行列式. 1.前置定义及性质 首先是Matrix-Tree定理相关的定义:对于一个无向图