正題

題目鏈接:https://www.luogu.com.cn/problem/P3317

---

題目大意

$n$個(gè)點(diǎn)若干條邊。告訴你每條邊出現(xiàn)的概率，求剛好出現(xiàn)一顆生成樹的概率是多少。

---

解題思路

矩陣樹定理是計(jì)算每個(gè)生成樹的每條邊乘積之和。

我們考慮將答案轉(zhuǎn)換為那個(gè)形式，$a_{i,j}$表示$i?>j$的邊出現(xiàn)的概率，那么對(duì)于每棵生成樹有答案
$$\prod _{x?>y}{a}_{i,j}?\prod _{x?/>y}(1?a_{i,j})$$
也就是
$$\prod _{i=1}^n\prod _{j=1}^n(1?a_{i,j})?\prod _{x?>y}\frac{a_{i,j}}{1?a_{i,j}}$$
所以我們先計(jì)算出所有的$1?a_{i,j}$的乘積，然后再將邊權(quán)變?yōu)?span id="ze8trgl8bvbq" class="katex--inline">$\frac{a_{i,j}}{1?a_{i,j}}$用矩陣樹計(jì)算即可。

時(shí)間復(fù)雜度$O(n^3)$

---

$code$

#include<cstdio> #include<cstring> #include<algorithm> #include<cmath> using namespace std; const int N=51; const double eps=1e-8; int n; double a[N][N],ans; void Gauss(){for(int i=1;i<n;i++){int mx=i;for(int j=i+1;j<n;j++)if(fabs(a[j][i])>fabs(a[mx][i]))mx=j;if(mx!=i)for(int j=1;j<n;j++)swap(a[i][j],a[mx][j]);for(int j=i+1;j<n;j++){double mul=a[j][i]/a[i][i];for(int k=i;k<n;k++)a[j][k]-=a[i][k]*mul;}if(fabs(a[i][i])<eps){ans=0;return;}}for(int i=1;i<n;i++)ans=ans*a[i][i];ans=fabs(ans); } int main() {scanf("%d",&n);ans=1;for(int i=1;i<=n;i++)for(int j=1;j<=n;j++){scanf("%lf",&a[i][j]);if(a[i][j]<eps) a[i][j]=eps;if(1-a[i][j]<eps) a[i][j]=1-eps;if(i<j) ans*=1-a[i][j];a[i][j]=a[i][j]/(1-a[i][j]);}for(int i=1;i<=n;i++){a[i][i]=0;for(int j=1;j<=n;j++)if(i!=j)a[i][i]-=a[i][j];}Gauss();printf("%.10lf",ans); }

總結(jié)

以上是生活随笔為你收集整理的P3317-[SDOI2014]重建【矩阵树定理,数学期望】的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。