Пример #1
0
        /// <inheritdoc />
        public OptimizationResult Solve(long timeLimitMilliseconds, EventHandler <ProgressReport> progress, EventHandler <string> consoleProgress)
        {
            var sw     = Stopwatch.StartNew();
            var result = new OptimizationResult
            {
                OptimizationInput = input
            };

            var numberOfFakeSantas = config.MaxNumberOfAdditionalSantas;
            var numberOfRoutes     = (input.Santas.Length + numberOfFakeSantas) * input.Days.Length;
            var numberOfDays       = input.Days.Length;
            var visits             = input.Visits.ToList();
            var breakDictionary    = new Dictionary <(int day, int santa), int>();

            for (var i = 1; i < numberOfDays; i++)
            {
                foreach (var breakVisit in input.Visits.Where(v => v.IsBreak))
                {
                    visits.Add(breakVisit);
                    breakDictionary.Add((day: i, santa: breakVisit.SantaId), visits.Count - 1);
                }
            }

            var vrptwBreaksOnWayTimeLimitFactor = 1 - config.VrpTimeLimitFactor - config.VrptwTimeLimitFactor;

            if (Math.Abs(config.VrpTimeLimitFactor + config.VrptwTimeLimitFactor + vrptwBreaksOnWayTimeLimitFactor - 1d) > 0.00001)
            {
                throw new Exception("programmer was stupid :). Time limit factors don't add up to 1.");
            }

            using (var localSolver = new localsolver.LocalSolver())
            {
                var model = localSolver.GetModel();

                var solverVariables = new SolverVariables(model, numberOfRoutes, visits, input.RouteCosts, input, numberOfFakeSantas);
                var modelBuilder    = new ModelBuilder(solverVariables);

                consoleProgress?.Invoke(this, "Starting to model");

                #region VRP

                modelBuilder.AddPartitionConstraint();
                modelBuilder.AddVisitsNotInLastSequenceConstraint();

                for (var day = 0; day < numberOfDays; day++)
                {
                    for (var santa = 0; santa < input.Santas.Length + numberOfFakeSantas; santa++)
                    {
                        modelBuilder.SetSantaUsed(day, santa);

                        // break
                        var breaks = visits.Where(v => v.IsBreak && (v.SantaId == santa)).ToList();
                        if (breaks.Any())
                        {
                            var breakIndex = day == 0 ? visits.IndexOf(breaks.First()) : breakDictionary[(day, santa)];
Пример #2
0
        public static void Execute()
        {
            const int numberOfSantas = 2;
            const int numberOfVisits = 6;
            const int m            = int.MaxValue;
            var       distanceHome = new int[numberOfVisits]
            {
                0, 0, 0, 0, 0, 0
            };
            var distance = new int[numberOfVisits][]
            {
                new int[numberOfVisits] {
                    m, 3, 6, 2, 8, 1
                },
                new int[numberOfVisits] {
                    4, m, 3, 4, 4, 5
                },
                new int[numberOfVisits] {
                    3, 2, m, 6, 3, 5
                },
                new int[numberOfVisits] {
                    4, 2, 5, m, 4, 4
                },
                new int[numberOfVisits] {
                    3, 3, 2, 6, m, 4
                },
                new int[numberOfVisits] {
                    7, 4, 5, 7, 6, m
                },
            };

            using (var localsolver = new localsolver.LocalSolver())
            {
                // Declares the optimization model.
                var model = localsolver.GetModel();

                var santaUsed      = new LSExpression[numberOfSantas];
                var visitSequences = new LSExpression[numberOfSantas];
                var routeDistances = new LSExpression[numberOfSantas];


                // Sequence of customers visited by each truck.
                for (int k = 0; k < numberOfSantas; k++)
                {
                    visitSequences[k] = model.List(numberOfVisits);
                }

                model.Constraint(model.Partition(visitSequences));

                // Create demands and distances as arrays to be able to access it with an "at" operator


                var distanceArray     = model.Array(distance);
                var distanceHomeArray = model.Array(distanceHome);

                for (int s = 0; s < numberOfSantas; s++)
                {
                    var sequence = visitSequences[s];
                    var c        = model.Count(sequence);

                    santaUsed[s] = c > 0;


                    var distSelector = model.Function(i => distanceArray[sequence[i - 1], sequence[i]]);
                    routeDistances[s] = model.Sum(model.Range(1, c), distSelector)
                                        + model.If(c > 0, distanceHomeArray[sequence[0]] + distanceHomeArray[sequence[c - 1]], 0);
                }


                var totalDistance = model.Sum(routeDistances);

                // Objective: minimize the number of trucks used, then minimize the distance traveled
                //model.Minimize(nbTrucksUsed);
                model.Minimize(totalDistance);

                model.Close();

                // Parameterizes the solver.
                var phase = localsolver.CreatePhase();
                phase.SetTimeLimit(3);

                localsolver.Solve();

                // output
                for (int i = 0; i < numberOfSantas; i++)
                {
                    System.Console.WriteLine($"santa {i+1}: ");
                    System.Console.WriteLine(string.Join("->", visitSequences[i].GetCollectionValue()));
                }
            }
        }
Пример #3
0
        public static void Execute()
        {
            const int numberOfSantas = 6;
            const int numberOfVisits = 30;
            const int m            = int.MaxValue;
            var       distanceHome = new int[numberOfVisits]
            {
                106, 378, 305, 337, 483, 337, 316, 142, 306, 242, 130, 0, 381, 358, 255, 167, 34, 345, 972, 1116, 1160, 1113, 467, 534, 1073, 956, 937, 1020, 254, 389
            };
            var distance = new int[numberOfVisits][]
            {
                new int[numberOfVisits] {
                    0, 273, 274, 366, 383, 429, 408, 234, 398, 334, 222, 92, 302, 253, 150, 156, 63, 240, 1064, 1208, 1252, 1205, 496, 563, 1165, 1048, 1029, 1112, 264, 418
                },
                new int[numberOfVisits] {
                    246, 0, 203, 219, 181, 675, 654, 480, 644, 580, 468, 338, 100, 428, 107, 331, 309, 120, 1310, 1454, 1498, 1451, 416, 483, 1411, 1294, 1275, 1358, 439, 504
                },
                new int[numberOfVisits] {
                    268, 215, 0, 180, 207, 621, 600, 426, 590, 526, 340, 284, 105, 450, 184, 353, 255, 264, 1256, 1400, 1444, 1397, 377, 444, 1357, 1240, 1221, 1304, 461, 348
                },
                new int[numberOfVisits] {
                    414, 275, 224, 0, 267, 709, 688, 514, 678, 614, 426, 372, 165, 632, 311, 539, 343, 324, 1344, 1488, 1532, 1485, 197, 264, 1445, 1328, 1309, 1392, 626, 311
                },
                new int[numberOfVisits] {
                    345, 170, 184, 200, 0, 774, 753, 579, 743, 679, 495, 437, 81, 527, 206, 430, 408, 219, 1409, 1553, 1597, 1550, 293, 360, 1510, 1393, 1374, 1457, 538, 485
                },
                new int[numberOfVisits] {
                    480, 753, 680, 712, 858, 0, 186, 236, 176, 355, 443, 375, 756, 733, 630, 525, 409, 720, 1085, 1229, 977, 1008, 842, 909, 1094, 1069, 1050, 1041, 612, 743
                },
                new int[numberOfVisits] {
                    441, 714, 641, 673, 819, 161, 0, 224, 8, 262, 404, 336, 717, 694, 591, 486, 370, 681, 992, 1136, 1069, 1100, 803, 870, 1093, 976, 957, 1040, 573, 704
                },
                new int[numberOfVisits] {
                    244, 517, 444, 476, 622, 195, 208, 0, 198, 119, 207, 139, 520, 497, 394, 289, 173, 484, 849, 993, 1037, 990, 606, 673, 950, 833, 814, 897, 376, 507
                },
                new int[numberOfVisits] {
                    433, 706, 633, 665, 811, 153, 11, 216, 0, 254, 396, 328, 709, 686, 583, 478, 362, 673, 984, 1128, 1061, 1092, 795, 862, 1085, 968, 949, 1032, 565, 696
                },
                new int[numberOfVisits] {
                    329, 602, 529, 561, 707, 299, 224, 104, 214, 0, 292, 224, 605, 582, 479, 374, 258, 569, 746, 874, 918, 871, 691, 758, 831, 714, 695, 778, 338, 592
                },
                new int[numberOfVisits] {
                    253, 526, 379, 409, 557, 423, 402, 228, 392, 328, 0, 148, 455, 506, 403, 298, 182, 493, 1058, 1202, 1246, 1199, 539, 606, 1159, 1042, 1023, 1106, 385, 393
                },
                new int[numberOfVisits] {
                    106, 378, 305, 337, 483, 337, 316, 142, 306, 242, 130, 0, 381, 358, 255, 167, 34, 345, 972, 1116, 1160, 1113, 467, 534, 1073, 956, 937, 1020, 254, 389
                },                                                                                                                                                                                                                                                            // depot
                new int[numberOfVisits] {
                    285, 110, 103, 119, 102, 695, 674, 500, 664, 600, 414, 358, 0, 467, 146, 370, 329, 159, 1330, 1474, 1518, 1471, 316, 383, 1431, 1314, 1295, 1378, 478, 404
                },
                new int[numberOfVisits] {
                    203, 405, 406, 569, 515, 632, 611, 437, 601, 537, 425, 295, 434, 0, 282, 207, 266, 372, 1267, 1411, 1455, 1408, 699, 766, 1368, 1251, 1232, 1315, 315, 621
                },
                new int[numberOfVisits] {
                    139, 123, 179, 271, 233, 568, 547, 373, 537, 473, 361, 231, 152, 321, 0, 224, 202, 90, 1203, 1347, 1391, 1344, 468, 535, 1304, 1187, 1168, 1251, 332, 421
                },
                new int[numberOfVisits] {
                    128, 330, 331, 465, 440, 448, 427, 253, 417, 353, 241, 128, 359, 229, 207, 0, 162, 297, 1083, 1227, 1271, 1224, 595, 662, 1184, 1067, 1048, 1131, 108, 517
                },
                new int[numberOfVisits] {
                    70, 343, 266, 298, 444, 365, 344, 170, 334, 270, 158, 28, 342, 323, 220, 195, 0, 310, 954, 1098, 1008, 1039, 711, 778, 1055, 938, 919, 1002, 481, 612
                },
                new int[numberOfVisits] {
                    217, 124, 256, 272, 234, 646, 625, 451, 615, 551, 439, 309, 153, 399, 78, 302, 280, 0, 1281, 1425, 1469, 1422, 469, 536, 1382, 1265, 1246, 1329, 410, 499
                },
                new int[numberOfVisits] {
                    1114, 1387, 1314, 1346, 1522, 1084, 1009, 889, 999, 794, 1077, 1009, 1420, 1367, 1264, 1159, 970, 1354, 0, 245, 316, 269, 1476, 1543, 229, 112, 276, 176, 1052, 1377
                },
                new int[numberOfVisits] {
                    1280, 1553, 1480, 1512, 1688, 1250, 1175, 1055, 1165, 951, 1243, 1175, 1586, 1533, 1430, 1325, 1136, 1520, 259, 0, 353, 306, 1642, 1709, 139, 252, 315, 184, 1289, 1543
                },
                new int[numberOfVisits] {
                    1355, 1628, 1555, 1587, 1763, 1034, 1145, 1130, 1135, 1026, 1318, 1250, 1661, 1608, 1505, 1400, 1082, 1595, 361, 391, 0, 149, 1532, 1599, 235, 250, 282, 182, 1364, 1540
                },
                new int[numberOfVisits] {
                    1311, 1584, 1511, 1543, 1719, 1068, 1179, 1086, 1169, 982, 1274, 1206, 1617, 1564, 1461, 1356, 1116, 1551, 317, 347, 152, 0, 1566, 1633, 191, 206, 238, 138, 1320, 1574
                },
                new int[numberOfVisits] {
                    577, 505, 454, 230, 387, 872, 851, 677, 841, 777, 589, 535, 395, 830, 541, 702, 758, 554, 1507, 1651, 1508, 1539, 0, 67, 1608, 1491, 1472, 1555, 789, 474
                },
                new int[numberOfVisits] {
                    646, 574, 523, 299, 456, 941, 920, 746, 910, 846, 658, 604, 464, 899, 610, 771, 827, 623, 1576, 1720, 1577, 1608, 69, 0, 1677, 1560, 1541, 1624, 858, 543
                },
                new int[numberOfVisits] {
                    1247, 1520, 1447, 1479, 1655, 1130, 1142, 1022, 1132, 918, 1210, 1142, 1553, 1500, 1397, 1292, 1103, 1487, 253, 156, 214, 167, 1609, 1676, 0, 142, 176, 45, 1256, 1510
                },
                new int[numberOfVisits] {
                    1125, 1398, 1325, 1357, 1533, 1095, 1020, 900, 1010, 796, 1088, 1020, 1431, 1378, 1275, 1170, 981, 1365, 131, 257, 224, 177, 1487, 1554, 137, 0, 184, 84, 1134, 1388
                },
                new int[numberOfVisits] {
                    1109, 1382, 1309, 1341, 1517, 1079, 1004, 884, 994, 780, 1072, 1004, 1415, 1362, 1259, 1154, 965, 1349, 290, 322, 251, 204, 1471, 1538, 166, 179, 0, 113, 1118, 1372
                },
                new int[numberOfVisits] {
                    1202, 1475, 1402, 1434, 1610, 1085, 1097, 977, 1087, 873, 1165, 1097, 1508, 1455, 1352, 1247, 1058, 1442, 208, 209, 169, 122, 1564, 1631, 54, 97, 131, 0, 1211, 1465
                },
                new int[numberOfVisits] {
                    241, 443, 444, 557, 553, 540, 519, 345, 509, 337, 333, 220, 472, 342, 320, 113, 426, 410, 996, 1211, 1255, 1208, 687, 754, 1168, 1051, 1032, 1115, 0, 609
                },
                new int[numberOfVisits] {
                    505, 606, 443, 350, 598, 775, 754, 580, 744, 680, 445, 463, 496, 758, 521, 630, 661, 611, 1410, 1554, 1525, 1551, 480, 547, 1511, 1394, 1375, 1458, 717, 0
                },
            };

            var visitDuration = new int[numberOfVisits]
            {
                2700,
                1500,
                2100,
                3600,
                1500,
                1200,
                1800,
                1500,
                1500,
                1800,
                2100,
                0,
                1800,
                1200,
                1500,
                1500,
                1800,
                1500,
                1500,
                1500,
                1500,
                1200,
                1200,
                1500,
                1800,
                1500,
                1500,
                1500,
                2400,
                1500
            };

            using (var localsolver = new localsolver.LocalSolver())
            {
                // Declares the optimization model.
                var model = localsolver.GetModel();

                var santaUsed           = new LSExpression[numberOfSantas];
                var visitSequences      = new LSExpression[numberOfSantas];
                var routeDistances      = new LSExpression[numberOfSantas];
                var santaVisitDurations = new LSExpression[numberOfSantas];


                // Sequence of customers visited by each truck.
                for (int k = 0; k < numberOfSantas; k++)
                {
                    visitSequences[k] = model.List(numberOfVisits);
                }

                model.Constraint(model.Partition(visitSequences));


                var distanceArray      = model.Array(distance);
                var distanceHomeArray  = model.Array(distanceHome);
                var visitDurationArray = model.Array(visitDuration);

                for (int s = 0; s < numberOfSantas; s++)
                {
                    var sequence = visitSequences[s];
                    var c        = model.Count(sequence);

                    santaUsed[s] = c > 0;

                    var distSelector  = model.Function(i => distanceArray[sequence[i - 1], sequence[i]] + visitDurationArray[i]);
                    var visitSelector = model.Function(i => visitDurationArray[i]);
                    routeDistances[s] = model.Sum(model.Range(1, c), distSelector)
                                        + model.If(santaUsed[s], distanceHomeArray[sequence[0]] + distanceHomeArray[sequence[c - 1]], 0);
                    santaVisitDurations[s] = model.Sum(model.Range(1, c), visitSelector);

                    model.Constraint(routeDistances[s] <= (21 - 17) * 60 * 60); // working hours from 17:00 to 21:00
                }


                var totalDistance = model.Sum(routeDistances);


                // model.Minimize(totalDistance);

                // every santa same visit duration:

                var totalDuration      = visitDuration.Sum();
                var numberOfSantasUsed = model.Sum(santaUsed);
                var avgDists           = new LSExpression[numberOfSantas];
                for (int s = 0; s < numberOfSantas; s++)
                {
                    avgDists[s] = model.Abs(santaVisitDurations[s] -
                                            (totalDuration / numberOfSantasUsed));
                }

                var avgTime = model.Sum(avgDists);
                model.Minimize(avgTime * 5 + totalDistance);

                model.Close();

                // Parameterizes the solver.
                var phase = localsolver.CreatePhase();
                phase.SetTimeLimit(500);

                localsolver.Solve();

                // output
                for (int i = 0; i < numberOfSantas; i++)
                {
                    System.Console.WriteLine($"santa {i + 1}: ");
                    System.Console.WriteLine(string.Join("->", visitSequences[i].GetCollectionValue()));
                    System.Console.WriteLine($"route distance: {routeDistances[i].GetValue()}");
                }
            }
        }