标签:等价 正整数 列合并 inpu time stream long 博客 www
本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作。
本文作者:ljh2000
作者博客:http://www.cnblogs.com/ljh2000-jump/
转载请注明出处,侵权必究,保留最终解释权!
输出仅一行,包含一个正整数,表示合法质量序列的数目对 10^9+7取模的结果。
不同的合法序列共5个,如下所示:
正解:DP+组合数学
解题报告:
考虑先把相等的变量用并查集合起来。
题目中说只存在最多一个质量不比i差的数,这不就是一棵树吗QAQ
新建一个虚拟根,向入度为$0$的点连边,跑一遍树形$DP$。
用$f[i][j]$表示以$i$为根的子树中分为$j$个不等的段的方案数。
转移的话不妨考虑一下:如何把长为$a$、$b$的两个序列合并。
假设合并后长度为$c$,那么$max(a,b)<=c<=a+b$,而且$f[i][c]=f[v1][a]*f[v2][b]*C_{c}^{a}*C_{a}^{b-(c-a)}$,
这个的含义是先在$c$个位置中选$a$个则剩下的$c-a$个显然只能放第二个序列中的数,下一步就是从$b$个数中选出一些数与之前的$a$个数中的一些合并,一共有$b-(c-a)$这样的需要合并的数。
所以上式就不难解释了。
//It is made by ljh2000 #include <iostream> #include <cstdlib> #include <cstring> #include <cstdio> #include <cmath> #include <algorithm> #include <ctime> #include <vector> #include <queue> #include <map> #include <set> #include <string> #include <complex> using namespace std; typedef long long LL; const int MOD = 1000000007; const int MAXN = 120; int n,m,father[MAXN],cnt,first[MAXN],ecnt,to[MAXN],next[MAXN],in[MAXN],size[MAXN]; LL C[MAXN][MAXN],f[MAXN][MAXN],p[MAXN],ans; char ch[12]; //只存在最多一个质量不比i差的数,这不就是一棵树吗QAQ struct edge{ int x,y; }e[MAXN]; inline int find(int x){ if(father[x]!=x) father[x]=find(father[x]); return father[x]; } inline void link(int x,int y){ next[++ecnt]=first[x]; first[x]=ecnt; to[ecnt]=y; } inline int getint(){ int w=0,q=0; char c=getchar(); while((c<‘0‘||c>‘9‘) && c!=‘-‘) c=getchar(); if(c==‘-‘) q=1,c=getchar(); while (c>=‘0‘&&c<=‘9‘) w=w*10+c-‘0‘,c=getchar(); return q?-w:w; } inline void dfs(int x,int fa){ f[x][0]=1; int lim; LL now; for(int i=first[x];i;i=next[i]) { int v=to[i]; if(v==fa) continue; dfs(v,x); size[x]+=size[v]; for(int j=0;j<=size[x];j++) { for(int k=1;k<=size[v];k++) { lim=min(n,j+k); for(int l=max(j,k);l<=lim;l++) { now=C[l][j]*C[j][k-(l-j)]; now%=MOD; now*=f[x][j]; now%=MOD; now*=f[v][k-1]; now%=MOD; p[l]+=now; p[l]%=MOD; } } } for(int j=0;j<=size[x];j++) f[x][j]=p[j],p[j]=0; } size[x]++; } inline void work(){ n=getint(); m=getint(); int x,y; for(int i=0;i<=n;i++) father[i]=i,C[i][0]=1; for(int i=1;i<=n;i++) for(int j=1;j<=i;j++) C[i][j]=C[i-1][j-1]+C[i-1][j],C[i][j]%=MOD; for(int i=1;i<=m;i++) { x=getint(); scanf("%s",ch); y=getint(); if(ch[0]==‘=‘) father[find(y)]=find(x); else e[++cnt].x=x,e[cnt].y=y; } for(int i=1;i<=cnt;i++) link(find(e[i].x),find(e[i].y)),in[find(e[i].y)]++; cnt=0; for(int i=1;i<=n;i++) if(find(i)==i) { if(in[find(i)]==0) link(n+1,i); cnt++; } dfs(n+1,0); for(int i=1;i<=n;i++) ans+=f[n+1][i],ans%=MOD; printf("%lld",ans); } int main() { work(); return 0; }
标签:等价 正整数 列合并 inpu time stream long 博客 www
原文地址:http://www.cnblogs.com/ljh2000-jump/p/6431328.html