可看周冬的论文《生成树的计数及其应用》,利用Matrix-Tree定理解决生成树计数的问题,复杂度是矩阵乘法的复杂度O(n^3)。
总结:
无向图,允许有重边。
四个重要矩阵A(邻接矩阵),D(度数矩阵),C(KirchHoff矩阵,C=D-A),B(关联矩阵,B其实是用来证明和理解的)。
构造出C矩阵后,C的任何一个n-1阶主子式的行列式的绝对值就是答案。
ps:注意double矩阵求行列式的精度问题。
SPOJ104是论文里的例题。
#include<iostream> #include<cstdio> #include<cstring> #include<cmath> #include<algorithm> using namespace std; const double eps=1e-8; const double inf=1e10; const int maxn=15; int Sign(double x) { return x<-eps?-1:x>eps; } double Fabs(double x) { return (Sign(x)<0)?-x:(x>eps?x:0); } //*****************************Work******************* int Cnt[maxn]; double C[maxn][maxn]; void Print(int n) { int i,j; cout<<"*******************************************"<<endl; for(i=0;i<n;++i) { for(j=0;j<n;++j) { cout<<" "<<C[i][j]; } cout<<endl; } } double Det(int n) //化成下三角形式 { double ret=1.0,tmp; int i,j,k,sign=0; for(i=0;i<n;++i) { for(j=i;j<n;++j) if(Sign(C[j][i])!=0) break; if(j==n) return 0.0; if(j!=i) sign++; for(k=0;k<n;++k) swap(C[i][k],C[j][k]); // Print(n); for(j=i+1;j<n;++j) { tmp=C[i][j]/C[i][i]; for(k=i+1;k<n;++k) { C[k][j]-=tmp*C[k][i]; } } // Print(n); } for(i=0;i<n;++i) ret*=C[i][i]; if(sign&1) ret=-ret; return ret; } int main() { int t,i,j,a,b,cnt,n,m; scanf("%d",&t); while(t--) { memset(Cnt,0,sizeof(Cnt)); memset(C,0,sizeof(C)); scanf("%d %d",&n,&m); while(m--) { scanf("%d %d",&a,&b); C[a-1][b-1]=C[b-1][a-1]=-1; Cnt[a-1]++,Cnt[b-1]++; } for(i=0;i<n;++i) C[i][i]=Cnt[i]; // Print(n); printf("%.0lf\n",Det(n-1)); } return 0; }