动态规划：背包问题（用C/C++实现）

一、0-1背包

题意：N个物品，容量是V的背包，每个物品有体积vi和价值wi，每个物品最多用1次。问题：在背包能装得下的情况下，我们能选择物品的最大总价值。

解法：状态表示： f(i, j)  集合（表示所有选法的集合，满足两个条件：只从前i个物品中选、选出来物品总体积 <= j）和属性（选出价值最大值）

状态计算： 集合的划分 （需要求的是f(N, V)）将当前的集合f(i, j)划分为两个子集，分别为含第i个物品的集合f(i-1, j)与不含第i个物品的集合f(i-1, j-vi) + wi，两集合取最大即为结果。划分原则：两个集合不能出现重复元素、不能有元素遗漏。

优化：从二维数组改为一维数组。f(i)只与f(i-1)有关，j与j – vi都只在j的一侧。

#include<iostream>
#include<algorithm> 
/*N个物品，容量是V的背包，每个物品有体积vi和价值wi，
每个物品最多用1次。
问题：在背包能装得下的情况下，我们能选择物品的最大总价值。*/ 
using namespace std;

const int N = 1010;
/*二维 
//n是物品个数，m是背包容量
int n, m;
int v[N], w[N];
int f[N][N];

int main()
{
	cin >> n >> m;
	
	for(int i = 1; i <= n; i ++) cin >> v[i] >> w[i];
	
	//f[0][0]全局变量定义时已经初始化
	for(int i = 1; i <= n; i ++)
		for(int j = 0; j <= m; j ++)
		{
			//含i的集合
			f[i][j] = f[i - 1][j];
			//不含i的集合，要保证背包能够放下i
			if(j >= v[i]) f[i][j] = max(f[i][j], f[i - 1][j - v[i]] + w[i]); 
			 
		}
		
	cout << f[n][m] << endl;
	return 0;
}
*/
/*一维*/ 
int n, m;
int v[N], w[N];
int f[N];

int main()
{
	cin >> n >> m;
	
	for(int i = 1; i <= n; i ++) cin >> v[i] >> w[i];
	
	//f[0][0]全局变量定义时已经初始化
	for(int i = 1; i <= n; i ++)
		for(int j = m; j >= v[i]; j --)//j从0~v[i -1]的部分是没有意义的 
		{
			
			f[j] = max(f[j], f[j - v[i]] + w[i]); 
			 
		}
		
		cout << f[m] << endl;
		return 0;
}

二、完全背包

题意：每件物品有无限个。

解法：状态表示：f(i, j)所有只考虑前i个物品且总体积不大于j的所有选法的集合。属性为集合中总价值最大的选法。

状态计算： 集合的划分，分成若干组，分别是第i个物品选0个f(i – 1, j)，第i个物品选1个，第i个物品选2个……第i个物品选k个(f(i - 1, j – k*vi) + k*wi)……。

优化：优化为二维。

#include<iostream>
#include<algorithm> 
/*每件物品有无限个*/
using namespace std;

const int N = 1010;

int n, m;
int v[N], w[N];
int f[N][N];

int main()
{
	cin >> n >> m;
	
	for(int i = 1; i <= n; i ++) cin >> v[i] >> w[i];
	
	for(int i = 1; i <= n; i ++)
		for(int j = 0; j <= m; j ++)
		/*朴素做法 
			for(int k = 0; k * v[i] <= j; k ++)
			{
				f[i][j] = max(f[i][j], f[i - 1][j - k * v[i]] + k * w[i]);
			}
			*/
		{
			f[i][j] = f[i - 1][j];
			if(j >= v[i]) f[i][j] = max(f[i][j], f[i][j - v[i]] + w[i]);
		}	
			
	cout << f[n][m] << endl;
			
	return 0;
}
/*一维 
//n是物品个数，m是背包容量
int n, m;
int v[N], w[N];
int f[N];

int main()
{
	cin >> n >> m;
	
	for(int i = 1; i <= n; i ++) cin >> v[i] >> w[i];
	
	//f[0][0]全局变量定义时已经初始化
	for(int i = 1; i <= n; i ++)
		for(int j = v[i]; j <= m; j ++)
			f[j] = max(f[j], f[j - v[i]] + w[i]); 
		
	cout << f[m] << endl;
	return 0;
}
*/

三、多重背包

题意：每个物品最多有si个。

解法：状态表示：f(i, j) 所有只考虑前i个物品且总体积不大于j的所有选法的集合。属性集合中总价值最大的选法。

状态计算：集合划分，分别是第i个物品选0个f(i – 1, j)，第i个物品选1个，第i个物品选2个……第i个物品选k个(f(i - 1, j – k*vi) + k*wi)……第i个物品选si个。

f(i, j) = max(f(i – 1, j – k * vi) + wi * k) (k = 0, 1, 2, ……，si)

#include<iostream>
#include<algorithm> 
/*每件物品有si个*/
using namespace std;

const int N = 1010;

int n, m;
int v[N], w[N], s[N];
int f[N][N];

int main()
{
	cin >> n >> m;
	
	for(int i = 1; i <= n; i ++) cin >> v[i] >> w[i] >> s[i];
	
	for(int i = 1; i <= n; i ++)
		for(int j = 0; j <= m; j ++)
			for(int k = 0; k  <= s[i] && k * v[i] <= j; k ++)
			{
				f[i][j] = max(f[i][j], f[i - 1][j - k * v[i]] + k * w[i]);
			}
			
	
			
	cout << f[n][m] << endl;
			
	return 0;
}

优化：

由于多出一项，故不能使用完全背包优化方法来优化该问题。

对于s，列出数列1 2 4 8 16 32 …… 2k c（c < 2k+1），使用这些数一定能凑出s，对于1~2k能够凑出0~2k+1-1中的任何数，对于1~c能够凑出0~2k+1-1 + c中的任何数，2k+1-1 + c即s。故我们能将背包中的若干物品分成多组，有效降低时间复杂度。

s被优化成了logs，时间复杂度NVs被优化成了NVlogs。

#include<iostream>
#include<algorithm> 
/*每件物品有si个*/
using namespace std;

//si -> logsi O(NVs)-> O(NVlogs)
const int N = 25000, M = 2000;//1000*log2000

int n, m;
int v[N], w[N];
int f[N];

int main()
{
	cin >> n >> m;
	
	int cnt = 0;//用于存储所有物品 
	for(int i = 1; i <= n; i ++)
	{
		//当前物品的体积，价值，数量 
		int a, b, s; 
		cin >> a >> b >> s;
		//开始列数，从1开始
		int k = 1;
		while(k <= s)
		{
			cnt ++;
			v[cnt] = a * k;
			w[cnt] = b * k;
			s -= k;
			k *= 2;
		}
		if(s > 0) 
		{
			cnt ++;
			v[cnt] = a * s;
			w[cnt] = b * s;
		 } 
	} 
	
	n = cnt;
	
	for(int i = 1; i <= n; i ++)
		for(int j = m; j >= v[i]; j --)
				f[j] = max(f[j], f[j - v[i]] + w[i]);
		
			
	
			
	cout << f[m] << endl;
			
	return 0;
}

四、分组背包

题意：物品有N组，每组中有若干物品，每组中只能选一个。

解法：状态表示：f(i, j) 所有只考虑前i个物品且总体积不大于j的所有选法的集合。属性集合中总价值最大的选法。

状态计算：集合划分，分为第i组物品选哪一个。对于第i组物品，一个都不选f(i – 1, j)，选第1个，选第2个，……选第k个f(i-1, j – v(i, k)) + w(i, k)，……。

#include<iostream>
#include<algorithm> 
/*物品有N组，每组中有若干物品，每组中只能选一个*/
using namespace std;

const int N = 110;

int n, m;
int v[N][N], w[N][N], s[N];
int f[N];

int main()
{
	cin >> n >> m;
	
	for(int i = 1; i <= n; i ++)
	{
		cin >> s[i];
	
		for(int j = 0; j < s[i]; j ++)
				cin >> v[i][j] >> w[i][j];
	
	} 
		
	for(int i = 1; i <= n; i ++)
		for(int j = m; j >= 0; j --)
			for(int k = 0; k < s[i]; k ++)
				if(v[i][k] <= j)
					f[j] = max(f[j], f[j - v[i][k]] + w[i][k]);
	
	cout << f[m] << endl;
	
	return 0;
}