[BZOJ5006][LOJ#2290][THUWC2017]随机二分图(概率+状压DP)

https://loj.ac/problem/2290

题解：https://blog.csdn.net/Vectorxj/article/details/78905660

不是很好理解，对于边(x1,y1)和(x2,y2)，可以分“x1或y1已匹配”，“x2或y2已匹配”，“x1,x2,y1,y2均未匹配”三种情况考虑拆边的正确性。

状压的时候，对于当前左边已经匹配的集合，只需要枚举左边已匹配的最后一个是用哪条边匹配的即可，也就是程序里的S<(1<<T)。

不要用顺推，记忆化搜索会忽略一些用不到的状态，所以会快一些。

 1 #include<map>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #define rep(i,l,r) for (int i=l; i<=r; i++)
 5 using namespace std;
 6
 7 const int N=16,mod=1000000007,inv2=(mod+1)>>1,inv4=(mod+1)>>2;
 8 int n,m,cnt,t,x,y;
 9 struct E{ int S,p,c; E(int _S=0,int _p=0,int _c=0):S(_S),p(_p),c(_c){} }G[N*N<<2];
10 map<int,int>f[1<<N];
11
12 inline void add(int &x,int y){ x+=y; if (x>=mod) x-=mod; }
13 int F(int S){
14     if (!S) return 1;
15     int T0=S>>n,S0=S^(T0<<n);
16     if (f[S0].count(T0)) return f[S0][T0];
17     int res=0;
18     rep(i,1,cnt){
19         int T=G[i].S;
20         if ((T&S)==T && S<(T<<1)) add(res,1ll*F(S^T)*G[i].p%mod);
21     }
22     return f[S0][T0]=res;
23 }
24
25 int main(){
26     freopen("bzoj5006.in","r",stdin);
27     freopen("bzoj5006.out","w",stdout);
28     scanf("%d%d",&n,&m);
29     rep(i,1,m){
30         scanf("%d%d%d",&t,&x,&y);
31         int S1=(1<<(x-1))|(1<<(y+n-1)); G[++cnt]=E(S1,inv2,1);
32         if (t){
33             scanf("%d%d",&x,&y);
34             int S2=(1<<(x-1))|(1<<(y+n-1)); G[++cnt]=E(S2,inv2,1);
35             if (S1 & S2) continue;
36             G[++cnt]=E(S1|S2,((t==1)?inv4:mod-inv4),2);
37         }
38     }
39     printf("%lld\n",(1ll<<n)*F((1<<(n*2))-1)%mod);
40     return 0;
41 }

转载于:https://www.cnblogs.com/HocRiser/p/8941925.html

[BZOJ5006][LOJ#2290][THUWC2017]随机二分图(概率+状压DP)相关推荐

【PKUWC2018】随机算法【状压dp】【组合计数】
题意:一张nnn个点mmm条边的无向无权图,求以下算法计算最大独立集的正确率:随机一个排列,依次考虑排列中每一个点,如果不与任何一个独立集中的点相邻则将其加入独立集.模9982443539982443 ...
【Gym - 101915D】Largest Group（二分图最大团，状压dp）
题干: 大黑山上有小小民和小小涛两种物种,山东人小李想要研究这两种物种的关系奇怪的是大黑山上有相同数量的小小民和小小涛.小李数了数一共有 P 个,小李分别给P个小小民和小小涛编号 1 - P 号,已 ...
loj 1316(spfa预处理+状压dp)
题目链接:http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=27024 题意:求0-(n-1)的经过最多的标记的点的最短路. 思路 ...
Educational Codeforces Round 13 E. Another Sith Tournament 状压dp
E. Another Sith Tournament 题目连接: http://www.codeforces.com/contest/678/problem/E Description The rul ...
POJ3163 King of Fighters 状压DP/费用流
回归后的第一篇blog,康复打板中这种题我做了2个多小时..(反正我现在已经失去对题目难度的辨识了) 很容易想到O(mn2n)O(mn2n)O(mn2^n)的状压做法,但是超时,考虑优化,想到不正确 ...
POJ 1038 Bugs Integrated Inc （复杂的状压DP）
\(POJ~1038~~*Bugs~Integrated~Inc:\) (复杂的状压DP) \(solution:\) 很纠结的一道题目,写了大半天,就想练练手,结果这手生的.其实根据之前那道炮兵阵地 ...
codeforces 8C. Looking for Order 状压dp
题目链接给n个物品的坐标, 和一个包裹的位置, 包裹不能移动. 每次最多可以拿两个物品, 然后将它们放到包里, 求将所有物品放到包里所需走的最小路程. 直接状压dp就好了. #include < ...
UVA10296 Jogging Trails（中国邮递员问题）（欧拉回路、一般图最大权匹配 / 状压DP）
整理的算法模板合集: ACM模板目录思路 UVA10296 Jogging Trails 题目翻译: 给你n个点,m条无向边,每条边有一定的距离数值,构造成一个连通图.问从任意一点出发,遍历所有的 ...
POJ 2411 Mondriaan‘s Dream（最清楚好懂的状压DP讲解）（连通性状态压缩DP）
poj 2411 Mondriaan's Dream(最清晰的状压DP解析) 闫氏DP大法好我们这里是一列一列地来,因为是一个棋盘性的状态压缩DP,从哪个方向都一样摆放的小方格总方案数等价于横 ...
【每日DP】day2、P1879 [USACO06NOV]Corn Fields G玉米地（状压DP模板题）难度⭐⭐⭐★
昨天的每日DP我还在写01背包,今天就到状压DP了,真刺激. P1879 [USACO06NOV]Corn Fields G 题目链接输入 2 3 1 1 1 0 1 0 输出 9 一道简单的状压D ...

[BZOJ5006][LOJ#2290][THUWC2017]随机二分图(概率+状压DP)

[BZOJ5006][LOJ#2290][THUWC2017]随机二分图(概率+状压DP)相关推荐

最新文章

热门文章