凸壳优化dp/spole trick

最新推荐文章于 2023-08-29 15:32:30 发布

Flame♡

最新推荐文章于 2023-08-29 15:32:30 发布

阅读量554

点赞数

分类专栏：学习总结文章标签：动态规划算法机器学习

本文链接：https://blog.csdn.net/m0_50170681/article/details/121984732

版权

学习总结专栏收录该内容

19 篇文章 0 订阅

订阅专栏

写给自己的一点复习与总结

（主要是今天高强度盯着东西罚坐以及晚上和d*h battle了很久，所以来写点东西）

（话说，最近用onenote记东西习惯了，一下子有点不适应）

spole trick

当一些dp直接进行复杂度不可接受的时候，如果这个dp的斜率都是单调递增（或递减的），并且都是整数，那么我们可以考虑这么优化

如果我们将某个函数 $g (x)$ 的所有斜率改变1的点维护下来,(如果在某个点 $x$ 斜率的变化大于1，那么这个位置要记录多次)，那么我们就维护出来了这个函数

那么问题是？我们怎么求想要的函数值呢？

对于最小值dp问题来说，斜率最接近0的部分，应该就是所求答案。

对于函数图像，我们可以这么理解，斜率为0的一段上的dp值都相同（但我们此时还不知道是多少），因此，在所有能给这一段上的点更新答案的位置，使得更新出来的答案最小的，就是我们想要找的点（同时，因为要使得答案在最小范围内，用于更新答案的位置应该在这段上）

例1

CF713C Sonya and Problem Wihtout a Legend

首先先把每个位置的值减去 $i$ ，就转单调递增为单调不降

我们定义 $f_{i,j}$ 为确定到第 $i$ 个数，且第 $i$ 个数小于等于 $j$

我们定义 $g_{i,j}$ 为确定到第 $i$ 个数，且第 $i$ 个数等于 $j$

这样的定义上 $f_{i,j}=\min\limits_{k<=j}g_{i,k}$

$g_{i,j}=f{i,j}+|j-a_i|$

两个函数同时具有凹性，所以考虑用这个技巧维护

但问题来了？我们怎么统计答案呢？

贴个代码

#include<bits/stdc++.h>
using namespace std;
int n;
long long ans;
priority_queue<int> q;
int main(){
	int x; scanf("%d",&n);
	for(int i=1;i<=n;i++){
		scanf("%d",&x); x-=i;
		q.push(x); q.push(x);
		ans+=q.top()-x; q.pop();
	}
	printf("%lld",ans);
	return 0;
}