BZOJ 3812 主旋律 (状压DP+容斥) + NOIP模拟赛巨神兵(obelisk)(状压DP)

这道题跟另一道题很像，先看看那道题吧

巨神兵(obelisk)

题面

欧贝利斯克的巨神兵很喜欢有向图，有一天他找到了一张nnn个点mmm条边的有向图。欧贝利斯克认为一个没有环的有向图是优美的，请问这张图有多少个子图（即选定一个边集）是优美的？答案对 1,000,000,0071,000,000,0071,000,000,007 取模。
n<=17n<=17n<=17

分析

这道题就是枚举拓扑序最后的点集来转移

#include <bits/stdc++.h>
using namespace std;
#define LL long long
const int MAXN = 17, MAXS = 1<<17, MAXM = 250, mod = 1e9+7;
int n, m, flag[MAXS], f[MAXS], mul[MAXM+1], sum[MAXS], in[MAXS];
bool g[MAXN][MAXN];
int main ()
{scanf("%d%d", &n, &m);for(int x, y, i = 1; i <= m; i++)scanf("%d%d", &x, &y), g[x-1][y-1] = true; //邻接表存图flag[0] = -1;for(int s = 1; s < (1<<n); s++) flag[s] = flag[s>>1] * (s & 1 ? -1 : 1); //求容斥系数 奇数个为1 偶数个为-1mul[0] = 1;for(int i = 1; i <= m; i++) mul[i] = mul[i-1] * 2 % mod; //预处理2^kf[0] = 1;for(int i = 0; i < (1<<n)-1; i++){for(int k = 0; k < n; k++) in[1<<k] = 0; //计算入度 把k点的入度存在 1<<k上for(int j = 0; j < n; j++)if(i & (1<<j))for(int k = 0; k < n; k++)in[1<<k] += g[j][k];int t = (1<<n)-1-i; //t为当前状态的补集,即剩下的点sum[0] = 0;for(int s = (t-1)&t; ; s = (s-1)&t) //枚举本次选取的点对于t的补集{int now = t ^ s, last = now & -now; //now 是本次选取的点sum[now] = sum[now-last] + in[last];f[i+now] = ((LL)f[i+now] + (LL)flag[now] * mul[sum[now]] * f[i]) % mod;if(!s) break;}}printf("%d\n", (f[(1<<n)-1]+mod)%mod);
}

BZOJ 3812 主旋律

这道题做法差不多，不过是枚举拓扑序最后的强连通分量来进行转移
详见大佬博客 Miskcoo’s Space

#include <bits/stdc++.h>
using namespace std;
const int mod = 1e9 + 7;
const int MAXS = 1<<15;
const int MAXN = 15;
int n, m, Out[MAXS], In[MAXS], mul[MAXN*MAXN];
int f[MAXS], g[MAXS], bitcnt[MAXS], h[MAXS], p[MAXS];
int main () {scanf("%d%d", &n, &m);for(int i = 0, x, y; i < m; ++i) {scanf("%d%d", &x, &y);x = 1 << (x-1);y = 1 << (y-1);Out[x] |= y;In[y] |= x;}mul[0] = 1;for(int i = 1; i < n*n; ++i)mul[i] = 2ll * mul[i-1] % mod;bitcnt[0] = 0;for(int i = 1; i < (1<<n); ++i)bitcnt[i] = bitcnt[i>>1] + (i&1);for(int state = 1; state < (1<<n); ++state) {int one = state & -state, Outside = state ^ one;for(int i = Outside; i; i = (i-1)&Outside)g[state] = (g[state] - 1ll * f[state^i] * g[i] % mod) % mod;h[state] = h[Outside] + bitcnt[In[one]&Outside] + bitcnt[Out[one]&Outside];f[state] = mul[h[state]];for(int sub = state; sub; sub = (sub-1)&state) {if(sub != state) {int del = (sub^state) & -(sub^state);p[sub] = p[sub^del] + bitcnt[Out[del]&sub] - bitcnt[In[del]&(sub^state)];}else p[sub] = 0;f[state] = (f[state] - 1ll * mul[h[state^sub]+p[sub]] * g[sub] % mod) % mod;}g[state] = (g[state] + f[state]) % mod;}printf("%d\n", (f[(1<<n)-1]+mod)%mod);
}

BZOJ 3812 主旋律 (状压DP+容斥) + NOIP模拟赛巨神兵(obelisk)(状压DP)相关推荐

P5664 [CSP-S2019] Emiya 家今天的饭（dp+容斥）
题目传送门题意这个题的意思就是说,每一行只能选一个(你可以选你也可以不选,但是你只能选一个),你选出来的种类数nnn,其中n/2n/2n/2不能是同一列的东西,请问你有多少种选法. 思路这道题, ...
[CQOI2012] 局部极小值（状压DP + 容斥 + 搜索）
problem luogu-P3160 solution 这么小的数据范围,非暴力不状压.暴力 O(28!)O(28!)O(28!) 呵呵呵可以拉走了. 我们不妨从小到大填数字,这样如果局部极小值点还 ...
bzoj 1042: [HAOI2008]硬币购物（dp+容斥）
1042: [HAOI2008]硬币购物 Time Limit: 10 Sec Memory Limit: 162 MB Submit: 2555 Solved: 1537 [Submit][St ...
bzoj 4455 [Zjoi2016]小星星树形dp容斥
4455: [Zjoi2016]小星星 Time Limit: 10 Sec Memory Limit: 512 MB Submit: 643 Solved: 391 [Submit][Statu ...
BZOJ2339: [HNOI2011]卡农(dp 容斥)
题意从$1 - n$中任意选择一些数,选$m$次构成$m$个集合保证: 集合不为空任意两个集合不相同集合内各个元素xor起来等于0 Sol 神仙题Orz 我看到两种做法,一种是洛谷题解上的直接 ...
CF296B dp\容斥
传送门文章目录题意: 思路: 题意: n≤1e5n\le1e5n≤1e5 思路: 求方案数基本就是考虑dpdpdp了,看到nnn这么大可以考虑一下分情况讨论的dpdpdp状态. 设f[i][j]f ...
P5405-[CTS2019]氪金手游【树形dp,容斥,数学期望】
前言话说在LojLojLoj下了个数据发现这题的名字叫fgofgofgo 正题题目链接:https://www.luogu.com.cn/problem/P5405 题目大意 nnn张卡的权值为1 ...
AT4352-[ARC101C] Ribbons on Tree【dp,容斥】
正题题目链接: https://www.luogu.com.cn/problem/AT4352 https://atcoder.jp/contests/arc101/tasks/arc101_c 题 ...
小奇遐想树状数组实现+容斥思想
问题 M: 小奇遐想时间限制: 1 Sec 内存限制: 128 MB 提交: 165 解决: 21 [提交] [状态] [讨论版] [命题人:admin] 题目描述撷来一缕清风飘渺方知今日书 ...
[LOJ#3119][Luogu5405][CTS2019]氪金手游(DP+容斥)
先考虑外向树的做法,显然一个点在其子树内第一个出现的概率等于它的权值除以它子树的权值和.于是f[i][j]表示i的子树的权值和为j时,i子树内所有数的相互顺序都满足条件的概率,转移直接做一个背包卷积即 ...

BZOJ 3812 主旋律 (状压DP+容斥) + NOIP模拟赛巨神兵(obelisk)(状压DP)

巨神兵(obelisk)

题面

分析

BZOJ 3812 主旋律

BZOJ 3812 主旋律 (状压DP+容斥) + NOIP模拟赛巨神兵(obelisk)(状压DP)相关推荐

最新文章

热门文章

BZOJ 3812 主旋律 (状压DP+容斥) + NOIP模拟赛 巨神兵(obelisk)(状压DP)

巨神兵(obelisk)

题面

分析

BZOJ 3812 主旋律

BZOJ 3812 主旋律 (状压DP+容斥) + NOIP模拟赛 巨神兵(obelisk)(状压DP)相关推荐

最新文章

热门文章

BZOJ 3812 主旋律 (状压DP+容斥) + NOIP模拟赛巨神兵(obelisk)(状压DP)

BZOJ 3812 主旋律 (状压DP+容斥) + NOIP模拟赛巨神兵(obelisk)(状压DP)相关推荐