链接:
Fruit
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 2632 Accepted Submission(s): 1466
于是,很多人们慕名而来,找Lele买水果。
甚至连大名鼎鼎的HDU ACM总教头 lcy 也来了。lcy抛出一打百元大钞,"我要买由M个水果组成的水果拼盘,不过我有个小小的要求,对于每种水果,个数上我有限制,既不能少于某个特定值,也不能大于某个特定值。而且我不要两份一样的拼盘。你随意搭配,你能组出多少种不同的方案,我就买多少份!"
现在就请你帮帮Lele,帮他算一算到底能够卖出多少份水果拼盘给lcy了。
注意,水果是以个为基本单位,不能够再分。对于两种方案,如果各种水果的数目都相同,则认为这两种方案是相同的。
最终Lele拿了这笔钱,又可以继续他的学业了~
每组测试第一行包括两个正整数N和M(含义见题目描述,0<N,M<=100)
接下来有N行水果的信息,每行两个整数A,B(0<=A<=B<=100),表示至少要买该水果A个,至多只能买该水果B个。
题目数据保证这个答案小于10^9
2 3 1 2 1 2 3 5 0 3 0 3 0 3
2 12
算法:分组背包
思路:坚持一切母函数皆背包原则
P06: 分组的背包问题
问题
有N件物品和一个容量为V的背包。第i件物品的费用是c[i],价值是w[i]。这些物品被划分为若干组,每组中的物品互相冲突,最多选一件。求解将哪些物品装入背包可使这些物品的费用总和不超过背包容量,且价值总和最大。
算法
这个问题变成了每组物品有若干种策略:是选择本组的某一件,还是一件都不选。也就是说设f[k][v]表示前k组物品花费费用v能取得的最大权值,则有:
f[k][v]=max{f[k-1][v],f[k-1][v-c[i]]+w[i]|物品i属于第k组}
使用一维数组的伪代码如下:
for 所有的组k
for v=V..0
for 所有的i属于组k
f[v]=max{f[v],f[v-c[i]]+w[i]}
注意这里的三层循环的顺序,甚至在本文的beta版中我自己都写错了。“for v=V..0”这一层循环必须在“for 所有的i属于组k”之外。这样才能保证每一组内的物品最多只有一个会被添加到背包中。
另外,显然可以对每组内的物品应用P02中“一个简单有效的优化”。
小结
分组的背包问题将彼此互斥的若干物品称为一个组,这建立了一个很好的模型。不少背包问题的变形都可以转化为分组的背包问题(例如P07),由分组的背包问题进一步可定义“泛化物品”的概念,十分有利于解题。
来自《背包九讲》优化出错的根本原因是滚动数组的不完全覆盖!
大多数的背包问题最高维(一般就是物品编号)的每一次循环求解的第二维区间(一般就是[0,背包容量])都是相同的。
如果用滚动数组优化成一维,每一次循环之后当前数据会完全覆盖之前的数据,下一次循环时调用的数据就全部是更新了的数据。
但是本题的最高维(i)的每一次循环求解的第二维区间([A[i],M])却不尽相同!
如果用滚动数组优化成一维,每一次循环之后当前数据仅会部分覆盖之前的数据!这就造成了数据残留,下一次循环时调用的数据就可能是残留数据!
正确的做法应该是用滚动数组优化成二维
code:
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int maxn = 110;
int dp[maxn][maxn];
int a[maxn];
int b[maxn];
int nValue;
int main()
{
int n,m;
while(scanf("%d%d", &n,&m) != EOF)
{
memset(dp,0,sizeof(dp));
nValue = m;
dp[0][0] = 1;
for(int i = 1; i <= n; i++)
scanf("%d%d", &a[i], &b[i]);
for(int i = 1; i <= n; i++) // n 种对应于 n 组物品, 每组中的物品相互冲突
{
for(int j = m; j >= a[i]; j--) //对应背包容量
{
for(int k = a[i]; k <= b[i] && k <= j; k++)
dp[i][j] += dp[i-1][j-k]; // 前 i 组物品花费 j
}
}
printf("%d\n", dp[n][m]);
}
return 0;
}
看学弟总结后滚动数组优化后的代码:
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int maxn = 110;
int dp[2][maxn];
int a[maxn];
int b[maxn];
int nValue;
int main()
{
int n,m;
while(scanf("%d%d", &n,&m) != EOF)
{
memset(dp,0,sizeof(dp));
nValue = m;
dp[0][0] = 1;
for(int i = 1; i <= n; i++)
scanf("%d%d", &a[i], &b[i]);
for(int i = 1; i <= n; i++) // n 种对应于 n 组物品, 每组中的物品相互冲突
{
memset(dp[i%2], 0, sizeof(dp[i%2]));
for(int j = m; j >= a[i]; j--) //对应背包容量
{
for(int k = a[i]; k <= b[i] && k <= j; k++) //所有的 k 属于组数 i
dp[i%2][j] += dp[(i-1)%2][j-k];
}
}
printf("%d\n", dp[n%2][m]);
}
return 0;
}