【NOIP2017提高A组模拟7.13】abcd

题目大意:

【NOIP2017提高A组模拟7.13】abcd_第1张图片
n200,25a[i]<b[i]25,1c[i]20,0d[i]100000

题解:

这其实是个多重背包问题。
只是这个物品的个数是a[i]..b[i],并且可以是负数。
其实是一样的,只需要把一开始的队列处理一下,队列出队的条件是一样的。
我们发现这样的常数有点大,所以需要卡卡。
Code:

#include 
#include 
#include 
#define fo(i, x, y) for(int i = x; i <= y; i ++)
#define max(a, b) ((a) > (b) ? (a) : (b))
using namespace std;

const int N = 100000, M = N + 600;

int o, n, a, b, w, v, l, r, f[2][N * 2 + 5000], d[2 * N + 5000];

void insert(int o, int x) {
    while(l <= r && f[o][x + M] >= f[o][d[r] + M] + (x - d[r]) / w * v) r --;
    d[++ r] = x;
}

int main() {
    freopen("a.in", "r", stdin);
    freopen("a2.out", "w", stdout);
    fo(i, 0, 1) fo(j, 0, N * 2 + 4000) f[i][j] = -1e9;
    f[o][M] = 0;
    scanf("%d", &n);
    fo(i, 1, n) {
        o = !o;
        scanf("%d %d %d %d", &a, &b, &w, &v);
        fo(mo, 0, w - 1) {
            l = 1, r = 0; int u = (-N - mo) / w, k = u * w + mo;
            fo(j, u - b, u - a) insert(!o, j * w + mo);
            f[o][k + M] = f[!o][d[l] + M] + (k - d[l]) / w * v;
            fo(j, u + 1, (N - mo) / w) {
                k += w;
                if(l <= r && k - d[l] > b * w) l ++;
                insert(!o, k - a * w);
                f[o][k + M] = f[!o][d[l] + M] + (k - d[l]) / w * v;
            }
        }
    }
    printf("%d", f[o][M]);
}

你可能感兴趣的:(动态规划)