lingo 程序

k/1..5/

z/1..8/:m

link(k,z):a,c,l,s

endsets

min=@sum(link:a)

@sum(z(j):a(4,j)*c(4,j))<=150

@sum(z(j):a(5,j)*c(5,j))<=400

@sum(z(j):a(1,j)*c(1,j))<=360+40*t

@sum(z(j):a(2,j)*c(2,j))<=600+30*t

@sum(z(j):a(3,j)*c(3,j))<=500+20*t

@for(z(j):@sum(k(i):a(i,j))=m(j))

t=@max(link(i,j):c(i,j)*(2*s(i,j)+2))

t<=10

data:

m=300 330 120 170 110 100 1000 700

s=3.08 2.5 2.6 5.74 4.28 6.2 2 3.86

1.16 3.14 4.12 5.06 2.36 5.52 1.96 2.96

4.48 6.64 6.74 2.9 3.28 1.86 3.34 2.04

4.8 4.5 5.4 8.1 8.2 6.55 5.74 4.32

1.84 2.58 3.92 3.74 3.84 3.98 3.82 2.36

enddata

程序如迟尘下(Lingo)

min=a1+a2+a3+a4+a5+a6+a7+a8

99*a1+139*a2+153*a3+167*a4+173*(a5+1)+185*a6+195*a7+237*a7+265*a8=753

@gin(a1)

@gin(a2)

@gin(a3)

@gin(a4)

@gin(a5)

@gin(a6)

@gin(a7)

@gin(a8)

结果的确没解。。雹旦搜。

但不限制173时源历，有组解如下：

A11.0000001.000000

A20.0000001.000000

A31.0000001.000000

A43.0000001.000000

A50.0000001.000000

A60.0000001.000000

A70.0000001.000000

A80.0000001.000000