（生成树计数）Organising the Organisation

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