回溯法实验报告
实验内容
-
利用回溯法编程求解0-1背包问题,并对算法进行时间复杂性分析;
-
利用回溯法编程求解TSP问题,并对算法进行时间复杂性分析。
【注】输出结果要求显示每一步的搜索过程。
实验目的
- 掌握回溯法的设计思想;
- 掌握解空间树的构造方法,以及在求解过程中如何存储求解路径;
- 考察回溯法求解问题的有效程度。
程序清单
/*
* @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
*/
运行结果
分析与思考
遇到的主要问题就是解空间树的构造方法,以及在求解过程中如何存储求解路径;
在递归过程中错误的改变了一些参数,导致结果错误,经过调试完善了代码。