AGC 016 F - Games on DAG(状压dp)
题意
给你一个有 \(n\) 个点 \(m\) 条边 DAG 图,点的标号和拓扑序一致。
现在有两个人进行博弈,有两个棋子分别在 \(1, 2\) 号点上,需要不断移动到它指向的点上。
如果当前两个点都无法移动,那么就视为当前操作的人失败。
问有多少边集满足先手必胜。
\(\displaystyle 2 \le n \le 15, m \le \frac{n \times (n+1)}{2}\)
题解
参考了 wxh010910 大佬的博客 。
首先利用博弈的 SG 函数易得,如果 \(1\) 号点和 \(2\) 号点的 SG 值异或不为 \(0\) 则先手必败。
但这样不好计数,我们考虑它的反面,求 \(1,2\) 号点 SG 值异或为 \(0\) 的方案数。
也就是 \(1,2\) 号点 SG 值为 \(0\) 的方案数。
考虑边集显然是不现实的,改成考虑点集,令 \(f(S)\) 为 \(S\) 这个集合满足 \(SG(1) = SG(2)\) 的方案数。
我们考虑枚举它的一个子集 \(T\) 假设满足点集中所有的点 \(SG\) 都不为 \(0\) ,并且它对于 \(S\) 的补集 \(U\) 都满足 \(SG\) 都为 \(0\) 。
这是连边后的情况,连边前 \(T\) 集合会存在 \(SG\) 为 \(0\) 的点。
我们考虑这些点的连边方案数。
- \(U\) 集合内部,显然不能连边。(不然不能保证 SG 为 \(0\) )
- \(U \to T\) 随意连边。
- \(T \to U\) 要保证 \(T\) 中每个点都需要有一条向 \(U\) 连的边(要保证 \(T\) 集合中的 SG 不为 \(0\) )
- \(T\) 内部的方案数,正好就是 \(f(T)\) 。这是为什么呢?把 \(T\) 中每个点都新连向一个 \(SG = 0\) 的点,那么它的 \(SG\) 都会增加 \(1\) 。
然后为了保证 \(SG(1)=SG(2)\) 枚举 \(S, T, U\) 的时候需要保证 \(1, 2\) 同时选或同时不选。
为什么这样是对的呢?
因为我们保证了 \(T\) 集合的方案本身满足 \(SG(1) = SG(2)\) 那么之后同时 \(+1\) 也会满足这个性质。
\(U\) 集合显然也是成立的( \(SG(1) = SG(2) = 0\) )。
然后枚举的时候有细节,我们考虑枚举 \(U\) 为 \(S\) 的子集,因为 \(U\) 必然不为 \(\varnothing\) ,而 \(T\) 可以为 \(\varnothing\) 。
总结
对于一类关于边集的状压 \(dp\) ,可以考虑转化成点集,然后计算连边的方案数。
代码
#include <bits/stdc++.h>
#define For(i, l, r) for(register int i = (l), i##end = (int)(r); i <= i##end; ++i)
#define Fordown(i, r, l) for(register int i = (r), i##end = (int)(l); i >= i##end; --i)
#define Set(a, v) memset(a, v, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define debug(x) cout << #x << ": " << x << endl
#define DEBTG(...) fprintf(stderr, __VA_ARGS__)
using namespace std;
inline bool chkmin(int &a, int b) {return b < a ? a = b, 1 : 0;}
inline bool chkmax(int &a, int b) {return b > a ? a = b, 1 : 0;}
inline int read() {
int x = 0, fh = 1; char ch = getchar();
for (; !isdigit(ch); ch = getchar()) if (ch == '-') fh = -1;
for (; isdigit(ch); ch = getchar()) x = (x << 1) + (x << 3) + (ch ^ 48);
return x * fh;
}
void File() {
#ifdef zjp_shadow
freopen ("F.in", "r", stdin);
freopen ("F.out", "w", stdout);
#endif
}
const int N = 15, Mod = 1e9 + 7;
int n, m, G[N], Table[1 << N], dp[1 << N];
int main () {
File();
n = read(); m = read();
For (i, 1, m) { int u = read() - 1, v = read() - 1; G[u] |= (1 << v); }
int maxsta = (1 << n) - 1; dp[0] = 1;
For (i, 0, maxsta) Table[i] = 1 << __builtin_popcount(i);
For (S, 2, maxsta) if ((S & 1) == ((S >> 1) & 1))
for (int U = S; U; U = (U - 1) & S) if ((U & 1) == ((U >> 1) & 1)) {
int Ways = 1, T = S ^ U;
For (p, 0, n - 1) if ((S >> p) & 1) {
if ((U >> p) & 1) Ways = 1ll * Table[G[p] & T] * Ways % Mod;
else Ways = 1ll * (Table[G[p] & U] - 1) * Ways % Mod;
}
dp[S] = (dp[S] + 1ll * Ways * dp[T]) % Mod;
}
int ans = 1; For (i, 1, m) ans = (ans << 1) % Mod;
printf ("%d\n", (ans + Mod - dp[maxsta]) % Mod);
return 0;
}
AGC 016 F - Games on DAG(状压dp)的更多相关文章
- AtCoder Grand Contest 016 F - Games on DAG
题目传送门:https://agc016.contest.atcoder.jp/tasks/agc016_f 题目大意: 给定一个\(N\)点\(M\)边的DAG,\(x_i\)有边连向\(y_i\) ...
- Codeforces Round #531 (Div. 3) F. Elongated Matrix(状压DP)
F. Elongated Matrix 题目链接:https://codeforces.com/contest/1102/problem/F 题意: 给出一个n*m的矩阵,现在可以随意交换任意的两行, ...
- 2019牛客多校第五场 F maximum clique 1 状压dp+最大独立集
maximum clique 1 题意 给出一个集合s,求每个子集的最大独立集的权值和(权值是独立集的点个数) 分析 n比较小,一股浓浓的暴力枚举每一个子集的感觉,但是暴力枚举模拟肯定会T,那么想一想 ...
- Atcoder Grand Contest 016 F - Games on DAG(状压 dp)
洛谷题面传送门 & Atcoder 题面传送门 如何看待 tzc 补他一个月前做的题目的题解 首先根据 SG 定理先手必输当且仅当 \(\text{SG}(1)=\text{SG}(2)\). ...
- DAG求最短路--TSP变形--状压dp
DAG状压dp的一种 题目: $m$个城市,$n$张车票,第i张车票上的时间是$t_i$, 求从$a$到$b$的最短时间,如果无法到达则输出“impossible” 解法: 考虑状态:“现在在城市$v ...
- 【XSY3042】石像 拓扑排序 状压DP 洲阁筛
题目大意 有 \(n\) 个整数 \(a_1,a_2,\ldots,a_n\),每个数的范围是 \([1,m]\).还有 \(k\) 个限制,每个限制 \(x_i,y_i\) 表示 \(a_{x_i} ...
- 【状压DP】poj2686 Traveling by Stagecoach
状压DP裸题,将({当前车票集合},当前顶点)这样一个二元组当成状态,然后 边权/马匹 当成边长,跑最短路或者DAG上的DP即可. #include<cstdio> #include< ...
- 状压dp专题复习
状压dp专题复习 (有些题过于水,我直接跳了) 技巧总结 : 1.矩阵状压上一行的选择情况 \(n * 2^n\) D [BZOJ2734][HNOI2012]集合选数 蒻得不行的我觉得这是一道比较难 ...
- BZOJ 1087: [SCOI2005]互不侵犯King [状压DP]
1087: [SCOI2005]互不侵犯King Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 3336 Solved: 1936[Submit][ ...
随机推荐
- OpenStack报错:MessagingTimeout: Timed out waiting for a reply to message ID
L3.agent中出现大量消息超时错误,对网络的操作各种异常. 报错如下: -- :: ERROR neutron.agent.l3.agent [req-db9207e6--4f23-8c19-0d ...
- CF954I Yet Another String Matching Problem 并查集、FFT
传送门 题意:给出两个由小写$a$到$f$组成的字符串$S$和$T$($|S| \geq |T|$),给出变换$c1\,c2$表示将两个字符串中所有$c1$字符变为$c2$,求$S$的每一个长度为$T ...
- socket、tcp、udp、http 的认识及区别
一.先来一个讲TCP.UDP和HTTP关系的 1.TCP/IP是个协议组,可分为三个层次:网络层.传输层和应用层. 在网络层有IP协议.ICMP协议.ARP协议.RARP协议和BOOTP协议. 在传输 ...
- BodeAbp服务端介绍
BodeAbp服务端只提供api,绝大部分api通过abp的动态WebApi机制提供,原理可以参考这篇文章:http://www.cnblogs.com/1zhk/p/5418694.html 与业务 ...
- checkpoint-BLCR部署和测试(源码)
1. 概述2. 部署过程2.1 源码下载2.2 解压安装2.3 添加库环境2.4 插入内核模块3. 测试3.1 创建测试程序3.2 功能测试4. 参考博客 1. 概述 checkpoint 2. 部署 ...
- CSS文本实例
CSS 文本属性可定义文本的外观. 通过文本属性,您可以改变文本的颜色.字符间距,对齐文本,装饰文本,对文本进行缩进,等等.#############################CSS 文本属性属 ...
- 常用rsync命令操作梳理
作为一个运维工程师,经常可能会面对几十台.几百台甚至上千台服务器,除了批量操作外,环境同步.数据同步也是必不可少的技能.说到“同步”,不得不提的利器就是rsync.rsync不但可以在本机进行文件同步 ...
- 2017-2018-2 1723《程序设计与数据结构》第八周作业 & 实验二 & 第一周结对编程 总结
作业地址 第八周作业:https://edu.cnblogs.com/campus/besti/CS-IMIS-1723/homework/1847 (作业界面已评分,可随时查看,如果对自己的评分有意 ...
- Beta版会议总结
开会时间:2015年06月12日 开会地点:基教601 开会人员:李想,王颖瑞,朱少辉,陈晨,侯涛亮. 开会内容:对于6月10日,大一同学的投票情况进行讨论和反思. 讨论结果如下: 一.目前存在的问题 ...
- logstash 解析日志文件
input { file { path => "/usr/local/test/log.log" } } filter { grok { match => { &quo ...