SPOJ.104.Highways([模板]Matrix Tree定理 生成树计数)

题目链接

(Description)

一个国家有1~n座城市，其中一些城市之间可以修建高速公路(无自环和重边)。
求有多少种方案，选择修建一些高速公路，组成一个交通网络，使得任意两座城市之间恰好只有一条路径。

(Solution)

生成树计数 直接上Matrix Tree
无解情况别忘了判
MatrixTree定理大体见这吧，证明别的应用什么的先不管了。
基尔霍夫矩阵=度数矩阵-边矩阵。

#include <cmath>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
const int N=15;

int n;
double K[N][N];

inline int read()
{
	int now=0,f=1;register char c=gc();
	for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now*f;
}
double Gauss()
{
	for(int j=1; j<n; ++j)//去掉第一行第一列 
	{
		int mxrow=j;
		for(int i=j+1; i<n; ++i)
			if(fabs(K[i][j])>fabs(K[mxrow][j])) mxrow=i;
		if(!K[mxrow][j]) return 0.0;//别忘判 
		if(mxrow!=j) std::swap(K[mxrow],K[j]);
		for(int i=j+1; i<n; ++i)
			if(K[i][j])
			{
				double t=K[i][j]/K[j][j];
				for(int k=j; k<n; ++k)
					K[i][k]-=t*K[j][k];
			}
	}
	double res=1;
	for(int i=1; i<n/*&& res*/; ++i) res*=K[i][i];
	return res>0?res:-res;
}

int main()
{
	int t=read(),m,u,v;
	while(t--)
	{
		memset(K,0,sizeof K);
		n=read(),m=read();
		while(m--)
			u=read()-1,v=read()-1, ++K[u][u], ++K[v][v], --K[u][v], --K[v][u];
		printf("%.0lf
",Gauss());
	}
	return 0;
}