[容斥][排列组合] BZOJ4710 && Jsoi2011.分特产

把M个物品分给N个人的方案数,可以用插板法得到,为 (M+N1N1)
不考虑每个人至少分到一个,这些特产的总方案数为 (Ai+N1N1)
发现每个人至少分到一个这个限制很麻烦,反过来的话就比较好做,容斥就可以了

#include <cstdio>
#include <iostream>
#include <algorithm>

using namespace std;

const int N=3010,P=1e9+7;

int n,m,a[N],fac[N],inv[N];

inline void Pre(){
  fac[0]=inv[0]=inv[1]=1;
  for(int i=1;i<=3000;i++) fac[i]=1LL*fac[i-1]*i%P;
  for(int i=2;i<=3000;i++) inv[i]=1LL*(P-P/i)*inv[P%i]%P;
  for(int i=1;i<=3000;i++) inv[i]=1LL*inv[i]*inv[i-1]%P;
}

inline void add(int &x,int y){
  if((x+=y)>=P) x-=P;
}

inline int C(int x,int y){
  return 1LL*fac[x]*inv[y]%P*inv[x-y]%P;
}

int main(){
  scanf("%d%d",&n,&m); Pre();
  for(int i=1;i<=m;i++)
    scanf("%d",&a[i]);
  int ans=0;
  for(int i=0,f=1;i<=n;i++,f=-f){
    int prod=1;
    for(int j=1;j<=m;j++)
      prod=1LL*prod*C(a[j]+n-i-1,n-i-1)%P;
    add(ans,(P+1LL*f*prod*C(n,i)%P)%P);
  }
  printf("%d\n",ans);
  return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值