九大背包问题专题--多重背包问题(二进制优化方法;单调队列问题)

3.多重背包问题1

题目：
有N件物品和一个容量是V的背包。
第i种物品最多有si件，每件的体积是vi，价值是wi。
求解将哪些物品装入背包，可使这些物品的总体积不超过背包的容量，且价值总和最大。

输入格式
第一行有两个整数，N,V用空格隔开，分别表示物品种数和背包容积。
接下来有N行，每行三个整数vi，wi，si，用空格隔开，分别表示第i种物品的体积、价值和数量。

输出格式
输出一个整数，表示最大价值

数据范围
0<N.V<=100
0<vi,wi，si<=100

输入样例
4 5
1 2 3
2 4 1
3 4 3
4 5 2

输出样例
10

分析思路：

时间复杂度O(n^3)

f[i]:表示总体积是i的情况下，最大价值是多少

每件物品选法，假设有s个，则选法有s+1种选法。可以选0,1,2,3…m个。

01背包：

从前往后考虑
for(int i=0;i<n;i++){
    for(int j=m,j>=v[i];j--)
       f[j]=max(f[j]，f[j-v[i]]+w[i])

选的话为f[j-v[i]]，不选的话为f[j]不变

多重背包
可以看做01背包的扩展

状态转移加一种循环

选法有很多种，0,1,2，…s

f[j]=max(f[j]，f[j-v[i]]+w[i]，f[j-2*v[i]]+w[i]...);

第一种情况，在初始化的时候，把所有f[i]都初始化为0，答案就是f[m]
f[i]=0;

第二种情况：在初始化时把f[0]=0,其余f[i]都初始化为负无穷，枚举0~m求最大值。
f[0]=0;
f[i]=-INF；（i!=0）
max{f[0....m]}

代码：

#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;

const int N=110;
int n,m; 
int f[N];
int main(){
	cin>>n>>m;
	for(int i=1;i<n;i++){  //枚举所有物品 
	int v,w,s; 
	cin>>v>>w>>s;
	 for(int j=m;j>=0;j--)  // 从大到小枚举 
	    for(int k=1;k<=s&&k*v<=j;k++)  //状态转移，s+1种选法，判断枚举个数不能大于j 
	   f[j]=max(f[j],f[j-k*v]+k*w);  
	   cout<<f[m]<<endl;
}
	   return 0;
} 

---

二进制优化方法

多重背包问题2

题目：
问题：
有N件物品和一个容量是V的背包。
第i种物品最多有si件，每件的体积是vi，价值是wi。
求解将哪些物品装入背包，可使这些物品的总体积不超过背包的容量，且价值总和最大。

输入格式
第一行有两个整数，N,V用空格隔开，分别表示物品种数和背包容积。
接下来有N行，每行三个整数vi，wi，si，用空格隔开，分别表示第i种物品的体积、价值和数量。

输出格式
输出一个整数，表示最大价值

数据范围
0<N<=1000
0<V<=1000
0<vi,wi，si<=2000

输入样例
4 5
1 2 3
2 4 1
3 4 3
4 5 2

输出样例
10

分析思路:

把多重背包转化为01背包
假设v,w，s
假设物品有s个，把物品拆成s份(重复s份)，转化后变为01物品，每个物品最多只能用一次，
其复杂度很高 10^9

二进制优化方法

假设质量为，用最少的物品的个数情况下表示出来
每个数两种选法，把0~7八种方案表示出来，至少需要3个数log以2为底8的对数，
乘法原理:8种方案，选3个数，2^3=8
下界一定为3，往上取整，注意3个数能不能满足要求
0
1
3=1+2
4=4
5=1+4
6=2+4
7=1+2+4

给定任意数s，最少把s可以分成多少个数（划分的数，两种选择选或者不选）拼成<=s.
log2(s)往上取整
但是遇到不是2的整数幂，无法算出整数。

举例:10
log2(10)上取整，4
1 2 4 8
（这样取的话不止可以表示0~10，还可以表示0 ~15，不能选，一共只有10个物品）
让s一直减，先减1再减2，减4，减8，一直减到负数为止，最后剩下的数直接放入
上式变为：
1 2 4 3（1+2+4=7，剩余3直接放入）
这4个数就一定可以表示0~10，
原因:1-4可以表示0~7，
3和0~7可以表示0 ~10

常规情况：
s-1 -2 -4 -8
减到为负为止，把剩余数，直接取出来，这个数和前面2的整数次幂组成的状态可以凑出0~s，这样对每个s，分成log(s)份

时间复杂度：2*10^7

把每个物品拆成log份，1个，2个，4个…(2的整数幂)

代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#include<algorithm>
using namespace std;

const int N=110;
int n,m; //n个物品，m总容量 
int f[N];
struct good
{
	int v,w; //体积，价值 
 };
 int main(){
 	vector<good>goods; //表示物品 
 	cin>>n>>m;
 	for(int i=0;i<n;i++){  //枚举所有物品 
	int v,w,s; //体积，价值，个数 
	cin>>v>>w>>s;
	for(int k=1;k<=s;k*=2) //2的整数幂 
	{
		s-=k; //s一直减，直到为负数
		goods.push_back({v*s,w*s}); //把拆完后的物品放到物品组 
	}
	if(s>0) goods.push_back({v*s,w*s}); //s有剩余，剩余的物品加入物品组 
}
   //01背包，从大到小，枚举所有 
	for(auto good: goods) //遍历vector,good
	   for(int j=m;j>=good.v;j--)
	      f[j]=max(f[j],f[j-good.v]+good.w);
		  cout<<f[m]<<endl; 
		  return 0;
 }

---

多重背包问题3—究极版

题目：
问题：
有N件物品和一个容量是V的背包。
第i种物品最多有si件，每件的体积是vi，价值是wi。
求解将哪些物品装入背包，可使这些物品的总体积不超过背包的容量，且价值总和最大。

输入格式
第一行有两个整数，N,V用空格隔开，分别表示物品种数和背包容积。
接下来有N行，每行三个整数vi，wi，si，用空格隔开，分别表示第i种物品的体积、价值和数量。

输出格式
输出一个整数，表示最大价值

数据范围
0<N<=1000
0<V<=20000
0<vi,wi，si<=20000

输入样例
4 5
1 2 3
2 4 1
3 4 3
4 5 2

输出样例
10

分析思路：

时间复杂度:
按上一道题时间复杂度算
1000*log（20000）20000=310^8
超出时限上限

单调队列问题
考虑算f[j]，枚举一共有多少个第i个物品，以此类推，看k个数，在所有和j归在一类余数相同的体积中
也就是说算j的时候，看前k-1个数中挑最大值
f[j]=f[j-v]+w,f[j-2v]+2w,…f[j-kv]+kw

算f[j+v]每个数都会变化,每个数都会+w
f[j+v]-f[j]+w,f[j+v]+2w
因为求最大值，所有数都+w，并不影响其他数之间的关系，只需求哪个数最大
可把f[j]变为越小的加的越大
f[j]变为f[j]-kw，余数是j模v的第k个数
第一个数减0w，f[0]=0
f[v]变为：f[v]-1w
f[2v]-2*w
以此类推，在比较的时候，越小（前）的数，加上之后会变得幅度比后面的数都大

算f[j]+v时，把k个数往后移1位，但总容量中仍然有k个数，变为经典单调队列问题

给定一个序列，动态求所有长度为k的背包里面的最大值

举例:
假设一共一串数，先求前3个数的最大值，再求下一个3个数的最大值。

多重背包单调队列优化，需要理解两点

1.容量为什么分类（分完类，同一类里的容量最大价值就可以用单调队列了）

2.怎么用单调队列处理。

#include <iostream>
using namespace std;

const int N=20010;
int n,m;
int f[N],g[N],q[N];
int main()
{
cin>>n>>m;
for(int i=0;i<n;i++){   //物品种类枚举O(n)
int v,w,s;
cin>>v>>w>>s;
for(int j=0;j<v;j++){  //容量类别枚举O(v)
int le=0,ri=-1;
for(int k=j;k<=m;k+=v){// O(每个类别里容量的个数)，
//为什么正序，因为单调队列里保存了q[N]容量最大价值的大小关系
//以下为单调队列的典型写法，q保存区间大小关系，g保存下标（这里是容量）用来判断s个物品最大的价值
while(le<=ri&&q[ri]<f[k]-k/v*w) ri--;
ri++;
q[ri]=f[k]-k/v*w;
id[ri]=k;
if(g[le]+s*v<k) le++;
f[k]=q[le]+g[ri]/v*w;
}
}
}
cout<<f[m]<<endl;
return 0;
}

代码：

#include<iostream>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;

const int N=20010;
int n,m; 
int f[N],g[N],q[N];
int main(){
	cin>>n>>m;
	for(int i=1;i<n;i++){  //枚举所有物品 
	int v,w,s; 
	cin>>v>>w>>s;
	memcpy(g,f,sizeof f);
	for(int j=0;j<v;j++){  //从0~c枚举余数 
	//余数可单独考虑 
		int hh=0,tt=-1;
		for(int k=j;k<=m;k+=v)  //枚举余数j里面的所有数
		{  //单调队列优化
			f[k]=g[k];
			if(hh<=tt&&k-s*v>q[hh])  //每次把队首取出来，它是最大的
			hh++;
			if(hh<=tt)
			f[k]=max(f[k],g[q[hh]]+(k-q[hh])/v*w); //用最大数更新当前数
			while(hh<=tt && g[q[tt]]-(q[tt]-j)/v*w<=g[k]-(k-j)/v*w)  //把当前数插入队列，把队列中不用元素剔除
			 tt--;
			 q[++tt]=k;把当前数加入队列
		}
	}
}
    cout<<f[m]<<endl;
	return 0;
}