算法：动态规划3

提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档

---

一、

1.0-1背包问题

给定n种物品和一背包。物品n的重量是wi，其价值为vi，背包的容量为C。问:应该如何选择装入背包的物品，使得装入背包中物品的总价值最大?

在选择装人背包的物品时,对每种物品i只有两种选择，即装入背包或不装人背包。不能将物品i装入背包多次，也不能只装入部分的物品i。因此,该问题称为0-1背包问题。

公式

价值最大：max(sum(xi*vi))
重量小于最大值：sum(xi*wi)<=c;
是否装入背包：m(i, j) = max(m(i + 1, j), m(i + 1, j - wi) + vi);

2. 代码示例

递归关系式 & 非递归关系式

//递归表示
int Knapsack(vector<int>& w, vector<int>& v, int i, int j, int n)
{
	if (i == n)//只有一个物品
	{
		return j >= w[i] ? v[i] : 0;
	}
	else
	{
		if (j < w[i])//不装，因为物品重量大于背包容量
		{
			return Knapsack(w, v, i + 1, j, n);
		}
		else
		{
			int maxv1= Knapsack(w, v, i + 1, j, n);
			int maxv2 = Knapsack(w, v, i + 1, j - w[i], n) + v[i];
			return maxv1 > maxv2 ? maxv1 : maxv2;
		}
	}
}
//非递归表示(逆序填充)
int Knapsack2(vector<int>& w, vector<int>& v, int n, int c, vector<vector<int>>& m)
{
	if (n == 0) return 0;
	for (int j = 0; j <= c; ++j)
	{
		m[n][j] = j >= w[n] ? v[n] : 0;
	}
	for (int i = n - 1; i >= 1; --i)
	{
		for (int j = 1; j <= c; ++j)
		{
			if (j >= w[i])
			{
				m[i][j] = max(m[i - 1][j], m[i - 1][j - w[i]] + v[i]);

			}
			else
			{
				m[i][j] = m[i - 1][j];
			}
		}
	}
	return m[1][c];
}
//打印二维数组
void print_vect(vector < vector<int>>& m)
{
	for (int i = 0; i < m.size(); ++i)
	{
		for (int j = 0; j < m[i].size(); ++j)
		{
			printf("%4d", m[i][j]);
		}
		printf("\n");
	}
	printf("\n");
}
//非递归表示(正序填充)
int Knapsack3(vector<int>& w, vector<int>& v, int n, int c, vector<vector<int>>& m)
{
	for (int i = 1; i <= n; ++i)
	{
		for (int j = 1; j <= c; ++j)
		{
			if (j < w[i])//物品不放入
			{
				m[i][j] = m[i - 1][j];
			}
			else//物品放入
			{
				//最大价值
				m[i][j] = std::max(m[i - 1][j], m[i - 1][j - w[i]] + v[i]);
			}
		}
		print_vect(m);
	}
	return m[n][c];
}

int main()
{
	const int n = 5;//物品个数
	const int c = 10;//最大重量
	vector<int> w = { 0,2,2,6,5,4 };//重量
	vector<int> v = { 0,6,3,5,4,6 };//价值

	vector<vector<int>> m(n + 1, vector<int>(c + 1, 0));//初始化
	/*vector<vector<int>> m(n + 1);
	for (int i = 0; i < n + 1; ++i)
	{
		m[i].resize(c + 1);
	}*/
	int maxv = Knapsack3(w, v, n, c, m);
	print_vect(m);
	cout << maxv << endl;
}

3. 打印序列

//递归表示
int Knapsack(vector<int>& w, vector<int>& v, int i, int j, int n)
{
	if (i == n)//只有一个物品
	{
		return j >= w[i] ? v[i] : 0;
	}
	else
	{
		if (j < w[i])//不装，因为物品重量大于背包容量
		{
			return Knapsack(w, v, i + 1, j, n);
		}
		else
		{
			int maxv1= Knapsack(w, v, i + 1, j, n);
			int maxv2 = Knapsack(w, v, i + 1, j - w[i], n) + v[i];
			return maxv1 > maxv2 ? maxv1 : maxv2;
		}
	}
}
//非递归表示(逆序填充)
int Knapsack2(vector<int>& w, vector<int>& v, int n, int c, vector<vector<int>>& m)
{
	if (n == 0) return 0;
	for (int j = 0; j <= c; ++j)
	{
		m[n][j] = j >= w[n] ? v[n] : 0;
	}
	for (int i = n - 1; i >= 1; --i)
	{
		for (int j = 1; j <= c; ++j)
		{
			if (j >= w[i])
			{
				m[i][j] = max(m[i - 1][j], m[i - 1][j - w[i]] + v[i]);

			}
			else
			{
				m[i][j] = m[i - 1][j];
			}
		}
	}
	return m[1][c];
}
//打印二维数组
void print_vect(vector < vector<int>>& m)
{
	for (int i = 0; i < m.size(); ++i)
	{
		for (int j = 0; j < m[i].size(); ++j)
		{
			printf("%4d", m[i][j]);
		}
		printf("\n");
	}
	printf("\n");
}
//非递归表示(正序填充)
int Knapsack3(vector<int>& w, vector<int>& v, int n, int c, vector<vector<int>>& m)
{
	for (int i = 1; i <= n; ++i)
	{
		for (int j = 1; j <= c; ++j)
		{
			if (j < w[i])//物品不放入
			{
				m[i][j] = m[i - 1][j];
			}
			else//物品放入
			{
				//最大价值
				m[i][j] = std::max(m[i - 1][j], m[i - 1][j - w[i]] + v[i]);
			}
		}
		print_vect(m);
	}
	return m[n][c];
}
//打印序列
void backx(vector<int>& w, vector<vector<int>>& m, int n, int c, vector<bool>& X)
{
	for (int i = n; i >= 1; --i)
	{
		if (m[i][c] != m[i - 1][c])
		{
			X[i] = true;
			c = c - w[i];
		}
	}
}
int main()
{
	const int n = 5;//物品个数
	const int c = 10;//最大重量
	vector<int> w = { 0,2,2,6,5,4 };//重量
	vector<int> v = { 0,6,3,5,4,6 };//价值

	vector<vector<int>> m(n + 1, vector<int>(c + 1, 0));//初始化
	vector<bool> X(n + 1, false);
	/*vector<vector<int>> m(n + 1);
	for (int i = 0; i < n + 1; ++i)
	{
		m[i].resize(c + 1);
	}*/
	int maxv = Knapsack3(w, v, n, c, m);
	print_vect(m);
	cout << maxv << endl;
	backx(w, m, n, c, X);

	for (int i = 1; i <= n; ++i)
	{
		if (X[i])
		{
			cout << i << endl;
		}
	}
}

4. 矩阵连乘问题

代码示例

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

int matrixchain(int i, int j, int p[])
{
	if (i == j)
	{
		return 0;
	}
	else
	{
		int k = i;
		int t = matrixchain(i, k, p) + matrixchain(k + 1, j, p) + p[i - 1] * p[k] * p[j];
		for (k = i + 1; k < j; ++k)
		{
			int mt = matrixchain(i, k, p) + matrixchain(k + 1, j, p) + p[i - 1] * p[k] * p[j];
			if (t > mt)
			{
				t = mt;
			}
		}
		return t;
	}
}
void Print_Vec(vector<vector<int>>& m)
{
	int n = m.size();
	for (int i = 0; i < n; ++i)
	{
		for (int j = 0; j < n; ++j)
		{
			printf("%7d", m[i][j]);
		}
		printf("\n");
	}
	printf("\n");
}
int matrixchain_1(int p[], int n, vector<vector<int>>& m, vector<vector<int>>& s)
{
	for (int i = 1; i <= n; ++i)
	{
		m[i][i] = 0;
	}
	for (int num = 2; num <= 6; ++num)
	{
		Print_Vec(m);
		for (int i = 1; i < n - num + 1; ++i)
		{
			int j = i + num - 1;
			m[i][j] = m[i][i] + m[i + 1][j] + p[i - 1] * p[i] * p[j];
			s[i][j] = i;
			for (int k = i + 1; k < j; ++k)
			{
				int t = m[i][k] + m[k + 1][j] + p[i - 1] * p[k] * p[j];
				if (m[i][j] > t)
				{
					m[i][j] = t;
					s[i][j] = k;
				}
			}
		}
	}
	return m[1][n];
}

int main()
{
	const int n = 6;
	int p[n + 1] = { 30,35,15,5,10,20,25 };
	vector < vector<int>>m(n + 1, vector<int>(n + 1, 0));
	vector < vector<int>>s(n + 1, vector<int>(n + 1, 0));


	int smin = matrixchain_1(p, n, m, s);
	Print_Vec(m);
	cout << smin << endl;
	return 0;
}

---