P2569 [SCOI2010]股票交易

题目描述

最近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能赚到的最多的钱数。

 

输入输出样例

输入样例#1：                   输出样例#1：

5 2 0                3
2 1 1 1                
2 1 1 1
3 2 1 1
4 3 1 1
5 4 1 1

 

说明

对于30%的数据，0<=W<T<=50,1<=MaxP<=50

对于50%的数据，0<=W<T<=2000,1<=MaxP<=50

对于100%的数据，0<=W<T<=2000,1<=MaxP<=2000

对于所有的数据，1<=BPi<=APi<=1000,1<=ASi,BSi<=MaxP

 

 

 

 

 

 

 

#include<iostream>
#include<cstdio>
#include<stdio.h>
#include<algorithm>
#include<cstring>
#include<string>
using namespace std;
const int maxn = 2010;
//const int inf = -2323232;
int t, maxp, w;
int f[maxn][maxn];
int head, till, que[10001];
struct days {
    int ap, bp, as, bs;
}day[maxn];
/*
int dp() {
    for (int i = 0; i <= t; i++) {
        for (int j = 0; j <= maxp; j++) {
            f[i][j] = inf;
        }
    }            
    memset(f, 128, sizeof(f));
    for (int i = 0; i <= t; i++) f[i][0] = 0;
    for (int i = 1; i <= t; i++) {
        for (int j = 0; j <= maxp; j++) {
            if (j <= day[i].as) {
                f[i][j] = max(f[i][j], -j * day[i].ap);
            }
            f[i][j] = max(f[i][j], f[i - 1][j]);
        }
        head = 1,till = 0;
        for (int j = 0; j <= maxp; j++) {
            while (head <= till && j - que[head] > day[i].as) {
                head++;
            }
            while (head <= till && f[i - w - 1][que[till]] + que[till] * day[i].ap <= f[i - w - 1][j] + j * day[i].ap) {
                till--;
            }
            que[++till] = j;
            if (head <= till) {
                f[i][j] = max(f[i][j], f[i - w - 1][que[head]] + que[head] * day[i].ap - j * day[i].ap);
            }
        }
        head = 1, till = 0;
        for (int j = maxp; j >= 0; j--) {
            while (head <= till && que[head] - j > day[i].bs) {
                head++;
            }
            while (head <= till && f[i - w - 1][que[till]] + que[till] * day[i].bp <= f[i - w - 1][j] + j * day[i].bp) {
                till--;
            }
            que[++till] = j;
            if (head <= till) {
                f[i][j] = max(f[i][j], f[i - w - 1][que[head]] + que[head] * day[i].bp - j * day[i].bp);
            }
        }
    }
    int ans = f[t][0];
    return ans;
}  */
int main() {
    scanf("%d %d %d", &t, &maxp, &w);
    for (int i = 1; i <= t; i++) {
        scanf("%d %d %d %d", &day[i].ap, &day[i].bp, &day[i].as, &day[i].bs);
    }
    //int aaa = dp();
    memset(f, 128, sizeof(f));
    for (int i = 0; i <= t; i++) f[i][0] = 0;
    for (int i = 1; i <= t; i++) {
        /*for (int j = 0; j <= maxp; j++) {
            if (j <= day[i].as) {
                f[i][j] = max(f[i][j], -j * day[i].ap);
            }
            f[i][j] = max(f[i][j], f[i - 1][j]);
        } */
        for (int j = 0; j <= day[i].as; j++) f[i][j] = -day[i].ap * j;
        for (int j = maxp; j >= 0; j--)  f[i][j] = max(f[i][j], f[i - 1][j]);
        if (i - w - 1 >= 0) {
            head = 1, till = 0;
            for (int j = 0; j <= maxp; j++) {
                while (head <= till && j - que[head] > day[i].as) {
                    head++;
                }
                while (head <= till && f[i - w - 1][que[till]] + que[till] * day[i].ap <= f[i - w - 1][j] + j * day[i].ap) {
                    till--;
                }
                que[++till] = j;
                if (head <= till) {
                    f[i][j] = max(f[i][j], f[i - w - 1][que[head]] + que[head] * day[i].ap - j * day[i].ap);
                }
            }
            head = 1, till = 0;
            for (int j = maxp; j >= 0; j--) {
                while (head <= till && que[head] - j > day[i].bs) {
                    head++;
                }
                while (head <= till && f[i - w - 1][que[till]] + que[till] * day[i].bp <= f[i - w - 1][j] + j * day[i].bp) {
                    till--;
                }
                que[++till] = j;
                if (head <= till) {
                    f[i][j] = max(f[i][j], f[i - w - 1][que[head]] + que[head] * day[i].bp - j * day[i].bp);    
                }
            }
        }
    }
    int ans = f[t][0];
    printf("%d\n", ans);
    return 0;
}