矩阵树定理

矩阵树定理是一个计数定理，用于解决连通图的生成树计数问题。

---

无向图

若连通图G的邻接矩阵为A，将A的对角线(i,i)元素依次换为节点i的度，其余元素(i,j) (j!=i) 取Aij的相反数，所得矩阵记为M，则M的每个代数余子式相等，且等于G的生成树的数目.这就是矩阵-树定理.

一个特例:完全图的生成树个数是

---

有向图

若连通图G的邻接矩阵为A，

外向树，对角线(i,i)元素为节点i的入度；

内向树，对角线(i,i)元素为节点的出度；

其余元素，均为Aij的相反数。

删去根所在的行列，剩余元素组成的行列式即为生成树的数目。

---

应用矩阵树定理解连通图的生成树计数问题，也就是求解这个联通图的基尔霍夫（Kirchhoff）矩阵的某个代数余子式的值。对于此类行列式求值的问题，常用的方法是高斯消元法，将原行列式变成一个等价的上三角行列式，进而快速得到其值。其时间复杂度是O()的，

---

例题

Luogu P4336 [SHOI2016]黑暗前的幻想乡

思路

矩阵树定理+容斥原理

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int Mod=1e9+7;
const int N=18;
inline int q_pow(int a,int b,int p){
	int ans=1;
	while(b){
		if(b&1)ans=1ll*ans*a%p;
		a=1ll*a*a%p;  b>>=1;
	}
	return ans%p;
}
int kir[N][N];//Kirchhoff矩阵 
int m[N],u[N][N*N],v[N][N*N],siz[1<<N];
int det(int n){
	int ans=1;  bool f=0;
	for(int i=1;i<=n;i++){
		if(!kir[i][i]){
			for(int j=i+1;j<=n;j++){
				if(kir[j][i]){
					swap(kir[i],kir[j]);
					f^=1;  break;
				}
			}
		}
		if(!kir[i][i])return 0;
		for(int j=i;j<=n;j++){
			if(!kir[j][i])continue;
			int div=q_pow(kir[j][i],Mod-2,Mod);  ans=1ll*ans*kir[j][i]%Mod;
			for(int k=i;k<=n;k++)kir[j][k]=1ll*kir[j][k]*div%Mod;
		}
		for(int j=i+1;j<=n;j++){
			if(!kir[j][i])continue;
			for(int k=i;k<=n;k++){
				kir[j][k]-=kir[i][k];
				if(kir[j][k]<0)kir[j][k]+=Mod;
			}
		}
	}
	return f?(ans?Mod-ans:0):ans;
}
int main(){
	int n,t,ans=0;
	scanf("%d",&n);  t=(1<<n-1)-1;
	for(int i=1;i<n;i++){
		scanf("%d",&m[i]);
		for(int j=1;j<=m[i];j++)
			scanf("%d%d",&u[i][j],&v[i][j]);
	}
	for(int i=1;i<=t;i++)siz[i]=siz[i>>1]+(i&1);
	for(register int i=1;i<=t;i++){
		memset(kir,0,sizeof(kir));
		for(register int j=1,p=i;p;p>>=1,j++){
			if(!(p&1))continue;
			for(register int k=1;k<=m[j];k++){
				int U=u[j][k],V=v[j][k];
				kir[U][U]++;  kir[V][V]++;
				kir[U][V]--;  kir[V][U]--;
				if(kir[U][V]<0)kir[U][V]+=Mod;
				if(kir[V][U]<0)kir[V][U]+=Mod;
			}
		}
		ans+=n-siz[i]&1?det(n-1):-det(n-1);//det(n-1) :去掉一行一列
		if(ans>=Mod)ans-=Mod;
		else if(ans<0)ans+=Mod;
	}
	printf("%d\n",ans);
	return 0;
}

---

拓展：带权图

求所有生成树的权值和问题，权值的定义：生成树的所有边权之积。

问题转化，将边权看做多条边，求解此时生成树的个数，即为题目所求。

例题：

P6178 【模板】Matrix-Tree 定理

思路：

Aij存的是边权，同理可求

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cctype> 
using namespace std;
const int Mod=1e9+7;
const int N=305;
int kir[N][N];
inline int read(){
	char c=getchar();  int x=0;
	while(!isdigit(c))c=getchar();
	while(isdigit(c))x=(x<<3)+(x<<1)+(c^48),c=getchar();
	return x;
}
inline int q_pow(int a,int b,int p){
	int ans=1;
	while(b){
		if(b&1)ans=1ll*ans*a%p;
		a=1ll*a*a%p;  b>>=1;
	}
	return ans%p;
}
int det(int l,int r){
	int ans=1;  bool f=0;
	for(int i=l;i<=r;i++){
		if(!kir[i][i]){
			for(int j=i+1;j<=r;j++)
				if(kir[j][i]){
					swap(kir[j],kir[i]);
					f^=1;  break;
				}
		}
		if(!kir[i][i])return 0;
		for(int j=i;j<=r;j++){
			if(!kir[j][i])continue;
			int div=q_pow(kir[j][i],Mod-2,Mod);  ans=1ll*ans*kir[j][i]%Mod;
			for(int k=i;k<=r;k++)kir[j][k]=1ll*kir[j][k]*div%Mod;
		}
		for(int j=i+1;j<=r;j++){
			if(!kir[j][i])continue;
			for(int k=i;k<=r;k++){
				kir[j][k]-=kir[i][k];
				if(kir[j][k]<0)kir[j][k]+=Mod;
			}
		}
	}
	return f?(ans?Mod-ans:0):ans;
}
int main(){
	int n=read(),m=read(),t=read();
	for(int i=1;i<=m;i++){
		int u=read(),v=read(),w=read();
		if(u==v)continue;
		if(t){
			kir[v][v]+=w;
			kir[u][v]-=w;
			if(kir[v][v]>=Mod)kir[v][v]-=Mod;
			if(kir[u][v]<0)kir[u][v]+=Mod;
		}
		else {
			kir[u][u]+=w;  kir[v][v]+=w;
			kir[u][v]-=w;  kir[v][u]-=w;
			if(kir[u][u]>=Mod)kir[u][u]-=Mod;
			if(kir[v][v]>=Mod)kir[v][v]-=Mod;
			if(kir[u][v]<0)kir[u][v]+=Mod;
			if(kir[v][u]<0)kir[v][u]+=Mod;
		}
	}
	printf("%d\n",det(2,n));
	return 0;
}
//https://www.luogu.com.cn/problem/P6178 

---

矩阵-树定理_百度百科