标签:等价 正整数 列合并 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