BZOJ5006 THUWC2017随机二分图（概率期望+状压dp）

　　下称0类为单边，1类为互生边，2类为互斥边。对于一种匹配方案，考虑其出现的概率*2ⁿ后对答案的贡献，初始为1，如果有互斥边显然变为0，否则每有一对互生边其贡献*2。于是有一个显然的dp，即设f[S1][S2]为左边选取S1右边选取S2对答案的贡献。转移时考虑S1中编号最小的点x与右边的点y匹配。首先将f[S1-(1<<x)][S2-(1<<y)]统计进去。然后若(x,y)是单边，或者虽存在互生互斥关系，但其对应边的左端点还不在S1中或就是x，或右端点还不在S2中或就是y，就不管了；否则若互斥将f[S1-(1<<x)-(1<<x')][S2-(1<<y)-(1<<y')]减掉，若互生将f[S1-(1<<x)-(1<<x')][S2-(1<<y)-(1<<y')]加上，其中(x',y')是(x,y)的对应边。这样大概就是C(30,15)*15的。

　　但是这只有暴力20分，甚至连空间都开不下。然而满分做法和他是一个复杂度的。对上面的做法改为记忆化搜索，map存储状态就可以了。多了log也多了80分。不是非常理解意义何在。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;
#define ll long long
#define N 15
#define P 1000000007
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<''||c>'')) c=getchar();return c;}
int gcd(int n,int m){return m==?n:gcd(m,n%m);}
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
int n,m,typ[N][N],match[N][N][],s[<<N];
map<int,int> f;
void inc(int &x,int y){x+=y;if (x>=P) x-=P;}
inline int trans(int n,int m){return n<<N|m;}
int solve(int i,int j)
{
    if (f.find(trans(i,j))!=f.end()) return f[trans(i,j)];
    int x=i&-i,ans=;
    for (int t=j,k=t&-t;t;t^=k,k=t&-t)
    if ((k&j)&&typ[s[x]][s[k]]>=)
    {
        inc(ans,solve(i^x,j^k));
        int u=match[s[x]][s[k]][],v=match[s[x]][s[k]][];
        if ((u&i)&&(v&j)&&x!=u&&k!=v)
        {
            if (typ[s[x]][s[k]]==) inc(ans,solve(i^x^u,j^k^v));
            if (typ[s[x]][s[k]]==) inc(ans,P-solve(i^x^u,j^k^v));
        }
    }
    f[trans(i,j)]=ans;return ans;
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj5006.in","r",stdin);
    freopen("bzoj5006.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),m=read();
    memset(typ,,sizeof(typ));
    for (int i=;i<=m;i++)
    {
        int op=read(),x=read()-,y=read()-;
        if (op==) typ[x][y]=;
        else
        {
            int p=read()-,q=read()-;
            typ[x][y]=typ[p][q]=op;
            match[x][y][]=<<p,match[x][y][]=<<q;
            match[p][q][]=<<x,match[p][q][]=<<y;
        }
    }
    for (int i=;i<n;i++) s[<<i]=i;
    f[]=;
    cout<<solve((<<n)-,(<<n)-);
    return ;
}