POJ-2754 Similarity of necklaces 2 区间取下界操作+DP

这题一个很巧妙的操作就是对一个给定的区间同意将下界拿出来, 然后转化为多重背包问题. 使用二进制优化.

代码如下:

#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <iostream>
#define ER 0x80808080
using namespace std;

/*
解题:给定一个Multi数组, 一个Pairs数组, 要求出最好的Table的数组. 其中应满足这个几
     个要求, 首先sum{ Table[i]*Multi[i] } = 0, 并且要求sum{ Table[i]*Pairs }最大
     
解法:将下界单独拿出来作为一部分进行计算, 因此[L[i], U[i]] 就转化为[1, U[i]-L[i]]
     的一个多重背包. M[i]和P[i]均进行单独的下界计算.
     计算出 T = L[1]*M[1]+L[2]*M[2].... 之后, 就是一个关于容量T刚好放满的多重背包 
*/

int N, P[205], M[205], L[205], U[205];
int LIM, dp[100005];

void zobag(int w, int p) {
    for (int i = LIM; i >= w; --i) {
        if (dp[i-w] != ER)
            dp[i] = max(dp[i], dp[i-w] + p);    
    }
}

int DP() {
    int q, n, k;
    dp[0] = 0;
    for (int i = 1; i <= N; ++i) { // 枚举的是物品的编号
        k = 1;
        while (U[i]-k > 0) { // 如果还能够分出这么多件出来  
            zobag(k*M[i], k*P[i]);
            U[i] -= k;
            k <<= 1;
        }
        if (U[i])
            zobag(U[i]*M[i], U[i]*P[i]);
    }
    return dp[LIM];
}

int main() {
    int ret;
    while (scanf("%d", &N) != EOF) {
        LIM = ret = 0;
        memset(dp, 0x80, sizeof (dp));
        for (int i = 1; i <= N; ++i) {
            scanf("%d %d %d %d", &P[i], &M[i], &L[i], &U[i]);
            LIM += L[i] * M[i], U[i] -= L[i];
            ret += L[i] * P[i]; // P[i]同样要对下界进行一个计算, 最后再把其加回去 
        }
        LIM *= -1; // 一定会是一个不大于0的数, 因为这已经是取得最小值, 而题目要求最后结果为0 
        printf("%d\n", DP() + ret);
    }
    return 0;
}