2020 动物园 · 方案一 子集枚举
逐位预处理 + 枚举 2^k 子集计数
正文
// CSP-S 2020 复赛 T2 · 动物园(方案一:预处理 + 2^k 子集计数,k ≤ 20)
// 原题:https://oj.yecheng.tv/p/CSPS2020B
// 题意:n 只动物各有 k 位特征掩码;m 条规则 (p, q):动物若有第 p 位
// 特征,就必须买 q 号饲料。已买饲料 = 原有动物触发规则的并集。
// 合法新动物 S:S 中每位触发的饲料都已买。求合法 S 总数 − n。
// 思路(逐位预处理 + 子集枚举):
// 1. appear = 所有原有动物掩码的 OR;已买饲料 = appear 中位触发的 q。
// 2. 对每个特征位 b(0..k-1):b 是否"可自由使用"?——若 b ∈ appear,
// 它触发的饲料必然已买 → 安全;若 b ∉ appear,检查它触发的每条
// 规则的 q 是否都在已买集合中(q 用排序数组 + 二分查找)。
// 不可用的位记 bad[b] = true。
// 3. 合法 S = 只由"安全位"组成的掩码 → 个数 = 2^(安全位数)。
// 4. 答案 = 2^(安全位数) − n。
// 本方案用 2^k 枚举来"数"合法 S(k ≤ 20 才可行),逻辑与方案二
// 完全一致,方便对照理解;方案二用公式直接算。
// 复杂度:O(n + m·d + 2^k)(d 为每位的平均规则数)。
// 易错点:
// 1. 饲料编号 q ≤ 1e8,不能位压,用数组 + 排序 + 二分;
// 2. appear 中的 1 位永远安全,别误标 bad;
// 3. 答案与总数用 unsigned long long。
#include <cstdio>
#include <vector>
#include <algorithm>
using namespace std;
int n, m, k, c;
unsigned long long a[1000005];
vector<unsigned long long> feedOf[65]; // 位 b 触发的饲料列表
unsigned long long bought[1000005];
int boughtN = 0;
int main() {
freopen("zoo.in", "r", stdin);
freopen("zoo.out", "w", stdout);
scanf("%d%d%d%d", &n, &m, &k, &c);
unsigned long long appear = 0;
for (int i = 0; i < n; i++) {
scanf("%llu", &a[i]);
appear |= a[i];
}
for (int j = 0; j < m; j++) {
unsigned long long p, q;
scanf("%llu%llu", &p, &q);
feedOf[p].push_back(q); // 题面保证 p < k
}
// 已买饲料:appear 中位触发的 q
for (int b = 0; b < k; b++) {
if (!((appear >> b) & 1ULL)) continue;
for (unsigned long long q : feedOf[b])
bought[boughtN++] = q;
}
sort(bought, bought + boughtN);
boughtN = unique(bought, bought + boughtN) - bought;
// 逐位判安全
int safeCnt = 0;
for (int b = 0; b < k; b++) {
if ((appear >> b) & 1ULL) { safeCnt++; continue; }
bool ok = true;
for (unsigned long long q : feedOf[b])
if (!binary_search(bought, bought + boughtN, q)) { ok = false; break; }
if (ok) safeCnt++;
}
// 2^k 枚举计数(与公式 2^safeCnt 对照)
long long total = 0;
for (unsigned long long S = 0; S < (1ULL << k); S++) {
bool ok = true;
for (int b = 0; b < k && ok; b++)
if ((S >> b) & 1ULL) {
if (!((appear >> b) & 1ULL)) {
for (unsigned long long q : feedOf[b])
if (!binary_search(bought, bought + boughtN, q)) { ok = false; break; }
}
}
if (ok) total++;
}
printf("%lld\n", total - (long long)n);
return 0;
}
CSP-S 标程 · 复赛真题的其它内容
- 2019 格雷码 · 方案一 递归构造
- 2019 格雷码 · 方案二 异或公式
- 2020 儒略日 · 方案一 逐天模拟
- 2020 儒略日 · 方案二 分段整块跳
- 2021 廊桥分配 · 方案一 枚举分配数模拟
- 2021 廊桥分配 · 方案二 预处理归属加前缀和
- 2022 假期计划 · 方案一 BFS 加平方枚举
- 2022 假期计划 · 方案二 预处理最佳中转
- 2023 密码锁 · 方案一 全域枚举
- 2023 密码锁 · 方案二 基准候选收敛
- 2024 决斗 · 方案一 排序贪心模拟
- 2024 决斗 · 方案二 桶计数线性扫描
- 2025 社团招新 · 方案一 状态计数DP
- 2025 社团招新 · 方案二 超额排序移人
- 2019 括号树 · 方案一 逐点重算
- 2019 括号树 · 方案二 栈加 DFS 递推
- 2020 动物园 · 方案二 位或统计加计数公式
- 2021 括号序列 · 方案一 立方区间 DP
- 2021 括号序列 · 方案二 平方递推
- 2022 策略游戏 · 方案一 暴力扫描
- 2022 策略游戏 · 方案二 ST 表区间极值
- 2023 消消乐 · 方案一 枚举区间加栈
- 2023 消消乐 · 方案二 记忆化递归
- 2024 超速检测 · 方案一 暴力判定
- 2024 超速检测 · 方案二 区间转化加贪心选点
- 2025 道路修复 · 方案一 逐子集重建MST
- 2025 道路修复 · 方案二 预筛MST全局排序
- 2019 树上的数 · 方案一 全排列暴力
- 2019 树上的数 · 方案二 贪心定序加时刻链
- 2020 函数调用 · 方案一 直接模拟
- 2020 函数调用 · 方案二 拓扑序乘子回推
- 2021 回文 · 方案一 环形配对逆向
- 2021 回文 · 方案二 位置表逆向构造
- 2022 星战 · 方案一 重建判定
- 2022 星战 · 方案二 出度计数维护
- 2023 结构体 · 方案一 顺序模拟
- 2023 结构体 · 方案二 统一类型表封装
- 2024 染色 · 方案一 平方 DP
- 2024 染色 · 方案二 last 指针线性 DP
- 2025 谐音替换 · 方案一 逐对逐位置暴力
- 2025 谐音替换 · 方案二 叠串+AC自动机
- 2019 Emiya 家今天的饭 · 方案一 逐列容斥 DP
- 2019 Emiya 家今天的饭 · 方案二 状态折叠
- 2020 贪吃蛇 · 方案一 multiset 模拟
- 2020 贪吃蛇 · 方案二 双端队列停时规律
- 2021 交通规划 · 方案一 Dinic 最小割
- 2021 交通规划 · 方案二 对偶图最短路
- 2022 数据传输 · 方案一 k 为 1 前缀和
- 2022 数据传输 · 方案二 倍增加矩阵
- 2023 种树 · 方案一 按深度贪心
- 2023 种树 · 方案二 堆加合并贪心
- 2024 擂台游戏 · 方案一 逐 K 模拟
- 2024 擂台游戏 · 方案二 倍增分层预处理
- 2025 员工招聘 · 方案一 集合记忆化搜索
- 2025 员工招聘 · 方案二 三维计数DP