AtCoder Grand Contest 016 F

题目传送门：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
",(mlt(2,m)-dp[(1<<n)-1]+p)%p);
	return 0;
}