【BZOJ3534】重建（矩阵树定理）

最新推荐文章于 2021-09-09 22:50:27 发布

小蒟蒻yyb

最新推荐文章于 2021-09-09 22:50:27 发布

阅读量230

点赞数 1

CC 4.0 BY-SA版权

分类专栏： BZOJ 各省省选

本文链接：https://blog.youkuaiyun.com/qq_30974369/article/details/81489150

BZOJ 同时被 2 个专栏收录

382 篇文章

订阅专栏

各省省选

232 篇文章

订阅专栏

本文介绍了一个使用矩阵树定理解决图论中特定问题的例子。通过C++代码实现了一个计算过程，该过程利用高斯消元法求解特定矩阵的行列式，以此来计算与图相关的一个特殊值。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题面

BZOJ
洛谷

题解

这。。。。
矩阵树定理神仙用法？？？？

#include<iostream>
#include<cmath>
using namespace std;
#define ll long long
#define eps 1e-8
int n;
double g[55][55],s=1;
double Gauss()
{
    double ret=1;
    for(int i=1;i<n;++i)
    {
        int nw=i;
        for(int j=i+1;j<n;++j)if(fabs(g[j][i])>fabs(g[nw][i]))nw=j;
        swap(g[nw],g[i]);
        for(int j=i+1;j<n;++j)
        {
            double t=g[j][i]/g[i][i];
            for(int k=i;k<n;++k)g[j][k]-=t*g[i][k];
        }
        ret*=g[i][i];
    }
    return fabs(ret);
}
int main()
{
    ios::sync_with_stdio(false);
    cin>>n;
    for(int i=1;i<=n;++i)
        for(int j=1;j<=n;++j)
        {
            cin>>g[i][j];
            if(g[i][j]<eps)g[i][j]=eps;
            if(g[i][j]>1-eps)g[i][j]=1-eps;
            if(i<j)s*=1-g[i][j];
            g[i][j]/=1-g[i][j];
        }
    for(int i=1;i<=n;++i)
        for(int j=1;j<=n;++j)
            if(i!=j)g[i][i]-=g[i][j];
    printf("%.8lf\n",Gauss()*s);
}