[THUWC2017]随机二分图

题目大意

给一张二分图，有左部点和右部点。

有三种边，第一种是直接从左部点连向右部点，出现概率为50%。

第二种边一组里有两条边，这两条边同时出现或者不出现，概率都是50%。

第三种边一组里有两条边，这两条边只能出现一条，概率都是50%。

求这张图完美匹配数的期望

题解

一条边能够带来贡献的条件不是它出现了，而是它出现在了匹配中。所以我们应当直接计算每条边出现在最大匹配中的概率。

第一种边不用研究。

第二种边每一条条边出现在最大匹配中的概率都是50%。

两条边出现在最大匹配中的概率也是50%。

如果我们直接连两条50%的边的话，两条边同时出现的概率就是25%。

所以我们补一条两组的边，概率为25%。

第三种边同理，两条边出现在最大匹配中的概率是0%。

所以我们补一组-25%的边就好了。

dp的话，拿map记搜一下就好了。

代码

#include<iostream>
#include<cstdio>
#include<map>
#define N 16
using namespace std;
typedef long long ll;
const int mod=1e9+;
map<int,int>mp;
map<int,int>::iterator it;
int lo[<<N],tot,head[N],n,m;
inline int rd(){
    int x=;char c=getchar();bool f=;
    while(!isdigit(c)){if(c=='-')f=;c=getchar();}
    while(isdigit(c)){x=(x<<)+(x<<)+(c^);c=getchar();}
    return f?-x:x;
}
inline ll power(ll x,ll y){
    ll ans=;
    while(y){
        if(y&)ans=ans*x%mod;x=x*x%mod;y>>=;
    }
    return ans;
}
inline void MOD(ll &x){while(x>=mod)x-=mod;}
struct edge{int n,to,l;}e[N*N*];
inline void add(int u,int v){
    int loo=lo[u&(-u)];
    e[++tot].n=head[loo];e[tot].to=u;head[loo]=tot;e[tot].l=v;
}
int DP(int s){
    if(!s)return ;
    it=mp.find(s);
    if(it!=mp.end())return it->second;
    int loo=lo[s&(-s)];ll ans=;
    for(int i=head[loo];i;i=e[i].n){
        int v=e[i].to;if((v&s)!=v)continue;
        ans+=1ll*DP(s^v)*e[i].l%mod;MOD(ans);
    }
    return mp[s]=ans;
}
int main(){
    n=rd();m=rd();int u1,v1,u2,v2,ma=(<<n)-,opt;
    ll n2=power(,mod-),n4=power(,mod-);
    lo[]=;int cc=;
    for(int i=;i<=n;++i)cc<<=,lo[cc]=i;
    for(int i=;i<=m;++i){
        opt=rd();
        if(opt==){
            u1=rd();v1=rd();u1--;v1--;
            int s1=(<<u1)|(<<v1+n);
            add(s1,n2);
        }
        else if(opt==){
            u1=rd()-;v1=rd()-;u2=rd()-;v2=rd()-;
            int s1=(<<u1)|(<<v1+n),s2=(<<u2)|(<<v2+n);
            add(s1,n2);add(s2,n2);
            if(!(s1&s2)){
                add(s1|s2,n4);
            }
        }
        else{
            u1=rd()-;v1=rd()-;u2=rd()-;v2=rd()-;
            int s1=(<<u1)|(<<v1+n),s2=(<<u2)|(<<v2+n);
            add(s1,n2);add(s2,n2);
                if(!(s1&s2)){
                add(s1|s2,mod-n4);
            }
        }
    }
    cout<<1ll*DP(ma|(ma<<n))*power(,n)%mod;
    return ;
}