TSP問題是指旅行家要旅行n個城市，要求各個城市經歷且僅經歷一次然後回到出發城市，並要求所走的路程最短。各個城市間的距離可以用代價矩陣來表示。

d(i,V')=min{cik+d(k,V－{k})}(k∈V')

d(k,{})=cki(k≠i)

    虛擬碼：

for (i=1; i<n; i++)     //初始化第0列
	d[i][0]=c[i][0]; 
 for (j=1; j<2n-1-1; j++)  
 	for (i=1; i<n; i++)   //依次進行第i次迭代
		if (子集V[j]中不包含i) 
 			對V[j]中的每個元素k，計算[i][j]=min(c[i][k]+d[k][j-1]);
對V[2n-1-1]中的每一個元素k，計算d[0][2n-1-1]=min(c[0][k]+d[k][2n-1-2]);
輸出最短路徑長度d[0][2n-1-1];
 

    筆者在寫此程式時卡在瞭如何判斷：子集V[j]中不包含i。解決方法是對n個頂點分別用0~n-1的數字編號，按個數為1,2,……,n-1的順序生成1~n-1個元素的子集存放在陣列V[2^n-1]中。這樣，在判斷時就可以將當前點與陣列編號按位與運算（如下圖所示），為0則為需要進行計算（這時才能品味按位與運算的美）。

    完整程式碼：

#include<iostream.h>
#include<math.h>
int main()
{
	int i,j,k,min,brief,n;
	int D[20][20];
	cout<<"頂點個數：";
	cin>>n;
	int b=(int)pow(2,n-1);
  	for(i=0;i<n;i++)
    	for(j=0;j<n;j++)
      		cin>>D[i][j];
  	int ** Route = (int **)calloc(n, sizeof(int *));
  	int ** Mat = (int **)calloc(n, sizeof(int *));
  	for(i=0;i<n;i++)
  	{
    	Route[i] = (int *)calloc(b*b, sizeof(int))+i*b;
    	Mat[i] = (int *)calloc(b*b, sizeof(int))+i*b;
  	}
  	for(i=0;i<b;i++)
   		for(j=0;j<n;j++)
	    {
			Route[j][i] = -1;
   			Mat[j][i] = -1;
	    }
  	for(i=0;i<n;i++)//初始化第0列 
    	Route[i][0] = D[i][0];    
  	for(i=1;i<b-1;i++)
	    for(j=1;j<n;j++)//依次進行第i次迭代 
			if( ((int)pow(2,j-1) & i) == 0)//子集V[j不包含i 
			{
				min=999;
				for(k=1;k<n;k++)
				if( (int)pow(2,k-1) & i )
				{
					brief = D[j][k] + Route[k][i-(int)pow(2,k-1)]; 
					if(brief < min)
					{
						min = brief;
						Route[j][i] = min;
						Mat[j][i] = k;//區域性最優決策
					}
				}    
			}
	min=999;
  	for(k=1;k<n;k++)
	{
		brief = D[0][k] + Route[k][b-1 - (int)pow(2,k-1)];
		if(brief < min)
		{
			min = brief;
			Route[0][b-1] = min;//最優解
			Mat[0][b-1] = k;
		}
	}
	cout<<"最短路徑長度："<<Route[0][b-1]<<endl;//最短路徑長度
	cout<<"最短路徑："<<"1";
	for(i=b-1,j=0; i>0; )
	{
		j = Mat[j][i];
		i = i - (int)pow(2,j-1);
		cout<<"->"<<j+1;
	}
	cout<<"->1"<<endl;
	for(i=0;i<n;i++)
	{
		for(j=0;j<b;j++)
  			cout<<Route[i][j]<<" ";
		cout<<endl;
	}
	free(Route);
	free(Mat);
	return 0;
}