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

来源：互联网发布：安卓放弃java 编辑：程序博客网时间：2024/05/20 16:14

看了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;}

阅读全文

0 0