传送门
题解:
首先算期望就是对于每一种完美匹配算出现概率然后求个和。
如果只有第 0 0 0类边的话我们可以考虑如下的DP,把现在左侧还没有匹配的点压起来,右侧还没有匹配的点压起来,用一个int表示状态。
那么显然我们需要知道的就是这些点自己期望能够产生多少种合法匹配使得自己匹配完。
由于边是无序的,转移的时候直接强行枚举一边把集合中编号最小的点搞掉即可。
考虑 1 1 1类边组和 2 2 2类边组。
这两类的限制太强了,考虑先拆成两个 0 0 0类边。
然后 1 1 1类再加一个概率为 1 4 \frac{1}{4} 41的边组, 2 2 2类再加一个概率为 − 1 4 -\frac{1}{4} −41的边组。
然后转移的时候边和边组都考虑一下即可。
为什么这样是对的。
实际上我们只关心我们要的匹配里面出现了哪些边,对于以上两类边,两条边在匹配中单独出现的概率都是 1 2 \frac{1}{2} 21,同时出现的概率是 1 4 \frac{1}{4} 41,所以我们加一个边组把同时出现的概率补回去即可。对于 1 1 1类,同时出现的概率需要 + 1 4 +\frac{1}{4} +41,对于 2 2 2类同时出现的概率需要 − 1 4 -\frac{1}{4} −41。
记忆化一下,开个哈希表,然后就做完了。
代码:
#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const
cs int mod=1e9+7;
inline int add(int a,int b){a+=b-mod;return a+(a>>31&mod);}
inline int mul(int a,int b){ll r=(ll)a*b;return r>=mod?r%mod:r;}
inline void Inc(int &a,int b){a+=b-mod;a+=a>>31&mod;}
inline void Mul(int &a,int b){a=mul(a,b);}
int n,m;
int to[20],e[20][20],val[20][20];
std::unordered_map<int,int> f;
inline int dp(int s){
if(s==0)return 1;if(f.count(s))return f[s];
int mi=-1,res=0;
for(int re i=0;i<n&&mi==-1;++i)if((s>>i)&1)mi=i;
int p=(s>>15)&to[mi],del=s^(1<<mi);
for(int re i=0;i<n;++i)if(p>>i&1){
int k=del^(1<<(i+15));
Inc(res,dp(k));int ne=e[mi][i];
if(ne&&(ne&k)==ne)Inc(res,mul(val[mi][i],dp(k^ne)));
}
return f[s]=res;
}
signed main(){
#ifdef zxyoi
freopen("match.in","r",stdin);
#endif
scanf("%d%d",&n,&m);
for(int re i=1;i<=m;++i){
int t,u1,v1,u2,v2;scanf("%d%d%d",&t,&u1,&v1);--u1,--v1;
switch(t){
case 0:to[u1]|=1<<v1;break;
case 1:
scanf("%d%d",&u2,&v2);--u2,--v2;
if(u1>u2)std::swap(u1,u2),std::swap(v1,v2);
to[u1]|=1<<v1;to[u2]|=1<<v2;
e[u1][v1]=1<<u2|(1<<v2+15);
val[u1][v1]=1;break;
case 2:
scanf("%d%d",&u2,&v2);--u2,--v2;
if(u1>u2)std::swap(u1,u2),std::swap(v1,v2);
to[u1]|=1<<v1;to[u2]|=1<<v2;
e[u1][v1]=1<<u2|(1<<v2+15);
val[u1][v1]=mod-1;break;
}
}
int s=(1<<n)-1;
std::cout<<dp(s<<15|s)<<"\n";
return 0;
}