首页 > 其他分享 >dp的优化(单队,斜率)

dp的优化(单队,斜率)

时间:2024-02-17 12:33:57浏览次数:31  
标签:单队 frac max times 斜率 dp sumf 单调

1.单调队列优化\(dp\)

维护最小值:\(x \leqslant q.tail()\)

维护最大值:\(x \geqslant q.tail()\)

其实原理不难,当\(dp\)的转移源头是一个区间时,往往使用单调队列来维护区间最值(一般队列里装下标以方便维护区间大小,但也只是一般情况),节省了处理区间的时间(甚至噶掉一维),重点是对区间的处理方式和细节,以及处理什么区间都因题而异,因此直接结合相关情景来进行学习

P2564 [SCOI2009] 生日礼物(例题)

(正片开始)在某两个城市之间有 $ n $座烽火台,每个烽火台都有一个代价。要求连续 $ m $个烽火台中至少要有一个发出信号。求总共最少的代价

设\(dp[i]\)表示\(1 \sim i\)中选了\(i\)的方案

那么为了合法,只能从\([i - m ,i - 1]\)中转移

用单调队列优化,最简单的一类

P1725 琪露诺 很像

类似的还有最大子序和,就是拿前缀和搞(记得提前往队列里塞一个\(0\))

P2627 [USACO11OPEN] Mowing the Lawn G

FJ有\(N\)只排成一排的奶牛,奶牛\(i\)的效率为\(E_i\)

计算没有连续的超过\(K\)只奶牛时可以得到的最大效率。

我们可以考虑从反面下手,找出不选的奶牛,将原来的序列分割成若干长度\(\leqslant k\)的区间

因此定义\(dp[i][1/0]\)表示选/不选第\(i\)头奶牛时的最大效率

如果不选,那么继承前一位的方案

\[dp[i][0] = \max(dp[i - 1][0],dp[i - 1][1]) \]

如果选了,那么在\([i - k,i]\)中必须有一个不选

\[dp[i][1] = \max(dp[j][0] + \sum\limits_{k = j + 1}^{i}E_k),j \in[i - k,i - 1] \]

拿前缀和优化:

\[dp[i][1] = \max(dp[j][0] +sum[i] - sum[j]) \]

对于枚举的\(j\),\(sum[i]\)是固定的,那么把他弄出去

\[dp[i][1] = \max(dp[j][0] - sum[j]) + sum[i] \]

用单调队列维护\(dp[j][0] - sum[j]\)即可

单调队列优化多重背包

一个初学背包时的禁忌领域

先摆出多重背包的基本递推式

\[dp[n] = \max(dp[n],dp[n - k \times v] + k \times w),k \in[1,s],k \times v\leqslant n \]

发现每次都是多减去一个\(v\),因此设\(n = pv + q\),那么

\[dp[n] = \max(dp[n],dp[i \times v + q]),i \in [0,\min(s,p)] \]

也就是说,\(dp[pv + q]\)只会从\(dp[i \times v + q]\)转移而来

后面的部分可以考虑使用单调队列,但别着急,还有\(w\)呢,对于不同的\(i\),\(w\)的系数也不一样

举个例子:

\[dp[3v + j] = \max(dp[j] + 3w,dp[j + v] + 2w,dp[j + 3w]) \]

考虑使用上一道题的方法,提出一个\(3w\):

\[dp[3v + j] = \max(dp[j],dp[j + v] - w,dp[j + 2v] - 2w,dp[j + 3v] - 3w) \]

这样队列维护的就是\(dp[iv + j] - iw,i \in [0,\min(s,p)]\)

wc太牛逼了

P3423 [POI2005] BAN-Bank Notes

共 \(n\) 种面值的硬币,为 \(b_1,b_2,\cdots,b_n\)。每种硬币有数量限制\(c_i\),现在我们想要凑出面值 \(m\),求最少要用多少个硬币

设\(dp[i]\)表示凑出金额\(i\)所需的最小硬币数

转移:\(dp[i] = \min(dp[i - j * b[l]] + j)\)

依照上述思路:维护\(dp[r + kb[i]] - k\)

一定要初始化啊啊啊啊啊啊啊

P2254 [NOI2005] 瑰丽华尔兹

描述略

第一想到理想的正方形

用两个单调队列维护两个维度

但不需要

这道题大抵是要用到\(dp\)的

\(dp[i][j]\)表示走到该格子的最大长

设一个方向上的时间区间长度为\(k\)

则(以向右为例)

\(dp[i][j] = \max(dp[i][l] + j - l + 1),l \in [j - k + 1,j - 1]\)

\(dp[i][j] = \max(dp[i][l] - l) + j + 1,l \in [j - k + 1,j - 1]\)

其他方向从略(向左/向上时,\(l \in[j + 1,j + k - 1]\),所以是\(dp[i][l]/dp[l][i] + l\),最后\(- j +1\))

维护:\(dp[i][l] +or- l\)

那就对于每次滑动,用单调队列这个窗口在每一行/列按这个方向滑动,当然会有无效操作,但没关系,因为顺序是对的,最终答案涉及的方格总能按照正确的顺序被修改和记录,所以无所谓,搞就完了

但是:

对于区间\([l,r]\),\(r - l + 1\)的计算方式包含两个端点,而观察样例可知:每次滑动的起点是不计入答案的

所以操作区间没有加一

方程:

\(dp[i][j] = \max(dp[i][l] + j - l),l \in [j - k + 1,j - 1]\)

\(dp[i][j] = \max(dp[i][l] - l) + j,l \in [j - k + 1,j - 1]\)

踢队头时的区间也是\(|q.front() - l|\)

md坑害我一天

P3089 [USACO13NOV] Pogo-Cow S

有\(N\)个目标点,目标点\(i\)在\(x_i\),该点得分为\(p_i\)。开始时可以选择站在一个目标点上,只允许朝一个方向跳跃,且每次跳跃的距离大于等于上一次跳跃的距离,并且必须跳到一个目标点,求所经目标点的最大得分和

首先,排序,方向未知,所以向左向各跑一遍

定义\(dp[i][j]\)表示跳到了目标点\(i\),并且是从\(j\)跳到\(i\)的最大和,得到关系:

\(dp[i][j] = \max(dp[j][k]) + p_i,|x_j - x_k| \leqslant |x_i - x_j|\)

\(dp[i][i] = p_i\)

先以向左为例,向右同理

那么\(0 \leqslant k < j < i \leqslant n\)

暴力是\(O(n^3)\)的,需要优化

肯定优化掉\(k\)那一维

我们采用单调队列的思想,既然已经按位置排了序,那么标号越小的点离\(j\)越远。利用这个天然的单调,从\(j\)开始倒着走找转移点,距离一旦不符合就不找了,节省时间

还有一个大坑:先枚举\(j\),再枚举\(i\)

如果不这样,其实相当于默认了起点是\(1\)

观察方程,可知每次转移相当于该方案继承了上一次跳跃的起点,所以必须先锁死起点,再去枚举终点

可见,很多时候单调队列的应用就是单调性的应用,从单调性入手遍历,可以少遍历很多非法情况,有时甚至不用建队,用单调性枚举即可

P2569 [SCOI2010] 股票交易

看到\(AP_i \geqslant BP_i\),知道了同一天内要么买,要么卖,不可能同时进行(否则必亏)

定义\(dp_{i,j}\)表示第\(i\)天有\(j\)张股票时最多能赚多少

分类讨论:

    • 在上一次交易的基础上买进
    • 白手起家(这天手头啥也没有)
    • 在上一次交易的基础上卖出
  • 啥也不干

四大情况,我们来一一分析

  • 在上一次交易的基础上买进

\[dp_{i,j} = \max(dp_{i,j},dp_{i - W - 1,j - k} - k \times AP_i),k \in [1,AS_i] \]

  • 白手起家:相当于直接赋值

\[dp_{i,j} = -AP_i \times k,k \in[1,AS_i] \]

  • 在上一次的基础上卖出

\[dp_{i,j} = \max(dp_{i,j},dp_{i - W - 1,j + k} + k \times BP_i) ,k\in [1,BS_i] \]

  • 啥也不干:\(dp_{i,j} = \max(dp_{i,j},dp_{i - 1,j}\))

暴力\(O(T\times MAXP \times \sum (AS_i + BS_i)),\),三次方级别,70\(pts\)

考虑使用单调队列优化

根据经验,肯定是维护一个\(dp_i + x_i\)的结构,但当前的状态定义没法子这样做(\(j + k,k\)和\(j - k,k\)都不好合并),我们改一下状态,以买进为例:

\[dp_{i,j} = \max(dp_{i,j},dp_{i - W - 1,k} - (j - k) \times AP_i),0 \leqslant j - k \leqslant AS_i \]

\[dp_{i,j} = \max(dp_{i - W - 1,k} + k \times AP_i) - j \times AP_i, \]

此时\(k,k\)对应,括号内的东西在遍历\(j\)时可以使用单调队列维护,以单调性维护队尾,以\(j - k\)的范围维护队头,能压到\(O(T \times MAXP)\)平方级别

可见状态的定义也非常重要 (服了)

2.斜率优化\(dp\)

可以认为是单调队列的拓展方法

观察发现,能用单调队列维护的\(dp\)方程,一定能分离出一个与当前枚举的\(i\)无关的量(单调队列的维护对象),但如果出现了形如\((i - j)^2 = i ^ 2 - 2ij + j ^ 2\)的东西时,\(2ij\)这个东西脚踏双船,\(i\)和\(j\)有一个变了他就变,无法使用普通队列来维护了

那么就针对他来下手

有\(n\)个数,每次取出一段数\([l,r]\)的花费是\((s[r] - s[l - 1])^2 + M\),\(M\)一定,求最小花费

设\(dp[i]\)表示取完\(1 \sim i\)的最小花费

可得:\(dp[i] = \min(dp[j] + (s[i] - s[j])^2 + M),j \in [1,i)\)

出现了麻烦项,如何处理

我们把式子这样写:

\(dp[i] = \min(dp[j] + s[j]^2 - 2s[i]s[j]) + s[i]^2 + M\)

\(dp[i] = \min(-2s[i](s[j]) + dp[j] + s[j]^2) + s[i]^2 + M\)

发现括号里是一个形如\(kx + b\)的式子,设\(k_i = -2s[i],x_j = s[j],b_j = dp[j] + s[j]^2\)

对于固定的\(i\),\(k_i\)是不变的,变得只有\(x_j,b_j\),而我们发现这两项只与\(j\)有关,似乎能用队列了

接下来看如何单调

如果想搞单调,即从\(j\)转移比从\(k\)转移更优,那么

\[dp[j] + (s[i] - s[j])^2 + M < dp[k] + (s[i] - s[k])^2 + M \]

化简得:

\[\frac{dp[j] + s[j]^2 - dp[k] - s[k]^2}{s[j] - s[k]} < 2s[i] = -k_i \]

设\(f[x] = dp[x] + s[x]^2\),则

\[\frac{f[j] - f[k]}{s[j] - s[k]} < -k_i \]

看到左边的形式想到小学的直线斜率

那么我们就以\((s[i],f[i])\)为点,代表该处状态

设现在枚举到了\(t\),有\(i,j,k\)三个备选状态,且状态如下

显然

\[\frac{f[j] - f[k]}{s[j] - s[k]} > \frac{f[i] - f[j]}{s[i] - s[j]} \]

对于\(2s[t] = -k_t\),有以下几种情况

  • \(\frac{f[j] - f[k]}{s[j] - s[k]}> \frac{f[i] - f[j]}{s[i] - s[j]}> -k_t\),那么\(j\)比\(i\)优,\(k\)比\(j\)优

  • \(\frac{f[j] - f[k]}{s[j] - s[k]}> -k_t > \frac{f[i] - f[j]}{s[i] - s[j]}\),那么\(i\)比\(j\)优,\(k\)比\(j\)优

  • \(-k_t > \frac{f[j] - f[k]}{s[j] - s[k]}>\frac{f[i] - f[j]}{s[i] - s[j]}\),那么\(i\)比\(j\)优,\(j\)比\(k\)优

也就是说,上凸的点一定不会是最优解,删掉上凸点后,剩下的点一定都往下凸,也就是一个下凸包,即所有线斜率单调不降

接下来,\(-k_t\)已知了,也就说待求直线的斜率已知,那么这根线与凸包最外侧的交点将是转移点(相关详情参考线性规划)

来看看交点的特征

显然,若\(k_a < k_b < \cdots < k_s < -k_t < \cdots\),则\(s\)点的状态就是转移点

那么用单调队列维护斜率,找到每个\(i\)对应的\(s\),转移即可

找\(s\):取出队首的两个元素\(a,b\),拿推出的式子比,小的为答案,大的就滚蛋

注:先维护队头更答案,再维护队尾塞进\(i\)

P3195 [HNOI2008] 玩具装箱

题目描述

P 教授有编号为 \(1 \cdots n\) 的 \(n\) 件玩具,第 \(i\) 件玩具经过压缩后的一维长度为 \(C_i\)。P 教授不关心容器的数目,他可以制作出任意长度的容器

要求:

  • 在一个一维容器中的玩具编号是连续的。

  • 如果将第 \(i\) 件玩具到第 \(j\) 个玩具放到一个容器中,那么容器的长度将为 \(x=j-i+\sum\limits_{k=i}^{j}C_k\)。

如果容器长度为 \(x\),其制作费用为 \((x-L)^2\)。其中 \(L\) 是一个常量。但他希望所有容器的总费用最小。

设\(dp[i]\)表示前\(i\)个玩具装箱的最小代价,那么

\[dp[i] = \min_{0 \leqslant j < i}(dp[j] + (\sum\limits_{k = j + 1}^{i}C_k + i - j - 1 - L)^2) \]

前缀和助阵

\[dp[i] = \min_{0 \leqslant j < i}(dp[j] + (s[i] - s[j] + i - j - 1 - L)^2) \]

再设\(f[i] = s[i] + i\)(可预处理)

\[dp[i] = \min_{0 \leqslant j < i}(dp[j] + (f[i] - f[j] - 1 - L)^2) \]

按鞋油推:

\[\frac{dp_j + (f_j + L + 1)^2 - dp_k - (f_k + L + 1)^2}{f_j - f_k} < 2f_i \]

设\(h_i = dp_i + (f_i + L + 1)^2\)

\[\frac{h_j - h_k}{f_j - f_k} < 2f_i \]

上板子就好了

[APIO2010] 特别行动队

设\(dp_i\)表示\(1 \sim i\)的士兵经修正后的最大战力

那么

\[dp[i] = \max_{0 \leqslant j < i}(dp[j] + a(s[i] - s[j])^2 + b(s[i] - s[j]) + c) \]

再套:

\[\frac{dp_j + as_j^2 - bs_j - (dp_k + as_k^2 -bs_k)}{s_j - s_k} > 2as_i \]

设\(f_i = dp_i + as_i^2 - bs_i\)

\[\frac{f_j - f_k}{s_j - s_k} > 2as_i \]

发现这次是大于号,采用类似的方法,假设是下凸包,发现下凸点是无法成为最优解的,所以这次维护的是上凸包

上板子

P2120 [ZJOI2007] 仓库建设

设\(dp[i]\)表示在\(i\)处建厂的最小费用

对于\(0 \leqslant j < i\),由定义的状态,可知\(j\)处已经建厂,再结合题目所说,只有\(j + 1 \sim i\)的货物会流入仓库\(i\),所以

\[dp_i = \min\limits_{0 \leqslant j < i}(dp_j + \sum\limits_{k = j + 1}^{i}p_k(x_i - x_k)) + c_i \]

考虑处理\(\sum\limits_{k = j + 1}^{i}p_k(x_i - x_k)\)

\[\sum\limits_{k = j + 1}^{i}p_k(x_i - x_k) = x_i\sum\limits_{k = j + 1}^{i}p_k - \sum\limits_{k = j + 1}^{i}p_kx_k \]

\(\sum\limits_{k = j + 1}^{i}p_k\),\(\sum\limits_{k = j + 1}^{i}p_kx_k\)用前缀和处理,记为\(s_i - s_j,t_i - t_j\)

\(dp_i = \min\limits_{0 \leqslant j < i}^{}(dp_j + x_i(s_i - s_j) + t_i - t_j) + c_i\)

\(dp_i = \min\limits_{0 \leqslant j < i}(dp_j - x_is_j - t_j) + x_is_i + t_i + c_i\)

维护\(dp_j - x_is_j-t_j\)

套板子

\[\frac{dp_j - t_j - (dp_k - t_k)}{s_j - s_k} < x_i \]

设\(f_i = dp_i - t_i\)

\[\frac{f_j - f_k}{s_j - s_k} < x_i \]

搞就完了

但这题有一堆细(da)节(bian):

  • 若序列中间存在\(p_i = 0\),会导致除数为零,要\(return\) 极大/极小值来代替

  • 中间的工厂没货的话就不可能成为仓库,仓库应该建在最后一个有货的工厂

  • 末尾可能会有若干没货的工厂,要从中取一个最小值

P2900 [USACO08MAR] Land Acquisition G

设长为\(x_i,\)宽为\(y_i\)

我们先按\(x\)从小到大排序,再让\(y\)从大到小(保证每次向前枚举的时候符合题意),同时去掉那些被大土地完全包含的小土地

\(dp_i\)表示前\(i\)块土地的最小花费,则

\(dp_i = \min\limits_{1 \leqslant j < i}(dp_j+ x_i \times y_{j + 1})\)

\[\frac{dp_j - dp_k}{y_{j + 1} - y_{k + 1}} < - x_i \]

小于号,下凸包?

发现这次的右边成为了一个负数,是单调减的,就好像是这样:

上凸包 (因为6翻了)

可见任何部分的单调性都对凸包是有影响的

还有一点就是队列内元素个数的问题,考虑到一个点就动用了\(k,k + 1\)两个,所以至少留三个(重一个)

P2365 任务安排

设\(dp_i\)表示前\(i\)个任务的最短完成时间

\[dp_i = \min(dp_j + (sumt_i + num \times s) \times (sumf_i - sumf_j)),j \in [0,i - 1] \]

\(sumt,sumf\)为时间,费用的前缀和,\(num\)表示分组数

看到式子里有个\(num\),考虑把它搞掉(事实上,这个东西不能随着方程转移)

化简式子

\[dp_i = \min(dp_j + sumt_i \times (sumf_i - sumf_j) + num \times s \times (sumf_i - sumf_j)) \]

考虑搞掉\(num\times s \times (sumf_i - sumf_j)\)

根据\(dp\)的尿性,这个\(num\)也不符合\(dp\)的风格(累计),所以考虑对于已知的分组,把该状态分组的影响全部并到此时的方程里一并计算,就可以秒掉\(num\)了

方程中已知的分组就是\(j + 1 \sim i\)


如图明显可得:分组后时间上多了一个\(s\)(\(i+ 1\)处)

代价增加\((sumf_n - sumf_i) \times s\)

但考虑到\(num\)还包含了这一次(\(j + 1\)处)的启动,所以影响其实是

\((sumf_i - sumf_j) * s + (sumf_n - sumf_i) * s = (sumf_n - sumf_j) * s\)

其实就是相对于\(j + 1 \sim n\)不分组(一次也不启动)的增加值

所以

\[dp_i = \min(dp_j + sumt_i \times (sumf_i - sumf_j) + (sumf_n - sumf_j) \times s) \]

\(O(n^2)\),可过弱化版

接下来到了激动人心的优化环节

\[dp_i = \min(dp_j -(sumt_i+s) \times sumf_j) + sumt_i \times sumf_i + sumf_n \times s \]

\(\cdots\)

\[\frac{dp_j - dp_k}{sumf_j - sumf_k} < sumt_i+s \]

下凸包

然后就来了个不要脸的

P5785 [SDOI2012] 任务安排ProMax

\(|T_i| ?\)

wc有负数?

\(sumt\)单调性不存在了

\(dp\)不变,下凸包不变,麻烦就麻烦在没有单调性

也就是说:

  • 弹掉对头的依据没有了

对于第一种情况,我们就不能弹掉队头了,答案也不一定是队头了,要二分查找了

物理学不存在了!(悲伤)(扭曲)(猴叫)(爬行)(蛇了)

标签:单队,frac,max,times,斜率,dp,sumf,单调
From: https://www.cnblogs.com/MLP123/p/18017878

相关文章

  • 树状DP心得
    说一下近日所学的主要两种题型,一个是分叉情况问题,一种是树上背包问题。分叉情况问题具体的题可以参考小胖守皇宫和三色二叉树。点击查看代码小胖守皇宫#include<bits/stdc++.h>usingnamespacestd;constintN=2000;vector<int>son[N];intfa[N],h[N],f[N][3];//f[0]......
  • DPInst64.exe difxapi64.dll 是什么 为什么
    DPInst64.exe:安装和卸载驱动程序包。默认情况下,该工具将搜索当前目录,并尝试安装找到的所有驱动程序包。用法:C:\Users\Administrator\Desktop\dp\dp\DPInst64.exe[/UINF文件][/S|/Q][/LM][/P][/F][/SH][/SA][/A][/PATH路径][/EL][/L语言ID][/C][/D][/LogTitle标题][/SW][/?......
  • 空间转移(dp+倍增+Floyd)
    第2题   空间转移 查看测评数据信息给定一个n个点,m条边的有向图,编号从1到n,所有边权值是1,现在有一个动点从点1开始一动,其每秒钟可以移动2的k次方千米(k是任意自然数,且2的k次方不超过1000000000)。最少需要几秒才能到达n号点。数据保证1到n至少有一条路径。n⩽50,m⩽1000......
  • ThreadPoolTaskExecutor以及通过注解实现异步任务
    ThreadPoolTaskExecutor是Spring框架的线程池,实现方式如下:1//声明一个name为asyncTaskExecutor的线程池bean到容器中2@Bean("asyncTaskExecutor")3publicExecutorgetAsyncExecutor(){4ThreadPoolTaskExecutorthreadPoolExecutor=newThreadPoolTaskExecuto......
  • 状压 dp 与 插头 dp
    矩阵上的dp:按格dp/轮廓线dp设\(f[i,j,S]\)表示考虑到第\(i\)行第\(j\)个格子,轮廓线上所有的格子的状态。复杂度为\(\Theta(nm|S|)\)。按行dp\(f[i,S]\)表示选完前\(i\)行的合法方案数。总复杂度为\(\Theta(n|S|^2)\)。P3226[HNOI2012]集合选数给定集合......
  • [学习笔记]换根 DP 的常用处理方式
    [学习笔记]换根DP的常用处理方式换根DP,又称作二次扫描法,通常用于“对每个根求出树形DP的答案”。以每个点作为根节点进行一次树形DP的代价通常无法承受,因此我们会使用两次DFS:第一次DFS指定一个点为根节点,运行一次常规的树形DP。第二遍DFS进行换根DP,得到将根转移......
  • DP 做题记录
    复杂DP各种巨大DP。CF123CBrackets题意:括号数组是一个只有“(”或“)”两类字符的二维数组。括号数组中的合法路径只能从任意位置开始,向右或向下移动。如果一个n×m括号数组中从(1,1)到(n,m)的所有路径经过的字符构成的字符串均为可以完全匹配的括号序列,则这个括号......
  • DP 优化
    DP优化一些典型或者不典型的DP优化实例。CF354DTransferringPyramid题意:有一个\(n\)层的金字塔,现在能进行两种操作,一是给某个点染色,代价是3,或者画一个底边是金字塔底边的子三角形,给每个点,代价是点数+2,要求用最小代价覆盖所有要染色的\(m\)个点。思路:定义\(h_i\)表......
  • 【无评级杂题】DP/贪心
    题目这题后来看了看网上的思路,发现贪心就能做,亏我还写了个O(2*N)的DP...浪费时间了属于是my-code#include<iostream>#include<cstring>#include<string>#include<cstdio>#include<cstdlib>#include<algorithm>#include<stack>#include<queue>......
  • 树形 DP
    简单的树上dp其实已经在普及组涉及过:自上而下和自下而上传递的性质。现在我们需要研究更复杂的树上dp,比如换根dp等等。【树上dp】最大子树和给出一棵带点权的树,求这棵树中的最大权连通块。因为是无根树,我们人为规定1号结点为根。\(dp[i]\)表示以\(i\)为根的子树......