请教一段lingo程序，多谢各位大侠帮忙啦~~帮顶也好~~~

hyltt7636

试试这个吧
sets:
supply/1..7/:p,s;
need/1..15/:y,t,l;
link(supply,need):x,c;
endsets
data:
s=800,800,1000,2000,2000,2000,3000;
p=160,155,155,160,155,150,160;
c=264.3 160.3 140.2 98.6 38.0 20.5 3.1 21.2 64.2 92.0 96.0 106.0 121.2 128.0 142.0
  309.3 205.3 190.2 171.6 111.0 95.5 86.0 71.2 114.2 142.0 146.0 156.0 171.2 178.0 192.0
  324.3 220.3 200.2 181.6 121.0 105.5 96.0 86.2 48.2 82.0 86.0 96.0 111.2 118.0 132.0
  354.3 250.3 235.2 216.6 156.0 140.5 131.0 116.2 84.2 62.0 51.0 61.0 76.2 83.0 97.0
  349.3 245.3 225.2 206.6 146.0 130.5 121.0 111.2 79.2 57.0 33.0 51.0 71.2 73.0 87.0
  359.3 255.3 235.2 216.6 156.0 140.5 131.0 121.2 84.2 62.0 51.0 45.0 26.2 11.0 28.0
  369.3 265.3 245.2 226.6 166.0 150.5 141.0 131.2 99.2 77.0 66.0 56.0 38.2 26.0 2.0;
enddata
init:
endinit
min=@sum(link(i,j)p(i)+c(i,j))*x(i,j))+0.05*@sum(need(i):y*(y+1)+(l-y)*(l-y-1));
@for(supply(i)sum(need(j):x(i,j))>=500*t(i));
      @for(supply(i)sum(need(j):x(i,j))<=s(i)*t(i));
      @for(supply(i)bin(t(i)));
@for(need(j)|j#NE#1 #AND# j#NE#15sum(supply(i):x(i,j))=y(j)+(l(j-1)-y(j-1)));
@sum(link(i,j)|j#EQ#15:x(i,j))=l(14)-y(14);
y(1)=0;
@for(need(i):y(i)<=l(i));
@for(need(i):y(i)>=0);
@for(link(i,j):x(i,j)>=0);
@sum(link(i,j):x(i,j))=5171;
end

hyltt7636

试试这个
sets:
supply/1..7/:p,s;
need/1..15/:y,t,l;
link(supply,need):x,c;
endsets
data:
s=800,800,1000,2000,2000,2000,3000;
p=160,155,155,160,155,150,160;
c=264.3 160.3 140.2 98.6 38.0 20.5 3.1 21.2 64.2 92.0 96.0 106.0 121.2 128.0 142.0
  309.3 205.3 190.2 171.6 111.0 95.5 86.0 71.2 114.2 142.0 146.0 156.0 171.2 178.0 192.0
  324.3 220.3 200.2 181.6 121.0 105.5 96.0 86.2 48.2 82.0 86.0 96.0 111.2 118.0 132.0
  354.3 250.3 235.2 216.6 156.0 140.5 131.0 116.2 84.2 62.0 51.0 61.0 76.2 83.0 97.0
  349.3 245.3 225.2 206.6 146.0 130.5 121.0 111.2 79.2 57.0 33.0 51.0 71.2 73.0 87.0
  359.3 255.3 235.2 216.6 156.0 140.5 131.0 121.2 84.2 62.0 51.0 45.0 26.2 11.0 28.0
  369.3 265.3 245.2 226.6 166.0 150.5 141.0 131.2 99.2 77.0 66.0 56.0 38.2 26.0 2.0;
enddata
init:
endinit
min=@sum(links(i,j)p(i)+c(i,j))*x(i,j))+0.05*@sum(need(i):y*(y+1)+(l-y)*(l-y-1))
@for(supply(i)sum(need(j):x(i,j))>=500*t(i));
      @for(supply(i)sum(need(j):x(i,j))<=s(i)*t(i));
      @for(supply(i)bin(t(i)));
@for(need(j)|j#NE#1 #AND# j#NE#15sum(supply(i):x(i,j))=y(j)+(l(j-1)-y(j-1)));
@sum(link(i,j)|j#EQ#15:x(i,j))=l(14)-y(14);
y(1)=0;
@for(need(i):y(i)<=l(i));
@for(need(i):y(i)>=0);
@for(link(i,j):x(i,j)>=0);
@sum(link(i,j):x(i,j))=5171;
end