【洛谷】[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段的方案数,应该为
除空集以外,共有个子集,在其中选个
如果最后一段是空集,则方案数就是f[i-1]
如果最后一段重复了,则在f[i-2]的基础上,在选出一个与前i-2个不重复的子集,有,
还要注意,这个与最后一段重复的子集有个位置可以放
故转移方程为
最后乘上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;
}
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列01:轻松3步本地部署deepseek,普通电脑可用
· 25岁的心里话
· 按钮权限的设计及实现