旅行商问题—蛮力法

旅行商问题之蛮力法求解最短路径

找到一个从任一顶点出发，经过所有顶点再回到此顶点的最短路径。

如图，假设从顶点0开始出发，经过顶点1，2，3再回到顶点0，求解路径最短的方法。

 暴力法求解，就是将所有情况列举出来，然后选择权值最小的路径为最优解。可以先对顶点进行全排列，然后用计算回路权值的函数将每种情况的权值之和计算出来，再选择其中最小的作为最优解。

具体c++代码如下：

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

//计算权值
double pathCost(const vector<vector<double>>& G, const vector<int>& path)
{
	double cost = 0;
	int k = path[0];
	for (int i = 1; i < path.size(); i++)
	{
		cost += G[k][path[i]];
		k = path[i];
	}
    //最后要构成回路，所以再加上一个回到起点的权值
	cost += G[k][path[0]];
	return cost;
}

//全排列求最短路径 G是权值数组 S是存放顶点  minPath存放最短路径
double TSP_BF(vector<vector<double>>& G, vector<int>& S, int k, vector<int> minPath)
{
	int n = S.size();
	if (k == n)   //全排列完毕
		return pathCost(G, S);

	double minCost = 100000;
	for (int i = k; i < n; i++)   //进行全排列
	{
		swap(S[i], S[k]);  //回溯法  交换位置之后为了不让原序列打乱，需要再换回来
		double cost = TSP_BF(G, S, k + 1, minPath);
		if (cost < minCost)  //比之前的花费更少，则记录此时的cost和此时的路径
		{
			minCost = cost;
			minPath = S;//此时S已经交换了 所以可以得到最小路径
		}
		swap(S[i], S[k]);  //换回来
	}
	return minCost;
}

//输出vector
void Pirnt_vec(vector<int>& vec)
{
	for (int i = 0; i < vec.size(); i++)
		cout << vec[i] << " ";
	cout << endl;
}

void main()
{
	double inf = 100000;//自己到自己的权值是无穷大
    //G用来存放顶点之间的权值，用了一个二维的集合，定义权值类型为double
	vector<vector<double>> G = { {inf,30,6,4},{30,inf,5,10},{6,5,inf,20},{4,10,20,inf} };
	//S用来存放顶点的集合
    vector<int> S = { 0,1,2,3 };
    //minPath用来存储最优解的路径
	vector<int> minPath;

	//设置0顶点作为初始位置，所以k从1开始
	double minCost = TSP_BF(G, S, 1, minPath);
	cout << minCost << endl;
	Pirnt_vec(minPath);
}