各位大师看一下,这个程序哪儿不对啊,谢谢了!!
Model:
Sets:
cities/1, 2, 3, 4, 5, 6, 7/;
mode/1, 2, 3/;
link1 (cities, cities, mode): poss, c, r, t, x, U;
link2 (cities, mode, mode): qoss, h, e, v, y;
Endsets
Min= @sum(link1:poss*c*x)+ @sum(link2:qoss*h*y);
@ sum(link1:poss*t*x)+ @ sum(link2:qoss*e*y) <= 96 ;
@for(link1:
@prod(link1(i, j, k)|poss(i, j, k)*x(i, j, k) #eq# 1: r(i, j, k))>=0.99);
@for(link2:
@prod(link2:1-qoss(j, k, l)*v(j, k, l)*y(j, k, l))>=0.975);
@for(link1:
Quantity<=poss(i, j, k)*U( i, j, k));
@for(cities(i):
@for(cities(j):
@sum(mode(k): poss(i, j, k)*x(i, j, k)) <=1));
@for(cities(i)| i# ne # 1 # and # i# ne # n:
@sum(mode(k):
@sum(mode(l):qoss(j, k, l)*y(j, k, l)))<=1);
@for(cities(i)| i #eq# 1:
@sum(cities(j):
@sum(mode(k):poss(i, j, k)*x(i, j, k)))=1);
@for(cities(j)| j # eq # 7:
@sum(cities(i):
@sum(mode(k): poss(i, j, k)*x(i, j, k)))=1);
@for(cities(i)| i # ne # 1 # and # i # ne# 7:
@sum(cities(j):
@sum(mode(k):poss(i, j, k)*x(i, j, k)))<=1);
@for(cities(i)| i # ne # 1 # and # i # ne# 7:
@sum(cities(j):
@sum(mode(k):poss(j, i, k)*x(i, j, k)))<=1);
@for(cities(i)| i # ne # 1 # and # i # ne# 7:
@sum(cities(j):
@sum(mode(k):poss(i, j, k)*x(i, j, k)))=
@sum(cities(j):
@sum(mode(k):poss(j, i, k)*x(i, j, k))));
@bin(x(i, j, k));
@bin(y(j, k, l));