题目

最近 \(\text{lxhgww}\) 又迷上了投资股票,通过一段时间的观察和学习,他总结出了股票行情的一些规律。

通过一段时间的观察,\(\text{lxhgww}\) 预测到了未来 \(T\) 天内某只股票的走势,第 \(i\) 天的股票买入价为每股 \(ap_i\),第 \(i\) 天的股票卖出价为每股 \(bp_i\)(数据保证对于每个 \(i\),都有 \(ap_i \geq bp_i\)),但是每天不能无限制地交易,于是股票交易所规定第 \(i\) 天的一次买入至多只能购买 \(as_i\) 股,一次卖出至多只能卖出 \(bs_i\) 股。

另外,股票交易所还制定了两个规定。为了避免大家疯狂交易,股票交易所规定在两次交易(某一天的买入或者卖出均算是一次交易)之间,至少要间隔 \(w\) 天,也就是说如果在第 \(i\) 天发生了交易,那么从第 \(i + 1\) 天到第 \(i + w\) 天,均不能发生交易。同时,为了避免垄断,股票交易所还规定在任何时间,一个人的手里的股票数不能超过 \(\text{MaxP}\)。

在第 \(1\) 天之前,\(\text{lxhgww}\) 手里有一大笔钱(可以认为钱的数目无限),但是没有任何股票,当然,\(T\) 天以后,\(\text{lxhgww}\) 想要赚到最多的钱,聪明的程序员们,你们能帮助他吗?

分析

设 \(f(i,j)\) 代表第 \(i\) 持有 \(j\) 股的最大收益。

\[f(i,j) = \max\left\{ \begin{array}{}
\max_{k = 1}^{as_i} f(i - w - 1,j - k) - ap_i\times k \\
\max_{k = 1}^{bs_i} f(i - w - 1,j + k) + bp_i\times k \\
f(i - 1,j)
\end{array} \right.
\]

很容易得出,该方程的时间复杂度为 \(\Theta (n^3)\)。显然会\(\rm T\)。

以只买入的 \(\max_{k = 1}^{as_i} f(i - w - 1,j - k) - ap_i\times k\) 看,我们发现我们要查询的 \(f\) 总是一段连续的区间,但是若使用单调队列维护 \(f\) ,\(ap_i\times k\) 又很恶心。

但是我们通过观察,对于每个 \(ap_i\times k\) 和 \(bp_i \times k\) 所影响到的要查询的 \(f\) ,若 \(j + 1\) ,则它们都会同时加上 \(ap_i\) 或 \(-bp_i\)。

比如原先的 \(f(i - w - 1, j) - ap_i,\ f(i - w - 1, j - 1) - 2ap_i,\cdots,f(i - w - 1, j + k) - ap_i\times k\)。

当区间右移一次时,它们的影响变成了:$f(i - w - 1, j - 1) - ap_i,\ f(i - w - 1, j - 2) - 2ap_i,\cdots $

所以,当要查询的区间右移一次时,它们的相对大小关系是不会改变的,我们可以对原 \(\rm dp\) 式进行变换来更好的运用这个性质。

以买入为例,\(\max_{k = 1}^{as_i} f(i - w - 1,j - k) - ap_i\times k\),若令 \(p = j - k\),则变为:

\[\begin{array}{}
& \max_{p = j - 1}^{j - as_i} f(i - w - 1,p) - ap_i\times (j - p) \\
= & \max_{p = j - 1}^{j - as_i} ( f(i - w - 1,p) + ap_i\times p ) - ap_i\times j
\end{array}
\]

这么一来,\(f\) 与 \(ap_i\) 就有了一一对应的关系,可以很容易的利用单调队列优化。

同理,卖出也可以这样处理,变式为 \(\max_{p = j + 1}^{j + bs_i} (f(i - w - 1,p) + bp_i\times p) - bp_i\times j\)。

接下来就是常规的单调队列优化,时间复杂度为 \(\Theta (n ^ 2)\)。

代码

#include <bits/stdc++.h>

using namespace std;

const int MAXN = 2000 + 5, inf = 0x3f3f3f3f;

int f[MAXN][MAXN], ap[MAXN], bp[MAXN], as[MAXN], bs[MAXN], w, t, maxp;

struct monQue {
pair<int, int> data[MAXN];
int head, tail; void clear() {head = 0, tail = 0;}
bool empty() {return tail <= head;}
void push(pair<int, int> x) {
while(!empty() && data[tail - 1] <= x) tail--;
data[tail++] = x;
}
void pop() {head++;}
pair<int, int> front() {return data[head];}
} sell, buy; int main() {
memset(f, ~0x3f, sizeof(f));
scanf("%d%d%d", &t, &maxp, &w);
for(int i = 1; i <= t; i++)
scanf("%d%d%d%d", &ap[i], &bp[i], &as[i], &bs[i]);
for(int i = 1; i <= t; i++) {
buy.clear(); sell.clear();
for(int j = 0; j <= maxp; j++) {
f[i][j] = max(f[i - 1][j], f[i][j]);
if(j <= as[i]) f[i][j] = max(-ap[i] * j, f[i][j]);
if(i - w - 1 <= 0) continue;
if(j) {
while(!buy.empty() && buy.front().second < j - as[i]) buy.pop();
buy.push(make_pair(f[i - w - 1][j - 1] + ap[i] * (j - 1), j - 1));
f[i][j] = max(buy.front().first - ap[i] * j, f[i][j]);
}
}
if(i - w - 1 <= 0) continue;
for(int j = maxp - 1; j >= 0; j--) {
while(!sell.empty() && sell.front().second > j + bs[i]) sell.pop();
sell.push(make_pair(f[i - w - 1][j + 1] + bp[i] * (j + 1), j + 1));
f[i][j] = max(sell.front().first - bp[i] * j, f[i][j]);
}
}
int ans = -inf;
for(int i = 0; i <= maxp; i++) {
ans = max(f[t][i], ans);
}
printf("%d\n", ans);
return 0;
}

后记

终于把这个坑填上了。。

【SCOI 2010】股票交易的更多相关文章

  1. [SCOI 2010] 股票交易

    [题目链接] https://www.lydsy.com/JudgeOnline/problem.php?id=1855 [算法] 单调队列优化动态规划 [代码] #include<bits/s ...

  2. SCOI 2010 连续攻击游戏(贪心,图论)

    SCOI 2010 连续攻击游戏 solution 直接就硬刚 我愿称贪心为暴力 因为题目中要求一定从小到大贪心,那么当前点的下标有能够选取的较大点,那么它一定可以和前面的一个较小点连接,所以可以直接 ...

  3. 【BZOJ 1857】【SCOI 2010】传送带

    三分套三分,虽然简单,但是也得掌握,,, 时间复杂度$O(log_{1.5}^2 n)$ 一开始WA好几次发现是快速读入里没有return,这样也能过样例?_(:3J∠)_ #include<c ...

  4. SCOI 2010 序列操作

    题目描述 lxhgww最近收到了一个01序列,序列里面包含了n个数,这些数要么是0,要么是1,现在对于这个序列有五种变换操作和询问操作: 0 a b 把[a, b]区间内的所有数全变成0 1 a b ...

  5. [SCOI 2010]传送带

    Description 题库链接 在一个 \(2\) 维平面上有两条传送带,每一条传送带可以看成是一条线段.两条传送带分别为线段 \(AB\) 和线段 \(CD\) .在 \(AB\) 上的移动速度为 ...

  6. [SCOI 2010]字符串

    Description lxhgww最近接到了一个生成字符串的任务,任务需要他把n个1和m个0组成字符串,但是任务还要求在组成的字符串中,在任意的前k个字符中,1的个数不能少于0的个数.现在lxhgw ...

  7. 解题:SCOI 2010 序列操作

    题面 线段树......模板题(雾? 然而两种标记会互相影响,必须保证每次只放一个(不然就不知道怎么放了),具体的影响就是: 翻转标记会使得覆盖标记一起翻转,下放的时候就是各种swap 覆盖标记会抹掉 ...

  8. Scoi 2010 幸运数字

    [题目描述]在中国,很多人都把6和8视为是幸运数字!lxhgww也这样认为,于是他定义自己的“幸运号码”是十进制表示中只包含数字6和8的那些号码,比如68,666,888都是“幸运号码”!但是这种“幸 ...

  9. 【SCOI 2010】传送带

    为了方便,我们不妨设$\rm P \lt Q,R$ 我们发现,有$\rm E$点在$\rm AB$上,$\rm F$点在$\rm CD$上,最优解一定是$\rm AE\rightarrow EF\ri ...

随机推荐

  1. Linux命令行环境与桌面环境护切换

    1.前言 在大部分情况下,我们在使用Linux时习惯使用命令行环境,但是有时候也还是会使用到安装桌面环境,所以在这里介绍一下如何给没有安装桌面环境的系统安装桌面环境.以Centos 6.5 为例演示一 ...

  2. System Center Configuration Manager 2016 域准备篇(Part2)

    对于" 服务器角色",请选择" Active Directory域服务",当系统提示您添加Active Directory域服务所需的功能时,请选择" ...

  3. Asp.net MVC 服务端验证多语言错误

    服务端验证用户提交信息时为了实现多语言使用了资源文件,如: using System.ComponentModel.DataAnnotations; public class UserModel { ...

  4. gitlab用户添加ssh免密钥认证后clone还是要求输入密码

    今天在centos 7公网服务器上安装gitlab在配置ssh免密钥时遇到一个奇怪的事,正确添加了本机的公钥到gitlab账户上,进行clone时死活都要你输入密码gitlab使用yum安装的,之前在 ...

  5. POJ 3734 Blocks (线性递推)

    定义ai表示红色和绿色方块中方块数为偶数的颜色有i个,i = 0,1,2. aij表示刷到第j个方块时的方案数,这是一个线性递推关系. 可以构造递推矩阵A,用矩阵快速幂求解. /*********** ...

  6. Why malloc+memset is slower than calloc?

    https://stackoverflow.com/questions/2688466/why-mallocmemset-is-slower-than-calloc/ The short versio ...

  7. object-detection-crowdai数据处理

    import os file=os.listdir('/home/xingyuzhou/object-detection-crowdai') file.sort(key= lambda x:int(x ...

  8. Bootstrap 弹出框(Popover)插件

    Bootstrap 弹出框(Popover)插件与Bootstrap 提示工具(Tooltip)插件类似,提供了一个扩展的视图,用户只需要把鼠标指针悬停到元素上面即可.弹出框的内容完全由Bootstr ...

  9. java算法面试题:写一个Singleton出来

    package com.swift; public class Singleton { public static void main(String[] args) { /* * 写一个Singlet ...

  10. 梁勇 java教材 编程练习题 第二章 2.6 键盘 读取一个在0 到 9999 之间的整数,并将该整数的各位数字相加。

    import java.util.Scanner; /** * 需求:从键盘输入任意0~9999之间的整数,输出这个整数各位数字之和. * 思路:方法一,使用扫描器Scanner类,扫描控制台输入流 ...