Solution -「AGC 016F」Games on DAG

(mathcal{Description})

  Link.

  给定一个含 (n) 个点 (m) 条边的 DAG，有两枚初始在 1 号点和 2 号点的棋子。两人博弈，轮流移动其中一枚棋子到邻接结点位置，无法移动者负。求 (2^m) 个边集中，加入图中能使先手必胜的方案数。答案对 (10^9+7) 取模。

  (nle15)。

(mathcal{Solution})

  先从博弈角度思考：两枚棋子显然独立，那么先手必胜等价于 (operatorname{sg}(1) ot=operatorname{sg}(2))，其中，(operatorname{sg}) 有

[operatorname{sg}(u)=operatorname{mex}_{vinoperatorname{adj}(u)}operatorname{sg}(v) ]

在不同边集的情况中，(operatorname{adj}(u)) 有所改变，所以情况不同。

  接下来，方便起见，尝试计数 (operatorname{sg}(1)=operatorname{sg}(2)) 的情况。令 (f(S)~(1,2in S)) 表示仅考虑点集 (S) 的导出子图时，(operatorname{sg}(1)=operatorname{sg}(2)) 的方案数。而 DP 中的“子问题”源于 (operatorname{sg}) 函数中 (operatorname{mex}) 的定义——如果我们挑除所有 (operatorname{sg}(u)=0) 的 (u)（一定存在）及其邻接边，剩下的子图仍然是同类问题，只是所有 (operatorname{sg}) 值减少 (1)，不影响 DP 定义，这就是转移的雏形。具体地，枚举点集 (Tsubseteq S)，钦定其中所有点的 (operatorname{sg}) 非 (0)，考虑：

• (1,2in T)，(T) 中每个点至少连向一个 (S-T) 中的点，其余边任意（注意 (T) 中的边是子问题，不考虑）。

• (1,2 otin T)，已经满足条件，在上种情况的基础上加上 (T) 中边任选的方案，直接贡献。

  综上，复杂度 (mathcal O(3^nn))。

(mathcal{Code})

/* Clearink */

#include <cstdio>

#define rep( i, l, r ) for ( int i = l, repEnd##i = r; i <= repEnd##i; ++i )
#define per( i, r, l ) for ( int i = r, repEnd##i = l; i >= repEnd##i; --i )

const int MAXN = 15, MAXM = MAXN * ( MAXN - 1 ) >> 1, MOD = 1e9 + 7;
int n, m, pwr[MAXM + 5], conc[1 << MAXN][MAXN], f[1 << MAXN];
bool adj[MAXN + 5][MAXN + 5];

inline int mul( const long long a, const int b ) { return a * b % MOD; }
inline int sub( int a, const int b ) { return ( a -= b ) < 0 ? a + MOD : a; }
inline void subeq( int& a, const int b ) { ( a -= b ) < 0 && ( a += MOD ); }
inline int add( int a, const int b ) { return ( a += b ) < MOD ? a : a - MOD; }
inline void addeq( int& a, const int b ) { ( a += b ) >= MOD && ( a -= MOD ); }

int main() {
	scanf( "%d %d", &n, &m ), pwr[0] = 1;
	rep ( i, 1, m ) {
		pwr[i] = add( pwr[i - 1], pwr[i - 1] );
		int u, v; scanf( "%d %d", &u, &v );
		adj[u - 1][v - 1] = true;
	}

	rep ( S, 1, ( 1 << n ) - 1 ) {
		int v = 0; for ( ; !( S >> v & 1 ); ++v );
		rep ( u, 0, n - 1 ) conc[S][u] = conc[S ^ 1 << v][u] + adj[u][v];
	}

	rep ( S, 0, ( 1 << n ) - 1 ) if ( ( S & 3 ) == 3 ) {
		f[S] = 1;
		for ( int T = S & ( S - 1 ); T; T = ( T - 1 ) & S ) {
			if ( T & 1 && T & 2 ) {
				int coe = 1;
				rep ( u, 0, n - 1 ) if ( S >> u & 1 ) {
					if ( T >> u & 1 ) {
						coe = mul( coe, pwr[conc[S ^ T][u]] - 1 );
					} else {
						coe = mul( coe, pwr[conc[T][u]] );
					}
				}
				addeq( f[S], mul( coe, f[T] ) );
			} else if ( !( T & 1 || T & 2 ) ) {
				int coe = 1;
				rep ( u, 0, n - 1 ) if ( S >> u & 1 ) {
					if ( T >> u & 1 ) {
						coe = mul( coe,
							mul( pwr[conc[S ^ T][u]] - 1, pwr[conc[T][u]] ) );
					} else {
						coe = mul( coe, pwr[conc[T][u]] );
					}
				}
				addeq( f[S], coe );
			}
		}
	}

	printf( "%d
", sub( pwr[m], f[( 1 << n ) - 1] ) );
	return 0;
}