洛谷 P4547 & bzoj 5006 随机二分图 —— 状压DP+期望

题目：https://www.luogu.org/problemnew/show/P4547

https://www.lydsy.com/JudgeOnline/problem.php?id=5006

参考博客：https://www.cnblogs.com/yanshannan/p/9452802.html

注意同一个点连出去的两条边本来就不能一起选！

每次调用 map 会很慢！所以修改的时候新定义一个 &tmp，就能过了。

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#define il inline
using namespace std;
typedef long long ll;
int rd()
{
  int ret=,f=; char ch=getchar();
  while(ch<''||ch>''){if(ch=='-')f=; ch=getchar();}
  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
  return f?ret:-ret;
}
int const xn=,xm=(<<)+,mod=1e9+;
int n,m,bin[],ct;
map<int,int>f[xm];
struct N{
  int s,t,w;
  N(int s=,int t=,int w=):s(s),t(t),w(w) {}
}ed[];
il ll pw(ll a,int b){ll ret=; for(;b;b>>=,a=a*a%mod)if(b&)ret=ret*a%mod; return ret;}
il int upt(int x){while(x>=mod)x-=mod; while(x<)x+=mod; return x;}
il int dfs(int s,int t)
{
  if(!s&&!t)return ;
  if(f[s][t])return f[s][t];
  int &tmp=f[s][t];//
  for(int i=;i<=ct;i++)
    {
      if((ed[i].s<<)<=s||(ed[i].s|s)!=s||(ed[i].t|t)!=t)continue;//t
      tmp=(tmp+(ll)dfs(s^ed[i].s,t^ed[i].t)*ed[i].w)%mod;
    }
  return tmp=upt(tmp);//-
}
int main()
{
  n=rd(); m=rd();
  bin[]=; for(int i=;i<=;i++)bin[i]=bin[i-]*;
  int inv2=pw(,mod-),inv4=pw(,mod-);
  for(int i=,tp,a,b,c,d;i<=m;i++)
    {
      tp=rd(); a=rd(); b=rd();
      ed[++ct]=N(bin[a-],bin[b-],inv2);
      if(!tp)continue; c=rd(); d=rd();
      ed[++ct]=N(bin[c-],bin[d-],inv2);
      if(a==c||b==d)continue;//!
      if(tp==)ed[++ct]=N(bin[a-]|bin[c-],bin[b-]|bin[d-],inv4);
      if(tp==)ed[++ct]=N(bin[a-]|bin[c-],bin[b-]|bin[d-],-inv4);
    }
  printf("%lld\n",(ll)dfs(bin[n]-,bin[n]-)*bin[n]%mod);
  return ;
}