[SDOI2014] 重建 - 矩阵树定理,概率期望

#include 
#define eps 1e-6
using namespace std;

const int N = 55;

namespace mat {
double a[N][N];
int n,p=1;

double gauss_jordan() {
    double ans = 1;
    for(int i=1;i<=n;i++) {
        int r=i;
        for(int j=i+1;j<=n;j++) if(fabs(a[j][i])>fabs(a[r][i])) r=j;
        if(r-i) {for(int j=1;j<=n+1;j++) swap(a[i][j],a[r][j]); ans*=-1;}
        //if(fabs(a[i][i])>n;
    for(int i=1;i<=n;i++) {
        for(int j=1;j<=n;j++) {
            cin>>p[i][j];
            if(abs(1-p[i][j])