题目链接

题目大意:

对于集合 \(\{1,2,\dots,n\}\) ,求它的子集族中,有多少个满足:

  • 任意两个子集互不相同;
  • \(1,2,\dots,n\) 都在其中至少出现了 \(2\) 次。

答案对 \(M\) 取模。

看到这种东西就要想到容斥。

设 \(F_i\) 表示至少有 \(i\) 个数字只出现了一次。

更具体的,就是 \(F_i\) 个数只出现一次,其他的数出现次数随便。

由容斥我们可以知道:

\[Ans = \sum_{i=0}^n(-1)^iF_i
\]

我们来考虑这时 \(n-i\) 个随便的部分构成的方案数。

首先,这 \(n - i\) 个数随便定是否只出现了一次,可以看出方案数是 \(2^{n-i}\) 。

然后对于每一个得出的情况,我们都可以选或不选,所以就是 \(2^{2^{n-i}}\) 。

再来看从 \(n\) 个数选出 \(i\) 个一定出现一次的方案数 \(C_n^i\) 很简单。

如果我们设 \(f_i\) 表示恰好有 \(i\) 个数只出现了一次,那么可以得到:

\[F_i=f_i\times2^{2^{n-i}}\times C_n^i
\]

现在着手考虑 \(f_i\) 是怎么得到的。

我们设 \(g_{i,j}\) 表示 \(i\) 个只出现一次的数放到 \(j\) 个集合的方案数。

  • 当此时第 \(j\) 个集合没有不合法的数,此时 \(i\) 只能填在此处 \(\rightarrow g_{i-1,j-1}\)
  • 意味着 \(j\) 个集合里面都有不合法的数字了,这样的话第 \(i\) 个不合法的数可以选择加入到 \(j\) 个集合中任意一个,也可以不加入任何集合。(因为不是强制的) \(\rightarrow (j+1)g_{i-1,j}\)

非常显然的 \(g_{i,j}\) 自然是这两种情况的和。

\[g_{i,j}=g_{i-1, j-1}+(j+1)g_{i-1,j}
\]

最后我们枚举有几个集合里有非法的元素就可以了。

注意:除了不合法的数字一个集合中还可以有其他的元素,及 \((2^{n-i})^j\)

\[f_i=\sum_{j=0}^ig_{i,j}\times (2^{n-i})^j\\
Ans = \sum_{i=0}^n(-1)^i\times\sum_{j=0}^ig_{i,j}\times (2^{n-i})^j\times2^{2^{n-i}}\times C_n^i
\]

呃呃呃,这个 \(Ans\) 的式子可能有一点点乱。。。

Code

#include <cstdio>
#include <iostream>
#include <algorithm> #define file(a) freopen(a".in", "r", stdin), freopen(a".out", "w", stdout) #define Enter putchar('\n')
#define quad putchar(' ') #define int long long
const int N = 3005; int n, mod, fac[N], g[N][N], ans; inline int power(int a, int n, int mod);
inline int C(int n, int m); signed main(void) {
std::cin >> n >> mod;
fac[0] = 1;
for (int i = 1; i <= n; i++)
fac[i] = fac[i - 1] * i % mod;
for (int i = 0; i <= n; i++) {
g[i][0] = 1;
for (int j = 1; j <= i; j++)
g[i][j] = (g[i - 1][j - 1] + (j + 1) * g[i - 1][j] % mod) % mod;
}
for (int i = 0, lala; i <= n; i++) {
int two = power(2, n - i, mod - 1);
two = power(2, two, mod);
int num = power(2, n - i, mod), F = 0, mul = 1;
for (int j = 0; j <= i; j++) {
F = (F + g[i][j] * mul) % mod;
mul = mul * num % mod;
}
if (i % 2 == 1) lala = mod - C(n, i);
else lala = C(n, i);
ans = (ans + F * lala % mod * two % mod) % mod;
ans = (ans % mod + mod) % mod;
}
std::cout << ans << std::endl;
return 0;
} inline int power(int a, int n, int mod) {
int ret = 1;
while (n) {
if (n & 1) ret = ret * a % mod;
a = a * a % mod;
n /= 2;
}
return ret;
}
inline int C(int n, int m) {
if (n < m) return 0;
int ret = fac[n];
ret = ret * power(fac[m], mod - 2, mod) % mod;
ret = ret * power(fac[n - m], mod - 2, mod) % mod;
return ret;
}

[ARC096C] Everything on It 补题记录的更多相关文章

  1. 【补题记录】ZJU-ICPC Summer Training 2020 部分补题记录

    补题地址:https://zjusummer.contest.codeforces.com/ Contents ZJU-ICPC Summer 2020 Contest 1 by Group A Pr ...

  2. 【JOISC 2020 补题记录】

    目录 Day 1 Building 4 Hamburg Steak Sweeping Day 2 Chameleon's Love Making Friends on Joitter is Fun R ...

  3. 【cf补题记录】Codeforces Round #608 (Div. 2)

    比赛传送门 再次改下写博客的格式,以锻炼自己码字能力 A. Suits 题意:有四种材料,第一套西装需要 \(a\).\(d\) 各一件,卖 \(e\) 块:第二套西装需要 \(b\).\(c\).\ ...

  4. 【cf补题记录】Codeforces Round #607 (Div. 2)

    比赛传送门 这里推荐一位dalao的博客-- https://www.cnblogs.com/KisekiPurin2019/ A:字符串 B:贪心 A // https://codeforces.c ...

  5. Codeforces 1214 F G H 补题记录

    翻开以前打的 #583,水平不够场上只过了五题.最近来补一下题,来记录我sb的调试过程. 估计我这个水平现场也过不了,因为前面的题已经zz调了好久-- F:就是给你环上一些点,两两配对求距离最小值. ...

  6. Yahoo Programming Contest 2019 补题记录(DEF)

    D - Ears 题目链接:D - Ears 大意:你在一个\(0-L\)的数轴上行走,从整数格出发,在整数格结束,可以在整数格转弯.每当你经过坐标为\(i-0.5\)的位置时(\(i\)是整数),在 ...

  7. Codeforces 补题记录

    首先总结一下前段时间遇到过的一些有意思的题. Round #474 (Div. 1 + Div. 2, combined)   Problem G 其实关键就是n这个数在排列中的位置. 这样对于一个排 ...

  8. 【补题记录】NOIp-提高/CSP-S 刷题记录

    Intro 众所周知原题没写过是很吃亏的,突然发现自己许多联赛题未补,故开此坑. 在基本补完前会持续更新,希望在 NOIp2020 前填完. 虽然是"联赛题",但不少题目还是富有思 ...

  9. ZJUT11 多校赛补题记录

    牛客第一场 (通过)Integration (https://ac.nowcoder.com/acm/contest/881/B) (未补)Euclidean Distance (https://ac ...

随机推荐

  1. RecyclerView + SQLite 简易备忘录-----中(2)

    (3)RecyclerView的实现 ---中间的内容 RecyclerView是一个比ListView更加强大的滚动控件.要使用这个控件需要先在项目的build.gradle中添加RecyclerV ...

  2. python基础练习题(题目 递归求等差数列)

    day20 --------------------------------------------------------------- 实例028:递归求等差数列 题目 有5个人坐在一起,问第5个 ...

  3. Ubuntu 百度飞桨和 CUDA 的安装

    Ubuntu 百度飞桨 和 CUDA 的安装 1.简介 本文主要是 Ubuntu 百度飞桨 和 CUDA 的安装 系统:Ubuntu 20.04 百度飞桨:2.2 为例 2.百度飞桨安装 访问百度飞桨 ...

  4. XCTF练习题---MISC---Banmabanma

    XCTF练习题---MISC---Banmabanma flag:flag{TENSHINE} 解题步骤: 1.观察题目,打开附件 2.发现是一张斑马图片,这是典型的条形码啊,直接开网站识别 网站:h ...

  5. 3D打印CLI文件格式的读取

    CLI 英文全称:COMMON LAYER INTERFACE 中文名称:通用层接口. 推荐从官方网址看一下其完整的内容: (https://www.hmilch.net/downloads/cli_ ...

  6. pycham的安装、优化、使用

    一.下载与安装 下载地址:https://confluence.jetbrains.com/display/PYH/Previous+PyCharm+Releases 注册:参考地址:http://i ...

  7. django-rest-framework 基础二 序列化器和路由

    django-rest-framework 基础二 序列化器和路由 目录 django-rest-framework 基础二 序列化器和路由 1. 序列化器 1.1 Serializer的使用 1.2 ...

  8. JWT 访问令牌

    JWT 访问令牌 更为详细的介绍jwt 在学习jwt之前我们首先了解一下用户身份验证 1 单一服务器认证模式 一般过程如下: 用户向服务器发送用户名和密码. 验证服务器后,相关数据(如用户名,用户角色 ...

  9. GDB的简单使用一

    GDB的简单使用一 一.概念 二.GDB的基本使用方法一 调试前预备知识 获取进程的内核转储 启动gdb调试 1.启动 2.设置断点 3.运行程序 4.显示栈帧 5.显示变量 6.显示寄存器 7.单步 ...

  10. Java变量, 常量和作用域

    目录 变量 作用域 局部变量 实例变量 类变量 常量 命名规范 视频课程 变量 变量就是可以变化的量 Java是一种强类型的语言, 每个变量都必须声明其类型 Java变量是程序中最基本的存储单元, 其 ...