[SDOI2014]重建

题意

每条边都有一定概率$p\in[0,1]$出现在图中,求生成一棵树的概率是多少

---

题解

可以想一个简单的暴力枚举所有树$T$,然后

$$Ans=\sum_T\prod_{e\in T}p_e\prod_{e\notin T}(1-p_e)$$

题目中要问的生成树的概率,既然是生成树那么想一想应该要用到矩阵树定理来求

根据基尔霍夫矩阵的定义我们可以知道他任意一个代数余子式是所有生成树的边权积的和

说人话就是求的

$$\sum_T\prod_{e\in T}w_e$$

和上面要求的长得好像,但是后面那一坨怎么办?

考虑到

$$\prod_{e\notin T}(1-p_e)=\frac{\prod_e(1-p_e)}{\prod_{e\in T}(1-p_e)}$$

把这个带回去就得到了

$$Ans=\prod_e(1-p_e)\sum_T\prod_{e\in T}\frac{p_e}{1-p_e}$$

于是我们就可以令$w_e=\frac{p_e}{1-p_e}$然后用高斯消元求一遍行列式再乘上前面那一坨就好了

要是$p_e=1$怎么办$?$

考虑到$\frac1{1-p_e}\to\infty$又有$\frac1\epsilon\approx\infty$所以近似令$p_e=1-\epsilon$就好了

#include<bits/stdc++.h>
#define fp(i,a,b) for(register int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(register int i=a,I=b-1;i>I;--i)
#define go(u) for(register int i=fi[u],v=e[i].to;i;v=e[i=e[i].nx].to)
#define file(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}
using namespace std;
const int N=55;
const double eps=1e-8;
typedef double d;
int n;d prod=1,G[N][N];
inline int cmp(double x){return fabs(x)<eps?0:(x<0?-1:1);}
inline d det(int n){
    int mx;d t,tp=1;
    fp(i,1,n){mx=i;
        fp(j,i+1,n)if(cmp(G[mx][i]-G[j][i])<0)mx=i;
        if(mx^i)fp(j,i,n)swap(G[i][j],G[mx][j]);
        if(!G[i][i])return 0;
        fp(j,i+1,n){
            t=G[j][i]/G[i][i];
            fp(k,i,n)G[j][k]-=G[i][k]*t;
        }tp*=G[i][i];
    }
    return fabs(tp);
}
int main(){
    #ifndef ONLINE_JUDGE
        file("s");
    #endif
    scanf("%d",&n);
    fp(i,1,n)fp(j,1,n){
        scanf("%lf",G[i]+j);
        if(i==j)continue;
        if(G[i][j]>1-eps)G[i][j]-=eps;
        if(i<j)prod*=(1-G[i][j]);
        G[i][j]/=(1-G[i][j]);
    }
    fp(i,1,n)fp(j,1,n)if(i^j)
        G[i][i]+=G[i][j],G[i][j]=-G[i][j];
    printf("%.8lf\n",det(n-1)*prod);
return 0;
}