spoj104 Highways（矩阵树定理）

Matrix-tree定理裸题。
对于一张无向图G，我们有度数矩阵$D[G]$

$$d_{ij}= \begin{cases} du[i], &\text{$i=j$}\\ 0, &\text{$i \neq j$} \end{cases}$$
我们还有邻接矩阵 $A[G]$
$$a_{ij}= \begin{cases} 1, &\text{i,j有边}\\ 0, &\text{i,j无边} \end{cases}$$
我们的基尔霍夫矩阵 $C[G]=D[G]-A[G]$
无向图G的生成树个数就是基尔霍夫矩阵的任一n-1阶主子式的行列式的绝对值。
证明可以参见周冬2007集训队论文:《生成树计数及其应用》

而求解行列式的值，我们可以高斯消元得到上三角矩阵。
然后把主对角线上的数都乘起来取绝对值就是行列式的绝对值了。
复杂度$O(n^3)$

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 20
#define eps 1e-8
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(T==S){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m;
double a[N][N];
inline void Gauss(){
    n--;
    for(int i=1;i<=n;++i){
        int r=i;
        for(int j=i+1;j<=n;++j) if(fabs(a[j][i])>fabs(a[r][i])) r=j;
        if(fabs(a[r][i])<eps){puts("0");return;}
        if(r!=i) for(int j=i;j<=n;++j) swap(a[r][j],a[i][j]);
        for(int j=i+1;j<=n;++j){
            double t=a[j][i]/a[i][i];
            for(int k=i;k<=n;++k) a[j][k]-=t*a[i][k];
        }
    }double ans=1;
    for(int i=1;i<=n;++i) ans*=a[i][i];
    printf("%.0lf\n",ans);
}
int main(){
//  freopen("a.in","r",stdin);
    int tst=read();
    while(tst--){
        n=read();m=read();memset(a,0,sizeof(a));
        while(m--){
            int x=read(),y=read();
            a[x][x]++;a[y][y]++;a[x][y]--;a[y][x]--;
        }Gauss();
    }
}