分支限界法求解旅行商售货员（TSP）问题 C语言

问题描述

某售货员要到若干城市去推销商品，已知各城市之间的路程（或旅费）。他要选定一条从驻地出发，经过每个城市一次，最后回到驻地的路线，使总的路程（或总旅费）最小。

分支限界法解题步骤

（1）先不考虑原问题的整数约束，求解相应的松弛问题。用图解法或单纯形法求得最优解，记为  。

（2）若求得的最优解  刚好就是整数解，则该整数解就是原整数规划问题的最优解；否则，对原问题进行分枝寻求整数最优解。

（3）分枝。根据对变量重要性的了解，在最优解中选择一个不符合整数约束条件的xj ，其值为bj ，以[bj]表示小于bj 的最大整数。构造两个约束条件: x≤ [bj]和 x≥[bj]+1分别加入原LP问题形成两个子问题，因为[bj] 与[bj]+1之间无整数，故这两个子集内的整数解必定与原可行解集合整数解一致，这一步称为分枝。

（4）定界。首先判断各个子问题是否存在整数解。若存在，找出目标函数值最大对应的整数解，设为Z*，则A问题的整数解目标函数Z≥Z*，这就是定界。而且分枝过程中，一旦有某个子问题Z≥Z*，则令Z*＝Z。

（5）若存在大于Z*的子问题则需分枝。第（4）步中若不存在整数解，也需继续分枝寻找整数解，并从目标函数值最大对应的子问题先分枝。

（6）若所有子问题的目标值都小于等于Z*，则不需继续分枝，Z*所对应的整数解即为最优解。

代码

#include<iostream>
using namespace std;
#define NoEdge -1
#define MAX 20
int G[MAX] [MAX];
int ans[MAX],x[MAX];
int bestc, cc;
 
void init(int n)
{
int i,j, len;
 
memset (G, NoEdge, sizeof (G)) ;
while (cin>>i>>j)
{
if (i==0&&j==0) break;cin>> len;
 
G[i] [j]=len;
 
G[j] [i]=len;
}
for(i=1;i<=n;i++) x[i]=i;bestc=0x7fffff;
 
cc=0;
}
 
void Swap(int &i,int &j)
{
  int t=i;
  i=j;
  j=t;
}
void Traveling(int i,int n)
{
int j;
 
if(i==n+1)
{
if(G[x[n-1]][x[n]]!=NoEdge && G[x[n]][1]!= NoEdge &&(cc+ G[x[n]][1]<bestc))
{
	for(j=1;j<=n;j++)ans[j]=x[j];
 
bestc=cc+= G[x[n]][1];
}}
else{
 
 
 
	for(j=i;j<=n;j++){
 
if(G[x[i-1]] [x[j]]!=NoEdge && (cc+ G[x[i-1]] [x[j]]<bestc))
{
       Swap (x[i],x[j]);
 
	   cc+= G[x[i-1]] [x[i]];Traveling ( i+1,n) ;cc-= G[x[i-1]][x[i]];Swap(x[i],x[j]);}
	}
}
}
 
void print(int n)
{
cout<<"最小的旅行费用为: "<<bestc<<endl;
cout<<"最佳路径是: ";
for(int i=1;i<=n;i++)
 cout<<ans[i] <<"->";
cout<<ans [1]<<endl;}
 
int main()
{
int n;
 
cout<<"请输人需要旅行多少个城市: "<<endl;
while (cin>>n&&n){
 
cout<<"輸人丙个城市之同的距高，例如1 2 20,輸人00結束"<<endl;
init(n) ;
 
Traveling(2,n) ;
 
print(n) ;}
return 0;
}

结果

     ’