1. reset;
2.  
3. option solver CPLEX;
4.  
5. var s1 >=0 , binary;
6. var s2 >=0 , binary;
7. var s3 >=0 , binary;
8. var s4 >=0 , binary;
9. var s5 >=0 , binary;
10. var s6 >=0 , binary;
11.  
12. var xA >=0 , binary;
13. var xB >=0 , binary;
14. var xC >=0 , binary;
15. var xD >=0 , binary;
16. var xE >=0 , binary;
17. var xF >=0 , binary;
18.  
19. var x1dB >=0 , binary;
20. var x1dD >=0 , binary;
21. var x2dA >=0 , binary;
22. var x2dE >=0 , binary;
23. var x3dC >=0 , binary;
24. var x3dD >=0 , binary;
25. var x3dF >=0 , binary;
26. var x4dA >=0 , binary;
27. var x4dC >=0 , binary;
28. var x4dE >=0 , binary;
29. var x5dB >=0 , binary;
30. var x5dD >=0 , binary;
31. var x6dB >=0 , binary;
32. var x6dC >=0 , binary;
33. var x6dF >=0 , binary;
34. maximize fcelu:
35.     s1+s2+s3+s4+s5+s6;
36.    
37. subject to ogr1:
38.     s1<=1;
39. subject to ogr2:
40.     s2<=1;
41. subject to ogr3:
42.     s3<=1;
43. subject to ogr4:
44.     s4<=1;
45. subject to ogr5:
46.     s5<=1;    
47. subject to ogr6:
48.     s6<=1;
49. subject to ogr7:
50.     xA<=1;
51. subject to ogr8:
52.     xB<=1;
53. subject to ogr9:
54.     xC<=1;
55. subject to ogr10:
56.     xD<=1;
57. subject to ogr11:
58.     xE<=1;
59. subject to ogr12:
60.     xF<=1;
61. subject to ogrA:
62.     x1dB<=1;
63. subject to ogrB:
64.     x1dD<=1;
65. subject to ogrC:
66.     x2dA<=1;
67. subject to ogrD:
68.     x2dE<=1;
69. subject to ogrE:
70.     x3dC<=1;      
71. subject to ogrF:
72.     x3dD<=1;
73. subject to ogrG:
74.     x3dF<=1;
75. subject to ogrI:
76.     x4dA<=1;
77. subject to ogrJ:
78.     x4dC<=1;
79. subject to ogrK:
80.     x4dE<=1;
81. subject to ogrL:
82.     x5dB<=1;
83. subject to ogrM:
84.     x5dD<=1;
85. subject to ogrO:
86.     x6dB<=1;
87. subject to ogrP:
88.     x6dC<=1;
89. subject to ogrR:
90.     x6dF<=1;
91. subject to ogr20:
92.     s1-x1dB-x1dD = 0;
93. subject to ogr21:
94.     s2-x2dA-x2dE = 0;
95. subject to ogr22:
96.     s3-x3dC-x3dD-x3dF = 0;
97. subject to ogr23:
98.     s4-x4dA-x4dC-x4dE = 0;
99. subject to ogr24:
100.     s5-x5dB-x5dD = 0;
101. subject to ogr25:
102.     s6-x6dB-x6dC-x6dF = 0;
103.    
104. subject to ogr26:
105.     xA-x2dA-x4dA = 0;
106. subject to ogr27:
107.     xB-x1dB-x5dB-x6dB = 0;
108. subject to ogr28:
109.     xC-x3dC-x4dC-x6dC = 0;
110. subject to ogr29:
111.     xD-x1dD-x3dD-x5dD = 0;
112. subject to ogr30:
113.     xE-x2dE-x4dE = 0;
114. subject to ogr31:
115.     xF-x3dF-x6dF = 0;
116. subject to ogrSpecjalne:
117.     3*x3dC+13*x3dD+17*x3dF-30*x5dB-7*x5dD-18*x6dB-8*x6dC-21*x6dF >= 0;
118.  
119. solve;
120. display fcelu, x1dB, x1dD, x2dA, x2dE, x3dC, x3dD, x3dF, x4dA, x4dC, x4dE, x5dB, x5dD, x6dB, x6dC, x6dF;