Maxtrix tree 生成树计数

/*度数矩阵D：是一个 N×N 的矩阵，其中
D[i][j]=0(i≠j)，D[i][i]=i号点的度数
邻接矩阵A：是一个 N×N 的矩阵，其中
A[i][i]=0,A[i][j]=A[j][i]=i,j之间的边数
然后基尔霍夫Kirchhoff矩阵K=D-A
对于一个无向图 G ，它的生成树个数等于其基尔霍夫Kirchhoff矩阵任何一个N-1阶主子式的行列式的绝对值。*/

//a邻接矩阵，d度数矩阵，g基尔霍夫矩阵//
#include 
#include 
#include 
#include 
#include 
using namespace std;
typedef long long ll;
const int maxn = 25;
ll a[maxn][maxn],d[maxn][maxn],g[maxn][maxn];
int t,n,m;

ll solve(ll G[][maxn],int N)//求G的n-1阶主子式的值
{
    ll res = 1;
    for(int i = 1;i <= N;i++)
    {
        for(int j = i + 1;j <= N;j++)
        {
            while(G[j][i])
            {
                ll o = G[i][i] / G[j][i];
                for(int k = 1;k <= N;k++)
                {
                    G[i][k] = G[i][k] - G[j][k] * o;
                    swap(G[i][k],G[j][k]);
                }
            }
            res *= -1;
        }
        res *= G[i][i];
    }
    return res;
}

int main()
{
    scanf("%d",&t);
    while(t--)
    {
        memset(a,0,sizeof(a));
        memset(d,0,sizeof(d));
        memset(g,0,sizeof(g));
        int u,v;
        scanf("%d%d",&n,&m);
        while(m--)
        {
            scanf("%d%d",&u,&v);
            d[u][u]++,d[v][v]++;
            a[u][v] = a[v][u] = 1;
        }
        for(int i = 1;i <= n;i++)
        {
            for(int j = 1;j <= n;j++)
            {
                g[i][j] = d[i][j] - a[i][j];
            }
        }
        ll ans = solve(g,n - 1);
        printf("%lld\n",abs(ans));
    }
    return 0;
}