最簡(jiǎn)單的那個(gè)模板題,這個(gè)題spoj有兩個(gè)同名的,該代碼為輸入第一行為4那個(gè)
兩個(gè)模板
//這個(gè)det函數(shù)有點(diǎn)不好用
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int MAX=20;
double deg[MAX];
double g[MAX][MAX];
double b[MAX][MAX];
double Det(int n){
int i,j,k;
int sign=0;
double ret=1;
for(i=1;i<=n;i++){
if(!b[i][i]){ //如果該點(diǎn)未被連接
j=i+1;
while(!b[j][i]&&j<=n){//找不為0的一行替換
j++;
}
if(j>n)return 0;//如果大于n,說明對(duì)角線一定有一個(gè)0钉寝,
for(k=i;k<=n;k++){
swap(b[i][k],b[j][k]); //換列
}
sign++;
}
if(!b[i][i]) return 0;
ret*=b[i][i];
for(j=i+1;j<=n;j++){ //處理矩陣
double tmp=b[j][i]/b[i][i];
for(k=i+1;k<=n;k++){
b[j][k]-=tmp*b[i][k];
}
}
}
if(sign&1) ret=-ret;
// cout << sign << endl;
return ret;
}
int main()
{
int T,n,m;
int u,v;
scanf("%d",&T);
while(T--){
scanf("%d%d",&n,&m);
memset(deg,0,sizeof(deg));
memset(g,0,sizeof(g));
memset(b,0,sizeof(b));
while(m--){
scanf("%d%d",&u,&v);
g[u][v]=1;g[v][u]=1; //存入鄰接矩陣
deg[u]++; deg[v]++; //計(jì)算度數(shù)
}
for(int i=1;i<=n;i++){
b[i][i]=deg[i]; //建立度數(shù)矩陣
// cout << b[i][i] << endl;
}
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
b[i][j]-=g[i][j]; //把鄰接矩陣減到度數(shù)矩陣中,負(fù)數(shù)方便分辨是鄰接矩陣還是度數(shù)矩陣
}
} //下面注釋的部分可以打印出矩陣
/*for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
printf(" %.0lf",b[i][j]);
}puts("");
}*/
printf("%.2lf\n",Det(n-1));//傳入的是n-1,任意n-1階就行
}
return 0;
}
下面這個(gè)轉(zhuǎn)自http://blog.csdn.net/jarily/article/details/8901363
#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)//生成樹計(jì)數(shù):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()
{
//freopen("C:\\Users\\Administrator\\Desktop\\kd.txt","r",stdin);
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;
}