六、 参考文献
薛毅，优化模型与lindo/lingo软件，清华大学出版社，2005.7。 ?1?谢金星，
?2?火车网，http://www.huoche.com.cn/，2014年5月25日。 ?3?携程网，http://www.ctrip.com/，2014年5月25日。
11
七、 附录
附表1： 南昌 九江 景德镇 鹰潭 上饶 抚州 南昌 0 20 40 24 40 30 九江 20 0 45 40 55 40 景德镇 40 45 0 50 70 67 鹰潭 24 40 20 0 17 50 上饶 40 55 70 17 0 70 抚州 30 40 67 50 70 0 Lingo程序
运用LINGO，输入程序：
MODEL:
!Traveling Sales Problem for the cities of six city; SETS:
CITY / 1..6/: U; ! U( I) = sequence no. of city; LINK( CITY, CITY): COST, ! The cost matrix; X; ! X( I, J) = 1 if we use link I, J; ENDSETS
DATA: !Cost matrix, it need not be symmetric; COST= 0
20 40 24 40 30 20 0 45 40 55 40 40 45 0 20 70 67 24 40 50 0 17 50 40 55 70 17 0 70 30 40 67 50 70 0;
ENDDATA
N= @SIZE( CITY);
MIN= @SUM( LINK: COST * X);
@FOR( CITY( K): !It must be entered; @SUM( CITY( I)| I #NE# K: X( I, K)) = 1; ! It must be departed;
@SUM( CITY( J)| J #NE# K: X( K, J)) = 1; !Weak form of the subtour breaking constrains; !These are not very powerful for large problems; @FOR( CITY( J)| J #GT# 1 #AND# J #NE# K:
U(J)>=U(K)+X(K,J)-(N-2)*(1-X(K,J))+(N-3)*X(J,K));); @FOR( LINK: @BIN( X)); !Make the X's 0/1; ! For the first and last stop we know...; @FOR( CITY( K)| K #GT# 1:
12
U( K) <= N - 1 - (N - 2) * X( 1, K); U( K) >= 1 + ( N - 2) * X( K, 1)); END
13