前言

斜率优化\(DP\)是难倒我很久的一个算法,我花了很长时间都难以理解。后来,经过无数次的研究加以对一些例题的理解,总算啃下了这根硬骨头。

基本式子

斜率优化\(DP\)的式子略有些复杂,大致可以表示成这样:

\[f_i=min_{j=1}^{i-1}(A(j)-B(j)*S(i)+C(i))\]

其中\(A(j)\)和\(B(j)\)是两个只与\(j\)有关的函数,\(S(i)\)和\(C(i)\)是两个只与\(i\)有关的函数,式子中的\(min\)其实也可以替换成\(max\),但这里以\(min\)为例。

不难发现,如果只有\(A(j)\)和\(C(i)\)两项,就是单调队列优化\(DP\)的基本式子了。

但是,由于式子中含有\(B(j)*S(i)\)这一项既与\(i\)相关,又与\(j\)相关的式子,就不能直接用单调队列,而要进行一定转化了。

考虑将\(A(j)\)移到等号左边,并将\(f_i\)移到等号右边,则原式可以转化成这样:

\[A(j)=B(j)*S(i)+(f_i-C(i))\]

注意,在\(i\)不变的时候,我们可以将只与\(i\)有关的项看成常数项。

于是,这个函数就可以看作一条直线,其中\(S(i)\)就相当于这条直线的斜率,而\(f_i-C(i)\)就相当于这条直线的截距

而\(C(i)\)是固定的,因此,如果要让\(f_i\)最小,则应让\(f_i-C(i)\)最小,对应到图像中就是让截距最小。

那么应如何让截距最小呢?

大致思路

首先,我们可以想象有一条斜率固定的直线(我太懒,不想画图... ...),然后图上有若干个点,现在要用这条直线从图的最下方往上慢慢移动,直至碰到第一个点,而这个点就是我们要找的最优点。

则不难发现,如果连续的三个点呈上凸状,则无论该直线斜率取多少,碰到的第一个点都不可能是中间这个点。

换句话说,就是中间这个点对答案没有任何贡献了。

于是就有一个策略:当我们要加入一个新的点时,比较当前点与前一个插入的点\(S_1\)、前一个插入的点与倒数第二个插入的点的斜率\(S_2\),如果\(S_1\le S_2\),则可将前一个插入的点弹出。

重复此操作,直至\(S_1>S_2\)或图上只剩一个点,然后将当前点插入。

如何求最优解

但是,这样一来,我们好像还是没能求出最优解。

此时又有两种操作方式:在凸包上二分单调队列维护最优解

对于某些问题,它可以确保决策单调性,即一个点如果当前不是最优解,则以后都不可能是最优解了。这样的问题可以直接用单调队列来维护最优解。

但有些问题却不一定满足这种性质,此时就需要在凸包上二分最优解了,但依然需要用单调队列来维护点集。

所以,如果你不会单调队列,最好赶紧去研究一下,然后再学习斜率优化\(DP\)。

几道例题

第一道例题: 【BZOJ2726】[SDOI2012] 任务安排

听说是入门题?洛谷上的弱化版可以直接单调队列维护,但是\(BZOJ\)上存在负数,需要在凸包上二分。

第二道例题: 【CF311B】Cats Transport

一眼看上去像是\(WQS\)二分,其实题意转换后可以直接斜率优化。

第三道例题: 【洛谷3648】[APIO2014] 序列分割

一开始觉得是区间\(DP\),结果一看数据范围,推了波性质,才发现其实可以用斜率优化\(DP\)来做。

动态规划专题(五)——斜率优化DP的更多相关文章

  1. 【学习笔记】动态规划—斜率优化DP(超详细)

    [学习笔记]动态规划-斜率优化DP(超详细) [前言] 第一次写这么长的文章. 写完后感觉对斜优的理解又加深了一些. 斜优通常与决策单调性同时出现.可以说决策单调性是斜率优化的前提. 斜率优化 \(D ...

  2. 斜率优化dp(POJ1180 Uva1451)

    学这个斜率优化dp却找到这个真心容易出错的题目,其中要从n倒过来到1的确实没有想到,另外斜率优化dp的算法一开始看网上各种大牛博客自以为懂了,最后才发现是错了. 不过觉得看那些博客中都是用文字来描述, ...

  3. HDU3507 print article【斜率优化dp】

    打印文章 时间限制:9000/3000 MS(Java / Others)内存限制:131072/65536 K(Java / Others) 总共提交:14521已接受提交:4531 问题描述 零有 ...

  4. bzoj-4518 4518: [Sdoi2016]征途(斜率优化dp)

    题目链接: 4518: [Sdoi2016]征途 Description Pine开始了从S地到T地的征途. 从S地到T地的路可以划分成n段,相邻两段路的分界点设有休息站. Pine计划用m天到达T地 ...

  5. bzoj-1096 1096: [ZJOI2007]仓库建设(斜率优化dp)

    题目链接: 1096: [ZJOI2007]仓库建设 Description L公司有N个工厂,由高到底分布在一座山上.如图所示,工厂1在山顶,工厂N在山脚.由于这座山处于高原内陆地区(干燥少雨),L ...

  6. [BZOJ3156]防御准备(斜率优化DP)

    题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3156 分析: 简单的斜率优化DP

  7. 【BZOJ-1096】仓库建设 斜率优化DP

    1096: [ZJOI2007]仓库建设 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 3719  Solved: 1633[Submit][Stat ...

  8. BZOJ 1010: [HNOI2008]玩具装箱toy 斜率优化DP

    1010: [HNOI2008]玩具装箱toy Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再 ...

  9. BZOJ 3156: 防御准备 斜率优化DP

    3156: 防御准备 Description   Input 第一行为一个整数N表示战线的总长度. 第二行N个整数,第i个整数表示在位置i放置守卫塔的花费Ai. Output 共一个整数,表示最小的战 ...

随机推荐

  1. C++ List的用法(转载)

    Lists将元素按顺序储存在链表中. 与 向量(vectors)相比, 它允许快速的插入和删除,但是随机访问却比较慢. assign() 给list赋值  back() 返回最后一个元素  begin ...

  2. 使用 PHPMailer 发送邮件出现诡异bug,间歇性发送失败

    场景 使用PHPMailer的SMTP发送邮件,用的是腾讯企业邮箱 smtp.exmail.qq.com 在邮箱设置里看到配置smtp方法 问题描述 本地windows开发环境发送邮件100%成功 远 ...

  3. Hadoop eclipse plugin

    我的eclipse是在win7上,hadoop在win7里的虚拟机里的ubuntu上,为了方便起见,想在eclipse上安装hadoop的插件,主要参考 https://my.oschina.net/ ...

  4. CF622F The Sum of the k-th Powers(拉格朗日插值)

    题意 给出 \(n,k\) , \(n\le10^9,k\le10^6\) ,求 \(\sum_{i=1}^n i^k(mod\;10^9+7)\) 题解 自然数幂次和,是一个\(k+1\)次多项式, ...

  5. Hadoop3.0 DataNode启动不成功——java.net.BindException: Port in use: localhost:0 Caused by: java.net.BindException: Cannot assign requested address解决办法

    一.问题出现的原因 启动Hadoop分布式环境时出现主节点的namenode.secondarynamenode启动成功,但是Worker节点datenode启动不成功. hadoop@master$ ...

  6. 分层图最短路【bzoj2662】[BeiJing wc2012]冻结

    分层图最短路[bzoj2662][BeiJing wc2012]冻结 Description "我要成为魔法少女!" "那么,以灵魂为代价,你希望得到什么?" ...

  7. PHP下载远程图片的几种方法总结

    1. 使用file_get_contents function dlfile($file_url, $save_to) { $content = file_get_contents($file_url ...

  8. 微信小程序 —— 微信小程序解析html富文本插件wxParse

    下载并把wxParse放到小程序的目录中 github下载地址:https://github.com/icindy/wxParse 一.基本使用方法 1. Copy文件夹wxParse,把wxPars ...

  9. STP-18-Port-Channl上的负载均衡

    Ether Channel通过在多条链路上传输多个数据帧,增加了可用带宽.一个以太网帧总是通过一个Ether Channel中的一条链路传输.针对数据帧地址字段执行散列计算能够产生一个编号,标识这个数 ...

  10. Java中如何实现代理机制(JDK动态代理和cglib动态代理)

    http://blog.csdn.net/skiof007/article/details/52806714 JDK动态代理:代理类和目标类实现了共同的接口,用到InvocationHandler接口 ...