生成树的计数 模板 Matrix-Tree定理

SPOJ104(Highways) 
题目地址： 
http://www.spoj.com/problems/HIGH/ 
题目大意： 
 *一个有n座城市的组成国家,城市1至n编号,其中一些城市之间可以修建高速公路; 
 *需要有选择的修建一些高速公路,从而组成一个交通网络; 
 *计算有多少种方案,使得任意两座城市之间恰好只有一条路径; 

代码：

#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>using namespace std;const int N=15;typedef long long LL;int degree[N];LL C[N][N];LL det(LL a[][N],int n)//生成树计数:Matrix-Tree定理{    LL ret=1;    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]);                ret=-ret;            }        if(a[i][i]==0)            return 0;        ret=ret*a[i][i];    }    if(ret<0)        ret=-ret;    return ret;}int main(){    int tcase;    scanf("%d",&tcase);    while(tcase--)    {        memset(degree,0,sizeof(degree));        memset(C,0,sizeof(C));        int n,m;        scanf("%d%d",&n,&m);        int u,v;        while(m--)        {            scanf("%d%d",&u,&v);            u--;            v--;            C[u][v]=C[v][u]=-1;            degree[u]++;            degree[v]++;        }        for(int i=0; i<n; ++i)            C[i][i]=degree[i];        printf("%lld\n",det(C,n));    }    return 0;}