【无标题】

对于分特产的容斥以及容斥原因

---

终于透彻地看懂一篇透彻的题解了! ! ! ! !

解决方案

• 若使g[i]表示正好i人分不到特产，则g[0]则为所求.
• 寻求包括关系，不难设f[i]表示至少i人分不到特产.

思考如何求解f数组

根据f的定义，强制规定i个人没有特产，则剩下的人可使用隔板求得方案数，且不保证全部分到，所以对于特产i，分配方式为：

${\text{C}}_{a[i]+n-i-1}^{n-i-1}$

故而

$f[i]$=${\text{C}}_n^i\times {\prod }_{j=0}^m{\text{C}}_{a[i]+n-i-1}^{n-i-1}$

在n个人中选取i个作为怨种，所以额外乘${\text{C}}_n^i$

此时可以发现，若规定x作为 怨种 ，在处理剩下$n-1$个人时正好存在单独y也作为 怨种 的情况，
则运算y时，x也会被考虑。
所以在f[1]中，g[2]被计算两次。
以样例n=5为例

$f[0]=g[0]+g[1]+g[2]+g[3]+g[4]$

$f[1]=g[1]+2\ast g[2]+g[3]\ast 3+g[4]\ast 4$

$f[2]=g[2]+g[3]\ast 3+g[4]\ast 6$

$f[3]=g[3]+g[4]\ast 4$

$f[4]=g[4]$

故而采用容斥原理

$ans={\sum }_{i=0}^{n-1}(-1{)}^{i}f[i];$

#include <bits/stdc++.h>
using namespace std;
const long long modd=1e9+7;
inline long long qr()
{
    long long x=0,f=0;char ch=0;
    while(!isdigit(ch)) {f|=ch=='-';ch=getchar();}
    while(isdigit(ch)) {x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
    return f ? -x :x;
}
long long n,m;
long long num[1010],C[2020][2020];
void firstset()
{
    C[0][0]=1;
    for(int i=1;i<=2005;i++)
    {
        C[i][0]=1;
        for(int j=1;j<=i;j++)
        C[i][j]=(C[i-1][j]+C[i-1][j-1])%modd;
    }
}
int main()
{
    freopen("speciall.in","r",stdin);
    freopen("speciall.out","w",stdout);
    firstset();
    n=qr(),m=qr();
    for(int i=1;i<=m;i++)
        num[i]=qr();
    long long ans=0;
    for(int i=0;i<n;i++)
    {
        long long t=C[n][i];
        for(int j=1;j<=m;j++)
        t=t*C[num[j]+n-i-1][n-i-1]%modd;
        if(i&1) ans=(ans-t+modd)%modd;
        else ans=(ans+t+modd)%modd;
    }
    printf("%lld",(ans+modd)%modd);
    return 0;
}