Lingo与最短路问题

Lingo与最短路问题

代码如下：

!最短路问题;
model:
data:
  n=10;
enddata
sets:
  cities/1..n/: F;  !10个城市;
  roads(cities,cities)/
    1,2  1,3
    2,4  2,5  2,6
    3,4  3,5  3,6
    4,7  4,8
    5,7  5,8  5,9
    6,8  6,9
    7,10
    8,10
    9,10
  /: D, P;
endsets
data:
  D=
    6  5
    3  6  9
    7  5  11
    9  1
    8  7  5
    4  10
    5
    7
    9;
enddata
  F(n)=0;
  @for(cities(i) | i #lt# n:
    F(i)=@min(roads(i,j): D(i,j)+F(j));
  );
  !显然，如果P(i,j)=1,则点i到点n的最短路径的第一步是i --> j，否则就不是。
   由此，我们就可方便的确定出最短路径;
  @for(roads(i,j):
    P(i,j)=@if(F(i) #eq# D(i,j)+F(j),1,0)
  );
end