uva10766生成树计数(矩阵树定理)

更正了我之前打错的地方，有边的话G[i][j]=-1;

WA了好多次，中间要转成long double才行。。这个晚点更新。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
using namespace std;

typedef long double ld;

const int N=60;
const ld eps=1e-9;
int map[N][N];
ld G[N][N];

ld myabs(ld x){return x>0 ? x:-x;}

ld guass(int n)
{
    ld ans=1;
    for(int i=1;i<=n;i++)
    {
        int r=i;
        for(int j=i+1;j<=n;j++)
            if(myabs(G[j][i]) > myabs(G[r][i])) r=j;
        if(r!=i)
        {
            for(int j=1;j<=n;j++) swap(G[i][j],G[r][j]);
            ans*=-1;
        }
        if(myabs(G[i][i])<eps) return 0;
        for(int j=i+1;j<=n;j++)
            for(int k=n;k>=i;k--)
                G[j][k]-=G[j][i]/G[i][i]*G[i][k];
    }
    for(int i=1;i<=n;i++) ans*=G[i][i];
    return myabs(ans);
}

int main()
{
    //freopen("a.in","r",stdin);
    int n,m,k;
    while(scanf("%d%d%d",&n,&m,&k)!=EOF)
    {
        int x,y;
        memset(map,0,sizeof(map));
        memset(G,0,sizeof(G));
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&x,&y);
            map[x][y]=map[y][x]=1;
        }
        for(int i=1;i<=n;i++)
            for(int j=1;j<i;j++)
                if(!map[i][j]) 
                {
                    G[i][j]=-1;G[j][i]=-1;
                    G[i][i]++;G[j][j]++;
                }
        printf("%.0lf
",(double)guass(n-1));
    }
    return 0;
}