【洛谷】[HNOI2011] 卡农

题目链接:https://www.luogu.com.cn/problem/P3214

题目大义:

给定集合S={1,2,...,n},选出m个子集,满足:
(1)所有选出的m个子集都不能为空。
(2)所有选出的m个子集中,不能存在两个完全一样的集合。
(3)所有选出的m个子集中,1到n每个元素出现的次数必须是偶数。

先不管重复,我们只需要在最后除以m!
记f[i]表示分成i段合法的方案数
对于限制3,其实在分好前i-1段时,最后一段就自动分好了,但是最后一段可能是空集或者跟前面重复,所以要用容斥。
对于前i-1段的方案数,应该为A2n1i1

除空集以外,共有2n1个子集,在其中选i1

如果最后一段是空集,则方案数就是f[i-1]
如果最后一段重复了,则在f[i-2]的基础上,在选出一个与前i-2个不重复的子集,有2n1(i2)=2ni+1
还要注意,这个与最后一段重复的子集有i1个位置可以放
故转移方程为f[i]=A2n1i1f[i1]f[i2]×(i1)×(2ni+1)
最后乘上m!的逆元即可

#include<bits/stdc++.h>
using namespace std;
const int N = 1e6 + 5, P = 1e8 + 7;
int n, m, p, fac[N], f[N];
inline int qpow(int a, int b){
	int ret = 1;
	for (; b; a = 1ll * a * a % P, b >>= 1)
		if (b & 1) ret = 1ll * ret * a % P;
	return ret;
}
int main(){
	scanf("%d%d", &n, &m);
	p = qpow(2, n);
	fac[0] = 1;
	for (int i = 1; i < m; ++i)
		fac[i] = 1ll * fac[i - 1] * (p - i + P) % P;
	
	f[0] = 1, f[1] = 0;
	for (int i = 2; i <= m; ++i)
		f[i] = ((fac[i - 1] - f[i - 1] - 1ll * f[i - 2] * (i - 1) % P * (p - i + 1) % P) % P + P) % P;
	
	int inv = 1;
	for (int i = 1; i <= m; ++i) inv = 1ll * inv * i % P;
	inv = qpow(inv, P - 2);
	printf("%d\n", 1ll * f[m] * inv % P);
	return 0;
}
posted @   浮游半笙  阅读(26)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列01:轻松3步本地部署deepseek,普通电脑可用
· 25岁的心里话
· 按钮权限的设计及实现
点击右上角即可分享
微信分享提示