矩阵树定理

最新推荐文章于 2022-12-10 15:55:10 发布

limboman

最新推荐文章于 2022-12-10 15:55:10 发布

阅读量435

点赞数

文章标签：矩阵树

本文链接：https://blog.csdn.net/limboman/article/details/54645652

版权

又看了一遍。。。。。
在周冬《生成树的计数及其应用》中有证明；
核心是利用Binet-Cauchy公式
但：
“如果图是一颗树，那么C的任一个n-1阶主子式的行列式均为1”
的结论
C好像是度数矩阵+邻接矩阵也没问题？

试着证明下又对了对拍也没问题？
这样就能推出错误结论啦
有点方。。。又略略看了其它证明平复一下心情
求路过神犇指正QAQ

#include<bits/stdc++.h>
#define rep(i,k,n) for(int i=k;i<=n;i++)
#define rep2(i,k,n) for(int i=k;i>=n;i--)
using namespace std;
typedef double db;
const int N=105;
const db eps=1e-7;
db a[N][N],deg[N],C[N][N];
int n,p[N];
db gauss(int m){
    db op=1.0;
    rep(i,1,m){int r=i;
        rep(j,i+1,m)if(fabs(a[j][i]) > fabs(a[r][i]))r=j;///
        if(r!=i){
            rep(j,1,m)swap(a[r][j],a[i][j]);
            op*=-1.0;//
        }
        if(fabs(a[i][i])<eps)return 0.0;
        rep(j,1,m)if(j!=i){
            rep2(k,m,i)
                a[j][k]-=a[i][k]*(a[j][i]/a[i][i]);
        }
    }db ans=1.0;
    rep(i,1,m)ans*=a[i][i];
    return op*ans;
}
int main(){
    int x,y,m;
    scanf("%d%d",&n,&m);
    rep(i,1,n)p[i]=i;
    random_shuffle(p+1,p+n+1);
    rep(i,1,m){
        scanf("%d%d",&x,&y);
        x=p[x],y=p[y];
        a[x][y]--;
        a[y][x]--;
        deg[x]++,deg[y]++;
    }rep(i,1,n)a[i][i]=deg[i];
    printf("%lf",gauss(n-1));
}