生产树计数

来源：互联网发布：彩票app源码下载编辑：程序博客网时间：2024/04/30 09:49

#include <iostream>#include <algorithm>#include <string.h>#include <queue>#include <stdio.h>#define mod 1000000000using 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;}

阅读全文

0 0