关于一类容斥原理设计 dp 状态的探讨
写在前面
为什么要写?因为自己学不明白希望日后能掌握。
大体思路大概是
- 设计一个容斥的方案,并使其贡献可以便于计算。
- 得出 dp 状态,然后优化以得出答案。
下列所有类似 \([l,r]\) 这样的都是离散的。
1.
\(n\) 个点,每个点有一个能选择的颜色 \(a_i\),左右相邻的点不能同色,求方案数。
如果我们使用容斥的思想,强制 \(k\) 段的颜色相同,这个限制下的方案数对答案的贡献的容斥系数就是 \((-1)^{n-k}\)。这应该是相邻颜色不同的方案数的一个非常平凡的trick。(但是我不会
可以设 \(f_i\) 表示统计到前 \(i\) 个点所容斥的答案和。枚举 \([j,i]\) 这一段强制颜色相等。
\]
这个东西可以用单调栈维护一下。
注意到这个东西可以拓展到环上。把 \(a_i\) 最小的位置轮换到最前面,然后你发现 \(f_i\) 其实就是强制了 \([i+1,n]\) 和 \(1\) 的颜色相同的答案。全部加起来就好了。
s[0] = f[0] = 1; int top = 0;
ll sum = 0;
fo(i, 1, c) {
while(top && b[stc[top]] > b[i])
sum = (sum + (ll)(s[stc[top] - 1] - (stc[top] == 1 ? 0 :
s[stc[top - 1] - 1]) + mod) * (b[i] - b[stc[top]] + mod)) % mod,
--top;
stc[++top] = i;
sum = (sum + (ll)f[i - 1] * b[i]) % mod;
f[i] = mod - sum;
s[i] = (f[i] + s[i - 1]) % mod;
}
2.
\(n\) 个点,一个区间可以覆盖 \([l_i,r_i]\) 这一段,每个区间有一个价值 \(v_i\) ,定义一种“覆盖”为每个点至少被一个区间所覆盖的方案,其价值为所有所选区间的价值积,求所有覆盖的价值之和。
考虑强制 \(k\) 个点不被覆盖,那么这种情况对答案的贡献的容斥系数就是 \((-1)^k\)。其贡献就是这些点之间的区间的乘积之和。
这样的话,设 \(f_i\) 表示 \(i\) 点被钦定,枚举 \(j\) 表示上一个钦定点,有
\]
这玩意可以线段树优化!考虑线段树的每一个位置记录的是它作为 \(j\) 造成的贡献,假设现在新加入一个区间 \(k\) ,它能使 \([0,l_k)\) 的位置的贡献发生变化,乘上 \((1+v_k)\)。
#include <cstdio>
#include <vector>
#include <cstring>
#include <algorithm>
#define ll long long
#define fo(i, a, b) for(int i = (a); i <= (b); ++i)
#define fd(i, a, b) for(int i = (a); i >= (b); --i)
using namespace std;
inline void read(int &x) {
x = 0; char ch = getchar();
while(ch < '0' || ch > '9') ch = getchar();
while(ch >= '0' && ch <= '9') x = (x << 1) + (x << 3) + (ch ^ 48), ch = getchar();
}
const int N = 2e5 + 10, mod = 1e9 + 7;
namespace Seg {
#define ls t << 1
#define rs ls | 1
#define mid ((l + r) >> 1)
int tr[N << 2], pro[N << 2];
inline void mul(int t, int v) {tr[t] = (ll)tr[t] * v % mod, pro[t] = (ll)pro[t] * v % mod;}
inline void push_down(int t) {
if(pro[t] > 1) {
mul(ls, pro[t]), mul(rs, pro[t]);
pro[t] = 1;
}
}
void build(int t, int l, int r) {
pro[t] = 1;
if(l == r) return;
build(ls, l, mid), build(rs, mid + 1, r);
}
void change(int t, int l, int r, int w, int v) {
tr[t] = (tr[t] + v) % mod;
if(l == r) return ;
push_down(t);
w <= mid ? change(ls, l, mid, w, v) : change(rs, mid + 1, r, w, v);
}
void update(int t, int l, int r, int fl, int fr, int v) {
if(fl <= l && r <= fr) return mul(t, v);
push_down(t);
fl <= mid && (update(ls, l, mid, fl, fr, v), 1);
fr > mid && (update(rs, mid + 1, r, fl, fr, v), 1);
tr[t] = (tr[ls] + tr[rs]) % mod;
}
int query(int t, int l, int r, int fl, int fr) {
if(fl <= l && r <= fr) return tr[t];
push_down(t);
int ret = 0;
fl <= mid && (ret = (ret + query(ls, l, mid, fl, fr)) % mod);
fr > mid && (ret = (ret + query(rs, mid + 1, r, fl, fr)) % mod);
return ret;
}
}
struct Op {
int l, r, v;
}p[N];
vector<int> q[N];
int n, m, f[N];
int main() {
freopen("gugugu.in", "r", stdin);
freopen("gugugu.out", "w", stdout);
read(n), read(m);
fo(i, 1, m) read(p[i].l), read(p[i].r), read(p[i].v), q[p[i].r].push_back(i);
Seg::build(1, 0, n);
Seg::change(1, 0, n, 0, 1);
fo(i, 1, n + 1) {
Seg::change(1, 0, n, i, f[i] = mod - Seg::query(1, 0, n, 0, i - 1));
for(auto k : q[i])
Seg::update(1, 0, n, 0, p[k].l - 1, (p[k].v + 1) % mod);
}
printf("%d\n", mod - f[n + 1]);
return 0;
}
To be continued..
关于一类容斥原理设计 dp 状态的探讨的更多相关文章
- HDU 4336 Card Collector (期望DP+状态压缩 或者 状态压缩+容斥)
题意:有N(1<=N<=20)张卡片,每包中含有这些卡片的概率,每包至多一张卡片,可能没有卡片.求需要买多少包才能拿到所以的N张卡片,求次数的期望. 析:期望DP,是很容易看出来的,然后由 ...
- [提升性选讲] 树形DP进阶:一类非线性的树形DP问题(例题 BZOJ4403 BZOJ3167)
转载请注明原文地址:http://www.cnblogs.com/LadyLex/p/7337179.html 树形DP是一种在树上进行的DP相对比较难的DP题型.由于状态的定义多种多样,因此解法也五 ...
- dp状态压缩
dp状态压缩 动态规划本来就很抽象,状态的设定和状态的转移都不好把握,而状态压缩的动态规划解决的就是那种状态很多,不容易用一般的方法表示的动态规划问题,这个就更加的难于把握了.难点在于以下几个方面:状 ...
- HDU 1074 Doing Homework (dp+状态压缩)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1074 题目大意:学生要完成各科作业, 给出各科老师给出交作业的期限和学生完成该科所需时间, 如果逾期一 ...
- hdu_4352_XHXJ's LIS(数位DP+状态压缩)
题目连接:hdu_4352_XHXJ's LIS 题意:这题花大篇篇幅来介绍电子科大的一个传奇学姐,最后几句话才是题意,这题意思就是给你一个LL范围内的区间,问你在这个区间内最长递增子序列长度恰为K的 ...
- hdu 4352 数位dp + 状态压缩
XHXJ's LIS Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total ...
- HDU 1074 Doing Homework(DP状态压缩)
题意:有n门功课需要完成,每一门功课都有时间期限以及你完成所需要的时间,如果完成的时间超出时间期限多少单位,就会被减多少学分,问以怎样的功课完成顺序,会使减掉的学分最少,有多个解时,输出功课名字典序最 ...
- 【bzoj1076】[SCOI2008]奖励关 期望dp+状态压缩dp
题目描述 你正在玩你最喜欢的电子游戏,并且刚刚进入一个奖励关.在这个奖励关里,系统将依次随机抛出k次宝物,每次你都可以选择吃或者不吃(必须在抛出下一个宝物之前做出选择,且现在决定不吃的宝物以后也不能再 ...
- hdu4336 Card Collector(概率DP,状态压缩)
In your childhood, do you crazy for collecting the beautiful cards in the snacks? They said that, fo ...
随机推荐
- RunLoop基础知识以及GCD
- 1.1 字面意思 a 运行循环 b 跑圈 - 1.2 基本作用(作用重大) a 保持程序的持续运行(ios程序因而能一直活着不会死) b 处理app中的各种事件(比如触摸事件 ...
- Dubbo服务分组
服务分组与多版本控制的使用方式几乎是相同的,只要将version替换为group即可.但使用目的不同.使用版本控制的目的是为了升级,将原有老版本替换掉,将来不再提供老版本的服务,所以不同版本间不能出现 ...
- 【JS】枚举类型
https://zhuanlan.zhihu.com/p/79137838 相当于用数字来代替一串字母 /** * 时间:2019年8月18日 * 前端教程: https://www.pipipi.n ...
- 如何查看Python的版本号
一.如何查看Python的版本号 win+r输入cmd在输入:python --version回车即可
- 删除其他列Table.SelectColumns(Power Query 之 M 语言)
数据源: "姓名""基数""个人比例""个人缴纳""公司比例""公司缴纳"&qu ...
- 创建项目文件(Project)
<Project2016 企业项目管理实践>张会斌 董方好 编著 按照张同学和董同学的说法,创建项目文件首选是利用可以参照的项目计划模板,如果找不到,那就利用现有的项目文件,实在这些都没有 ...
- 3、回溯算法解题套路框架——Go语言版
前情提示:Go语言学习者.本文参考https://labuladong.gitee.io/algo,代码自己参考抒写,若有不妥之处,感谢指正 关于golang算法文章,为了便于下载和整理,都已开源放在 ...
- 小迪安全 Web安全 基础入门 - 第一天 - 操作系统&名词&文件下载&反弹SHELL&防火墙绕过
一.专业名词 1.POC:(Proof of Concept),即概念验证.漏洞报告中的POC是一段说明或一个攻击的样例使读者能够确认这个漏洞是真实存在的. 2.EXP:exploit,即漏洞利用.对 ...
- .NET 云原生架构师训练营(建立系统观)--学习笔记
目录 目标 ASP .NET Core 什么是系统 什么是系统思维 系统分解 什么是复杂系统 作业 目标 通过整体定义去认识系统 通过分解去简化对系统的认识 ASP .NET Core ASP .NE ...
- CF570A Elections 题解
Content 有 \(n\) 个候选人和 \(m\) 个城市,每个城市可以给每个候选人投票,已知第 \(i\) 个城市给第 \(j\) 个人投的选票数是 \(a_{i,j}\).我们将第 \(i\) ...