• 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);
    }
    
  • 相关阅读:
    4.2. 入门案例
    4.1. Scrapy配置安装
    4. Scrapy框架
    2.6. 案例:使用BeautifuSoup4的爬虫
    2.03_01_Python网络爬虫urllib2库
    2.06_Python网络爬虫_正则表达式
    2.04_Python网络爬虫_Requests模块
    2.03_Python网络爬虫Http和Https
    TCP协议的三次握手
    从零开始搭建VUE项目
  • 原文地址:https://www.cnblogs.com/lcxer/p/10578415.html
Copyright © 2020-2023  润新知