linux - AMPL 脚本错误

Question

这里有一个 AMPL 脚本(TSP 成本问题)，但我不明白错误。

这里是代码:

param n, integer, >= 3;
/* number of nodes */

param MAX_TIME := 600;
param MAX_X :=20;
param maxspeed := 1;

set CITIES := 1..n;
/* set of nodes */

param p{CITIES};
param a{CITIES};
param b{CITIES};

set ARCS, within CITIES cross CITIES;
/* set of arcs */

param DIST{(i,j) in ARCS};
/* distance from node i to node j */
param NEXTC, integer, >= 0;
/* Next city after i in the solution */
var x{(i,j) in ARCS}, binary;
/* x[i,j] = 1 means that the salesman goes from node i to node j */


param bigM := 5000;
var tar{CITIES};         /*arrival */
var tlv{CITIES};         /* departure */
var tea{CITIES} >= 0;    /* early arrival (arrival before the designated   time window) */
var tla{CITIES} >= 0;    /* late arrival (arrival after the designated time window) */
var ted{CITIES} >= 0;    /* early departure (departure before the designated time window) */
var tld{CITIES} >= 0;    /* late departure (departure after the designated time window) */

set days := 1..5;
param root := 1;



var y{(i,j) in ARCS}, >= 0;
/* y[i,j] is a flow through arc (i,j) */


minimize total: sum{(i,j) in ARCS} DIST[i,j] * x[i,j];

s.t. t1 {i in CITIES} : tlv[i] >= tar[i]; 
s.t. t2 {i in CITIES, j in CITIES} : tar[j] >= tlv[i] + p[i] + DIST[i,j]/maxspeed - bigM*(1-x[i,j]);
s.t. t3 {i in CITIES} : tea[i] >= a[i] - tar[i];   /* early arrival */
s.t. t4 {i in CITIES} : tla[i] >= tar[i] - b[i];   /* late arrival */
s.t. t5 {i in CITIES} : ted[i] >= a[i] - tlv[i];   /* early departure */
s.t. t6 {i in CITIES} : tld[i] >= tlv[i] - b[i];   /* late departure */


data;

param n := 9 ;

param : ARCS : DIST :=
1  2   21
1  3   8
1  4   6
1  5   10
1  6   2
1  7   4
1  8   5
1  9   5
2  1   21
2  3   18
2  4   21
2  5   23
2  6   22
2  7   20
2  8   22
2  9   25
3  1   7
3  2   16
3  4   3
3  5   9
3  6   8
3  7   4
3  8   7
3  9   9
4  1   8
4  2   18
4  3   3
4  5   10
4  6   9
4  7   6
4  8   9
4  9   11
5  1   9
5  2   25
5  3   11
5  4   9
5  6   11
5  7   8
5  8   11
5  9   13
6  1   4
6  2   23
6  3   8
6  4   7
6  5   11
6  7   5
6  8   7
6  9   8
7  1   3
7  2   20
7  3   6
7  4   5
7  5   8
7  6   4
7  8   3
7  9   8
8  1   3
8  2   20
8  3   5
8  4   4
8  5   8
8  6   5
8  7   4
8  9   7
9  1   7
9  2   24
9  3   10
9  4   8
9  5   13
9  6   7
9  7   7
9  8   6
;

param : a :=
  1  705
  2  420
  3  420
  4  420
  5  420
  6  420
  7  420
  8  420
  9  420
;

param : b :=
  1  785
  2  795
  3  725
  4  500
  5  785
  6  785
  7  430
  8  785
  9  785
;

param : p :=
  1  65
  2  55
  3  125
  4  65
  5  65
  6  65
  7  65
  8  65
  9  65
;

option solver './cplex';

solve;

代码来自(神奇!)Stackoverflow。

在我的 Kubuntu-pc 中，使用 AMPL 和 cplex，我看到这个错误:

Error at _cmdno 2 executing "solve" command (file tsp_nodi_archi_costo.run, line 170, offset 2582): error processing constraint t2[1,1]: invalid subscript DIST[1,1]

我还有一个 .run 文件，它们运行良好，但这个不行!为什么？谢谢!

Answer 1

请注意，您设置的 ARCS 不包含 (1,1), (2,2), ....

因此 DIST[i,i] 未定义，并且 x[i,i] 是无效索引，因为您将 ARCS 声明为 x 的索引集。

但是，正如错误消息所建议的那样，您的约束 t2 在 {i in CITIES, j in CITIES} 上建立了索引，并且它将尝试创建的第一种情况是 i=1, j=1，这会导致无效下标。

您可以通过将 t2 中的索引更改为:{i in CITIES, j in CITIES: j <> i} 来解决此问题。

编辑:我注意到其他一些问题:

• “设置 ARCS”后有一个逗号，它不应该出现在此处。
• 您似乎缺少一些约束(如果您尝试运行它并检查解决方案，您就会明白我的意思)
• 您有一个参数“NEXTC”，它被标识为“解决方案中 i 之后的下一个城市”。大概这应该是一个 var 而不是 param。它似乎还缺少一个索引，并且程序的其余部分中没有任何内容引用它。这可能与缺少约束问题有关。