Hello,
I have an error in my code.
Could you help me?
In the objective function, it is supposed to get the cost according to the truck or lorry. But I am getting an error. I wonder why?
The error:
obj …
z =E= SUM((i,k,r), x(i,k,r) * (cost_truck(i)(k **EQ 'truck')** + cost_tir(i)(k EQ ‘tir’)));
Full code:
SETS
i distributors /D2, D23, D24, D25, D32, D33, D34, D35, D37, D38, D39, D41, D42, D44, D48, D51, D53, D57, D60, D61, D68, D71, D73, D74, D75, D76, D80, D85, D90, D104, D105, D107, D112, D125, D128, D130, D132, D141/
k vehicles /truck, tir/
r routes /R1*R100/;
ALIAS (i,ip);
PARAMETERS
dist(i,ip) distance between distributors i and ip
/D2.D23 588, D2.D24 611, D2.D25 586, D2.D32 332, D2.D33 915, D2.D34 440, D2.D35 89, D2.D37 357, D2.D38 681, D2.D39 395, D2.D41 458, D2.D42 225, D2.D44 164, D2.D48 816, D2.D51 402, D2.D53 369, D2.D57 294, D2.D60 340, D2.D61 508, D2.D68 119, D2.D71 799, D2.D73 601, D2.D74 433, D2.D75 916, D2.D76 219, D2.D80 554, D2.D85 127, D2.D90 364, D2.D104 383, D2.D105 466, D2.D107 590, D2.D112 266, D2.D125 296, D2.D128 412, D2.D130 824, D2.D132 427, D2.D141 379,
D23.D24 22, D23.D25 71, D23.D32 459, D23.D33 1436, D23.D34 941, D23.D35 541, D23.D37 841, D23.D38 946, D23.D39 331, D23.D41 206, D23.D42 387, D23.D44 714, D23.D48 1069, D23.D51 852, D23.D53 618, D23.D57 844, D23.D60 507, D23.D61 752, D23.D68 459, D23.D71 1349, D23.D73 952, D23.D74 939, D23.D75 1409, D23.D76 359, D23.D80 1067, D23.D85 472, D23.D90 804, D23.D104 495, D23.D105 254, D23.D107 1140, D23.D112 701, D23.D125 304, D23.D128 962, D23.D130 1374, D23.D132 742, D23.D141 409,
D24.D25 50, D24.D32 474, D24.D33 1471, D24.D34 978, D24.D35 576, D24.D37 877, D24.D38 962, D24.D39 347, D24.D41 228, D24.D42 446, D24.D44 749, D24.D48 1084, D24.D51 987, D24.D53 634, D24.D57 879, D24.D60 522, D24.D61 768, D24.D68 492, D24.D71 1384, D24.D73 967, D24.D74 979, D24.D75 1459, D24.D76 392, D24.D80 1103, D24.D85 505, D24.D90 882, D24.D104 510, D24.D105 270, D24.D107 1174, D24.D112 705, D24.D125 362, D24.D128 997, D24.D130 1409, D24.D132 758, D24.D141 425,
D25.D32 524, D25.D33 1448, D25.D34 969, D25.D35 553, D25.D37 866, D25.D38 1079, D25.D39 396, D25.D41 277, D25.D42 424, D25.D44 726, D25.D48 1202, D25.D51 963, D25.D53 751, D25.D57 856, D25.D60 572, D25.D61 817, D25.D68 469, D25.D71 1361, D25.D73 1110, D25.D74 966, D25.D75 1436, D25.D76 369, D25.D80 1094, D25.D85 477, D25.D90 871, D25.D104 627, D25.D105 319, D25.D107 1151, D25.D112 742, D25.D125 345, D25.D128 974, D25.D130 1386, D25.D132 903, D25.D141 474,
D32.D33 1043, D32.D34 443, D32.D35 274, D32.D37 343, D32.D38 490, D32.D39 133, D32.D41 291, D32.D42 248, D32.D44 418, D32.D48 613, D32.D51 655, D32.D53 162, D32.D57 548, D32.D60 52, D32.D61 296, D32.D68 311, D32.D71 956, D32.D73 495, D32.D74 445, D32.D75 910, D32.D76 310, D32.D80 569, D32.D85 318, D32.D90 347, D32.D104 38, D32.D105 205, D32.D107 763, D32.D112 240, D32.D125 248, D32.D128 648, D32.D130 985, D32.D132 286, D32.D141 55,
D33.D34 612, D33.D35 922, D33.D37 703, D33.D38 618, D33.D39 1148, D33.D41 1295, D33.D42 1056, D33.D44 723, D33.D48 582, D33.D51 650, D33.D53 867, D33.D57 617, D33.D60 1001, D33.D61 733, D33.D68 981, D33.D71 124, D33.D73 558, D33.D74 669, D33.D75 233, D33.D76 1078, D33.D80 540, D33.D85 989, D33.D90 696, D33.D104 980, D33.D105 1219, D33.D107 299, D33.D112 818, D33.D125 1132, D33.D128 471, D33.D130 210, D33.D132 751, D33.D141 1069,
D34.D35 320, D34.D37 100, D34.D38 267, D34.D39 567, D34.D41 776, D34.D42 537, D34.D44 299, D34.D48 415, D34.D51 407, D34.D53 309, D34.D57 320, D34.D60 398, D34.D61 298, D34.D68 493, D34.D71 528, D34.D73 176, D34.D74 65, D34.D75 482, D34.D76 589, D34.D80 140, D34.D85 500, D34.D90 93, D34.D104 399, D34.D105 638, D34.D107 335, D34.D112 215, D34.D125 627, D34.D128 274, D34.D130 557, D34.D132 205, D34.D141 489,
D35.D37 317, D35.D38 632, D35.D39 318, D35.D41 401, D35.D42 162, D35.D44 196, D35.D48 754, D35.D51 433, D35.D53 308, D35.D57 326, D35.D60 252, D35.D61 446, D35.D68 78, D35.D71 831, D35.D73 558, D35.D74 419, D35.D75 885, D35.D76 183, D35.D80 543, D35.D85 92, D35.D90 322, D35.D104 299, D35.D105 403, D35.D107 621, D35.D112 208, D35.D125 238, D35.D128 444, D35.D130 856, D35.D132 365, D35.D141 292,
D37.D38 336, D37.D39 470, D37.D41 627, D37.D42 440, D37.D44 234, D37.D48 484, D37.D51 412, D37.D53 211, D37.D57 304, D37.D60 300, D37.D61 224, D37.D68 403, D37.D71 617, D37.D73 245, D37.D74 105, D37.D75 571, D37.D76 499, D37.D80 229, D37.D85 410, D37.D90 8, D37.D104 302, D37.D105 541, D37.D107 424, D37.D112 118, D37.D125 515, D37.D128 317, D37.D130 646, D37.D132 108, D37.D141 391,
D38.D39 620, D38.D41 778, D38.D42 724, D38.D44 565, D38.D48 149, D38.D51 690, D38.D53 339, D38.D57 583, D38.D60 490, D38.D61 195, D38.D68 712, D38.D71 602, D38.D73 82, D38.D74 323, D38.D75 512, D38.D76 808, D38.D80 305, D38.D85 720, D38.D90 330, D38.D104 453, D38.D105 692, D38.D107 409, D38.D112 426, D38.D125 735, D38.D128 446, D38.D130 631, D38.D132 289, D38.D141 542,
D39.D41 165, D39.D42 208, D39.D44 507, D39.D48 744, D39.D51 745, D39.D53 293, D39.D57 637, D39.D60 182, D39.D61 427, D39.D68 309, D39.D71 1087, D39.D73 627, D39.D74 576, D39.D75 1042, D39.D76 270, D39.D80 700, D39.D85 311, D39.D90 479, D39.D104 170, D39.D105 79, D39.D107 895, D39.D112 365, D39.D125 188, D39.D128 755, D39.D130 1116, D39.D132 417, D39.D141 84,
D41.D42 243, D41.D44 574, D41.D48 902, D41.D51 812, D41.D53 451, D41.D57 704, D41.D60 339, D41.D61 585, D41.D68 345, D41.D71 1209, D41.D73 785, D41.D74 781, D41.D75 1200, D41.D76 247, D41.D80 858, D41.D85 351, D41.D90 637, D41.D104 327, D41.D105 87, D41.D107 1000, D41.D112 523, D41.D125 164, D41.D128 822, D41.D130 1234, D41.D132 575, D41.D141 242,
D42.D44 331, D42.D48 775, D42.D51 551, D42.D53 303, D42.D57 443, D42.D60 248, D42.D61 414, D42.D68 176, D42.D71 964, D42.D73 589, D42.D74 480, D42.D75 955, D42.D76 162, D42.D80 601, D42.D85 70, D42.D90 399, D42.D104 248, D42.D105 330, D42.D107 746, D42.D112 259, D42.D125 79, D42.D128 572, D42.D130 930, D42.D132 366, D42.D141 200,
D44.D48 499, D44.D51 228, D44.D53 192, D44.D57 91, D44.D60 370, D44.D61 243, D44.D68 384, D44.D71 574, D44.D73 234, D44.D74 130, D44.D75 514, D44.D76 449, D44.D80 97, D44.D85 394, D44.D90 242, D44.D104 370, D44.D105 602, D44.D107 311, D44.D112 123, D44.D125 591, D44.D128 245, D44.D130 517, D44.D132 214, D44.D141 375,
D48.D51 591, D48.D53 249, D48.D57 463, D48.D60 615, D48.D61 167, D48.D68 642, D48.D71 558, D48.D73 179, D48.D74 378, D48.D75 468, D48.D76 738, D48.D80 263, D48.D85 649, D48.D90 327, D48.D104 589, D48.D105 828, D48.D107 393, D48.D112 389, D48.D125 871, D48.D128 469, D48.D130 590, D48.D132 362, D48.D141 624,
D51.D53 429, D51.D57 319, D51.D60 682, D51.D61 412, D51.D68 437, D51.D71 793, D51.D73 422, D51.D74 293, D51.D75 733, D51.D76 646, D51.D80 131, D51.D85 456, D51.D90 318, D51.D104 680, D51.D105 764, D51.D107 494, D51.D112 289, D51.D125 732, D51.D128 40, D51.D130 773, D51.D132 243, D51.D141 599,
D53.D57 181, D53.D60 138, D53.D61 167, D53.D68 229, D53.D71 741, D53.D73 358, D53.D74 261, D53.D75 681, D53.D76 345, D53.D80 290, D53.D85 242, D53.D90 183, D53.D104 153, D53.D105 382, D53.D107 447, D53.D112 107, D53.D125 372, D53.D128 389, D53.D130 638, D53.D132 114, D53.D141 153,
D57.D60 361, D57.D61 174, D57.D68 437, D57.D71 486, D57.D73 143, D57.D74 91, D57.D75 436, D57.D76 504, D57.D80 59, D57.D85 444, D57.D90 91, D57.D104 361, D57.D105 584, D57.D107 241, D57.D112 77, D57.D125 573, D57.D128 193, D57.D130 442, D57.D132 168, D57.D141 358,
D60.D61 334, D60.D68 230, D60.D71 935, D60.D73 479, D60.D74 428, D60.D75 900, D60.D76 275, D60.D80 552, D60.D85 253, D60.D90 319, D60.D104 39, D60.D105 221, D60.D107 727, D60.D112 231, D60.D125 264, D60.D128 610, D60.D130 949, D60.D132 277, D60.D141 69,
D61.D68 420, D61.D71 671, D61.D73 226, D61.D74 270, D61.D75 636, D61.D76 504, D61.D80 196, D61.D85 417, D61.D90 211, D61.D104 343, D61.D105 574, D61.D107 391, D61.D112 204, D61.D125 563, D61.D128 309, D61.D130 667, D61.D132 185, D61.D141 341,
D68.D71 851, D68.D73 522, D68.D74 387, D68.D75 816, D68.D76 238, D68.D80 449, D68.D85 251, D68.D90 410, D68.D104 242, D68.D105 410, D68.D107 643, D68.D112 251, D68.D125 401, D68.D128 478, D68.D130 866, D68.D132 348, D68.D141 315,
D71.D73 501, D71.D74 612, D71.D75 176, D71.D76 924, D71.D80 665, D71.D85 860, D71.D90 618, D71.D104 923, D71.D105 1166, D71.D107 256, D71.D112 726, D71.D125 1079, D71.D128 424, D71.D130 258, D71.D132 698, D71.D141 1023,
D73.D74 168, D73.D75 476, D73.D76 628, D73.D80 221, D73.D85 484, D73.D90 252, D73.D104 478, D73.D105 717, D73.D107 326, D73.D112 250, D73.D125 745, D73.D128 409, D73.D130 522, D73.D132 242, D73.D141 500,
D74.D75 532, D74.D76 536, D74.D80 122, D74.D85 463, D74.D90 113, D74.D104 428, D74.D105 669, D74.D107 329, D74.D112 197, D74.D125 697, D74.D128 296, D74.D130 525, D74.D132 185, D74.D141 452,
D75.D76 897, D75.D80 654, D75.D85 886, D75.D90 621, D75.D104 899, D75.D105 1158, D75.D107 306, D75.D112 693, D75.D125 1068, D75.D128 422, D75.D130 282, D75.D132 678, D75.D141 1012,
D76.D80 482, D76.D85 113, D76.D90 491, D76.D104 268, D76.D105 430, D76.D107 584, D76.D112 180, D76.D125 429, D76.D128 473, D76.D130 817, D76.D132 341, D76.D141 243,
D80.D85 461, D80.D90 136, D80.D104 520, D80.D105 759, D80.D107 301, D80.D112 174, D80.D125 787, D80.D128 216, D80.D130 477, D80.D132 229, D80.D141 520,
D85.D90 417, D85.D104 253, D85.D105 366, D85.D107 628, D85.D112 267, D85.D125 254, D85.D128 535, D85.D130 874, D85.D132 378, D85.D141 219,
D90.D104 290, D90.D105 533, D90.D107 338, D90.D112 110, D90.D125 562, D90.D128 270, D90.D130 558, D90.D132 118, D90.D141 376,
D104.D105 210, D104.D107 712, D104.D112 217, D104.D125 243, D104.D128 580, D104.D130 919, D104.D132 274, D104.D141 95,
D105.D107 823, D105.D112 418, D105.D125 37, D105.D128 767, D105.D130 1130, D105.D132 487, D105.D141 259,
D107.D112 453, D107.D125 860, D107.D128 485, D107.D130 181, D107.D132 528, D107.D141 742,
D112.D125 401, D112.D128 289, D112.D130 667, D112.D132 129, D112.D141 302,
D125.D128 753, D125.D130 1123, D125.D132 480, D125.D141 253,
D128.D130 755, D128.D132 238, D128.D141 559,
D130.D132 522, D130.D141 916,
D132.D141 310
/
dep_dist(i) distance from depot to distributor i
/D2 368, D23 409, D24 427, D25 477, D32 42, D33 1053, D34 480, D35 287, D37 382, D38 529, D39 85, D41 243, D42 201, D44 422, D48 652, D51 660, D53 201, D57 552, D60 90, D61 335, D68 302, D71 1006, D73 516, D74 483, D75 874, D76 263, D80 608, D85 303, D90 387, D104 81, D105 157, D107 803, D112 275, D125 200, D128 670, D130 1035, D132 325, D141 7/
demand(i) demand of distributor i (kg)
/D2 12030, D23 54924, D24 19998, D25 101397, D32 10699, D33 7768, D34 36105, D35 47644, D37 28726, D38 7982, D39 90830, D41 14139, D42 19242, D44 17492, D48 10589, D51 21791, D53 80824, D57 4233, D60 12312, D61 13734, D68 9921, D71 6248, D73 40506, D74 60838, D75 26154, D76 18834, D80 11360, D85 9075, D90 4399, D104 46467, D105 40376, D107 14029, D112 4920, D125 17631, D128 7300, D130 5971, D132 12071, D141 99297/
cost_truck(i) cost of using a truck to distributor i
/D2 6320, D23 7195, D24 7513, D25 7703, D32 1656, D33 15630, D34 7645, D35 5091, D37 6440, D38 8149, D39 2200, D41 3943, D42 3417, D44 7874, D48 9578, D51 12284, D53 4458, D57 9233, D60 1751, D61 5732, D68 4467, D71 16044, D73 7908, D74 12622, D75 13453, D76 4521, D80 9133, D85 4467, D90 14281, D104 1792, D105 3417, D107 11974, D112 6188, D125 5091, D128 10281, D130 16486, D132 6440, D141 984/
cost_tir(i) cost of using a tir to distributor i
/D2 8553, D23 9363, D24 9413, D25 9622, D32 1803, D33 18490, D34 8990, D35 6827, D37 8631, D38 10994, D39 2642, D41 5139, D42 4312, D44 9015, D48 11322, D51 16428, D53 5235, D57 12972, D60 3183, D61 6598, D68 5806, D71 20882, D73 9612, D74 17050, D75 20209, D76 8384, D80 10512, D85 5806, D90 19304, D104 2213, D105 4086, D107 16078, D112 7638, D125 6827, D128 12889, D130 20705, D132 8631, D141 1117/;
SCALAR
max_capacity_truck /77120/
min_capacity_truck /72300/
max_capacity_tir /129176/
min_capacity_tir /126766/
max_distributors_per_route /3/;
VARIABLES
z total cost
x(i,k,r) binary variable indicating if vehicle j is used for distributor i in route k
BINARY VARIABLES x;
EQUATIONS
obj objective function
route_constraints1(r) ensure a route has at most 3 distributors
route_constraints2(i) ensure a distributor is assigned to exactly one route
capacity_constraints1(r,k) ensure vehicle capacity constraints
capacity_constraints2(r,k);
obj ..
z =E= SUM((i,k,r), x(i,k,r) * (cost_truck(i)$(k EQ 'truck') + cost_tir(i)$(k EQ 'tir')));
route_constraints1(r) ..
SUM((i,k), x(i,k,r)) =L= max_distributors_per_route;
route_constraints2(i) ..
SUM((k,r), x(i,k,r)) =E= 1;
capacity_constraints1(r,'truck') ..
SUM(i, x(i,'truck',r) * demand(i)) =L= max_capacity_truck;
capacity_constraints2(r,'tir') ..
SUM(i, x(i,'tir',r) * demand(i)) =L= max_capacity_tir;
MODEL distributor_routing /all/;
SOLVE distributor_routing USING mip MINIMIZING z;
DISPLAY x.l, z.l;