洛谷P4547 [THUWC2017]随机二分图

题目描述
题解

考虑一个完美匹配出现的概率 × 2 n \times 2^n ×2n 对答案的贡献,初始是 1 1 1 ,如果有出现一组 t = 2 t=2 t=2 的边的话,那贡献就是 0 0 0 ,否则每出现一组 t = 1 t=1 t=1 的边就要 × 2 \times 2 ×2 ,所以有个暴力的 dp \text{dp} dp f [ s 1 ] [ s 2 ] f[s1][s2] f[s1][s2] 表示左边出现点的状态为 s 1 s1 s1 ,右边为 s 2 s2 s2 的答案,为了避免重复我们每次挑 s 1 s1 s1 中编号最小的点进行转移,即 f [ s 1 ] [ s 2 ] = ∑ y ∈ s 2 f [ s 1 ∧ x ] [ s 2 ∧ y ] f[s1][s2]=\sum_{y \in s2}f[s1 \wedge x][s2 \wedge y] f[s1][s2]=ys2f[s1x][s2y] ,考虑如果 ( x , y ) (x,y) (x,y) 有一组 t = 1 t=1 t=1 的边 ( u , v ) (u,v) (u,v) 也在 s 1 , s 2 s1,s2 s1,s2 集合中的话,那要再加上 f [ s 1 ∧ x ∧ u ] [ s 2 ∧ y ∧ v ] f[s1 \wedge x \wedge u][s2 \wedge y \wedge v] f[s1xu][s2yv] ,反之如果有一组 t = 2 t=2 t=2 的在集合中的话就要扣掉,记忆化搜索即可。

代码
#include <bits/stdc++.h>
using namespace std;
const int P=1e9+7;map<int,int>f;
int n,m,k,a[15][15],c[1<<15],b[15][15][2];
int X(int x){return x>=P?x-P:x;}
int F(int i,int j){
	if (f.count(i<<n|j)) return f[i<<n|j];
	int w=0,x=i&-i,u,v;
	for (int y=j&-j,k=j;k;k^=y,y=k&-k){
		if (!~a[c[x]][c[y]]) continue;
		w=X(w+F(i^x,j^y));
		u=b[c[x]][c[y]][0];v=b[c[x]][c[y]][1];
		if (a[c[x]][c[y]] && (u&i) && (v&j) && u!=x && v!=y){
			if (a[c[x]][c[y]]&1) w=X(w+F(i^x^u,j^y^v));
			else w=X(w+P-F(i^x^u,j^y^v));
		}
	}
	return f[i<<n|j]=w;
}
int main(){
	cin>>n>>m;k=1<<n;
	for (int i=0;i<n;i++) c[1<<i]=i;
	memset(a,-1,sizeof a);
	for (int x,y,u,v,o,i=1;i<=m;i++){
		scanf("%d%d%d",&o,&x,&y);
		x--;y--;a[x][y]=o;
		if (o){
			scanf("%d%d",&u,&v);
			u--;v--;a[u][v]=o;
			b[x][y][0]=1<<u;b[x][y][1]=1<<v;
			b[u][v][0]=1<<x;b[u][v][1]=1<<y;
		}
	}
	f[0]=1;
	printf("%d\n",F(k-1,k-1));
	return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值