Joint assignment, scheduling and routing models to Home Care optimization: a pattern based approach Online Supplement Paola Cappanera Maria Grazia Scutellà March 17, 2014 1 Details on symmetry management efficiency Tables 1 and 2, respectively for the maxmin and the maxmin with EC cases, report for the Heur, the ImplSol and two flow based pattern generation policies, the node of the branch and bound tree where the best integer solution is found with respect to the total number of nodes explored: as an example, 0/20 reports that 20 nodes have been overall explored within the time limit, and that the best solution has been found at root node (node number 0). These information allow to measure the efficiency of each approach: the bigger the number of explored nodes the faster the computing at each branch and bound node; likewise, the smaller the node where the best solution is found the bigger the capability of a truncated branch and bound approach, based on the ILP models proposed in this paper, to provide solutions of good quality. In particular, Tables 1 and 2 reveal that the symmetry management generally allows one to explore more nodes, and that the very good solutions found by policy FB-0.50 are computed quite early within the branch and bound tree. Table 1: Maxmin - a measure of efficiency FlowBased Heur ImplSol 0.50 0.75 January 2006 NoSymm 0/20 0/1 821/63253 n.a. Symm 837/872 922/924 1360/2382 n.a. April 2007 NoSymm 180/189 0/3 824/5633 16083/16083 Symm 390/905 n.a. 5846/16112 2746/5473 Dipartimento di Ingegneria dell Informazione, Università degli Studi di Firenze, Firenze, Italy - email: paola.cappanera@unifi.it Dipartimento di Informatica, Università di Pisa, Pisa, Italy - email: scut@di.unipi.it 1
Table 2: Maxmin with EC - a measure of efficiency FlowBased Heur ImplSol 0.50 0.75 January 2006 NoSymm n.a. 0/17 1012/1884 n.a. Symm 824/824 n.a. 1729/14850 n.a. April 2007 NoSymm n.a. n.a. 4/26 1545/1545 Symm n.a. n.a. 2129/4557 806/21447 Table 3: Impact of symmetry management on LP dimension Heur ImplSol FlowBased 0.50 LPRows LPCols LPRows LPCols LPRows LPCols January 2006 NoSymm 85042 712375 90304 743486 39813 304141 Symm 80517 669563 85158 696822 37618 285417 April 2007 NoSymm 139962 1236919 142368 1214617 71353 593663 Symm 123209 1061379 125538 1044417 63248 510963 2 Details on symmetry management and LP dimension Here we provide some comments and explanations in terms of dimension of the solved problems. To this end, Table 3 presents the dimension of the LP problems in terms of constraints (columns LPRows) and variables (columns LPCols) for the three policies Heur, ImplSol and FB-0.50 when the maxmin objective function is used. The data refer to the LP dimension reported by CPLEX after the preprocessing phase. Similar results hold for the alternative objective function, i.e. minmax. Results relative to the symmetry management are given in rows labeled Symm while the rows NoSymm report data relative to the base model. The symmetry management allows to reduce the LP dimension and the three policies exhibit a quite similar behavior in terms of reduction. In fact, for the selected week in January 2006 the average reduction in the number of constraints when the model is equipped with the symmetry management is equal to 5.51% while the percentage reduction of the variables is equal to 6.15. For April 2007 the figures are respectively 11.72% for the rows and 14.04% for the columns. Despite the reduction in terms of constraints and variables, the time required to solve the LP may increase when the symmetries are treated. However, as observed in the paper, the symmetry management is crucial almost everywhere either to halve the ILP gap or to find feasible solutions in the most critical cases. 3 Details on minmax efficiency Table 4, referred to minmax, gives for each of the three pattern generation policies Heur, ImplSol and FB-0.50 the value of the best integer solution (IPValue) and information about the total number of nodes explored in the branch and bound tree with the node where the best solution is found (NodeInf). 2
Table 4: minmax - efficiency results Heur ImplSol FB-0.50 IPValue NodeInf IPValue NodeInf IPValue NodeInf 0106-5-40-4-0 0.2442 907/1220 0.2454 815/972 0.2271 40297/41268 0106-5-40-4-1 0.2883 4/1575 0.2883 0/1600 0.2883 0/3029 0106-5-40-4-2 0.2883 104/1602 0.2883 616/1640 0.2883 376/1670 0106-8-40-4-0 0.2733 17654/17654 0.2756 9114/9154 0.2733 17128/17128 0106-8-40-4-1 0.4200 1593/1593 0.4200 2075/12662 0.4200 0/1606 0106-8-40-4-2 0.3033 253/1434 0.3017 794/794 0.3017 0/1641 0106-5-60-5-0 0.4060 0/137 0.3983 3/107 0.3983 1610/1704 0106-5-60-5-1 0.4450 388/391 0.5461 0/68 0.4350 1575/1575 0106-5-60-5-2 0.3233 0/88 0.2971 288/448 0.2767 13912/14550 0106-8-60-5-0 0.4433 4/55 0.4000 0/44 0.3539 0/1617 0106-8-60-5-1 0.4050 0/90 0.3672 1258/1562 0.3583 774/1798 0106-8-60-5-2 0.4583 431/606 0.4706 244/282 0.4294 1695/1695 0407-5-50-5-0 0.3394 311/625 0.3394 192/1090 0.3394 239/323 0407-5-50-5-1 0.3529 3/9 0.3561 3/11 0.3150 1099/1125 0407-5-50-5-2 0.4128 951/951 0.4128 795/795 0.4128 4/4 0407-8-50-5-0 0.4000 408/869 0.4000 761/761 0.4094 2042/2042 0407-8-50-5-1 0.3844 392/3237 0.3844 1010/1010 0.3894 798/800 0407-8-50-5-2 0.4867 3/4 0.2671 1144/1792 0407-5-80-6-0 n.a. n.a. n.a. n.a. 0.5611 801/801 0407-5-80-6-1 n.a. n.a. n.a. n.a. 0.5033 520/685 0407-5-80-6-2 n.a. n.a. n.a. n.a. 0.5244 1553/1553 0407-8-80-6-0 n.a. n.a. n.a. n.a. n.a. n.a. 0407-8-80-6-1 n.a. n.a. n.a. n.a. 0.7306 375/785 0407-8-80-6-2 n.a. n.a. n.a. n.a. 0.6650 1616/1616 4 Details on care continuity constraints Columns in Table 5 report, for the minmax case, respectively the LPTime, the IPTime, the optimality gap, the IPValue and the node where the best integer solution is found with respect to the total number of nodes explored. For each instance, results relative to the case where care continuity is assured and to the case where care continuity is relaxed are reported in consecutive rows. In particular, a suffix CC is appended to the instance name in the former case whereas suffix NoCC is used in the latter case. In general, we observe that relaxing the care continuity constraints has no impact on the value of the solutions except in few cases; on the other hand, disregarding the care continuity, the model seems to be faster in finding feasible solutions with respect to the case where the care continuity is ensured. In regards to the IPTime, we observe a non monotonic behavior since the relaxation of the care continuity constraints can involve either a reduction or an increase of computational time with respect to the more constrained model, at least on instances with up to 60 patients. Results are a little bit clearer on the most critical instances in the test bed, i.e. on instances characterized by 80 patients. In such a case we report either a valuable reduction of the IPTime or of the optimality gap; indeed, the elimination of care continuity constraints allows to find a feasible solution also for instance 0407-8-80-6-0 whereas no feasible solution is found in the other case. For some instances, the IPValue is smaller for the NoCC variant than in the CC counterpart. Similar results are obtained in the maxmin case. Furthermore, information on the quality of the solutions under maxmin and minmax for the restricted instances related to April 2007 are reported in Figures 1 and 2. There, 3
Table 5: minmax - impact of care continuity on efficiency and efficacy FB-0.50 LPTime IPTime %Gap IPValue NodeInf 0106-5-40-4-0-CC 19.92 43200 6.41 0.2271 40297/41268 0106-5-40-4-0-NoCC 31.38 43200 7.56 0.2300 17448/43899 0106-5-40-4-1-CC 8.82 279.49 0.53 0.2883 0/3029 0106-5-40-4-1-NoCC 62.39 445.46 0.24 0.2883 0/1629 0106-5-40-4-2-CC 6.07 380.96 0.00 0.2883 376/1670 0106-5-40-4-2-NoCC 3.96 153.64 0.54 0.2883 0/497 0106-8-40-4-0-CC 4.15 1946.90 0.00 0.2733 17128/17128 0106-8-40-4-0-NoCC 66.60 1352.38 0.00 0.2733 17342/17342 0106-8-40-4-1-CC 0.87 113.01 0.00 0.4200 0/1606 0106-8-40-4-1-NoCC 0.55 69.38 0.00 0.4200 806/806 0106-8-40-4-2-CC 1.06 295.91 0.00 0.3017 0/1641 0106-8-40-4-2-NoCC 0.91 390.46 0.00 0.3017 0/1608 0106-5-60-5-0-CC 23.07 7454.13 0.38 0.3983 1610/1704 0106-5-60-5-0-NoCC 5.71 2748.38 0.29 0.3983 1659/1677 0106-5-60-5-1-CC 6.83 3169.17 0.32 0.4350 1575/1575 0106-5-60-5-1-NoCC 5.62 2660.03 0.00 0.4350 789/789 0106-5-60-5-2-CC 89.63 16959.97 0.50 0.2767 13912/14550 0106-5-60-5-2-NoCC 20.69 6797.70 0.00 0.2767 17557/17557 0106-8-60-5-0-CC 4.75 1086.50 0.00 0.3539 0/1617 0106-8-60-5-0-NoCC 4.08 2528.13 0.00 0.3539 1590/1590 0106-8-60-5-1-CC 33.61 5470.26 0.23 0.3583 774/1798 0106-8-60-5-1-NoCC 6.15 5869.57 0.17 0.3583 1191/1645 0106-8-60-5-2-CC 4.14 3570.36 0.00 0.4294 1695/1695 0106-8-60-5-2-NoCC 4.18 1462.03 0.11 0.4294 815/1286 0407-5-50-5-0-CC 9.97 685.17 0.07 0.3394 239/323 0407-5-50-5-0-NoCC 13.72 1146.12 0.45 0.3394 270/444 0407-5-50-5-1-CC 184.77 43200 11.01 0.3150 1099/1125 0407-5-50-5-1-NoCC 138.33 43200 8.26 0.3054 826/1198 0407-5-50-5-2-CC 235.35 787.64 0.31 0.4128 4/4 0407-5-50-5-2-NoCC 311.56 1889.92 0.37 0.4128 0/180 0407-8-50-5-0-CC 108.61 4306.24 0.00 0.4094 2042/2042 0407-8-50-5-0-NoCC 238.74 1412.99 0.10 0.4094 627/2062 0407-8-50-5-1-CC 379.89 2013.19 0.00 0.3894 800/800 0407-8-50-5-1-NoCC 304.00 2874.76 0.09 0.3894 67/3224 0407-8-50-5-2-CC 112.94 43200 9.70 0.2671 1144/1792 0407-8-50-5-2-NoCC 167.86 43200 9.56 0.2667 1084/2688 0407-5-80-6-0-CC 235.88 39240.70 0.00 0.5611 801/801 0407-5-80-6-0-NoCC 59.48 34469.97 0.00 0.5611 735/1555 0407-5-80-6-1-CC 164.30 43200 18.78 0.5033 520/685 0407-5-80-6-1-NoCC 4875.17 43200 9.45 0.4508 695/1123 0407-5-80-6-2-CC 91.10 32900.24 0.00 0.5244 1551/1553 0407-5-80-6-2-NoCC 56.39 23941.20 0.00 0.5244 773/775 0407-8-80-6-0-CC 38.02 43200 n.a. n.a. n.a. 0407-8-80-6-0-NoCC 53.65 33082.74 0.00 0.6417 561/1560 0407-8-80-6-1-CC 100.56 43200 6.62 0.7306 375/785 0407-8-80-6-1-NoCC 66.34 20276.95 0.00 0.6878 462/1609 0407-8-80-6-2-CC 65.11 43200 1.02 0.6650 1616/1616 0407-8-80-6-2-NoCC 48.06 43200 1.26 0.6606 458/556 4
as defined throughout the paper, UF refers to the operator utilization factor, while T T F gives the fraction of travel time with respect to the maximum weekly workload. Here and in the following figure as well, for each instance, information relative to the case where care continuity is assured and to the case where care continuity is relaxed are reported consecutively so as to facilitate their comparison. Finally, see Figures 3, 4, 5 and 6 for an analysis of the number of operators per patient in both scenarios (maxmin and minmax), under the care continuity constraints and their removal. In these figures the percentage of patients that have been assigned to one, two, three, four and five operators is given (in the legend Op stands for Operator ). We observe that the number of patients assigned to more than two operators is quite low when the care continuity constraints are ignored, for both the alternative objective functions. For restricted instances referring to January 2006 this behavior is more evident than for instances related to April 2007. Metrics 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 QoS Ind min UF max UF avg UF min TTF max TTF avg TTF 0.0 5-50-5-0 5-50-5-1 5-50-5-2 8-50-5-0 8-50-5-1 8-50-5-2 5-80-6-0 Instances 5-80-6-1 5-80-6-2 8-80-6-0 8-80-6-1 8-80-6-2 Figure 1: maxmin - solution quality analysis on 0407 instances 5
Metrics 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 QoS Ind min UF max UF avg UF min TTF max TTF avg TTF 5-50-5-0 5-50-5-1 5-50-5-2 8-50-5-0 8-50-5-1 8-50-5-2 5-80-6-0 Instances 5-80-6-1 5-80-6-2 8-80-6-0 8-80-6-1 8-80-6-2 Figure 2: minmax - solution quality analysis on 0407 instances 100% 90% 80% 70% 60% 5 Op 4 Op 3 Op 2 Op 1 Op 50% 40% 30% 20% 10% 0% Figure 3: maxmin - frequency analysis on 0106 instances 6
100% 90% 80% 70% 60% Op 5 Op 4 Op 3 Op 2 Op 1 50% 40% 30% 20% 10% 0% Figure 4: maxmin - frequency analysis on 0407 instances 100% 90% 80% 70% 60% Op 5 Op 4 Op 3 Op 2 Op 1 50% 40% 30% 20% 10% 0% Figure 5: minmax - frequency analysis on 0106 instances 7
100% 90% 80% 70% 60% Op 5 Op 4 Op 3 Op 2 Op 1 50% 40% 30% 20% 10% 0% Figure 6: minmax - frequency analysis on 0407 instances 8