TSP问题的Lingo求解

Lingo编程解决TSP问题

TSP问题即巡回旅行商问题，一个商人旅行经过所有城市一次最后回到原点，问什么走法使走的路程最短。接下来用Lingo编程实现，此方法具有普遍性，建议小本本记下。
假设有6个城市

思路

利用01矩阵表示一次巡回旅行的方案。如下面矩阵代表一种方案。
$$\left(\begin{array}{c}010000\\ 001000\\ 000100\\ 000010\\ 000001\\ 100000\end{array}\right)$$
a12=1表示从城市1走到城市2。那么这个矩阵的路线为1–>2–>3–>4–>5–>6–>1。
观察这个矩阵的特点为每一行只有一个1，每一列只有一个1。但并不是这样的矩阵都是合理的，例如如下矩阵：
$$\left(\begin{array}{c}010000\\ 001000\\ 100000\\ 000010\\ 000001\\ 000100\end{array}\right)$$
这个矩阵的路线为：1–>2–>3–>1，4–>5–>6–>4。说明路线矩阵不能出现子圈。
总距离=路线矩阵*距离矩阵。

建立线性规划模型

设城市之间的距离用矩阵d表示，dij表示城市i与城市j的距离。
设0–1矩阵X表示经过各城市之间的路线。xij=0表示城市i不到城市j，为1则反之。
路线矩阵每一行和为1：
$$\sum _{j=1}^n{x}_{ij}=1,i=1,...,n(i\ne j)$$
每一列和为1：
$$\sum _{i=1}^n{x}_{ij}=1,j=1,...,n(j\ne i)$$
路线矩阵不可以出现子圈（破圈），引入额外变量u：
$$u_i-u_j+n{x}_{ij}\le n-1,1<i\ne j\le n$$
对于破圈公式的解释：

1. 两个点出现子圈的情况：$x_{ij}=1,x_{ji}=1。$
   把$x_{ij},x_{ji}$代入破圈公式得：$$u_i-u_j\le -1,u_j-u_i\le -1$$
   两式相加得：$0\le -2$，显然不成立，保证了两个点不可能出现子圈。

2. 对于三个点的情况：$$x_{ij}=1,x_{jk}=1,x_{ki}=1$$ 代入破圈公式得：
   $$u_i-u_j\le -1,u_j-u_k\le -1,u_k-u_i\le -1$$
   相加得：$0\le -3$，不成立，保证了三个点不出现子圈。

3. 对于n个点以此类推。

此外:
$$x_{ij}=0或1，i,j=1..n$$$$u_i为实数，i=1..n$$
目标函数为：
$$minz=\sum _{i=1}^n{d}_{ij}\sum _{j=1}^n{x}_{ij}$$

Lingo编程代码如下

!TSP问题;
MODEL:
SETS:
city/1..6/:u;
link(city,city):d,x;
ENDSETS
DATA:
d=@FILE(E:\d.txt);!导入距离矩阵;
@text()=@writefor(link(i,j)|x(i,j)#GT#0:'x(',i,',',j,')=',x(i,j));!指定输出格式，输出为1的x;
ENDDATA
min=@sum(link:d*x);
@for(city(j):@sum(city(i)|j#ne#i:x(i,j))=1);!一列只有一个1;
@for(city(i):@sum(city(j)|j#ne#i:x(i,j))=1);!一行只有一个1;
@for(link(i,j)|i#ne#j#and#i#gt#1:u(i)-u(j)+6*x(i,j)<=5);!破圈;
@for(link:@bin(x));!x取0或1;
End

结果

解得最短路程是35，其路线为1–>4–>3–>2–>5–>6–>1。