【笔记篇】单调队列优化dp学习笔记&&luogu2569_bzoj1855股票交♂易
DP颂
DP之神 圣洁美丽 算法光芒照大地
我们怀着 崇高敬意 跪倒在DP神殿里
你的复杂 能让蒟蒻 试图入门却放弃
在你光辉 照耀下面 AC真心不容易
dp大概是最经久不衰 亘古不化的算法了吧.
而且有各种各样的类型 优化之类的.
一直dp都不怎么好. 而且也不太知道应该怎么提高.
基本见到不认识的dp方程就不大会推(但我会打表啊= =
所以dp还是很有的学的~
正好最近刚刚肝了计算几何, 所以就顺带搞一下斜率优化dp一类的...
单调队列优化dp
单调队列大家都会吧?
不会的先出去学一下, 这里不讲.
好的, 我们来看一下这个柿子
\]
其中\(\omega(j)\)是一个费用函数, 一般会根据题目的不同而变化.
这个dp能做到什么复杂度呢?
首先一眼\(O(n^2)\)...
然而我们可以用一个变量记录一下之前出现过的最大值.
这样转移是\(O(1)\)的了, 总复杂度就降到了\(O(n)\).
但是如果是这样呢?
\]
那就不能只维护一个变量了, 因为最大值如果出现在\(j\)的取值区间之外则转移是不合法的.
这样我们就考虑用单调队列来维护最大值, 这样转移依然可以做到\(O(1)\), 总复杂度\(O(n)\).
看道题: (woc辣鸡bzoj给的什么zz数据范围, T都没给怎么做...)
这题可以写出这么一个状态转移方程
f[i][j]=max\left\{
\begin{matrix}
f[i-w-1][k]-ap[i]*(j-k), (k\in[j-as[i],j)) //买入\\
f[i-w-1][k]+bp[i]*(k-j),(k\in(j,j+bs[i]]) //卖出\\
f[i-1][j]//不交♂易
\end{matrix}
\right.
\]
其中不交易的情况好处理, 但是如果前面两种枚举\(k\)的话就要做到\(O(n^*maxP^2)\), 显然是过不了的, 我们必须考虑优化.
我们以买入为栗化一波柿子(因为卖出同理) :
=max\{f[i-w-1][k]+ap[i]*k\}-ap[i]*j (k\in[j-as[i],j])
\]
我们令\(\omega(x)=ap[i]*x\), 而我们枚举\(i\), 就可以视为\(i\)是定值, 于是\(ap[i],as[i]\)都是定值.
我们就可以看出第二维形成了一个能用单调队列优化的柿子了.
这样优化之后复杂度成功降到了\(O(n*maxP)\), 就可以通过此题了.
根据贪心原则, 为了获得最多的现金, 手里不应该留股票, 所以用每个\(f[i][0]\)更新答案即可.
不过要注意一下边界条件... 挺扯淡的..
代码:
#include <cstdio>
#include <cstring>
const int N=2020;
const int INF=0x7fffffff;
int q[N<<1],d[N<<1],h=1,t=0;
int f[N][N],ap[N],bp[N],as[N],bs[N];
inline int gn(int a=0,char c=0){
for(;c<'0'||c>'9';c=getchar());
for(;c>47&&c<58;c=getchar())a=a*10+c-48;return a;
}
inline int max(const int& a,const int& b){
return a>b?a:b;
}
int main(){
int n=gn(),m=gn(),w=gn(),ans=0;
for(int i=1;i<=n;++i)
ap[i]=gn(),bp[i]=gn(),as[i]=gn(),bs[i]=gn();
memset(f,192,sizeof(f));
for(int i=1;i<=n;++i){
for(int j=0;j<=as[i];++j) f[i][j]=-ap[i]*j;
for(int j=0;j<=m;++j) f[i][j]=max(f[i][j],f[i-1][j]);
if(i>w){
h=1; t=0;
for(int j=0;j<=m;++j){
int val=f[i-w-1][j]+j*ap[i];
while(t>=h&&val>=q[t]) --t;
q[++t]=val; d[t]=j;
while(t>=h&&d[h]<j-as[i]) ++h;
f[i][j]=max(f[i][j],q[h]-ap[i]*j);
}
h=1; t=0;
for(int j=m;j>=0;--j){
int val=f[i-w-1][j]+bp[i]*j;
while(t>=h&&val>=q[t]) --t;
q[++t]=val; d[t]=j;
while(t>=h&&d[h]>j+bs[i]) ++h;
f[i][j]=max(f[i][j],q[h]-bp[i]*j);
}
}
ans=max(ans,f[i][0]);
}
printf("%d",ans);
}
反正差不多就这样吧....
【笔记篇】单调队列优化dp学习笔记&&luogu2569_bzoj1855股票交♂易的更多相关文章
- 「学习笔记」单调队列优化dp
目录 算法 例题 最大子段和 题意 思路 代码 修剪草坪 题意 思路 代码 瑰丽华尔兹 题意 思路 代码 股票交易 题意 思路 代码 算法 使用单调队列优化dp 废话 对与一些dp的转移方程,我们可以 ...
- POJ - 1821 单调队列优化DP + 部分笔记
题意:n个墙壁m个粉刷匠,每个墙壁至多能被刷一次,每个粉刷匠要么不刷,要么就粉刷包含第Si块的长度不超过Li的连续墙壁(中间可不刷),每一块被刷的墙壁都可获得Pi的利润,求最大利润 避免重复粉刷: 首 ...
- 算法笔记--单调队列优化dp
单调队列:队列中元素单调递增或递减,可以用双端队列实现(deque),队列的前面和后面都可以入队出队. 单调队列优化dp: 问题引入: dp[i] = min( a[j] ) ,i-m < j ...
- Parade(单调队列优化dp)
题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=2490 Parade Time Limit: 4000/2000 MS (Java/Others) ...
- 【bzoj1855】 [Scoi2010]股票交易 单调队列优化DP
上一篇blog已经讲了单调队列与单调栈的用法,本篇将讲述如何借助单调队列优化dp. 我先丢一道题:bzoj1855 此题不难想出O(n^4)做法,我们用f[i][j]表示第i天手中持有j只股票时,所赚 ...
- 1855: [Scoi2010]股票交易[单调队列优化DP]
1855: [Scoi2010]股票交易 Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1083 Solved: 519[Submit][Status] ...
- bzoj1499: [NOI2005]瑰丽华尔兹&&codevs1748 单调队列优化dp
这道题 网上题解还是很多很好的 强烈推荐黄学长 码风真的好看 神犇传送门 学习学习 算是道单调队列优化dp的裸题吧 #include<cstdio> #include<cstring ...
- 洛谷P3195 [HNOI2008]玩具装箱TOY(单调队列优化DP)
题目描述 P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再放到一种特殊的一维容器中.P教授有编号为1...N的N件玩具, ...
- 单调队列优化DP,多重背包
单调队列优化DP:http://www.cnblogs.com/ka200812/archive/2012/07/11/2585950.html 单调队列优化多重背包:http://blog.csdn ...
随机推荐
- 自动化监控系统(二)连接数据库,创建app,添加model,同步数据库
数据库我使用:mysql5.7 程序连接数据库的模块:pymysql 一.创建数据库: dbname:automatedmonitor 二.使用pip安装pymysql,这里我直接在pycharm上安 ...
- android Intent和IntentFilter
android的应用程序包含三种重要的组件:Activity.Service.BroadcastReceiver,应用程序采用一致的方式来启动他们——都是依靠Intent来进行启动.Intent就封装 ...
- C#实体类克隆
public static T Clone<T>(T source) { if (!typeof(T).IsSerializable) { throw new ArgumentExcept ...
- 2018-6-17-win10-UWP-全屏
title author date CreateTime categories win10 UWP 全屏 lindexi 2018-06-17 17:51:19 +0800 2018-2-13 17: ...
- 标准 IO fprintf 与 sprintf 函数使用
函数原型 fprintf int fprintf(FILE *stream, const char *format, ...); 把数据写到流中 int sprintf(char *str, con ...
- grunt中调用webpack
项目中一些很老的技术需要运用webpack管理,这时候必须使用grunt+webpack 配合使用 .做个笔记,不限多说直接上代码 grunt结合webpack 下面是 Gruntfile.js co ...
- delphi xe10 安卓设备信息
//引用单元 FMX.Android.DeviceInfo.GetInformation; Memo1.Lines.Add('ID:'+FMX.Android.DeviceInfo.ID); Memo ...
- 【NOI2019模拟2019.7.1】三格骨牌(轮廓线dp转杨图上钩子定理)
Description \(n,m<=1e4,mod ~1e9+7\) 题解: 显然右边那个图形只有旋转90°和270°后才能放置. 先考虑一个暴力的轮廓线dp: 假设已经放了编号前i的骨牌,那 ...
- TopCoder[TCO2016 Round 1A]:EllysTree(1000)
Problem Statement Elly has a graph with N+1 vertices, conveniently numbered from 0 to N. The gr ...
- 字符串dp——牛客多校第五场G
比赛的时候脑瘫了没想出来..打多校以来最自闭的一场 显然从s中选择大于m个数组成的数必然比t大,所以只要dp求出从s中选择m个数大于t的方案数 官方题解是反着往前推,想了下反着推的确简单,因为高位的数 ...