编程：Lingo旅行商问题

编程：Lingo旅行商问题
MODEL:
! 旅行商问题：从城市O出发经过城市A, B, C, D,E,F后又回到O 的过程;
SETS:
CITY / O A B C D E F/: P;
PRED( CITY,CITY)/A,B C,D E,F/旅行优先顺序;
STEP/1..7/STEP（I）表示第I个经过的城市;
LINK( CITY, CITY):
T, ! 时间矩阵;
X; ! X( I, J);
TXS( CITY,STEP):Y;
ENDSETS
DATA: !时间矩阵;
T = 0 1.7 1.5 1.4 3 1.9 2.8
1.7 0 3
2.6 4 3
3.5
1.5 3 0 1.5 3.7
2.6 2.8
1.4
2.6 1.5 0 11.5 3
3.2
3 4 3.7 11.5 0 3.6 3.3
1.9 3
2.6 3
3.6 0 15.6
2.8
3.5 2.8 3.2 3.3 15.6 0;
ENDDATA
!城市个数;
N = @SIZE( CITY);
MIN = @SUM( LINK:T * X);
!第K个经过城市I;
@FOR( CITY( I): @SUM( STEP( K): Y( I, K)) = 1);
!城市A和B经过的次序相邻，同理C和D，E和F;
@FOR( PRED( I, J): @SUM( STEP( K): P(K) * Y( J, K) - P(K) * Y( I, K)) =1#OR#-1);
@FOR( CITY( K):
!进入城区;
@SUM( CITY( I)| I #NE# K: X( I, K)) = 1;
!离开城区;
@SUM( CITY( J)| J #NE# K: X( K, J)) = 1;
@FOR( CITY( J)| J #GT# 1 #AND# J #NE# K:
P( J) >= P( K) + X ( K, J) -
( N - 2) * ( 1 - X( K, J)) +
( N - 3) * X( J, K)
);
);
! 使X、Y取0或1;
@FOR( LINK: @BIN( X));
@FOR( TXS: @BIN( Y));
!出发点和终点;
@FOR( CITY( K)| K #GT# 1:
P( K) <= N - 1 - ( N - 2) * X( 1, K);
P( K) >= 1 + ( N - 2) * X( K, 1)
);
END
用Lingo编的,条件是：从城市O出发经过城市A, B, C, D,E,F后又回到O的过程。

约束：第n个经过城市A，则第n+1或n-1个经过城市B；同理C和D，E和F皆是如此。

该程序可以运行，但结果不符合约束条件，求正解。