题目链接 http://www.rqnoj.cn/Problem_57.html
此题是二维费用背包
问题
二维费用的背包问题是指:对于每件物品,具有两种不同的费用;选择这件物品必须同时付出这两种代价;对于每种代价都有一个可付出的最大值(背包容量)。问怎样选择物品可以得到最大的价值。设这两种代价分别为代价1和代价2,第i件物品所需的两种代价分别为a[i]和b[i]。两种代价可付出的最大值(两种背包容量)分别为V和U。物品的价值为w[i]。
算法
费用加了一维,只需状态也加一维即可。设f[i][v][u]表示前i件物品付出两种代价分别为v和u时可获得的最大价值。状态转移方程就是:
f[i][v][u]=max{f[i-1][v][u],f[i-1][v-a[i]][u-b[i]]+w[i]}
如前述方法,可以只使用二维的数组:当每件物品只可以取一次时变量v和u采用逆序的循环,当物品有如完全背包问题时采用顺序的循环。当物品有如多重背包问题时拆分物品。这里就不再给出伪代码了,相信有了前面的基础,你能够自己实现出这个问题的程序。
物品总个数的限制
有时,“二维费用”的条件是以这样一种隐含的方式给出的:最多只能取M件物品。这事实上相当于每件物品多了一种“件数”的费用,每个物品的件数费用均为1,可以付出的最大件数费用为M。换句话说,设f[v][m]表示付出费用v、最多选m件时可得到的最大价值,则根据物品的类型(01、完全、多重)用不同的方法循环更新,最后在f[0..V][0..M]范围内寻找答案。
复数域上的背包问题
另一种看待二维背包问题的思路是:将它看待成复数域上的背包问题。也就是说,背包的容量以及每件物品的费用都是一个复数。而常见的一维背包问题则是实数域上的背包问题。(注意:上面的话其实不严谨,因为事实上我们处理的都只是整数而已。)所以说,一维背包的种种思想方法,往往可以应用于二位背包问题的求解中,因为只是数域扩大了而已。
小结
当发现由熟悉的动态规划题目变形得来的题目时,在原来的状态中加一纬以满足新的限制是一种比较通用的方法。希望你能从本讲中初步体会到这种方法。
当然此题还有一个时间最少的限制我们可以直接比较泡MM的个数,个数相同时在比较时间的多少。
#include <stdio.h>
#include <string.h>
#define MAX 1000
int n, rmb[MAX], rp[MAX], time[MAX], m, r;
int f[MAX][MAX], t[MAX][MAX];
int main()
{
memset(f, 0, sizeof(f));
memset(t, 0, sizeof(t));
scanf("%d", &n);
for(int i = 1; i <= n; i++)
{
scanf("%d%d%d", &rmb[i], &rp[i], &time[i]);
}
scanf("%d%d", &m, &r);
for(int i = 1; i <= n; i++)
{
for(int j = m; j >= rmb[i]; j--)
{
for(int k = r; k >= rp[i]; k--)
{
if(f[j][k] < f[j - rmb[i]][k - rp[i]] + 1)
{
f[j][k] = f[j - rmb[i]][k - rp[i]] + 1;
t[j][k] = t[j - rmb[i]][k - rp[i]] + time[i];
}
else
{
if(f[j][k] == f[j - rmb[i]][k - rp[i]] + 1)
{
if(t[j][k] > t[j - rmb[i]][k - rp[i]] + time[i])
{
t[j][k] = t[j - rmb[i]][k - rp[i]] + time[i];
}
}
}
}
}
}
printf("%d\n", t[m][r]);
}
我们还可以通过一个简单的转换不需要判断,把时间最少转化成时间最多。
#include <stdio.h>
#include <string.h>
#define MAX 105
int n, i, j, k, r, m, t;
int dp[MAX][MAX], rmb[MAX], rp[MAX], T[MAX];
int main()
{
scanf("%d", &n);
for(int i = 1; i <=n; i++)
{
scanf("%d%d%d", &rmb[i], &rp[i], &t);
T[i] = 100000 - t;
}
scanf("%d%d", &m, &r);
memset(dp, 0, sizeof(dp));
for(int i = 1; i <= n; i++)
{
for(j = m; j >= rmb[i]; j--)
{
for(k = r; k >= rp[i]; k--)
{
if(dp[j - rmb[i]][k - rp[i]] + T[i] > dp[j][k])
{
dp[j][k] = dp[j - rmb[i]][ k - rp[i]] + T[i];
}
}
}
}
printf("%d\n", 100000 - dp[m][r] % 100000);
return 0;
}
资料引自《背包九讲》