BZOJ3534: [Sdoi2014]重建（变元矩阵树）

看了xehoth大佬的博客,%%%.
https://blog.xehoth.cc/MatrixTree/

#include<bits/stdc++.h>
using namespace std;
inline int gcd(int x,int y){return (y?(gcd(y,x%y)):x);}
double H[150][150],*matrix[150],f=1;
const double eps=-1e-9;
int n;
inline double Det()
{
    for(int i=1;i<=n;i++)
    {
        int l=i;
        for(int j=i+1;j<=n;j++)if(fabs(matrix[j][i])>fabs(matrix[l][i]))l=j;
        if(l!=i)swap(matrix[i],matrix[l]);
        if(fabs(matrix[i][i])<eps)return 0;
        for(int j=i+1;j<=n;j++)
        {
            double t=matrix[j][i]/matrix[i][i];
            for(int k=i;k<=n;k++)matrix[j][k]-=t*matrix[i][k];
        }
    }
    double x=1;
    for(int i=1;i<=n;i++)x*=matrix[i][i];
    return x;
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(NULL);cout.tie(NULL);
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        matrix[i]=H[i];
        for(int j=1;j<=n;j++)
        {
            double x;cin>>x;
            if(j>=i)continue;
            f*=(1.0-x);
            double t=x/(1.0-x);
            H[i][i]+=t;H[j][j]+=t;
            H[i][j]-=t;H[j][i]-=t;
        }
    }
    n--;
    cout<<fixed<<setprecision(13)<<fabs(Det()*f)<<endl;
}