spoj highways(生成樹(shù)計(jì)數(shù))

最簡(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,說(shuō)明對(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)//生成樹(shù)計(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;
}
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書(shū)系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

推薦閱讀更多精彩內(nèi)容