数学建模社区-数学中国

标题: 求高手解决一个lingo程序 [打印本页]

作者: LWQGG 时间: 2012-12-16 04:10
标题: 求高手解决一个lingo程序
min=z;
SETS:
s/1..250/:k;
b/1..27/:l,T;
links(s,b):a;
ENDSETS
DATA:
T=30 41 40 70 10 80 32 24 50 80 10 44 90 32 21 80 36 45 70 70 90 42 44 50 70 90 12;
ENDDATA
z>200;
@for(s(i): @for(b(j):@bin(a(i,j))));
@FOR(b(j): @sum(s(i):i*a(i,j))+T(j)<=z);
@FOR(b(j): @sum(s(i):a(i,j))=1);
@for(s(i): @for(s(k): (i*a(i,17)-k*a(k,18)-45)*(k*a(k,18)-i*a(i,17)-36)<=0));
@for(s(i): @for(s(k): (i*a(i,16)-k*a(k,17)-36)*(k*a(k,17)-i*a(i,16)-80)<=0));
@for(s(i): @for(s(k): (i*a(i,16)-k*a(k,18)-45)*(k*a(k,18)-i*a(i,16)-80)<=0));
@for(s(i): @for(s(k): (i*a(i,16)-k*a(k,15)-21)*(k*a(k,15)-i*a(i,16)-80)<=0));
@for(s(i): @for(s(k): (i*a(i,14)-k*a(k,15)-21)*(k*a(k,15)-i*a(i,14)-32)<=0));
@for(s(i): @for(s(k): (i*a(i,14)-k*a(k,16)-80)*(k*a(k,16)-i*a(i,14)-32)<=0));
@for(s(i): @for(s(k): (i*a(i,14)-k*a(k,13)-90)*(k*a(k,13)-i*a(i,14)-32)<=0));
@for(s(i): @for(s(k): (i*a(i,14)-k*a(k,19)-70)*(k*a(k,19)-i*a(i,14)-32)<=0));
@for(s(i): @for(s(k): (i*a(i,19)-k*a(k,25)-70)*(k*a(k,25)-i*a(i,19)-70)<=0));
@for(s(i): @for(s(k): (i*a(i,19)-k*a(k,20)-70)*(k*a(k,20)-i*a(i,19)-70)<=0));
@for(s(i): @for(s(k): (i*a(i,25)-k*a(k,26)-90)*(k*a(k,26)-i*a(i,25)-70)<=0));
@for(s(i): @for(s(k): (i*a(i,25)-k*a(k,27)-12)*(k*a(k,27)-i*a(i,25)-70)<=0));
@for(s(i): @for(s(k): (i*a(i,26)-k*a(k,27)-12)*(k*a(k,27)-i*a(i,26)-90)<=0));
@for(s(i): @for(s(k): (i*a(i,25)-k*a(k,20)-70)*(k*a(k,20)-i*a(i,25)-70)<=0));
@for(s(i): @for(s(k): (i*a(i,20)-k*a(k,22)-42)*(k*a(k,22)-i*a(i,20)-70)<=0));
@for(s(i): @for(s(k): (i*a(i,20)-k*a(k,21)-90)*(k*a(k,21)-i*a(i,20)-70)<=0));
@for(s(i): @for(s(k): (i*a(i,21)-k*a(k,24)-50)*(k*a(k,24)-i*a(i,21)-90)<=0));
@for(s(i): @for(s(k): (i*a(i,21)-k*a(k,23)-44)*(k*a(k,23)-i*a(i,21)-90)<=0));
@for(s(i): @for(s(k): (i*a(i,23)-k*a(k,24)-50)*(k*a(k,24)-i*a(i,23)-44)<=0));
@for(s(i): @for(s(k): (i*a(i,21)-k*a(k,22)-42)*(k*a(k,22)-i*a(i,21)-90)<=0));
@for(s(i): @for(s(k): (i*a(i,13)-k*a(k,7)-32)*(k*a(k,7)-i*a(i,13)-90)<=0));
@for(s(i): @for(s(k): (i*a(i,13)-k*a(k,8)-24)*(k*a(k,8)-i*a(i,13)-90)<=0));
@for(s(i): @for(s(k): (i*a(i,8)-k*a(k,7)-32)*(k*a(k,7)-i*a(i,8)-24)<=0));
@for(s(i): @for(s(k): (i*a(i,8)-k*a(k,9)-50)*(k*a(k,9)-i*a(i,8)-24)<=0));
@for(s(i): @for(s(k): (i*a(i,8)-k*a(k,10)-80)*(k*a(k,10)-i*a(i,8)-24)<=0));
@for(s(i): @for(s(k): (i*a(i,9)-k*a(k,10)-80)*(k*a(k,10)-i*a(i,9)-50)<=0));
@for(s(i): @for(s(k): (i*a(i,10)-k*a(k,11)-10)*(k*a(k,11)-i*a(i,10)-80)<=0));
@for(s(i): @for(s(k): (i*a(i,10)-k*a(k,12)-44)*(k*a(k,12)-i*a(i,10)-80)<=0));
@for(s(i): @for(s(k): (i*a(i,11)-k*a(k,12)-44)*(k*a(k,12)-i*a(i,11)-10)<=0));
@for(s(i): @for(s(k): (i*a(i,7)-k*a(k,5)-10)*(k*a(k,5)-i*a(i,7)-32)<=0));
@for(s(i): @for(s(k): (i*a(i,7)-k*a(k,6)-80)*(k*a(k,6)-i*a(i,7)-32)<=0));
@for(s(i): @for(s(k): (i*a(i,5)-k*a(k,6)-80)*(k*a(k,6)-i*a(i,5)-10)<=0));
@for(s(i): @for(s(k): (i*a(i,5)-k*a(k,4)-70)*(k*a(k,4)-i*a(i,5)-10)<=0));
@for(s(i): @for(s(k): (i*a(i,5)-k*a(k,3)-40)*(k*a(k,3)-i*a(i,5)-10)<=0));
@for(s(i): @for(s(k): (i*a(i,3)-k*a(k,4)-70)*(k*a(k,4)-i*a(i,3)-40)<=0));
@for(s(i): @for(s(k): (i*a(i,6)-k*a(k,2)-41)*(k*a(k,2)-i*a(i,6)-80)<=0));
@for(s(i): @for(s(k): (i*a(i,6)-k*a(k,1)-30)*(k*a(k,1)-i*a(i,6)-80)<=0));
@for(s(i): @for(s(k): (i*a(i,1)-k*a(k,2)-41)*(k*a(k,2)-i*a(i,1)-30)<=0));
@for(s(i): @for(s(k): (i*a(i,13)-k*a(k,19)-70)*(k*a(k,19)-i*a(i,13)-90)<=0));
end

作者: madio 时间: 2012-12-16 12:21
你把这句s/1..250/:k;改成
s/1..250/:kk;
因为你后面用k做了指标变量，这里就不能再做为集的对象名字了，这样的重名是不允许的

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)