BZOJ 4013: [HNOI2015]实验比较

挺不错的题啊 然而本蒟蒻一开始一脸懵逼

然后在dalao的指导下 再看了看 注意到有“对每张图片 i，小 D 都最多只记住了某一张质量不比 i 差的另一张图片 Ki。”这么一句话
所以每个点入度就为1，就是一棵树了。

当然先要判断一下有没有环，环的话就输出0 否则就treedp 乘个组合数什么的就好了
如果要看详细题解的推荐这个：传送门

留个代码

#include<bits/stdc++.h>#define me(a,x) memset(a,x,sizeof a)using namespace std;typedef long long LL;const int N=105,mod=1e9+7;inline LL read(){    char ch=getchar(); LL x=0,f=1;    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();}    return x*f;}struct node{int y,next;}a[N]; int len,first[N];void ins(int x,int y){    a[++len]=(node){y,first[x]},first[x]=len;}int c[N][N],fa[N],ax[N],ay[N],al;int fi(int x){return fa[x]==x?x:fa[x]=fi(fa[x]);}bool v[N],ok;void fc(int x){    v[x]=1;    for(int k=first[x];k;k=a[k].next){        if(v[a[k].y]){ok=0; return;}        fc(a[k].y);    }}int siz[N],d[N],f[N][N],g[N];void dfs(int x){    bool o=1;    for(int k=first[x];k;k=a[k].next){        int y=a[k].y; dfs(y);        if(o){            siz[x]=siz[y],o=0;            for(int i=1;i<=siz[x];++i)f[x][i]=f[y][i];        }        else{            me(g,0);            for(int i=1;i<=siz[x];++i)if(f[x][i])            for(int j=1;j<=siz[y];++j)if(f[y][j])            for(int p=max(i,j);p<=i+j;++p)              g[p]=(g[p]+1ll*f[x][i]*f[y][j]%mod*c[p][i]%mod*c[i][j-p+i]%mod)%mod;            siz[x]+=siz[y];            for(int i=1;i<=siz[x];++i)f[x][i]=g[i];        }    }    if(x){        siz[x]++;        if(o)f[x][1]=1;        else for(int i=siz[x];i;--i)f[x][i]=f[x][i-1];    }}int main(){    int n=read(),m=read(),i,j;    for(i=1;i<=n;++i)fa[i]=i;    for(i=0;i<=n;++i)        for(c[i][0]=j=1;j<=i;++j) c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;    al=0;    for(int i=1;i<=m;++i){        int x=read(); char c; scanf("\n%c",&c);        int y=read();        if(c=='=') x=fi(x),y=fi(y),fa[x]=y;        else ax[++al]=x,ay[al]=y;    }    for(int i=1;i<=al;++i){        int x=fi(ax[i]),y=fi(ay[i]);        if(x==y) return printf("0\n"),0;        ins(x,y),d[y]++;    }    ok=1;    for(int i=1;i<=n;++i) if(!d[i])fc(i);    if(!ok) return printf("0\n"),0;    for(int i=1;i<=n;++i) if(!d[i] && fi(i)==i)ins(0,i);    dfs(0); int ans=0;    for(int i=1;i<=n;++i) ans=(ans+f[0][i])%mod;    printf("%d\n",ans);    return 0;}