[bzoj5006][THUWC2017]Bipartite 随机二分图——概率+状压DP

题面：

这里

思路：

先考虑只有t=0的时候该怎么做。
即枚举每一条边选不选入集合，枚举出所有的情况后对于每一种情况计算概率和二分图匹配数的乘积。
上述做法太暴力，于是考虑优化，不难发现一个完备匹配一定是n条匹配边构成并且不重复地覆盖了所有的点。于是设计一个类似于背包的DP，每一条边都是一种物品，依次加入每一种已经匹配好的点集：
如果点集中存在这条边的端点，那说明无论加不加入都不可能从这个点集转移出去，即这个点集在这一步的加边中可以选择加与不加，则贡献为$100\%$。
如果点集中不存在这条边的端点，则说明可以选择加入这条边来获得一个新的匹配，于是这个点集的贡献对于新的点集的贡献为$50\%$。
对于每一组边都重复上述的抉择，则可以得到最后的完备匹配出现的概率和=期望。
于是对于t=1和t=2的边组，也可以采用上述的做法去转移：
对于t=1的，两条边可以同时选择的时候，则新添加的匹配边有三种情况，即添加(a1,b1)，添加(a2,b2)，添加(a1,b2)和(a2,b2)。
如果只有一种可以添加，则只转移一种。
好的，按照上面的做法，在THUWC时可以得到20分的高分。
好像即是时这样复杂度还是无法接受，发现复杂度大部分花在枚举每一种物品，可对于这一种点集却又没有办法转移的过程上面。
于是不用背包的那种转移方式，即直接从0到all枚举匹配的点集并且依次加入两端都不在点集中的边转移即可，为了保证不重复转移，每次只转移一个标号最小的边。
等等，边组怎么保证不重复转移？好像很难处理的样子。
于是这题吊炸天的地方来了，把每个边组拆成互相独立的两条边，并且概率都是$50\%$，对于两条边都可以选择并且只转移一条边的情况当然时满足的，但是如何两条边都可以选择并且两条边都匹配，那么两条边一起加入边集的概率就变成了$25\%$,于是在第一条边独自加入匹配集合的时候便判断一下是否可以继续转移到两条边都加入集合，如果是这样那么提前对两条边都加入集合的状态加上那少算的$25\%$的贡献。
两条边只可以选一条的情况类似，多算了$25\%$便多减去$25\%$的多余贡献。

/**************************************************************
    Problem: 5006
    User: ylsoi
    Language: C++
    Result: Accepted
    Time:8208 ms
    Memory:27216 kb
****************************************************************/
 
#include<bits/stdc++.h>
 
#define REP(i,a,b) for(int i=a,i##_end_=b;i<=i##_end_;++i)
typedef long long ll;
 
using namespace std;
 
void File(){
    freopen("loj2290.in","r",stdin);
    freopen("loj2290.out","w",stdout);
}
 
const int maxn=15+10;
const int maxm=15*15+10;
const int maxw=(1<<17)+10;
const ll mod=1e9+7;
int n,m,all;
int beg[maxn],fr[maxm<<1],to[maxm<<1],las[maxm<<1],ex[maxm<<1],cnte=1;
ll w[maxm<<1],inv2=500000004,ans;
map<int,ll>dp[maxw];
map<int,ll>::iterator it;
 
void add(int u,int v){las[++cnte]=beg[u];beg[u]=cnte;to[cnte]=v;fr[cnte]=u;}
 
void init(){
    scanf("%d%d",&n,&m);
    int ty,a1,b1,a2,b2;
    REP(i,1,m){
        scanf("%d",&ty);
        if(ty==0){
            scanf("%d%d",&a1,&b1);
            add(a1,b1);
            ex[cnte]=-1;
        }
        else{
            scanf("%d%d%d%d",&a1,&b1,&a2,&b2);
            add(a1,b1); add(a2,b2);
            ex[cnte]=cnte-1;
            ex[cnte-1]=cnte;
            if(ty==1)w[cnte]=w[cnte-1]=inv2*inv2%mod;
            else w[cnte]=w[cnte-1]=-inv2*inv2%mod;
        }
    }
}
 
void add(ll &_,ll __){_=_+__;if(_>mod)_-=mod;if(_<0)_+=mod;}
 
void work(){
    dp[0][0]=1;
    all=(1<<n)-1;
    REP(S,0,all){
        REP(i,1,n){
            if((1<<(i-1))&S)continue;
            for(it=dp[S].begin();it!=dp[S].end();++it){
                int T=it->first; ll val=it->second;
                for(int j=beg[i];j;j=las[j]){
                    if((1<<(to[j]-1))&T)continue;
                    add(dp[S|(1<<(i-1))][T|(1<<(to[j]-1))],val*inv2%mod);
                    if(ex[j]==-1)continue;
                    int x=fr[ex[j]],y=to[ex[j]];
                    if(x==i || y==to[j])continue;
                    x=1<<(x-1); y=1<<(y-1);
                    if((x & S) || (y & T))continue;
                    add(dp[S|(1<<(i-1))|x][T|(1<<(to[j]-1))|y],val*w[j]%mod);
                }
            }
            break;
        }
    }
    ans=dp[all][all];
    REP(i,1,n)ans=ans*2%mod;
    printf("%lld\n",ans);
}
 
int main(){
    //File();
    init();
    work();
    return 0;
}