如何用lingo10编写这题的代码

Min z=∑∑CijXij (求和，i=1,…6, j=1,…5)

s.t. 2*yij<=xij<=aij*yij, i=1,…,6,j=1,…5
∑xij>=8,分别对i=1,…4 ,求j=1, …5的和
∑xij>=7, 分别对i=5,6时, j=1, …5的和
∑xij=14, j=1,…5, i=1,…,6的和
∑yij<=3, i=1,…6 , j=1, …5的和
∑yij<=3, j=1,…5, i=1,…,6的和
y5j+y6j>=1, j=1,…5
xij>=0, yij=0 或1, i=1,…6, j=1,…5

谢谢这位大哥了，不好意思我把模型弄错了。
原模型是：

Min z=∑∑CiXij (求和，i=1,…6, j=1,…5) (Ci=10,10,9.9,9.8,10.8,11.3)

s.t. 2*yij<=xij<=aij*yij, i=1,…,6,j=1,…5
∑xij>=8,分别对i=1,…4 ,求j=1, …5的和
∑xij>=7, 分别对i=5,6时, j=1, …5的和
∑xij=14, j=1,…5, i=1,…,6的和
∑yij<=3, i=1,…6 , j=1, …5的和
∑yij<=3, j=1,…5, i=1,…,6的和
y5j+y6j>=1, j=1,…5
xij>=0, yij=0 或1, i=1,…6, j=1,…5

真的很不好意思了,把一个不完整的模型麻烦大家.下面是完整的模型:

Min z=∑∑CiXij (求和，i=1,…6, j=1,…5) (Ci=10,10,9.9,9.8,10.8,11.3)

s.t. 2*yij<=xij<=aij*yij, i=1,…,6,j=1,…5
∑xij>=8,分别对i=1,…4 ,求j=1, …5的和
∑xij>=7, 分别对i=5,6时, j=1, …5的和
∑xij=14, j=1,…5, i=1,…,6的和
∑yij<=3, i=1,…6 , j=1, …5的和
∑yij<=3, j=1,…5, i=1,…,6的和
y5j+y6j>=1, j=1,…5
xij>=0, yij=0 或1, i=1,…6, j=1,…5

aij的数据由一个矩阵给出(i是行,j是列):
6 0 6 0 7
0 6 0 6 0
4 8 3 0 5
5 5 6 0 4
3 0 4 8 0
0 6 0 6 3

model:

sets:
Ci/1..6/ : c;
Cj/1..5/;
links(Ci,Cj) : x,y,a;!分别是系数a；变量x，y;
endsets

data:
c=10 10 9.9 9.8 10.8 11.3; !给c，a赋予初值，具体的自己输入;

a=6 0 6 0 7
0 6 0 6 0
4 8 3 0 5
5 5 6 0 4
3 0 4 8 0
0 6 0 6 3;
enddata

min=@sum(links(i,j) : c(i)*x);

@for(links : 2*y<=x);
@for(links : x<=a*y);
@for(Ci(i) | i #le# 4 : @sum(links(i,j) : x)>=8);
@for(Ci(i) | i #ge# 5 : @sum(links(i,j) : x)>=7);
@for(Cj(j) : @sum(links(i,j) : x)=14);
@for(Ci(i) : @sum(links(i,j) : y)<=3);
@for(Cj(j) : @sum(links(i,j) : y)<=3);
@sum(links(i,j) | i #ge# 5 : y)>=1;
@for(links : x>=0);
@for(links : @bin(y));

最优解是709.9 运行报告如下：
Global optimal solution found.
Objective value: 709.9000
Extended solver steps: 2
Total solver iterations: 49

Variable Value Reduced Cost
C( 1) 10.00000 0.000000
C( 2) 10.00000 0.000000
C( 3) 9.900000 0.000000
C( 4) 9.800000 0.000000
C( 5) 10.80000 0.000000
C( 6) 11.30000 0.000000
X( 1, 1) 5.000000 0.000000
X( 1, 2) 0.000000 0.1000000
X( 1, 3) 5.000000 0.000000
X( 1, 4) 0.000000 0.000000
X( 1, 5) 6.000000 0.000000
X( 2, 1) 0.000000 0.000000
X( 2, 2) 2.000000 0.000000
X( 2, 3) 0.000000 0.000000
X( 2, 4) 6.000000 0.000000
X( 2, 5) 0.000000 0.000000
X( 3, 1) 4.000000 0.000000
X( 3, 2) 7.000000 0.000000
X( 3, 3) 0.000000 0.000000
X( 3, 4) 0.000000 0.000000
X( 3, 5) 5.000000 0.000000
X( 4, 1) 5.000000 0.000000
X( 4, 2) 5.000000 0.000000
X( 4, 3) 6.000000 0.000000
X( 4, 4) 0.000000 0.000000
X( 4, 5) 0.000000 0.000000
X( 5, 1) 0.000000 0.000000
X( 5, 2) 0.000000 0.1000000
X( 5, 3) 3.000000 0.000000
X( 5, 4) 4.000000 0.000000
X( 5, 5) 0.000000 0.000000
X( 6, 1) 0.000000 0.000000
X( 6, 2) 0.000000 0.1000000
X( 6, 3) 0.000000 0.000000
X( 6, 4) 4.000000 0.000000
X( 6, 5) 3.000000 0.000000
Y( 1, 1) 1.000000 0.000000
Y( 1, 2) 0.000000 0.000000
Y( 1, 3) 1.000000 0.000000
Y( 1, 4) 0.000000 0.000000
Y( 1, 5) 1.000000 0.000000
Y( 2, 1) 0.000000 0.000000
Y( 2, 2) 1.000000 0.2000000
Y( 2, 3) 0.000000 0.000000
Y( 2, 4) 1.000000 0.000000
Y( 2, 5) 0.000000 0.000000
Y( 3, 1) 1.000000 -0.4000000
Y( 3, 2) 1.000000 0.000000
Y( 3, 3) 0.000000 -0.3000000
Y( 3, 4) 0.000000 0.000000
Y( 3, 5) 1.000000 -0.5000000
Y( 4, 1) 1.000000 -1.000000
Y( 4, 2) 1.000000 -0.5000000
Y( 4, 3) 1.000000 -1.200000
Y( 4, 4) 0.000000 0.000000
Y( 4, 5) 0.000000 -0.8000000
Y( 5, 1) 0.000000 0.000000
Y( 5, 2) 0.000000 0.000000
Y( 5, 3) 1.000000 0.000000
Y( 5, 4) 1.000000 0.000000
Y( 5, 5) 0.000000 0.000000
Y( 6, 1) 0.000000 0.000000
Y( 6, 2) 0.000000 0.000000
Y( 6, 3) 0.000000 0.000000
Y( 6, 4) 1.000000 0.000000
Y( 6, 5) 1.000000 0.000000
A( 1, 1) 6.000000 0.000000
A( 1, 2) 0.000000 0.000000
A( 1, 3) 6.000000 0.000000
A( 1, 4) 0.000000 0.000000
A( 1, 5) 7.000000 0.000000
A( 2, 1) 0.000000 0.000000
A( 2, 2) 6.000000 0.000000
A( 2, 3) 0.000000 0.000000
A( 2, 4) 6.000000 0.000000
A( 2, 5) 0.000000 0.000000
A( 3, 1) 4.000000 0.000000
A( 3, 2) 8.000000 0.000000
A( 3, 3) 3.000000 0.000000
A( 3, 4) 0.000000 0.000000
A( 3, 5) 5.000000 0.000000
A( 4, 1) 5.000000 0.000000
A( 4, 2) 5.000000 0.000000
A( 4, 3) 6.000000 0.000000
A( 4, 4) 0.000000 0.000000
A( 4, 5) 4.000000 0.000000
A( 5, 1) 3.000000 0.000000
A( 5, 2) 0.000000 0.000000
A( 5, 3) 4.000000 0.000000
A( 5, 4) 8.000000 0.000000
A( 5, 5) 0.000000 0.000000
A( 6, 1) 0.000000 0.000000
A( 6, 2) 6.000000 0.000000
A( 6, 3) 0.000000 0.000000
A( 6, 4) 6.000000 0.000000
A( 6, 5) 3.000000 0.000000

Row Slack or Surplus Dual Price
1 709.9000 -1.000000
2 3.000000 0.000000
3 0.000000 0.000000
4 3.000000 0.000000
5 0.000000 0.000000
6 4.000000 0.000000
7 0.000000 0.000000
8 0.000000 0.1000000
9 0.000000 0.000000
10 4.000000 0.000000
11 0.000000 0.000000
12 2.000000 0.000000
13 5.000000 0.000000
14 0.000000 0.000000
15 0.000000 0.000000
16 3.000000 0.000000
17 3.000000 0.000000
18 3.000000 0.000000
19 4.000000 0.000000
20 0.000000 0.000000
21 0.000000 0.000000
22 0.000000 0.000000
23 0.000000 0.000000
24 1.000000 0.000000
25 2.000000 0.000000
26 0.000000 0.000000
27 0.000000 0.000000
28 0.000000 0.000000
29 0.000000 0.000000
30 2.000000 0.000000
31 1.000000 0.000000
32 1.000000 0.000000
33 0.000000 0.000000
34 1.000000 0.000000
35 0.000000 0.000000
36 1.000000 0.000000
37 0.000000 0.000000
38 4.000000 0.000000
39 0.000000 0.000000
40 0.000000 0.000000
41 0.000000 0.000000
42 0.000000 0.1000000
43 1.000000 0.000000
44 0.000000 0.1000000
45 0.000000 0.1000000
46 0.000000 0.1000000
47 0.000000 0.2000000
48 0.000000 0.1000000
49 0.000000 0.2000000
50 0.000000 0.2000000
51 0.000000 0.2000000
52 0.000000 0.000000
53 0.000000 0.000000
54 1.000000 0.000000
55 4.000000 0.000000
56 0.000000 0.000000
57 0.000000 0.000000
58 0.000000 0.000000
59 0.000000 0.000000
60 2.000000 0.000000
61 0.000000 0.000000
62 8.000000 0.000000
63 0.000000 0.000000
64 8.000000 0.000000
65 8.000000 0.000000
66 0.000000 -0.8000000
67 0.000000 -1.300000
68 0.000000 -10.00000
69 0.000000 -9.900000
70 0.000000 -10.00000
71 0.000000 -10.00000
72 0.000000 -10.00000
73 0.000000 0.000000
74 1.000000 0.000000
75 0.000000 0.000000
76 0.000000 0.000000
77 1.000000 0.000000
78 1.000000 0.000000
79 0.000000 0.000000
80 0.000000 0.000000
81 0.000000 0.000000
82 0.000000 0.000000
83 0.000000 0.000000
84 3.000000 0.000000
85 5.000000 0.000000
86 0.000000 0.000000
87 5.000000 0.000000
88 0.000000 0.000000
89 6.000000 0.000000
90 0.000000 0.000000
91 2.000000 0.000000
92 0.000000 0.000000
93 6.000000 0.000000
94 0.000000 0.000000
95 4.000000 0.000000
96 7.000000 0.000000
97 0.000000 0.000000
98 0.000000 0.000000
99 5.000000 0.000000
100 5.000000 0.000000
101 5.000000 0.000000
102 6.000000 0.000000
103 0.000000 0.000000
104 0.000000 0.000000
105 0.000000 0.000000
106 0.000000 0.000000
107 3.000000 0.000000
108 4.000000 0.000000
109 0.000000 0.000000
110 0.000000 0.000000
111 0.000000 0.000000
112 0.000000 0.000000
113 4.000000 0.000000
114 3.000000 0.000000

1

f

你好！ sets: one/1..6/:c; two/1..5/:; link(one,two):x,y,a; endsets data: C=10,10,9.9,9.8,10.8,11.3; enddata min=@sum(link(i,j): c(i)*x(i,j)); @for(link(i,j): x(i,j)>2*y(i,j)); @for(link(i,j): x(i,j)8); @sum(two(j):x(5,j))>7; @sum(two(j):x(6,j))>7; @for(two(j): @sum(link(i,j):x(i,j))=14); @for(one(i): @sum(link(i,j):y(i,j))<3); @for(two(j): @sum(link(i,j):y(i,j))<3); @sum(two(j):y(5,j))+@sum(two(j):y(6,j))>1; @for(link(i,j): @bin(y(i,j))); ==================================================================== 上面说得对，a的数据都没是个问题。本身模型不知道对不对 如有疑问，请追问。