算法设计与分析实验——回溯法实验报告

回溯法实验报告

实验内容

1. 利用回溯法编程求解0-1背包问题，并对算法进行时间复杂性分析；

2. 利用回溯法编程求解TSP问题，并对算法进行时间复杂性分析。

   【注】输出结果要求显示每一步的搜索过程。

实验目的

1. 掌握回溯法的设计思想；
2. 掌握解空间树的构造方法，以及在求解过程中如何存储求解路径；
3. 考察回溯法求解问题的有效程度。

程序清单

/*
 * @Description: 01背包
 * @version:
 * @Author:
 * @Date: 2021-05-26 17:00:46
 * @LastEditors: Please set LastEditors
 * @LastEditTime: 2021-05-26 19:00:34
 */
#include <iostream>
using namespace std;
#define N 100      // 假设最大有99个物品
int w[N];          // 每个物品的重量
int v[N];          // 每个物品的价值
int x[N];          // x[i]=1：物品i放入背包，0代表不放入
int n, c;          // n：一共有多少物品，c：背包的最大容量
int CurWeight = 0; //当前放入背包的物品总重量
int CurValue = 0;  //当前放入背包的物品总价值
int BestValue = 0; //最优值；当前的最大价值，初始化为0
int BestX[N];      //最优解；BestX[i]=1代表物品i放入背包，0代表不放入
void input()
{
	cout << "请输入物品的个数：" << endl;
	cin >> n;
	cout << "请输入每个物品的重量及价值:" << endl;
	for (int i = 1; i <= n; i++)
		cin >> w[i] >> v[i];
	cout << "请输入背包的容量：" << endl;
	cin >> c;
}
void output()
{
	cout << "--------------------" << endl;
	cout << "当前的最优值是:" << BestValue << endl;
	cout << "(";
	for (int i = 1; i <= n; i++)
		cout << BestX[i] << " ";
	cout << ")" << endl;
}
// t=1时表示当前在决定是否将第一个物品放入背包
void backtrack(int t)
{
	// 完成了所有物品的决策
	if (t > n)
	{
		//如果找到了一个更优的解，记录、输出，否则退回上一层
		if (CurValue > BestValue)
		{
			//保存更优的值和解
			BestValue = CurValue;
			for (int i = 1; i <= n; ++i)
				BestX[i] = x[i];
			output();
		}
	}
	else
	{
		//遍历当前节点的子节点：0 不放入背包，1放入背包
		for (int i = 0; i <= 1; ++i)
		{
			// 记录当前决策，放入与否
			x[t] = i;
			// 如果这个物品不放入，进行下一物品的决策
			if (i == 0)
				backtrack(t + 1);
			// 如果这个物品放入的话
			else
			{
				// 判断是否能放下
				if ((CurWeight + w[t]) <= c)
				{
					// 放入第t个物品，记录当前的装入物品的容量与价值
					CurWeight += w[t];
					CurValue += v[t];
					backtrack(t + 1);
					// 进行回溯
					CurWeight -= w[t];
					CurValue -= v[t];
				}
			}
		}
	}
}
int main(void)
{
	input();
	backtrack(1);
	system("pause");
	return 0;
}
/*
6
1 1
2 2
3 3
4 4
5 5
6 6
15
*/

#include <iostream>
#include <cstring>
using namespace std;
#define N 100
int n;
int cost[N][N];      // 表示城市之间的旅行耗费，不存在路径则耗费等于-1
int bestTour[N];     // 记录最佳路径
int currentTour[N];  // 记录当前路径
int bestCost = -1;   // 当前最低耗费
int currentCost = 0; // 当前耗费
void output()
{
	cout << "--------------------" << endl;
	cout << "当前的最优值是:" << bestCost << endl;
	cout << "(";
	for (int i = 1; i <= n; i++)
		cout << bestTour[i] << " ";
	cout << ")" << endl;
}
void rTSP(int currentLevel)
{
	// 已经到达最后一个点
	if (currentLevel == n)
	{
		// 若前一点到这点可达且该点可以回到起点，且这个可行解的花费比当前的低，就保存
		if (cost[currentTour[n - 1]][n] != -1 && cost[currentTour[n]][1] != -1 && (bestCost == -1 || currentCost + cost[currentTour[n - 1]][n] + cost[currentTour[n]][1] < bestCost))
		{
			bestCost = currentCost + cost[currentTour[n - 1]][n] + cost[currentTour[n]][1];
			copy(currentTour + 1, currentTour + n + 1, bestTour + 1);
			output();
		}
	}
	else
	{
		// 从当前层开始循环判断该层的所有可能
		for (int j = currentLevel; j <= n; ++j)
			// 如果上一层到该点可达，且当前花费+这个边的花费小于当前保存的最小值，则进行处理，否则不处理，进行返回上层 ，即剪枝
			if (cost[currentTour[currentLevel - 1]][currentTour[j]] != -1 && (bestCost == -1 || currentCost + cost[currentTour[currentLevel - 1]][currentTour[j]] < bestCost))
			{
				// 交换当前路线的该层值与满足条件的结点值
				swap(currentTour[currentLevel], currentTour[j]);
				// 保存当前的花费
				currentCost += cost[currentTour[currentLevel - 1]][currentTour[currentLevel]];
				// 进行下一层的决策
				rTSP(currentLevel + 1);
				// 退出下一层的决策后，进行回溯，消除这层决策的影响，在这层寻找之后的可行结点，以进入兄弟子树
				currentCost -= cost[currentTour[currentLevel - 1]][currentTour[currentLevel]];
				swap(currentTour[currentLevel], currentTour[j]);
			}
	}
}
int main(void)
{
	memset(cost, -1, sizeof(cost));
	int i;
	cout << "请输入结点数：" << endl;
	cin >> n;
	for (i = 1; i <= n; i++)
		currentTour[i] = i;
	int numEdge;
	cout << "请输入边数：" << endl;
	cin >> numEdge;
	cout << "请输入边的顶点与边的权值：" << endl;
	for (int i = 0; i < numEdge; i++)
	{
		int a, b, p;
		cin >> a >> b >> p;
		cost[a][b] = p;
	}
	rTSP(2); //起点定为1，从第二层开始
	cout << "--------------------" << endl;
	cout << "最少的运费为:" << bestCost << endl;
	cout << "最佳路径为: ";
	for (i = 1; i <= n; i++)
		cout << bestTour[i] << "->";
	cout << bestTour[1] << endl;
	system("pause");
	return 0;
}
/*
4
12
1 2 100
1 3 6
1 4 100
2 1 5
2 3 2
2 4 3
3 1 6
3 2 1
3 4 2
4 1 3
4 2 7
4 3 5
 */

运行结果

分析与思考

遇到的主要问题就是解空间树的构造方法，以及在求解过程中如何存储求解路径；

在递归过程中错误的改变了一些参数，导致结果错误，经过调试完善了代码。