model:
sets:
inputs/ip1..ip3/;
badoutputs/bo1..bo2/;
provinces/p1..p5/;
links1(inputs,provinces):xt;
links2(badoutputs,provinces):bt;
goodoutputs(provinces):yt;
weighting(provinces):zt;
endsets
max=y_t;
@sum(provinces(K):zt(K)*yt(K))>=y_t;
@for(provinces(K):@for(badoutputs(J):@sum(provinces(K):zt(K)*bt(J,K))=bt(J,K)));
@for(provinces(K):@for(inputs(N):@sum(provinces(K):zt(K)*xt(N,K))<=xt(N,K)));
@for(provinces(K):zt(K)>=0);
data:
xt=337.33 240.41 475.86 309.73 189.76

566.50 455.98 2555.43 1154.10 856.60

2291.00 1680.00 4918.00 4395.00 1871.00
bt=30.97 21.99 63.46 47.45 37.99

6087.17 4254.68 13769.51 12185.54 6108.83
yt= 211.91
156.13
325.10
170.34
127.04
；

enddata
end