0

我有一个线性规划模型，用于最小路径问题。这是模型：

/* Min path problem

file: minPath.mod */

set V;
set E within V cross V;

param cost{E};
param S symbolic;
param T symbolic;

var flow{E} integer, >= 0;


minimize min_path: sum{(a,b) in E} cost[a,b] * flow[a,b];

s.t. conservazione{v in V: v != S and v != T}:
    sum{(a,b) in E: a == v} flow[a,b] == 
    sum{(a,b) in E: b == v} flow[a,b];
s.t. sorgente: sum{(a,b) in E: a == S} flow[a,b] == 1;
s.t. destinazione: sum{(a,b) in E: b == T} flow[a,b] == 1;

display {(a,b) in E} flow[a,b];

data;

set V := A B C D E;
set E := (A,B) (A,C) (B,D) (B,E) (C,D) (D,E);

param S := A;
param T := D;

param cost := [A,B] 2 [A,C] 1 [B,D] 3 [B,E] 1 [C,D] 1 [D,E] 1;

end;

我的示例目标值为 3，最小路径为：

A -> C -> D -> E

出于这个原因，向量流在边缘上必须为 1，我在上面写过。顺便说一句，当我用语句显示向量流时：

display {(a,b) in E} flow[a,b];

向量在所有位置都是 0。

flow[A,B].val = 0
flow[A,C].val = 0
flow[B,D].val = 0
flow[B,E].val = 0
flow[C,D].val = 0
flow[D,E].val = 0

我试图更改语法，但我无法强制 glpsol 打印实际值。

我错过了什么吗？

4

1 回答 1

1

你必须把声明：

solve;

前显示声明！

于 2013-11-04T15:11:38.873 回答