bzoj4710 [Jsoi2011]分特产（容斥原理）

题目链接

分析：
这道题的限制其实挺不明显的，应该是“每个人都至少有一个”
也就是说对于所有的物品，将其划分成n部分，每部分不能为空，问总的方案数

我们把每一种特产分开考虑
把$m$个特产分成$n$份，简单的我们可以利用插板法：$C(m+n-1,n-1)$
$ans=∏_{i=1}^{m}C(a_i+n-1,n-1)$
但是这样求出的方案允许有空集

$没有空集=至少有0个空集-至少有1个空集+至少有2个空集...$

至少有$i$个空集的方案：
我们可以从$n$个人中选出$i$个：$C(n,i)$，强制为空集
$ans_i=C(n,i)*∏_{i=1}^{m}C(a_i+n-i-1,n-i-1)$

#include<cstdio>
#include<cstring>
#include<iostream>
#define ll long long

using namespace std;

const int N=2005;
const ll p=1e9+7;
ll c[N][N];
int n,m,a[N];

void prepare() {
    c[0][0]=1;
    for (int i=1;i<=2000;i++) {
        c[i][0]=1;
        for (int j=1;j<=i;j++)
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%p;
    }
}

int main() {
    prepare();
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++) scanf("%d",&a[i]);
    ll ans=1;
    for (int i=1;i<=m;i++)
        ans=(ans*c[a[i]+n-1][n-1])%p;
    ll o=-1;
    for (int i=1;i<=n;i++) {
        ll t=c[n][i]%p;
        for (int j=1;j<=m;j++)
            t=(t*c[a[j]+n-i-1][n-i-1])%p;
        ans=(ans+o*t%p+p)%p;
        o*=-1;
    }
    printf("%lld\n",(ans%p+p)%p);   
    return 0;
}