生产树计数

#include <iostream>
#include <algorithm>
#include <string.h>
#include <queue>
#include <stdio.h>
#define mod 1000000000
using namespace std;
typedef long long ll;
typedef long long ll;
char map[120][120];
int num[120][120];
int c[120][120];
int n,m;
int dir[4][2]={ {1,0},{-1,0} ,{0,1},{0,-1} };
bool judge(int x,int y)
{
    if(x<0||x>=n||y<0||y>=m) return 0;
    return 1;
}
ll det(int a[][120],int n)
{
    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]-t*a[j][k]+mod)%mod;
                }
                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]) %mod;
    }
    return (ret+mod)%mod;
}

int main()
{
    memset(num,-1,sizeof(num));
    memset(c,0,sizeof(c));
    scanf("%d%d",&n,&m);
    int cnt=0;
    for(int i=0;i<n;++i)
    {
        for(int j=0;j<m;++j)
        {
            scanf(" %c",&map[i][j]);
            if(map[i][j]=='.')   num[i][j]= cnt++;
        }
    }
    for(int i=0;i<n;++i)
    {
        for(int j=0;j<m;++j)
        {
            if(num[i][j]!=-1)
            {
                for(int k=0;k<4;++k)
                {
                    int xx=i+dir[k][0],yy=j+dir[k][1];
                    if(judge(xx,yy)==0||num[xx][yy]==-1) continue;
                    int u=num[i][j],v=num[xx][yy];
                    c[u][u]++; c[u][v]=-1;
                }
            }
        }
    }
    printf("%I64d\n",det(c,cnt));
    return 0;
}

如果需要虚根的话，虚根到所有根的权值为总权值加1，最后再减去即可

模板 HDU 2121

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <queue>
typedef long long ll;
using namespace std;
const ll inf=0x3f3f3f3f;
const int maxn=1010;
struct node
{
     int u,v;
     ll w;
}edge[10010];

ll in[maxn];
int pre[maxn],vis[maxn],id[maxn];
int n,m;
int pos=0;
ll sum=0;
ll mst(int root,int nn,int mm)
{
    ll ret=0;
    while(1)
    {
        for(int i=0;i<nn;++i) in[i]=inf;
        for(int i=0;i<mm;++i)
        {
            int u=edge[i].u  ,  v=edge[i].v;
            if(edge[i].w<in[v] && u!=v)
            {
                in[v]=edge[i].w; pre[v]=u;
                if(u==root) pos=i;
            }
        }

        for(int i=0;i<nn;++i)
        {
            if(i==root) continue;
            if(in[i]==inf) return -1;
        }
        int cnt=0;
        memset(id,-1,sizeof(id));
        memset(vis,-1,sizeof(vis));
        in[root]=0;
        for(int i=0;i<nn;++i)
        {
            ret+=in[i];
            int v=i;
            while(vis[v]!=i && id[v]==-1 && v!=root)
            {
                  vis[v]=i;    v=pre[v];
            }
            if(v!=root && id[v]==-1)
            {
                for(int u=pre[v]; u!= v; u=pre[u]) {  id[u] =cnt; }
                id[v]=cnt++;
            }
        }
        if(cnt==0) break;
        for(int i=0;i<nn;++i)
        {
            if(id[i]==-1) id[i]=cnt++;
        }
        for(int i=0;i<mm;++i)
        {
            int u=edge[i].u, v=edge[i].v;
            edge[i].u=id[u]; edge[i].v=id[v];
            if(id[u]!=id[v]) edge[i].w-=in[v];
        }
        nn=cnt; root=id[root];
    }
    return ret;
}
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
         sum=0;
         for(int i=0;i<m;++i)
         {
             int a,b;
             ll c;
             scanf("%d %d %I64d",&a,&b,&c);
             a++; b++;   sum+=c;
             edge[i].u=a; edge[i].v=b; edge[i].w=c;
         }
         sum+=1;
         for(int i=m;i<n+m;++i)
         {
             edge[i].u=0; edge[i].v=i-m+1; edge[i].w=sum;
         }
         ll temp=mst(0,n+1,n+m);
         if(temp==-1||temp-sum>=sum) printf("impossible\n");
         else printf("%I64d %d\n",temp-sum,pos-m);
         puts("");
    }
    return 0;
}