【回溯法】-排列树-旅行商(TSP)问题C++实现

旅行商问题描述

有n个城市，任意两个城市之间的距离已知。一个旅行商从某个城市出发经过每个城市且任一城市只经过一次，最后回到出发城市，如何确定代价最小的线路

商人从西安出发，经过图中所示的所有城市回到西安，可以选择先到成都，然后到达广州，继续访问HK，访问完香港后不能回到广州，因为每个城市智能访问一次，离开香港后可以去上海，然后继续北上来到徐州，最后经过北京回到西安。上述路线是该问题的一个可行解，当然这个解不一定是最优解。
现在明确了什么是商旅问题。思考一下回溯法求解问题的基本步骤，

回溯法步骤

1、针对所给的问题，定义问题的解空间。
2、确定易于搜索的解空间结构（排列数，子集树）。
3、从根结点开始，采用深度优先的方式遍历解空间树。

在遍历的过程中为了避免无效搜索，采用剪枝策略提高搜索的速度。
1、利用约束函数减去不满足约束条件的子树。
2、利用限界函数减去得不到最优解的子树。

举例

假设有四个城市，编号分别为1,2,3,4.
该问题不是从4个城市中选择若干城市，而是求解对4个城市的访问顺序
因此需要使用排列数组织解空间。
从根到叶子的路径对应问题的一个可能解。采用深度优先的方法搜索这个树。
从跟结点A开始，到达结点B，继续深度优先，到达结点C，此时的含义是从城市1出发到城市2，继续深度优先到达结点F，此时的含义是，访问完城市2之后访问城市3，继续深度优先到达结点L，此时的含义是访问完城市3之后访问城市4，结点L是一个叶子结点，因此得到第一个解，即1,2,3,4。该可行解的距离代价是59。
到达叶子结点L后，从L开始回溯到结点F，此时结点C的左子树访问完毕，开始访问结点C的右子树，到达结点G，此时的含义是访问城市2之后，访问城市4，当前的距离是40,如果继续深度优先，到达点M，访问的代价是60，说明还没有回到原点，访问代价就已经超过了目前得到的最优解。因此使用限界函数将结点M剪枝，

源代码

/**
   回溯法-旅行商(TSP)问题
*/
#include<iostream>
#include<algorithm>
#define MAX 100
using  namespace std;
int n;                               //城市个数
int a[MAX][MAX];                   //城市间距离
int x[MAX];                       //记录路径
int bestx[MAX]  = {0};           //记录最优路径
int bestp = 63355;              //最短路径长
int cp = 0;                    //当前路径长
void backpack(int t){
     if(t>n){
        if((a[x[n]][1])&&(a[x[n]][1]+cp<bestp)){
              bestp = a[x[n]][1]+cp;
              for(int i = 1;i<=n;i++){
                 bestx[i] = x[i];
              }
        }
     }else{
         for(int i = t;i<=n;i++){
             /*约束为当前节点到下一节点的长度不为0,限界为走过的长度+当前要走的长度之和小于最优长度*/
            if((a[x[t-1]][x[i]])&&(cp+a[x[t-1]][x[i]]<bestp)){
                swap(x[t],x[i]);   
                cp+=a[x[t-1]][x[t]];
                backpack(t+1);
                cp-=a[x[t-1]][x[t]];
                swap(x[t],x[i]);
            }
         }
    }
}
int main(){
    cout<<"输入城市个数:"<<endl;
    cin>>n;      //顶点数
    for(int i = 1;i<=n;i++){
         x[i] = i;
    }
    cout<<"输入城市之间的距离(0表示城市间不通):"<<endl;
    for(int i = 1;i<=n;i++){
        for(int j = 1;j<=n;j++){
            cin>>a[i][j];
        }
    }
    backpack(2);
    cout<<"最少旅行费用为: "<<bestp<<endl;
    cout<<"旅行路径为:"<<endl;
    for(int i = 1;i<=n;i++){
       cout<<bestx[i]<<" ";
    }
    cout<<bestx[1];
    return 0;
}
/*
输出结果
4
0 30 6 4
30 0 5 10
6 5 0 20
4 10 20 0
*/