推荐的相关题目显示
题目描述
最近 lxhgww 又迷上了投资股票,通过一段时间的观察和学习,他总结出了股票行情的一些规律。
通过一段时间的观察,lxhgww 预测到了未来 T 天内某只股票的走势,第 i 天的股票买入价为每股 APi?,第 i 天的股票卖出价为每股 BPi?(数据保证对于每个 i,都有 APi?≥BPi?),但是每天不能无限制地交易,于是股票交易所规定第 i 天的一次买入至多只能购买 ASi? 股,一次卖出至多只能卖出 BSi? 股。
另外,股票交易所还制定了两个规定。为了避免大家疯狂交易,股票交易所规定在两次交易(某一天的买入或者卖出均算是一次交易)之间,至少要间隔 W 天,也就是说如果在第 i 天发生了交易,那么从第 i+1 天到第 i+W天,均不能发生交易。同时,为了避免垄断,股票交易所还规定在任何时间,
一个人的手里的股票数不能超过MaxP。
在第 1 天之前,lxhgww 手里有一大笔钱(可以认为钱的数目无限),但是没有任何股票,当然,T 天以后,lxhgww 想要赚到最多的钱,聪明的程序员们,你们能帮助他吗?
输入输出格式
输入格式:
输入数据第一行包括 3 个整数,分别是 T,MaxP,W。
接下来 T 行,第 i 行代表第 i?1 天的股票走势,每行 4 个整数,分别表示 APi?, BPi?, ASi?, BSi?。
输出格式:
输出数据为一行,包括 1 个数字,表示 lxhgww 能赚到的最多的钱数。
输入输出样例
/*
思路:
对于我这样的蒟蒻来说是一道难题
首先我们努力的列出一个DP式子
dp[i][j]表示第i天持有j支股票
那么我们便有四个决策
首先如果持有股量不超当日上界的情况直接购入(其实这种决策被第三种决策完全包括,反正你加上了也不影响正确性)
第二,我们不买也不卖,dp[i][j]=dp[i-1][j];
第三,我们买入x支股(x<=as[i]) dp[i][j]=max(dp[i][j],dp[i-m-1][j-y])(1<=y<=min(j,as[i]))
第四,我们出手x支股,与上种决策类似
这样的时间复杂度是O(TMaxPW) 期望50pts
然后,可以发现dp[i-m-1][j-y]的j-y为单调上升,这样我们可以用单调队列维护一下
然后,时间复杂度变成O(TMAXP)
应该就可以通过除了COJ以外一切测评机了
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<time.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
deque<pair<int,int> > Q,P;
int ap[2005],bp[2050],as[2020],bs[2050],dp[2050][2050],n,m,t;
int main(){
scanf("%d%d%d",&n,&t,&m);
m++;
rep(i,1,n) scanf("%d%d%d%d",&ap[i],&bp[i],&as[i],&bs[i]);
rep(i,1,t) dp[0][i]=-99999999;
rep(i,1,n){
if(i<=m){
rep(j,0,t){
dp[i][j]=dp[i-1][j];
if(j<=as[i] && dp[i][j]<-j*ap[i]) dp[i][j]=-j*ap[i];
//printf("dp[%d][%d]=%d\n",i,j,dp[i][j]);
}
continue;
}
while(!Q.empty()) Q.pop_back();
while(!P.empty()) P.pop_back();
rep(j,1,min(t,bs[i])){
while(!Q.empty() && Q.back().second<=dp[i-m][j]+j*bp[i]) Q.pop_back();
Q.push_back(make_pair(j,dp[i-m][j]+j*bp[i]));
}
rep(j,0,t){
dp[i][j]=dp[i-1][j];
if(j<=as[i] && dp[i][j]<=-j*ap[i]) dp[i][j]=-j*ap[i];
if(!Q.empty() && !P.empty()) dp[i][j]=max(dp[i][j],max(Q.front().second-j*bp[i],P.front().second-j*ap[i]));
else if(!Q.empty()) dp[i][j]=max(dp[i][j],Q.front().second-j*bp[i]);
else if(!P.empty()) dp[i][j]=max(dp[i][j],P.front().second-j*ap[i]);
while(!Q.empty() && Q.front().first==j+1) Q.pop_front();
while(!P.empty() && P.front().first+as[i]<j+1) P.pop_front();
while(!Q.empty() && j+bs[i]+1<=t && Q.back().second<=dp[i-m][j+bs[i]+1]+bp[i]*(j+bs[i]+1)) Q.pop_back();
while(!P.empty() && P.back().second<=dp[i-m][j]+ap[i]*j) P.pop_back();
P.push_back(make_pair(j,dp[i-m][j]+ap[i]*j));
if(j+bs[i]+1<=t) Q.push_back(make_pair(j+bs[i]+1,dp[i-m][j+bs[i]+1]+bp[i]*(j+bs[i]+1)));
}
}
printf("%d",dp[n][0]);
return 0;
}
