题目传送门:https://agc016.contest.atcoder.jp/tasks/agc016_f

题目大意:

给定一个\(N\)点\(M\)边的DAG,\(x_i\)有边连向\(y_i\),保证\(x_i<y_i\),原图有\(2^M\)个生成子图,对于每个子图\(G'\),\(A,B\)两人正在玩一个游戏:初始时点1,2上有棋子,每次操作可以把某个棋子沿有向边移动一步,最后不能操作的人为输。问有多少个子图\(G'\)满足先手必胜


这种神题一看就不会写……首先考虑博弈,先手必胜的话当且仅当\(sg[1]!=sg[2]\),这样不好求,我们考虑求\(sg[1]=sg[2]\)的方案

考虑状压dp,记\(f[S]\)表示只考虑\(S\)这个点集,使得\(sg[1]=sg[2]\)的方案数

枚举\(S\)的一个子集\(T\),其补集为\(U\),假设\(U\)集合的\(sg\)值都为0,而\(T\)集合都不为0,,考虑转移:

\(U\)内部的边:一条都不能连

\(U\)到\(T\)的边:随便连

\(T\)到\(U\)的边:\(T\)中的点至少有一条出边

如何保证1,2的\(sg\)相同,保证他们在同一个集合里即可

至少连边和随意连边我们可以预处理出来,一次转移为\(O(n)\),所以总复杂度为\(O(3^n*n)\)

/*program from Wolfycz*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
#define min(x,y) (x<y?x:y)
#define max(x,y) (x>y?x:y)
#define lowbit(x) ((x)&-(x))
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline char gc(){
static char buf[1000000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;
}
inline int frd(){
int x=0,f=1; char ch=gc();
for (;ch<'0'||ch>'9';ch=gc()) if (ch=='-') f=-1;
for (;ch>='0'&&ch<='9';ch=gc()) x=(x<<3)+(x<<1)+ch-'0';
return x*f;
}
inline int read(){
int x=0,f=1; char ch=getchar();
for (;ch<'0'||ch>'9';ch=getchar()) if (ch=='-') f=-1;
for (;ch>='0'&&ch<='9';ch=getchar()) x=(x<<3)+(x<<1)+ch-'0';
return x*f;
}
inline void print(int x){
if (x<0) putchar('-'),x=-x;
if (x>9) print(x/10);
putchar(x%10+'0');
}
const int N=15,p=1e9+7;
bool map[N+10][N+10];
int f[N+10][(1<<N)+10],g[N+10][(1<<N)+10];
int dp[(1<<N)+10],ID[(1<<N)+10];
int mlt(int a,int b){
int res=1;
for (;b;b>>=1,a=1ll*a*a%p) if (b&1) res=1ll*res*a%p;
return res;
}
int main(){
int n=read(),m=read();
for (int i=1;i<=n;i++) ID[1<<(i-1)]=i;
for (int i=1;i<=m;i++){
int x=read(),y=read();
map[x][y]=1;
}
for (int i=1;i<=n;i++)
for (int sta=1;sta<1<<n;sta++)
f[i][sta]=f[i][sta^lowbit(sta)]+map[i][ID[lowbit(sta)]];
for (int i=1;i<=n;i++){
for (int sta=1;sta<1<<n;sta++){
int k=ID[lowbit(sta)];
if (!map[i][k]) g[i][sta]=g[i][sta^lowbit(sta)];
else g[i][sta]=(2ll*g[i][sta^lowbit(sta)]%p+1)%p;
}
}
for (int sta=1;sta<1<<n;sta++){
dp[sta]=1;
for (int sub=(sta-1)&sta;sub;sub=(sub-1)&sta){
if ((sub&1)&&((sta^sub)&2)) continue;
if ((sub&2)&&((sta^sub)&1)) continue;
int res=1,tmp=0;
for (int i=1;i<=n;i++) if (sub&(1<<(i-1))) res=1ll*res*g[i][sta^sub]%p;
for (int i=1;i<=n;i++) if ((sta^sub)&(1<<(i-1))) tmp+=f[i][sub];
res=1ll*res*mlt(2,tmp)%p;
dp[sta]=(dp[sta]+1ll*res*dp[sub])%p;
}
}
printf("%d\n",(mlt(2,m)-dp[(1<<n)-1]+p)%p);
return 0;
}

AtCoder Grand Contest 016 F - Games on DAG的更多相关文章

  1. Atcoder Grand Contest 016 F - Games on DAG(状压 dp)

    洛谷题面传送门 & Atcoder 题面传送门 如何看待 tzc 补他一个月前做的题目的题解 首先根据 SG 定理先手必输当且仅当 \(\text{SG}(1)=\text{SG}(2)\). ...

  2. AtCoder Grand Contest 002 F:Leftmost Ball

    题目传送门:https://agc002.contest.atcoder.jp/tasks/agc002_f 题目翻译 你有\(n*k\)个球,这些球一共有\(n\)种颜色,每种颜色有\(k\)个,然 ...

  3. AtCoder Grand Contest 016 E - Poor Turkeys

    题目传送门:https://agc016.contest.atcoder.jp/tasks/agc016_e 题目大意: 有\(N\)只火鸡,现有\(M\)个人,每个人指定了两只火鸡\(x,y\),每 ...

  4. AtCoder Grand Contest 017 F - Zigzag

    题目传送门:https://agc017.contest.atcoder.jp/tasks/agc017_f 题目大意: 找出\(m\)个长度为\(n\)的二进制数,定义两个二进制数的大小关系如下:若 ...

  5. AtCoder Grand Contest 016 C - +/- Rectangle

    题目传送门:https://agc016.contest.atcoder.jp/tasks/agc016_c 题目大意: 给定整数\(H,W,h,w\),你需要判断是否存在满足如下条件的矩阵,如果存在 ...

  6. AtCoder Grand Contest 003 F - Fraction of Fractal

    题目传送门:https://agc003.contest.atcoder.jp/tasks/agc003_f 题目大意: 给定一个\(H×W\)的黑白网格,保证黑格四连通且至少有一个黑格 定义分形如下 ...

  7. AtCoder Grand Contest 011 F - Train Service Planning

    题目传送门:https://agc011.contest.atcoder.jp/tasks/agc011_f 题目大意: 现有一条铁路,铁路分为\(1\sim n\)个区间和\(0\sim n\)个站 ...

  8. AtCoder Grand Contest 016 B - Colorful Hats

    题目传送门:https://agc016.contest.atcoder.jp/tasks/agc016_b 题目大意: 有\(N\)只猫,每只猫头上带着一个帽子,帽子有颜色,现在告诉你每只猫能看到的 ...

  9. AtCoder Grand Contest 010 F - Tree Game

    题目传送门:https://agc010.contest.atcoder.jp/tasks/agc010_f 题目大意: 给定一棵树,每个节点上有\(a_i\)个石子,某个节点上有一个棋子,两人轮流操 ...

随机推荐

  1. Java Unit Testing - JUnit & TestNG

    转自https://www3.ntu.edu.sg/home/ehchua/programming/java/JavaUnitTesting.html yet another insignifican ...

  2. JAVA学习第十九课(java程序的异常处理 (二))

    异常处理的捕捉形式: 这是能够对异常进行针对性处理的方式 六.try.catch的理解 详细格式: try { //须要被检測异常的代码 } catch(异常类 变量)//改变量用于接受发生异常的对象 ...

  3. 设计模式学习笔记——Observer观察者模式

    观察者模式里面有两个东西:观察者(Observer)和目标(Subject).当目标发生变化的时候,观察者随之起舞,也作出相应的变化.此为观察者模式. 这是怎么做到的?主要是目标里面存有一份观察者的名 ...

  4. 序列化模块 json pickel shelve

    一.json  模块 1.定义 将字典.列表等内容转换成字符串的过程就是序列化. 操作的数据类型有限,但是可以支持所有编程语言操作. 2.为什么要有序列化? 1.以某种存储形式使自定义对象持久化. 2 ...

  5. using the flume-ng-sql-source plugin to push data from Mysql DB to Spark

    org.apache.flume.EventDeliveryException while running flume and sending data to spark · Issue #21 · ...

  6. 李雅普诺夫函数 LyapunovFunction 李雅普诺夫意义下的稳定性

    https://zh.wikipedia.org/zh-hans/李亞普諾夫函數 李雅普诺夫函数(Lyapunov function)是用来证明一动力系统或自治微分方程稳定性的函数.其名称来自俄罗斯数 ...

  7. CH 5402 选课(分组背包+树形DP)

    CH 5402 选课 \(solution:\) 很有讨论套路的一道题,利用树的结构来表示出不同课程之间的包含关系(这里还要建一个虚点将森林变成一颗打大树).然后用子树这个概念巧妙的消除了因为这些包含 ...

  8. field 属性操作

    首先必须明一点 Field类主要是用来辅助获取和操作类的属性的! 1.怎么通过反射获取类的属性 先来看JDK提供的方法有如下几种: a)Class.getDeclaredField(String na ...

  9. spring、spring MVC、spring Boot

    Spring 是一个“引擎” Spring MVC 是基于 Spring 的一个 MVC 框架 Spring Boot 是基于 Spring4 的条件注册的一套快速开发整合包 Spring 最初利用“ ...

  10. Delphi通过Get获取来自PHP的返回值

    Delphi代码 unit Unit1; interface uses Windows, Messages, SysUtils, Variants, Classes, Graphics, Contro ...