【数论】【矩阵树定理】BZOJ3534重建

分析：

首先，很显然我们要求的是：
$$\sum _{e\in T}{P}_e\ast \sum _{e^{\prime }\notin T}(1-P_{e^{\prime }})$$

前者可以用矩阵树定理算出来，问题就是转化后半部分。

$$\sum _{e^{\prime }\notin T}(1-P_{e^{\prime }})=\frac{{\sum }_e(1-P_e)}{{\sum }_{e\in T}(1-P_e)}$$

换言之，可以将边权设置为$\frac{P_e}{(1-P_e)}$，然后用矩阵树暴力求，最后再乘上${\sum }_e(1-P_e)$即可

注意，矩阵树定理中，$A_{i,i}=-{\sum }_j{A}_j$

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define SF scanf
#define PF printf
#define MAXN 60
#define EPS 1e-6
using namespace std;
int n;
double p[MAXN][MAXN],P=1;
double gauss(int row,int col){
	int r=1,c=1;
	for(;r<=row&&c<=col;r++,c++){
		int maxr=r;
		for(int i=r;i<=row;i++)
			if(p[i][c]>p[maxr][c])
				maxr=i;
		if(p[maxr][c]<EPS)
			return 0;
		swap(p[maxr],p[r]);
		for(int i=r+1;i<=row;i++){
			double k=p[i][c]/p[r][c];
			for(int j=c;j<=col;j++)
				p[i][j]=p[i][j]-p[r][j]*k;
		}
	}
	double res=1;
	for(int i=1;i<=row;i++)
		res*=p[i][i];
	return res;
}
int main(){
	SF("%d",&n);	
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++){
			SF("%lf",&p[i][j]);
			p[i][j]=min(p[i][j],1-EPS);
			p[i][j]=max(p[i][j],EPS);
			if(i>j)
				P*=(1.0-p[i][j]);
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			if(i!=j)
				p[i][j]=-p[i][j]/(1.0-p[i][j]);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			if(i!=j)
				p[i][i]-=p[i][j];
//	for(int i=1;i<=n;i++){
//		for(int j=1;j<=n;j++)
//			PF("%.2lf ",p[i][j]);
//		PF("\n");
//	}
	double res=gauss(n-1,n-1);
//	PF("[%lf]",res);
	res=res*P;
	PF("%lf",res);
}