https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1707
生成树计数即是统计一张图中一共有多少种构造生成树的方案,需要用到统计数学。
Matrix-Tree定理是解决生成树计数最有力的武器之一。
前置概念:
1.G的度数矩阵D[G]是一个nn的矩阵,并且满足i当i≠j时,dij=0;当i=j时,dij等于vi的度数。
2.G的邻接矩阵A[G]是也是一个nn的矩阵,并且满足如果vi、vj之间有边直接相连,则aij=1,否则为0。
我们定义G的Kirchhoff矩阵C[G]=D[G]-A[G],则Matrix-Tree定理可以描述为:
G的所有不同的生成树的个数等于其Kirchhoff矩阵C[G]任何一个n-1阶子式的行列式的绝对值。
所谓n-1阶子式,就是对于r(1<=r<=n),将C[G]的第r行、第r列同时去掉后得到的新矩阵,用Cr[G]表示。
#include <iostream>
#include <algorithm>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <vector>
#include <string>
#include <cmath>
using namespace std;
typedef long long ll;
const int maxn = 50 + 3;
int g[maxn][maxn], a[maxn][maxn];
ll Det(int n) {
ll ret = 1;
int sign = 0;//记录行交换的次数
//矩阵初等变换
for(int i = 1; i <= n; ++i) {
for(int j = i + 1; j <= n; ++j) {
while(a[j][i]) {
ll t = a[i][i] / a[j][i];
for(int k = i; k <= n; ++k)
a[i][k] = (a[i][k] - a[j][k] * t);
for(int k = i; k <= n; ++k)
swap(a[i][k], a[j][k]);
++sign;
}
}
if(a[i][i] == 0) return 0;
ret *= a[i][i];
/*cout << endl;
for(int i = 1; i <= n; ++i) {
for(int j = 1; j <= n; ++j) {
cout << a[i][j] << ' ';
}
cout << endl;
}
*/
}
if(sign & 1) ret = -ret;
return ret;
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);
int n, m, k;
while(cin >> n >> m >> k) {
memset(g, 0, sizeof(g));
memset(a, 0, sizeof(a));
int u, v;
for(int i = 0; i < m; ++i) {
cin >> u >> v;
g[u][v] = g[v][u] = 1;
}
for(int i = 1; i <= n; ++i) {
int d = 0;
for(int j = 1; j <= n; ++j) {
if(i != j && g[i][j] == 0) {
d++;
a[i][j] = -1;
}
}
a[i][i] = d;
}
cout << Det(n - 1) << endl;
}
}