SPOJ 104 Highways(生成树计数)

题目链接：
SPOJ 104 Highways
题意：
题意：
给$n$个点，$m$条无向边，无重边和自环，选择一些边连通使得任意两点之间恰好有一条路径。
求不同的选择方案数？$(n \leq 12)$
分析：
$Matrix-Tree$定理($Kirchhoff$矩阵-树定理)
1、$G$的度数矩阵$D[G]$是一个$n*n$的矩阵，并且满足：当$i≠j$时,$d[i][j]=0$；当$i=j$时，$d[i][j]$等于$v_i$的度数。
2、$G$的邻接矩阵$A[G]$也是一个$n*n$的矩阵， 并且满足：如果$v_i、v_j$之间有边直接相连，则$a[i][j]=1$，否则为$0\quad (a[i][i]=0)$。

定义$G$的$Kirchhoff$矩阵(也称为拉普拉斯算子)$C[G]$为$C[G]=D[G]-A[G]$，则$Matrix-Tree$定理可以描述为：$G$的所有不同的生成树的个数等于其$Kirchhoff$矩阵$C[G]$任何一个$n-1$阶主子式的行列式的绝对值。所谓$n-1$阶主子式，就是对于$r(1≤r≤n)$，将$C[G]$的第$r$行、第$r$列同时去掉后得到的新矩阵，用$C_r[G]$表示。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <climits>
#include <cmath>
#include <ctime>
#include <cassert>
#define IOS ios_base::sync_with_stdio(0); cin.tie(0);
using namespace std;
typedef long long ll;
const int MAX_N = 20;
const double eps = 1e-6;

int degree[MAX_N];
double C[MAX_N][MAX_N];

int sgn(double x)
{
    if(fabs(x) <= eps) return 0;
    else if(x > 0) return 1;
    else return -1;
}

double det(double mat[][MAX_N], int n)
{
    double res = 1;
    int cnt = 0;
    for(int i = 0; i < n; ++i) {
        if(sgn(mat[i][i]) == 0) { //正对角线上元素为0
            for(int j = i + 1; j < n; ++j) { //从后面的行中找到第i列不为0的行，并交换
                if(sgn(mat[j][i]) != 0) {
                    for(int k = i; k < n; ++k) {
                        swap(mat[i][k], mat[j][k]);
                    }
                    cnt++; //交换行次数
                    break;
                }
            }
        }
        if(sgn(mat[i][i]) == 0) return 0; 
        for(int j = i + 1; j < n; ++j) {
            mat[j][i] /= mat[i][i];  
            //消去倍数,double可以防止整数相除误差，必要时使用long double
        }
        for(int j = i + 1; j < n; ++j) {
            for(int k = i + 1; k < n; ++k) {
                mat[j][k] -= mat[j][i] * mat[i][k];
            }
        }
        res *= mat[i][i];
    }
    if(cnt & 1) res = -res;
    return  res;
}

int main()
{
    int T, n, m;
    scanf("%d", &T);
    while(T--) {
        memset(degree, 0, sizeof(degree));
        memset(C, 0, sizeof(C));
        scanf("%d%d", &n, &m);
        for(int i = 0; i < m; ++i){
            int a, b;
            scanf("%d%d", &a, &b);
            a--; b--;
            C[a][b] = C[b][a] = -1;
            degree[a]++; degree[b]++;
        }
        for(int i = 0; i < n; ++i) { 
            C[i][i] = degree[i];
        }
        printf("%.0lf\n", det(C, n - 1)); 
        //计算n - 1阶,相当于去除了最后一行和最后一列
    }
    return 0;
}