bzoj千题计划273:bzoj4710: [Jsoi2011]分特产
2018-03-11 20:04
513 查看
http://www.lydsy.com/JudgeOnline/problem.php?id=4710
答案=总方案数-不合法方案数
f[i][j] 前i种特产分给j个人(可能有人没有分到特产)的总方案数
考虑第i种特产的分配f[i][j]=f[i-1][j]*C(a[i]+j-1 , j-1)
g[i] 表示有i个人,每个人至少分到一种特产,其他人都没有分到的方案数
g[i]=f[m][i]-Σg[j]*C(i,j) j∈[1,i-1]
即有i个人分到特产=总方案数-只有1个人分到特产-只有2个人分到特产……-只有i-1个人分到特产
答案=总方案数-不合法方案数
f[i][j] 前i种特产分给j个人(可能有人没有分到特产)的总方案数
考虑第i种特产的分配f[i][j]=f[i-1][j]*C(a[i]+j-1 , j-1)
g[i] 表示有i个人,每个人至少分到一种特产,其他人都没有分到的方案数
g[i]=f[m][i]-Σg[j]*C(i,j) j∈[1,i-1]
即有i个人分到特产=总方案数-只有1个人分到特产-只有2个人分到特产……-只有i-1个人分到特产
#include<cstdio> #include<iostream> #define N 1001 using namespace std; const int mod=1e9+7; int f ,g ; int a ; int C[N<<1][N<<1]; void read(int &x) { x=0; char c=getchar(); while(!isdigit(c)) c=getchar(); while(isdigit(c)) { x=x*10+c-'0'; c=getchar(); } } int main() { int n,m; read(n); read(m); for(int i=1;i<=m;++i) read(a[i]); C[0][0]=1; int mm=1000+n; for(int i=1;i<=mm;++i) { C[i][0]=1; for(int j=1;j<=i;++j) C[i][j]=(C[i-1][j-1]+C[i-1][j])%mod; } for(int i=1;i<=n;++i) f[0][i]=1; for(int i=1;i<=m;++i) for(int j=1;j<=n;++j) f[i][j]=1LL*f[i-1][j]*C[a[i]+j-1][j-1]%mod; for(int i=1;i<=n;++i) { g[i]=f[m][i]; for(int j=1;j<i;++j) g[i]=(g[i]-1LL*C[i][j]*g[j]%mod+mod)%mod; } printf("%d",g ); }
相关文章推荐
- [BZOJ4710][JSOI2011]分特产(组合数+容斥原理)
- Bzoj4710--Jsoi2011分特产
- bzoj 4710: [Jsoi2011]分特产 排列组合+容斥原理
- bzoj 4710 [Jsoi2011]分特产
- 【bzoj4710】[Jsoi2011]分特产 容斥原理+组合数学
- [bzoj4710][JSOI2011]分特产
- BZOJ 4710: [Jsoi2011]分特产 [容斥原理]
- [BZOJ4710]4710: [Jsoi2011]分特产 容斥原理+组合数学
- bzoj 4710: [Jsoi2011]分特产 (容斥原理+DP)
- bzoj 4710: [Jsoi2011]分特产
- [Jsoi 2011] bzoj4710 分特产 [容斥原理+组合数学]
- ●BZOJ 4710 [Jsoi2011]分特产
- [BZOJ4710][JSOI2011]分特产(组合数学+容斥原理)
- bzoj 4710 : [Jsoi2011]分特产
- 【BZOJ 4710】【JSOI 2011】分特产【计数&容斥】
- Bzoj4710 [Jsoi2011]分特产
- bzoj 4710: [Jsoi2011]分特产
- bzoj4710 [Jsoi2011]分特产 容斥原理
- [BZOJ4710][Jsoi2011]分特产(容斥原理+组合数学)
- 【BZOJ 4710】 4710: [Jsoi2011]分特产 (容斥原理)