和Leo一起做爱数学的好孩子之容斥 [Jsoi2011]分特产

Description
JYY 带队参加了若干场ACM/ICPC 比赛，带回了许多土特产，要分给实验室的同学们。
JYY 想知道，把这些特产分给N 个同学，一共有多少种不同的分法？当然，JYY 不希望任
何一个同学因为没有拿到特产而感到失落，所以每个同学都必须至少分得一个特产。
例如，JYY 带来了2 袋麻花和1 袋包子，分给A 和B 两位同学，那么共有4 种不同的
分配方法：
A：麻花，B：麻花、包子
A：麻花、麻花，B：包子
A：包子，B：麻花、麻花
A：麻花、包子，B：麻花
Input
输入数据第一行是同学的数量N 和特产的数量M。
第二行包含M 个整数，表示每一种特产的数量。
N, M 不超过1000，每一种特产的数量不超过1000
Output
输出一行，不同分配方案的总数。由于输出结果可能非常巨大，你只需要输出最终结果
MOD 1,000,000,007 的数值就可以了。
Sample Input
5 4

1 3 3 5
Sample Output
384835

很不错的容斥题
本来想的是除一个排列实际上不行
发现特产本身是独立的
所以考虑单独考虑
计算随便分的减去一个不行。。。
∏i=1mCAi+n−1n−1\prod_{i=1}^mC_{A_i+n-1}^{n-1}∏i=1m​CAi​+n−1n−1​这是全部不考虑分不分
容斥后实际就是这个东西
∑i=0n(−1)iCni∏j=1mCAi+n−i−1n−i−1\sum_{i=0}^n(-1)^iC_n^i\prod_{j=1}^mC_{A_i+n-i-1}^{n-i-1}∑i=0n​(−1)iCni​∏j=1m​CAi​+n−i−1n−i−1​

/**************************************************************
    Problem: 4710
    User: Leo_JAM
    Language: C++
    Result: Accepted
    Time:676 ms
    Memory:50360 kb
****************************************************************/
 
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
typedef int INT;
#define int long long 
const int mod=1e9+7;
int C[2200][2200];
int F[1200][1200];
void Pre(){
    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-1]+C[i-1][j])%mod;
        }
    }
}
int ans=0;
int n,m;
int G[1200];
INT main(){
    Pre();
    scanf("%lld%lld",&n,&m);
    for(int i=1;i<=m;++i)scanf("%lld",&G[i]);
    int cur=-1;
    for(int i=0;i<n;++i){
        cur=cur*(-1);
        int now=cur;
        for(int j=1;j<=m;++j){
            now=(now*C[G[j]+n-i-1][n-i-1]%mod);
        }   
        now=now*C[n][i]%mod;
        ans=(ans+now)%mod;
    }
    cout<<(ans+mod)%mod;
    return 0;
} 


转载于:https://www.cnblogs.com/Leo-JAM/p/10079087.html