洛谷·[JSOI2011]分特产

最新推荐文章于 2024-05-17 10:00:00 发布

樱狸❀

最新推荐文章于 2024-05-17 10:00:00 发布

阅读量178

点赞数

分类专栏：数论文章标签：广义容斥定理

本文链接：https://blog.csdn.net/qq_43326267/article/details/103715871

版权

数论专栏收录该内容

47 篇文章 2 订阅

订阅专栏

初见安~这里是传送门：洛谷P5505 分特产

题解

比较明显吧，是个容斥的题。

题目要求每个同学至少分到一个特产，所以我们可以假设性质为：该同学没有分到特产。

对于 $\alpha(i)$ 表示有至少i个人没有分到特产，也就是至少满足i个性质。很好算，单独剥离i个人出来，其余特产分给剩下的n-i个人即可。看似受到特产种类的限制，其实算每个种类都分别插板法一下然后乘起来的方案数就可以了。也就是：

$\alpha[i]=C(n,i)*\sum_{j=0}^{m}C_{a[j]-(n-i-1)}^{n-i-1}$ 。

所以所求即为 $\beta(0)$ ，即没有一个人没有分到特产。套容斥定理即可。

上代码——

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#define maxn 2005
using namespace std;
typedef long long ll;
const int mod = 1e9 + 7, mx = 2000;
int read() {
	int x = 0, f = 1, ch = getchar();
	while(!isdigit(ch)) {if(ch == '-') f = -1; ch = getchar();}
	while(isdigit(ch)) x = (x << 1) + (x << 3) + ch - '0', ch = getchar();
	return x * f;
}

int n, m, a[maxn];
ll C[maxn][maxn], alpha[maxn], ans = 0;
signed main() {
	for(int i = 0; i <= mx; i++) C[i][0] = 1;
	for(int i = 1; i <= mx; i++) for(int j = 1; j <= i; j++) 
        C[i][j] = (C[i - 1][j - 1] + C[i - 1][j]) % mod;
	
	n = read(), m = read();//n个学生，m种食物
	for(int i = 1; i <= m; i++) a[i] = read();
	for(int i = 0; i < n; i++) {//至少i个人没有分到食物。不存在都没有食物的情况
		alpha[i] = 1;
		register int k = n - i;
		for(int j = 1; j <= m; j++) alpha[i] = alpha[i] * C[a[j] + k - 1][k - 1] % mod;
        alpha[i] = alpha[i] * C[n][i] % mod;
	}

	for(int i = 0, kd = 1; i < n; kd = -kd, i++) 
        ans = (ans + kd * alpha[i] % mod + mod) % mod;
	printf("%lld\n", ans);
	return 0;	
}

迎评：）
——End——