• bzoj4013:[HNOI2015]实验比较


    传送门

    好神的树形dp啊!
    首先考虑=的情况,显然可以将他们缩成一个点
    然后<号可以建出一棵树或者森林来,特判一下成环就无解
    由于有可能出现森林,考虑新建一个根,然后跑tree dp
    我们可以设用(f[i][j])表示以(i)为根的子树中分为(j)个不等的段的方案数。
    然后发现就是一个背包,转移的时候考虑合并
    假设合并后长度为(len),那么(f[i][len]=f[x][len_x]∗f[son][len_{son}]∗C^{len}_{len_x}∗C^{len_{son}−(len−len_x)}_{len_x})
    借用网上形象的比喻
    看做有(k)个盒子,(a)个白球,(b)个黑球,并且(a+b>=k)
    那么(a)个白球放的方案数就是(C^{a}_{k})
    然后(b)个黑球部分要与白球合并,那么合并的方案数就是(C^{b-(k-a)}_a)
    代码:

    #include<cstdio>
    #include<iostream>
    #include<algorithm>
    using namespace std;
    void read(int &x) {
    	char ch; bool ok;
    	for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch=='-') ok=1;
    	for(x=0; isdigit(ch); x=x*10+ch-'0',ch=getchar()); if(ok) x=-x;
    }
    #define rg register
    const int maxn=110,mod=1e9+7;char s[maxn];
    int x[maxn],y[maxn],n,m,fa[maxn],ans,size[maxn],rt,f[maxn],cnt;
    int gg[maxn],g[maxn][maxn],pre[maxn],nxt[maxn],h[maxn],id[maxn],c[maxn][maxn];
    int find(int x){return x==fa[x]?x:fa[x]=find(fa[x]);}
    void add(int x,int y){pre[++cnt]=y,nxt[cnt]=h[x],h[x]=cnt;}
    void prepare(int n)
    {
    	for(rg int i=0;i<=n;i++)c[i][0]=1;
    	for(rg int i=1;i<=n;i++)
    		for(rg int j=1;j<=i;j++)
    			c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
    }
    void dfs(int x)
    {
    	size[x]=1,g[x][1]=1;
    	for(rg int i=h[x];i;i=nxt[i])
    	{
    		dfs(pre[i]);
    		for(rg int j=1;j<=size[x]+size[pre[i]];j++)
    		{
    			gg[j]=0;
    			for(rg int k=1;k<=size[x];k++)
    				for(rg int t=1;t<=size[pre[i]];t++)
    					if(k+t>=j)(gg[j]+=1ll*g[x][k]*g[pre[i]][t]%mod*c[j-1][k-1]%mod*c[k-1][k+t-j]%mod)%=mod;
    		}
    		size[x]+=size[pre[i]];
    		for(rg int j=1;j<=size[x];j++)g[x][j]=gg[j];
    	}
    }
    int main()
    {
    	read(n),read(m),rt=n+1,prepare(n);
    	for(rg int i=1;i<=n;i++)fa[i]=i;
    	for(rg int i=1;i<=m;i++)scanf("%d %c %d",&x[i],&s[i],&y[i]);
    	for(rg int i=1;i<=m;i++)if(s[i]=='=')fa[find(x[i])]=find(y[i]);
    	for(rg int i=1;i<=n;i++)id[i]=find(i);
    	for(rg int i=1;i<=m;i++)
    		if(s[i]!='=')
    		{
    			add(id[x[i]],id[y[i]]),f[id[y[i]]]=id[x[i]];
    			int a=find(id[x[i]]),b=find(id[y[i]]);
    			if(a==b){printf("0
    ");return 0;}
    			else fa[a]=b;
    		}
    	for(rg int i=1;i<=n;i++)if(!f[i]&&id[i]==i)add(rt,i);
    	dfs(rt);for(rg int i=1;i<=rt;i++)(ans+=g[rt][i])%=mod;
    	printf("%d
    ",ans);
    }
    
  • 相关阅读:
    深入理解link和@import到底有什么区别?
    你应该知道的简单易用的CSS技巧
    META标签的设置
    前端webp图片
    PAT 1130 Infix Expression[难][dfs]
    PAT 1118 Birds in Forest [一般]
    生信-cufflinks输入输出文件分析
    PAT 1121 Damn Single[简单]
    PAT 1049 Counting Ones [难]
    NGS中的一些软件功能介绍
  • 原文地址:https://www.cnblogs.com/lcxer/p/10578415.html
Copyright © 2020-2023  润新知