《算法竞赛进阶指南》0x58 数据结构优化DP

最新推荐文章于 2024-05-21 23:16:53 发布
最新推荐文章于 2024-05-21 23:16:53 发布
阅读量323 收藏
点赞数
分类专栏： 算法竞赛进阶指南 文章标签： 算法 数据结构 c++ 动态规划
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_53899788/article/details/130157712
版权

0x58 数据结构优化DP

清理班次

题意：

给定区间 [ 1 , T ] [1,T] [1,T]，给定 n n n 条线段 [ l i , r i ] [l_i, r_i] [li​,ri​]，选择最少数量的线段，使区间每个点都被覆盖

解析：

将所有线段按右端点 r i r_i ri​ 升序排序。

令 f i f_i fi​ 为覆盖区间 [ 1 , i ] [1,i] [1,i] 所有点的最少线段数。对每条线段 [ l i , r i ] [l_i, r_i] [li​,ri​]，都有转移 f r i = min ⁡ ( f r i , min ⁡ j = l i − 1 r i f j + 1 ) f_{r_i} = \min(f_{r_i},\min\limits_{j = l_i-1}\limits^{r_i}f_j+1) fri​​=min(fri​​,j=li​−1minri​​fj​+1)

时间复杂度为 O ( n 2 ) O(n^2) O(n2)

但可以发现 min ⁡ j = l i − 1 r i f j \min\limits_{j = l_i-1}\limits^{r_i}f_j j=li​−1minri​​fj​ 是区间最小值，可以用线段树维护最小值，支持区间查询和单点修改。

初值 f 0 f_0 f0​ = 0，将所有点的坐标+1。因此初值为 f 1 = 0 f_1 = 0 f1​=0

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef double db;
#define fi first
#define se second
#define debug(x) cerr << #x << ": " << (x) << endl
#define rep(i, a, b) for(int i = (a); i <= (b); i++)
const int maxn = 1e6+10;
const int maxm = 1e5+10;
const int INF = 0x3f3f3f3f;
typedef pair<int, int> pii;

struct seg{
	int l, r;
	seg(){}
	seg(int l, int r) : l(l), r(r){}
	bool operator < (const seg &b) const{
		return r < b.r;
	}
}s[maxn];
int f[maxn];
int n, T, l, r;

inline int ls(int x){return x << 1;}
inline int rs(int x){return x << 1 | 1;}

int t[maxn << 2];
void pushup(int k){
	t[k] = min(t[ls(k)], t[rs(k)]);
}
void build(int k, int l, int r){
	if(l == r){
		t[k] = INF;
		return;
	}
	int mid = (l+r) >> 1;
	build(ls(k), l, mid);
	build(rs(k), mid+1, r);
	pushup(k);
}
void modify(int k, int l, int r, int pos, int v){
	if(l == r && l == pos){
		t[k] = v;
		return;
	}
	int mid = (l+r) >> 1;
	if(pos <= mid)
		modify(ls(k), l, mid, pos, v);
	else
		modify(rs(k), mid+1, r, pos, v);
	pushup(k);
}
int query(int k, int l, int r, int x, int y){
	if(x <= l && y >= r)
		return t[k];
	int mid = (l+r) >> 1;
	int res = INF;
	if(x <= mid)
		res = min(res, query(ls(k), l, mid, x, y));
	if(y > mid)
		res = min(res, query(rs(k), mid+1, r, x, y));
	return res;
}

int main(){
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	
	cin >> n >> T;
	T++;
	rep(i, 1, n){
		cin >> l >> r;
		l++, r++;
		l = max(2, l); r = min(T, r);
		s[i] = seg(l, r);
	}
	
	sort(s+1, s+1+n);
	build(1, 1, T);
	modify(1, 1, T, 1, 0);
	memset(f, INF, sizeof(f));
	f[1] = 0;
	
	rep(i, 1, n){
		f[s[i].r] = min(f[s[i].r], query(1, 1, T, s[i].l-1, s[i].r)+1);
		modify(1, 1, T, s[i].r, f[s[i].r]);
	}
	
	if(f[T] == INF)
		cout << -1 << endl;
	else
		cout << f[T] << endl;
	return 0;
}


清理班次2

题意：

给定区间 [ M , E ] [M,E] [M,E]，给定 n n n 条线段 [ l i , r i ] [l_i, r_i] [li​,ri​]，选择该线段的代价为 c i c_i ci​，询问使 [ M , E ] [M,E] [M,E] 被覆盖的最小代价。

解析：

基本同上题。

f i f_i fi​ 表示覆盖区间 [ M , i ] [M,i] [M,i] 的最小代价。

对于每条线段， f b i = min ⁡ ( f b i , min ⁡ j = l i − 1 r i f j + c i ) f_{b_i} = \min(f_{b_i},\min\limits_{j = l_i-1}\limits^{r_i}f_j+c_i) fbi​​=min(fbi​​,j=li​−1minri​​fj​+ci​)

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef double db;
#define fi first
#define se second
#define debug(x) cerr << #x << ": " << (x) << endl
#define rep(i, a, b) for(int i = (a); i <= (b); i++)
const int maxn = 1e5+10;
const int maxm = 1e5+10;
const ll INF = 0x3f3f3f3f3f3f3f3f;
typedef pair<int, int> pii;

struct seg{
	int l, r; 
	ll c;
	seg(){}
	seg(int l, int r, ll c) : l(l), r(r), c(c){}
	bool operator < (const seg &b) const{
		return r < b.r;
	}
}s[maxn];
ll f[maxn];
int n, L, R, l, r, w;

inline int ls(int x){return x << 1;}
inline int rs(int x){return x << 1 | 1;}

ll t[maxn << 2];
void pushup(int k){
	t[k] = min(t[ls(k)], t[rs(k)]);
}
void build(int k, int l, int r){
	if(l == r){
		t[k] = INF;
		return;
	}
	int mid = (l+r) >> 1;
	build(ls(k), l, mid);
	build(rs(k), mid+1, r);
	pushup(k);
}
void modify(int k, int l, int r, int pos, ll v){
	if(l == r && l == pos){
		t[k] = v;
		return;
	}
	int mid = (l+r) >> 1;
	if(pos <= mid)
		modify(ls(k), l, mid, pos, v);
	else
		modify(rs(k), mid+1, r, pos, v);
	pushup(k);
}
ll query(int k, int l, int r, int x, int y){
	if(x <= l && y >= r)
		return t[k];
	int mid = (l+r) >> 1;
	ll res = INF;
	if(x <= mid)
		res = min(res, query(ls(k), l, mid, x, y));
	if(y > mid)
		res = min(res, query(rs(k), mid+1, r, x, y));
	return res;
}


int main(){
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	
	cin >> n >> L >> R;
	L++, R++;
	rep(i, 1, n){
		cin >> l >> r >> w;
		l++, r++;
		l = max(l, L); r = min(r, R);
		s[i] = seg(l, r, w);
	}
	
	sort(s+1, s+1+n);
	build(1, L-1, R);
	modify(1, L-1, R, L-1, 0);
	memset(f, INF, sizeof(f));
	f[L-1] = 0;
	
	rep(i, 1, n){
		f[s[i].r] = min(f[s[i].r], query(1, L-1, R, s[i].l-1, s[i].r) + s[i].c);
		modify(1, L-1, R, s[i].r, f[s[i].r]);
	}
	
	if(f[R] == INF)
		cout << -1 << endl;
	else
		cout << f[R] << endl;
	return 0;
}


赤壁之战

题意：

给定长度为 n n n 的序列 a a a，询问有多少个长度为 m m m 的严格递增子序列。

解析：

令 f i , j f_{i,j} fi,j​ 为在以第 i i i 个数为结尾的，长度为 j j j 的严格递增子序列的个数

长度大的递增子序列是由长度小的递增子序列后面接上几个数形成的，所以按长度划分dp 的阶段。

f i , j = ∑ k = 1 i − 1 f k , j − 1 [ a k < a i ] f_{i,j} = \sum\limits_{k=1}\limits^{i-1} f_{k, j-1}[a_k < a_i] fi,j​=k=1∑i−1​fk,j−1​[ak​<ai​]

时间复杂度为 O ( n 3 ) O(n^3) O(n3)，需要优化。

对于长度为 j j j 的阶段，需要快速求和上一阶段满足条件的状态。对 i i i 而言，上一阶段有贡献的状态需要满足条件： k < i k<i k<i 且 a k < a i a_k<a_i ak​<ai​，二维偏序，可以用树状数组在 log ⁡ n \log n logn 的时间内解决。

所以，用树状数组维护上一阶段的状态，在本阶段，查询 a i − 1 a_i-1 ai​−1 的前缀和。

初始条件： f i , 1 = 1 f_{i,1} = 1 fi,1​=1，需要离散化，答案为 ∑ i = 1 n f i , m \sum\limits_{i=1}\limits^n f_{i, m} i=1∑n​fi,m​

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef double db;
#define fi first
#define se second
#define debug(x) cerr << #x << ": " << (x) << endl
#define rep(i, a, b) for(int i = (a); i <= (b); i++)
const int maxn = 5e3+10;
const int maxm = 5e3+10;
const ll mod = 1e9+7;
const int N = 5e3;
const int INF = 0x3f3f3f3f;
typedef pair<int, int> pii;

inline int lowbit(int x){return x & (-x);}
struct BIT{
	ll c[maxn];
	void add(int x, int v){
		for(; x <= N; x += lowbit(x))
			c[x] = (c[x] + v) % mod;
	}
	ll query(int x){
		ll res = 0;
		while(x){
			res = (res+c[x]) % mod;
			x -= lowbit(x); 
		}
		return res;
	}
	void init(){
		memset(c, 0, sizeof(c));
	}
}tr;
ll dp[maxn][maxn];
vector<ll> t;
ll a[maxn], n, m;
void solve(int T){
	
	cin >> n >> m;
	t.clear();
	tr.init();
	for(int i = 1; i <= n; i++){
		cin >> a[i];
		t.push_back(a[i]);
	}
	sort(t.begin(), t.end());
	t.erase(unique(t.begin(), t.end()), t.end());
	for(int i = 1; i <= n; i++)
		a[i] = lower_bound(t.begin(), t.end(), a[i])-t.begin()+1;
	
	for(int i = 1; i <= n; i++)
		dp[i][1] = 1;
	for(int len = 2; len <= m; len++){
		tr.init();
		for(int i = 1; i <= n; i++){
			dp[i][len] = tr.query(a[i]-1);
			tr.add(a[i], dp[i][len-1]);
		}
	}
	int res = 0;
	for(int i = 1; i <= n; i++)
		res = (res + dp[i][m]) % mod;
	printf("Case #%d: ", T);
	cout << res << endl;
	return;
}
int main(){
	
	int T;
	cin >> T;
	for(int i = 1; i <= T; i++){		
		solve(i);
	}
	return 0;
}
hydqiln
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
leetcode中国-ACAG-Code:《算法竞赛进阶指南》（AlgorithmCompetitionAdvancedGuide）中例题和习
06-29
《算法竞赛进阶指南》练习代码 《算法竞赛进阶指南》全部例题和习题的练习代码，希望能给大家带来帮助，如有错误欢迎批评指正。 总进度：0/323 关于作者：WHZ0325 QQ：1669135547 Blog：（书中题目的题解以后会整理...
算法竞赛进阶指南 0x58 数据结构优化DP
xjsc01的博客
07-31 149
表示覆盖1~i的区间，所需要的最少的奶牛数目，这样也是混乱，因为并不知道这些情况中区间覆盖的最右端究竟是在哪里。你需要帮助约翰排列出一个合理的奶牛的清理班次，使得每个班次都有奶牛在进行清理，而且动用的奶牛数量可以尽可能的少。在实际的情况下，并不需要完全枚举T，因为如果对于某一个i，不是以右端点结尾的，那么。输出一个整数，表示在每个班次都有奶牛清理的情况下，所需的奶牛最小数量。]时间段内都有奶牛在清理，并且所需支付给奶牛的报酬最少。如果无法做到在要求时间段内都有奶牛清理，则输出−1。...
算法竞赛进阶指南0x58 数据结构优化DP
m0_73935863的博客
04-12 58
算法竞赛进阶指南0x58 数据结构优化DP
【算法竞赛进阶指南】学习笔记
juruo_c的博客
01-16 1104
下面是个目录呀0xFF 前言0x00 基本算法0x01 位运算0x02 递推与递归0x03 前缀和与差分0x04 二分0x05 排序0x06 倍增0x07 贪心0x08 总结与习题0x10 基本数据结构 0xFF 前言 一时兴起想要来更新一下，希望能把这个系列做完吧！加油！！ 每一章我会把书里有难度的知识点总结一下，再附上书里的一些例题与习题。 平时遇到相关的习题也会加到相关习题部分。所以这是一个进阶版的博客，不适合小白阅读哦！ 0x00 基本算法 0x01 位运算 内容简介与例题习题 0x02
算法竞赛进阶指南学习目录
qq_61143965的博客
12-17 481
算法学习
五大基础算法——枚举，算法数据结构
04-07
五大基础算法——枚举，算法数据结构 五大常用算法
VL53L0X中文数据手册.pdf
04-06
VL53L0X中文数据手册，
数据结构常见算法代码实现.zip
04-25
常见的顺序表，单链表，双向链表，队列，栈，二叉树，哈希表等完整代码介绍
Cow Exhibition G的来龙去脉
m0_61150571的博客
05-15 454
从暴力方法，到初步方法，再优化最后得到正解，完整展现思考过程。
质数和约数
2302_77099705的博客
05-21 699
最大公约数（Greatest Common Divisor）指两个或多个整数共有约数中最大的一个。也称最大公因数、最大公因子，a， b的最大公约数记为（a，b），同样的，a，b，c的最大 公约数记为（a，b，c），多个 整数的最大公约数也有同样的记号。求最大公约数有多种 方法，常见的有 质因数分解法、 短除法、 辗转相除法、 更相减损法。
算法刷题笔记 逆序对的数量（C++实现）
hanmo22357的博客
05-18 311
AcWing上的算法基础课题目 逆序对的数量 刷题笔记记录
P1135 奇怪的电梯
loongnoy的博客
05-17 278
大楼的每一层楼都可以停电梯，而且第 𝑖i 层楼（1≤𝑖≤𝑁1≤i≤N）上有一个数字 𝐾𝑖Ki​（0≤𝐾𝑖≤𝑁0≤Ki​≤N）。例如： 3,3,1,2,53,3,1,2,5 代表了 𝐾𝑖Ki​（𝐾1=3K1​=3，𝐾2=3K2​=3，……），从 11 楼开始。第一行为三个用空格隔开的正整数，表示 𝑁,𝐴,𝐵N,A,B（1≤𝑁≤2001≤N≤200，1≤𝐴,𝐵≤𝑁1≤A,B≤N）。对于 100%100% 的数据，1≤𝑁≤2001≤N≤200，1≤𝐴,𝐵≤𝑁1≤A,B≤N，0≤𝐾𝑖≤𝑁0≤Ki​≤N。
第三十七天 | 860.柠檬水找零 406.根据身高重建队列 452.用最少数量的箭引爆气球
xxlyss的博客
05-18 300
其实你定义三个变量就好了呀，分别代表手上剩余的5，10，20面额纸币的张数，然后用for循环依次遍历bills，对三个变量进行运算，就可以模拟一次交易。当手上某一种纸币的数量不够找了，就false了。排序得到的结果：对于任意一个数据，它左边都是大于他的，右边都是小于他的。所以每一个元素前插，不会影响其他元素的 j 值，那就大胆前插，他的 j 是几就给他插到满意的位置。先将所有数据的身高（h）从大到小排序，使身高这个维度达到有序，然后再考虑前面有几个人的身高大于或者等于他。自己琢磨的七七八八了。
【算法】前缀和——二维前缀和模板题
最新发布
2302_79031646的博客
05-21 179
用“二维前缀和模板题”来介绍二位前缀和算法模板
【二叉树算法题记录】从中序与后序遍历序列构造二叉树，从前序与中序遍历序列构造二叉树
weixin_43627638的博客
05-15 382
这样我们就能确定整棵树的根节点在。
LeetCode 108.将有序数组转换为二叉搜索树
2301_80211119的博客
05-18 811
给你一个整数数组 nums ，其中元素已经按 升序 排列，请你将其转换为一棵 平衡 二叉搜索树。
LeetCode198：打家劫舍
m0_50127633的博客
05-18 150
你是一个专业的小偷，计划偷窃沿街的房屋。每间房内都藏有一定的现金，影响你偷窃的唯一制约因素就是相邻的房屋装有相互连通的防盗系统，如果两间相邻的房屋在同一晚上被小偷闯入，系统会自动报警。给定一个代表每个房屋存放金额的非负整数数组，计算你 不触动警报装置的情况下 ，一夜之内能够偷窃到的最高金额。
22. 括号生成
编程SHARE的博客
05-20 238
数字 n 代表生成括号的对数，请你设计一个函数，用于能够生成所有可能的并且 有效的 括号组合。
UUID: 0x2A01的数据结构
03-26
在蓝牙技术中，GATT（Generic Attribute Profile）用于定义设备之间的通信协议和数据格式。每个GATT特征都有一个唯一的UUID来标识它。 UUID: 0x2A01对应的是"Device Name"（设备名称）特征。它用于表示蓝牙设备的...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值