AtCoder Grand Contest 016 F - Games on DAG SG函数+状压dp

题意

给出一个n个点m条边的有向无环图G，可以删掉任意的边得到一个新的图G’，显然这样的图有$2^m$种。现在两个人在G’上玩游戏：初始在点1和点2上分别有一个标记。每次一个人可以把某个标记沿着某条出边移动。不能移动者输。问在这$2^m$个图里面有多少个满足先手必胜。
n<=15

分析

看到这个数据范围就知道是状压了。。。
首先不难发现若先手必胜则点1和点2的SG值不能相同。正难则反，我们考虑如何求SG值相同时的答案。
设f[S]表示只看S里面的点时有多少种情况使得1和2的SG值相同，其中S中要么a和b都在，要么都不在。
转移考虑枚举S的一个子集A，设A为所有SG值为0的点，其补集B表示所有SG值不为0的点。
那么A中的点两两不能连边，A到B可以乱连，B的每个点至少要有一条出边指向A，然后A中的方案是f[A]。
可以通过预处理一些东西使得总的复杂度达到$O(n*3^n)$

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;

const int MOD=1000000007;

int n,m,bin[205],g[20][40005],f[40005],q[40005];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

int lowbit(int x) {return x&(-x);}

int main()
{
    n=read();m=read();
    bin[0]=1;
    for (int i=1;i<=max(n,m);i++) bin[i]=bin[i-1]*2%MOD;
    for (int i=1;i<=m;i++)
    {
        int x=read()-1,y=read()-1;
        g[x][bin[y]]=1;
    }
    for (int i=0;i<n;i++)
        for (int j=1;j<bin[n];j++)
            g[i][j]=g[i][j-lowbit(j)]+g[i][lowbit(j)];
    for (int i=1;i<bin[n];i++)
    {
        if ((i&1)!=(i&2)/2) continue;
        f[i]=1;int h=0;
        for (int s=i;s;s=(s-1)&i) q[++h]=s;
        while (h--)
        {
            int s1=q[h+1],s2=s1^i,w=f[s1];
            if ((s1&1)!=(s1&2)/2) continue;
            for (int j=0;j<n;j++)
                if (s1&bin[j]) w=(LL)w*(bin[g[j][s2]]-1)%MOD;
                else if (s2&bin[j]) w=(LL)w*bin[g[j][s1]]%MOD;
            f[i]+=w;f[i]-=f[i]>=MOD?MOD:0;
        }
    }
    printf("%d",(bin[m]+MOD-f[bin[n]-1])%MOD);
    return 0;
}