Next Article in Journal
Machine Learning Model for Multiomics Biomarkers Identification for Menopause Status in Breast Cancer
Next Article in Special Issue
Ensemble Heuristic–Metaheuristic Feature Fusion Learning for Heart Disease Diagnosis Using Tabular Data
Previous Article in Journal
The Archimedean Origin of Modern Positional Number Systems
Previous Article in Special Issue
Enhancing Cryptocurrency Price Forecasting by Integrating Machine Learning with Social Media and Market Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compact Models to Solve the Precedence-Constrained Minimum-Cost Arborescence Problem with Waiting Times

by
Mauro Dell’Amico
1,2,
Jafar Jamal
1 and
Roberto Montemanni
1,2,*
1
Department of Sciences and Methods for Engineering, University of Modena and Reggio Emilia, Via Amendola 2, 42122 Reggio Emilia, RE, Italy
2
Interdepartmental Center En&Tech, University of Modena and Reggio Emilia, Capannone 19 Tecnopolo, Piazza Europa 1, 42122 Reggio Emilia, RE, Italy
*
Author to whom correspondence should be addressed.
Algorithms 2024, 17(1), 12; https://doi.org/10.3390/a17010012
Submission received: 2 December 2023 / Revised: 23 December 2023 / Accepted: 26 December 2023 / Published: 27 December 2023
(This article belongs to the Special Issue 2022 and 2023 Selected Papers from Algorithms Editorial Board Members)

Abstract

:
The minimum-cost arborescence problem is a well-studied problem. Polynomial-time algorithms for solving it exist. Recently, a new variation of the problem called the Precedence-Constrained Minimum-Cost Arborescence Problem with Waiting Times was presented and proven to be N P -hard. In this work, we propose new polynomial-size models for the problem that are considerably smaller in size compared to those previously proposed. We experimentally evaluate and compare each new model in terms of computation time and quality of the solutions. Several improvements to the best-known upper and lower bounds of optimal solution costs emerge from the study.

1. Introduction

The Minimum-Cost Arborescence (MCA) problem asks to find a directed minimum-cost spanning tree (i.e., an arborescence) rooted at vertex r—the root of the arborescence—in a directed graph. Chu and Liu [1] and Edmonds [2] proposed the same polynomial time algorithm to solve the problem, independently from each other. A faster implementation was later proposed by Gabow and Tarjan [3], and a different polynomial time algorithm [4] that operates directly on the cost matrix was introduced by Bock [5].
Several variations of the MCA problem were introduced in the literature since its introduction. Given a directed graph with a resource associated to each of its vertices, the Resource-Constrained Minimum-Weight Arborescence problem aims at retrieving an arborescence with minimum total cost, with the additional constraint that outgoing arcs from each vertex have to have a cost at most equal to that of the resource of the vertex itself (Fischetti and Vigo [6]). Given a weighted directed graph G = ( V , A ) with a vertex r V identified as the root and an integer p, the p-Arborescence Star problem asks to identify a minimum-cost arborescence rooted at r. The arborescence spans the set of vertices H V { r } of size p, and there must be an assignment between each vertex v V { H r } and one of the vertices in H (Pereira et al. [7], Morais et al. [8], Hakimi [9]). Given a directed graph with a number assigned to each vertex, the Restricted Fathers Tree problem seeks a minimum-cost arborescence, with the constraint that each path between a vertex and the root has to touch vertices with ranking not lower than the vertex (Guttmann-Beck and Hassin [10]). Given a directed graph with a vertex r designed as the root and a set A v V for every v V { r } such that r A v , the Restricted Ancestors Tree problem aims at finding a minimum-cost arborescence rooted at r, with the additional constraint that vertex i can be an ancestor of vertex j only when i A j (Guttmann-Beck and Hassin [10]). The Minimum Spanning Tree Problem with Conflict Pairs (Carrabs and Gaudioso [11]) is a variation of the classic Minimum-Spanning Tree problem (Kruskal [12]), characterized by an undirected graph and a set S containing conflicting pairs of conflicting edges. The objective is to retrieve a minimum-cost spanning tree with at most one edge from the pair in S. The Capacitated Minimum Spanning Tree problem, introduced in Gouveia and Lopez [13], is another variation, characterized by non-negative integer node demands q j for each node j V { r } and a budget Q for the sum of the weights in any root–leaf path. Further related problems can be found in Frieze and Tkocz [14], Fertin et al. [15], Eswaran and Tarjan [16], Li et al. [17], Kawatra and Bricker [18], Galbiati et al. [19], Bérczi et al. [20], Bang-Jensen [21], Yingshu et al. [22], Carrabs et al. [23], Darmann et al. [24] and Viana and Campêlo [25]. The Sequential Ordering Problem, introduced in Escudero [26], is relevant to the present study and can be described as follows. Given a weighted graph and a set of precedence constraints between vertex pairs and start and end vertices, the goal is to find a Minimum-Cost Hamiltonian Path that respects the precedence constraints. Solving algorithms can be found in Moon et al. [27], Balas et al. [28], Hernádvölgyi [29], Escudero et al. [30], Gambardella and Dorigo [31], Karan and Skorin-Kapov [32], Ascheuer et al. [33], Ascheuer et al. [34], Montemanni et al. [35] and Fiala Timlin and Pulleyblank [36]. The Precedence-Constrained Minimum-Cost Arborescence is an extension of the MCA problem first introduced in Dell’Amico et al. [37]. Precedence constraints have the following meaning. A precedence set R containing pairs of vertices is given. For each ( s , t ) R , if both vertices s and t are on a same path of the arborescence, then vertex s has to be visited before vertex t. The optimization seeks to find an arborescence of minimum total cost such that all the precedence constraints are satisfied. Several models for the problem, all based on Mixed Integer Linear Programming (MILP), were proposed. We refer the interested reader to Chou et al. [38].
The Precedence-Constrained Minimum-Cost Arborescence problem with Waiting Times (PCMCA-WT)—which is the object of the current paper—was first introduced in Chou et al. [38], where the problem was shown to be NP -hard through a reduction to the Rectilinear Steiner Arborescence problem (Shi and Su [39]), and different MILP models were proposed. The problem is about retrieving an arborescence where traveling times are present among vertices and temporal precedences relative to the time of visit have to be fulfilled among pairs of vertices. Waiting times at vertices are allowed to enforce such precedences, but these waiting times are accounted for in the objective function together with travel times. The optimization is to minimize such an objective function.
The organization of the paper is as follows. The PCMCA-WT is formally defined in Section 2. Section 3 describes a new family of compact models for the problem (characterized by a polynomial number of variables and constraints). Section 4 discusses the results of a vast experimental campaign where the new models are compared to those previously disclosed in the literature. Some conclusions are the content of Section 5.
The contributions of the paper can be summarized as follows:
  • New models for the PCMCA-WT that are polynomial in size and are characterized by a substantially smaller memory footprint compared to the known models are introduced. This result is achieved by exploiting some theoretical properties emerging from the current study and previously unobserved.
  • The new models are solved both with MILP and Constraint Programming (CP) solvers. The experimental results substantially improve the state of the art for the instances commonly adopted in the literature. Out of the 88 open instances from the literature, improved lower bounds are provided for 71 instances and improved upper bounds are provided for 80 instances. Finally, seven instances are closed for the first time.

2. Problem Description

The PCMCA-WT can be described according to the following definitions. A directed graph G = ( V , A , R ) is given, with  V = { 1 , , n } being a set of vertices and  A V × V a set of arcs, with a non-negative cost c i j associated with every arc ( i , j ) A . It represents the traversing time for that arc. The set R V × V contains precedence relationships. Let d j be the time step at which the flow enters vertex j V , with  d r = 0 . For any ( s , t ) R , we impose d t d s . This implies that the flow cannot enter vertex t before entering vertex s, but it can wait at any vertex before servicing it. We define w j as the waiting time at vertex j V , with  w r = 0 . The waiting time at a vertex j that is visited after another vertex i in the current solution is defined as w j = d j ( d i + c i j ) . Given a vertex r V being the root of the arborescence, the target of the optimization is to retrieve an arborescence T (with root r) that provides the lowest possible sum of the total cost plus the total waiting time.
An example of PCMCA-WT is provided in Figure 1. In the instance (top part of the figure), the precedence relationship ( 1 , 3 ) R is represented as a dashed arrow, while in the bottom part of the figure, an optimal solution for the given instance is shown. The corresponding values of d i and w i are exposed near each vertex. The cost of the solution is 8 (obtained by adding the cost of the traversed arcs and the waiting times payed at vertices). In the example, observe the waiting time of 1 unit at vertex 3 ( d 1 = 4 , d 3 = 3 and ( 1 , 3 ) R ).

3. New Compact Models

An MILP model for the PCMCA-WT that is polynomial in size was recently proposed in Chou et al. [38]. The model is based on a multicommodity flow formulation [40] that extends the flow conservation constraints in order to satisfy the precedence relationships between vertex pairs. However, the model suffers from computational limitations caused by the large number of variables and constraints, both in the order of O ( n 3 ) . In this section, we derive two new models for the PCMCA-WT that are also polynomial in size. Compared to the previous models, the new ones use less variables and constraints for the description of the precedence relationships among pairs of vertices. This solves the memory issues that were encountered when using the multicommodity flow model originally introduced in [38], making compact models competitive against other solutions.

3.1. The Complete Model

We first define V R as the set of vertices of V involved in at least one precedence relation as a head. Formally, V R = { t V | s V : ( s , t ) R } . Let x i j be a binary variable modeling if arc ( i , j ) A is visited: x i j = 1 if ( i , j ) T , and 0 otherwise. Let y i be an integer variable used to identify the position of vertex i V along the only path of the arborescence connecting the root r to vertex i itself. Let u j t be a binary variable with indices representing vertex j V and vertex t V R . This variable will be used to model precedences. Let d j be a continuous variable containing entering time of the flow at vertex j V . Finally, let w j be a continuous variable modeling the time waited at vertex j before letting the flow enter the node itself.
An MILP model for the PCMCA-WT is as follows.
minimize ( i , j ) A c i j x i j + i V w i
subject to : ( i , j ) A x i j = 1 j V { r }
y i y j + 1 n ( 1 x i j ) ( i , j ) A : j r
u s t = 0 ( s , t ) R
u t t = 1 t V R
u j t u i t x i j 1 t V R , ( i , j ) A
d r = 0
w r = 0
d j d i M + ( M + c i j ) x i j ( i , j ) A
w j d j d i M + ( M c i j ) x i j ( i , j ) A
d t d s ( s , t ) R
x i j { 0 , 1 } ( i , j ) A
y i { 0 , 1 , , n 1 } i V
u j t { 0 , 1 } t V R , j V
d i 0 i V
w i 0 i V
The objective function (1) minimizes the sum of the total travel and waiting times, as described in Section 2. The set of constraints (2) enforces that each vertex apart from the root needs to have one incoming arc. Constraints (3) model subtour elimination and dictate that any feasible solution cannot contain any cycles. The set of constraints (2) and (3) work together in order to enforce only solutions in the form of an arborescence rooted at vertex r. Constraints (4), (5) and (6) regulate precedence constraints among the vertices visited in a same branch of the tree. The logic behind these constraints will be explained in detail in Section 3.1.1, entirely devoted to this purpose. Constraints (7) and (8) initialize the distance traveled and the waiting time at the root r to 0. Constraints (9), activated once an arc ( i , j ) A is selected, force the arrival time at vertex j to be not lower than the arrival time vertex i plus the travel time c i j . Note that here and in the following set of constraints, M is an arbitrarily large constant. Constraints (10) push the waiting time at vertex j to be no smaller than the the service time at vertex j minus the service time at vertex i plus c i j . Constraints (11) impose that the service time at vertex t cannot be smaller than the service time at vertex s for all ( s , t ) R . This set of constraints, in conjunction with the previous ones, regulates the value of the waiting times. Finally, constraints (12)–(16) define the domain of the variables.
The value of the large constant M, appearing in constraints (9) and (10), is an approximation by excess of the optimal cost of the problem. In our case, the solution cost of solving the instance as a Sequential Ordering Problem [34] using a nearest neighbor algorithm [41] is taken as the value of M. This is a valid upper bound for the optimal cost of a PCMCA-WT instance, being a valid solution for the Sequential Ordering Problem a simple directed path that includes all the vertices of the graph, with  the constraint that t never precede s for all ( s , t ) R . This implies that d t d s for all ( s , t ) R .
The model proposed in this section has a considerably smaller memory footprint compared to the polynomial-size model proposed for the PCMCA-WT in [38]. In detail, the number of variables is reduced from O ( n 3 ) to O ( n 2 ) . The number of constraints remains instead in the order of O ( n 3 ) , although now the hidden multiplicative factor depends on the number of vertices involved in at least one precedence as a head, instead of all the vertices. This makes the number of constraints much smaller. All together, these improvements reduce substantially the memory footprint of the model, with great advantages for practical tractability.

3.1.1. The New Precedence Constraints

The approach proposed in this work to deal with precedences is similar to the idea originally introduced in Dell’Amico et al. [42] for the Precedence-Constrained Minimum-Cost Arborescence (PCMCA) problem. Constraints (4) and (5) impose the values of u s t and u t t to be 0 and 1, respectively, for all ( s , t ) R , and  t V : ( s , t ) R . On the other hand, the set of constraints (6) enforces that u j t u i t whenever arc ( i , j ) A is selected to be part of the solution. These concepts together forbid any violation of precedence constraints along a same path of the solution arborescence T.
Figure 2 shows an example of how a precedence-violating path is detected using the set of constraints (6). In the figure, the range/value of variable u j t is written on the left of each vertex, and black arcs show the arcs that are part of the solution, while the red arcs show a precedence relationship ( s , t ) R . In the figure, constraints (6) enforce that u 1 t and u 2 t have to be greater than or equal to 1. However, once u s t 1 is imposed through this logic, the constraint (4) relative to variable u s t is violated, rendering therefore the solution infeasible.

3.2. The Reduced Model

It can be observed that by removing the set of constraints (4)–(6) from the model described in Section 3.1, the set of constraints (11) in general still enforces the precedence relationships between s and t with ( s , t ) R , apart from the following special case. A directed path that visits t before visiting s implies that d t d s , which violates the set of constraints (11). However, the set of constraints (11) might fail to enforce a precedence relationship between s and t if a zero-cost path in G that reaches s from t exists. Therefore, variables u j t and constraints (4)–(6) need to be defined only for those t for which there exist at least an s for which ( s , t ) R and a zero-cost path that connects t to s is available in G. The identification of such pairs ( s , t ) R for which a zero-cost path exists has, however, to be performed in a preprocessing phase, adding some complexity to the overall process. The procedure we devised for the retrieval of such pairs will be detailed in Section 3.2.1.
The reduced set of constraints can be described formally as follows. Let G 0 = ( V 0 , A 0 , R ) be the graph obtained from G by considering only arcs with cost zero and the relevant nodes. Therefore, A 0 = { ( i , j ) A | c i j = 0 } and V 0 = { i V | j V : ( i , j ) A 0 or ( j , i ) A 0 } (a node is considered relevant for the residual precedence constraints if it has at least an outgoing or incoming zero-cost arc). Let S P i j A 0 be a shortest path that starts from i to j in G, and let c ( S P i j ) = ( k , l ) S P i j c k l be its cost. For each t V R , let V 0 ( t ) = { s V 0 | ( s , t ) R with c ( P t s ) = 0 } be the set of vertices involved in a precedence constraint with t as a head, and such that a zero-cost path from t to s exists. Finally, let V 0 R = { t V R | V 0 ( t ) } be the set of vertices involved as a head in at least one precedence constraint for which an inverse zero-cost path exists.
The Reduced model can be obtained from the Complete model described in Section 3.1 by substituting constraints (4)–(6) and (14) with the following specialized version of them, characterized by a reduced domain:
u s t = 0 t V 0 R , s V 0 ( t )
u t t = 1 t V 0 R
u j t u i t x i j 1 t V 0 R , ( i , j ) A 0
u j t 0 t V 0 R , j V 0
Notice that the domain of the u variables is reduced according to (20).
Compared to the model introduced in Section 3.1, and given that zero-cost paths are rare, the Reduced model uses less variables and constraints (although the theoretical complexity remains unchanged), thus further reducing the memory footprint. However, it might be characterized by a weaker linear relaxation due to the elimination of redundancy in the constraints, on top of having the burden of a preprocessing phase.

3.2.1. Selecting the Precedence Constraints involving a Zero-Cost Path

Zero-cost paths in G 0 that start from t and end in s for some ( s , t ) R can be retrieved by running the procedure described in Algorithm 1.
Algorithm 1 Retrieve the Relevant Zero-Cost Paths from an Instance
1:
Compute the shortest path S P i j for each pair of vertices i , j of G 0
2:
V 0 R =
3:
for all  t V 0  do
4:
       V 0 ( t ) =
5:
      for all  s V 0 : ( s , t ) R  do
6:
            if  C ( S P i j ) = 0 then
7:
                  V 0 ( t ) = V 0 ( t ) { s }
8:
            end if
9:
       end for
10:
      if  V 0 ( t )  then
11:
           V 0 R = V 0 R { t }
12:
      end if
13:
end for
Line 1 can be implemented by running the algorithm of Floyd-Warshall [43] to retrieve the shortest path between each pair of vertices of a graph. The algorithm has a computational complexity of O ( n 3 ) . Lines 2–13 scan the results to populate the sets V 0 ( t ) , containing vertices involved in relevant zero-cost paths for each vertex t V 0 , and the set V 0 R , with a total computational complexity of O ( n 2 ) . Therefore, the overall computational complexity of the procedure remains polynomial, in the order of O ( n 3 ) . This guarantees negligible computation times for the graphs considered in the study and in real applications of the problem, for which n does not exceed 700.

3.3. Solving the New Models via Constraint Programming

Modern Constraint Programming solvers such as Google OR-Tools CP-SAT [44]—the one adopted for the present work—are able to solve compact MILP models efficiently [45]. In particular, these solvers are very effective in treating logical inferences that can be expressed effectively without the use of big-M coefficients (that weaken linear relaxations and consequently worsen solving times) in their syntax. The two MILP models discussed in Section 3.1 and Section 3.2 use such a technique to describe the nonlinear relation between the variable x i j and the set of variables { y i , u j t , d j , w j } in order to turn the constraints off whenever the value of x i j is equal to zero. In this section, we will therefore manipulate the models previously introduced in order to transform the constraints involving big-Ms into logical inferences. Notice that this operation mentioned above is not strictly required, since CP-SAT is able to deal with big-M constraints natively, but according to some preliminary tests, using logical inferences enhances the performance of the solver. Conversely, MILP solvers such as CPLEX [46]—the one adopted for the present work—that are also able to treat logical inferences without big-M constraints present a strong degradation of the performances when big-M constraints are removed. For this reason, in the experiments reported in Section 4, we will use big-M constraints for the MILP solver and logical inferences for the CP solver.
Finally, it can be observed that CP solvers only accept integer-valued variables, which means that the value of the c i j s should be discretized before being passed to the model in case they are not integer. See Montemanni and Dell’Amico [45] for a deeper traction.
In detail, the Complete model can be adapted by modifying constraints (3), (6), (9) and (10), which are substituted by the following ones:
x i j y j = y i + 1 ( i , j ) A : j r
x i j u j t u i t t V R , ( i , j ) A
x i j d j = d i + w j + c i j ( i , j ) A
Constraints (21) implement subtour elimination. The nonlinear relationship y j = ( y i + 1 ) x i j is modeled by setting the value of y j to y i + 1 if x i j = 1 (true). Technically, the logical implication is implemented through the OnlyEnforceIf construct of the CP-SAT solver [44]. Constraints (22) are the precedence-enforcing constraints that set the value of u j t to be greater than or equal to u j t if x i j = 1 and model the nonlinear relationship u j t u i t x i j . Constraints (23) set the value of d j to d i + w j + c i j if x i j = 1 . The set of constraints (23) deals therefore with the nonlinear relationship ( d j d i w j c i j ) x i j = 0 . Notice that the two constraints (9) and (10) are now combined in a single set of constraints.
Analogously, it is possible to obtain a version of the Reduced model based on logical inferences by changing constraints (3) to (21) and substituting constraints (19), (9) and (10) with the following new ones:
x i j u j t u i t t V 0 R , ( i , j ) A 0
x i j d j = d i + w j + c i j ( i , j ) A 0
Notice that constraints (24) and (25) are the versions of (22) and (23) specialized to the reduced graph G 0 for what concerns zero-cost precedences and that constraints (25) cover again both the sets (9) and (10).

4. Computational Experiments

The experimental settings and conditions are described in Section 4.1 together with the instances adopted for the tests. The detailed results are instead presented and commented on in Section 4.2.

4.1. Experimental Settings

The computational experiments we present in order to position the proposed models within the existing literature are based on the benchmark instances of TSPLIB [47], SOPLIB [48] and COMPILERS [49]. All these datasets had originally been proposed for the Sequential Ordering Problem, and they are commonly adopted in the PCMCA-WT literature so far, see [38]. In total, the benchmark sets considered contain 116 instances with sizes ranging between 9 and 700 vertices (and an average of 248 vertices). Currently, the benchmark set has a total of 88 open instances (i.e., without a known optimal solution).
The MILP solver adopted is CPLEX v12.8 [46] with standard settings. The CP solver used is OR-Tools v9.5 [44] CP-SAT, also run with standard settings. The computation time of both solvers is limited to 1 h, while the preprocessing time required for the Reduced methods is not accounted for, being in the order of fractions of a second for all the instances considered.
The best-known solutions used as reference are those appearing in [38]. The value reported is for each instance the best of the results achieved by the three models introduced there. This biases the comparison in favor of the old methods, since—according to the No Free Lunch Theorem [50]—taking the best of different approaches might give a substantial advantage. Among the three ideas disclosed in [38], the one improved in the present work had shown some potential, however, it was the weakest of the set due to severe scalability issues (now solved, see Section 3). On the other hand, the results of [38] were obtained on an Intel i7-8550U processor running at 1.8 GHz and with 8 GB of RAM, while the new experiments are obtained on Intel Xeon Platinum 8375C running at 2.9 GHz and using up to 16 GB of RAM. This gives a hardware advantage to the new models, somehow balancing back the comparison. Notice that the newly proposed models are a direct improvement aiming at overcoming the crucial scalability of one of the three models of [38], making computational fairness considerations less central, in our view.

4.2. Results

An aggregated summary of the results is presented in Table 1. The average optimality gap across all the instances for which all the models were able to find a feasible/optimal solution is reported under Average optimality gap. The average solution time among all the instances that were solved to optimality by all the models can be found under Average solution time. The detailed results of each model can be found instead in Table 2, Table 3 and Table 4, where the following data are reported for each instance. The name and size can be found in the columns with these names. The best-known bounds found in [38] are reported in the column Best-Known [38], where LB shows the best lower bound found, and UB the best-known solution. For each model, the following columns are displayed. The lower and upper bound can be found in the column with these names. The optimality gap, computed as U B L B U B , is reported in the column Gap. The solution time in seconds, which is reported only for those instances that are closed in the given time, can be found in the column Time. Entries in bold indicate new best-known lower or upper bounds. Finally, the name of the instances for which optimality is proven for the first time in this work is highlighted in bold across the tables.
Comparing the average optimality gap of each model, it can be observed that the MILP solver run on the Complete model has an optimality gap of 0.206 on average (0.418 when all the instances solved by the model are considered) but fails to solve two instances, as it runs out of memory. The MILP solver on the Reduced model has an optimality gap of 0.153 on average (with a 25.7% improvement over the previous model) and an optimality gap of 0.340 on average (an 18.7% improvement) across all the instances. The MILP solver on the Reduced model also runs out of memory on one instance. When considering the CP solver, the Complete model shows an optimality gap of 0.159 on average (with a 22.8% improvement), with an optimality gap of 0.157 on average across all the instances. However, the largest instances, with size larger than 200 or with a very dense precedence graph, are not solved, since the model runs out of memory due to the large number of constraints (19). When solving the Reduced model, the CP solver achieves an optimality gap of 0.122 on average (with a 40.7% improvement), with an optimality gap of 0.286 averaged across all the instances.
In terms of solution time, and comparing the instances that are optimally solved by all models (27 instances), using the MILP solver takes the Complete model to a solution time of 690.8 s on average, while the Reduced takes 270.8 s on average (with a 60.8% improvement). When the CP solver is used, the Complete model has a solution time of 166.4 s on average (with a 75.9% improvement), while solving the Reduced model takes 36.4 s on average (with a 94.7% improvement). Cross-comparing a same model when treated by the two different solvers, it emerges that generally the CP models outperform the MILP models on instances with medium to high density precedence graphs.
In terms of solution costs, the Reduced model solved by an MILP solver finds new best-known lower bounds for 24 out of 88 instances (27.3%) compared to the 2 retrieved by the Complete model solved by the same solver. Furthermore, the Reduced model finds new best-known upper bounds for 15 (17.1%) compared to the 1 only found by the Complete model. This indicates that the strength of the linear relaxation of the model is not drastically affected after removing a subset of the variables and constraints from the model. Furthermore, this shows that the Reduced model is generally easier to solve by the MILP solver adopted, and therefore the solver is able to find new bounds more frequently compared to solving the Complete model. When considering the CP solver, the Reduced model finds new best-known lower bounds for 32 (36.4%), while the Complete model finds new lower bounds for 13 (14.8%). For new best-known upper bounds, solving the Reduced model leads to new 54 new bests (61.4%), while solving the Complete model leads to 10 new bests (11.4%). This indicates that the Reduced model is generally more effective to solve by the CP solver when compared to the Complete model. Moreover, the use of the CP solver led to seven newly proven optimal solutions. In general, using the MILP solver seems to produce better lower bounds, while the CP solver is better at finding lower cost solutions.
In summary, the computational results show that the Reduced model generally outperforms the Complete model independently of the solver adopted. This is due to the fact that the Reduced model has a substantially smaller number of variables and constraints, giving an advantage to the solvers. Moreover, the CP solver performs better than the MILP solver in terms of the quality of the solutions, the average solution time and the average optimality gap. Furthermore, the CP solver finds new best-known lower/upper bounds for some instances.

5. Conclusions

This work introduced new models for the Precedence-Constrained Minimum-Cost Arborescence Problem with Waiting-Times that are polynomial in size and are characterized by a smaller memory footprint with respect to the polynomial-sized models previously proposed in the literature. A first model is based on a new set of variables to model precedences. The number of variables and constraints are further reduced, at the price of a preprocessing phase, in a second model. The two models are solved both by Mixed Integer Linear Programs and Constraint Programming solvers. The computational results show that the model characterized by the need of preprocessing outperforms the other one. Furthermore, the Constraint Programming solver achieves the best overall results in terms of both optimality gap and solution time. However, the Mixed Integer Linear Programming solver generally finds better lower bound estimates on the instances. Finally, the models proposed were able to close 7 new instances that were previously open, to provide improved lower bounds for 71 instances, and to find improved upper bounds for 80 instances, out of a total of 88 open instances.
Future work should cover aspects such as robustness of the approaches and the addition to the models of other realistic constraints. Given the progress on the solvers, new instances should be also introduced in order to extend the study on scalability of the different models. Finally, a deeper analysis of the characteristics of the instances that mainly affect the different approaches presented should be in order.

Author Contributions

Conceptualization, J.J., R.M. and M.D.; methodology, J.J.; software, J.J.; validation, J.J. and R.M.; formal analysis, J.J., R.M. and M.D.; investigation, J.J. and R.M.; resources, R.M.; data curation, J.J.; writing—original draft preparation, J.J. and R.M.; writing—review and editing, R.M. and M.D.; visualization, R.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The instances used in the paper are available from the literature.

Conflicts of Interest

The authors declare no conflicts of interest.

Correction Statement

This article has been republished with a minor correction to the Document Type. This change does not affect the scientific content of the article.

References

  1. Chu, Y.J.; Liu, T. On the shortest arborescence of a directed graph. Sci. Sin. 1965, 14, 1396–1400. [Google Scholar]
  2. Edmonds, J. Optimum branchings. J. Res. Natl. Bur. Stand. 1967, 71, 233–240. [Google Scholar] [CrossRef]
  3. Gabow, H.N.; Galil, Z.; Spencer, T.; Tarjan, R.E. Efficient algorithms for finding minimum spanning trees in undirected and directed graphs. Combinatorica 1986, 6, 109–122. [Google Scholar] [CrossRef]
  4. Chou, X.; Gambardella, L.M.; Montemanni, R. A tabu search algorithm for the probabilistic orienteering problem. Comput. Oper. Res. 2021, 126, 105107. [Google Scholar] [CrossRef]
  5. Bock, F. An algorithm to construct a minimum directed spanning tree in a directed network. Dev. Oper. Res. 1971, 29–44. [Google Scholar]
  6. Fischetti, M.; Vigo, D. A branch-and-cut algorithm for the resource-constrained minimum-weight arborescence problem. Netw. Int. J. 1997, 29, 55–67. [Google Scholar] [CrossRef]
  7. Pereira, A.H.; Mateus, G.R.; Urrutia, S. Branch-and-cut algorithms for the p-arborescence star problem. Int. Trans. Oper. Res. 2021, 29, 2374–2400. [Google Scholar] [CrossRef]
  8. Morais, V.; Gendron, B.; Mateus, G.R. The p-arborescence star problem: Formulations and exact solution approaches. Comput. Oper. Res. 2019, 102, 91–101. [Google Scholar] [CrossRef]
  9. Hakimi, S.L. Optimum distribution of switching centers in a communication network and some related graph theoretic problems. Oper. Res. 1965, 13, 462–475. [Google Scholar] [CrossRef]
  10. Guttmann-Beck, N.; Hassin, R. On two restricted ancestors tree problems. Inf. Process. Lett. 2010, 110, 570–575. [Google Scholar] [CrossRef]
  11. Carrabs, F.; Gaudioso, M. A Lagrangian approach for the minimum spanning tree problem with conflicting edge pairs. Networks 2021, 78, 32–45. [Google Scholar] [CrossRef]
  12. Kruskal, J.B. On the shortest spanning subtree of a graph and the traveling salesman problem. Proc. Am. Math. Soc. 1956, 7, 48–50. [Google Scholar] [CrossRef]
  13. Gouveia, L.; Lopes, M.J. The capacitated minimum spanning tree problem: On improved multistar constraints. Eur. J. Oper. Res. 2005, 160, 47–62. [Google Scholar] [CrossRef]
  14. Frieze, A.M.; Tkocz, T. A Randomly Weighted Minimum Arborescence with a Random Cost Constraint. Math. Oper. Res. 2021, 47, 1664–1680. [Google Scholar] [CrossRef]
  15. Fertin, G.; Fradin, J.; Jean, G. Algorithmic Aspects of the Maximum Colorful Arborescence Problem. In Theory and Applications of Models of Computation; TAMC 2017; Springer: Cham, Switzerland, 2017; pp. 216–230. [Google Scholar]
  16. Eswaran, K.P.; Tarjan, R.E. Augmentation problems. SIAM J. Comput. 1976, 5, 653–665. [Google Scholar] [CrossRef]
  17. Li, J.; Liu, X.; Lichen, J. The constrained arborescence augmentation problem in digraphs. In Proceedings of the 2017 3rd IEEE International Conference on Computer and Communications (ICCC), Chengdu, China, 13–16 December 2017; pp. 1204–1209. [Google Scholar]
  18. Kawatra, R.; Bricker, D. Design of a degree-constrained minimal spanning tree with unreliable links and node outage costs. Eur. J. Oper. Res. 2004, 156, 73–82. [Google Scholar] [CrossRef]
  19. Galbiati, G.; Gualandi, S.; Maffioli, F. On minimum changeover cost arborescences. Lect. Notes Comput. Sci. 2011, 6630, 112–123. [Google Scholar]
  20. Bérczi, K.; Fujishige, S.; Kamiyama, N. A linear-time algorithm to find a pair of arc-disjoint spanning in-arborescence and out-arborescence in a directed acyclic graph. Inf. Process. Lett. 2009, 109, 1227–1231. [Google Scholar] [CrossRef]
  21. Bang-Jensen, J. Edge-disjoint in- and out-branchings in tournaments and related path problems. J. Comb. Theory—Ser. B 1991, 51, 1–23. [Google Scholar] [CrossRef]
  22. Li, Y.; Thai, M.; Wang, F.; Du, D.Z. On the construction of a strongly connected broadcast arborescence with bounded transmission delay. IEEE Trans. Mob. Comput. 2006, 5, 1460–1470. [Google Scholar]
  23. Carrabs, F.; Cerulli, R.; Pentangelo, R.; Raiconi, A. Minimum spanning tree with conflicting edge pairs: A branch-and-cut approach. Ann. Oper. Res. 2021, 298, 65–78. [Google Scholar] [CrossRef]
  24. Darmann, A.; Pferschy, U.; Schauer, J. Determining a Minimum Spanning Tree with Disjunctive Constraints. In Algorithmic Decision Theory; Springer: Berlin/Heidelberg, Germeny, 2009; pp. 414–423. [Google Scholar]
  25. Viana, L.A.d.C.; Campêlo, M. Two dependency constrained spanning tree problems. Int. Trans. Oper. Res. 2020, 27, 867–898. [Google Scholar] [CrossRef]
  26. Escudero, L. An inexact algorithm for the sequential ordering problem. Eur. J. Oper. Res. 1988, 37, 236–249. [Google Scholar] [CrossRef]
  27. Moon, C.; Kim, J.; Choi, G.; Seo, Y. An efficient genetic algorithm for the traveling salesman problem with precedence constraints. Eur. J. Oper. Res. 2002, 140, 606–617. [Google Scholar] [CrossRef]
  28. Balas, E.; Fischetti, M.; Pulleyblank, W. The precedence-constrained asymmetric traveling salesman polytope. Math. Program. 1995, 68, 241–265. [Google Scholar] [CrossRef]
  29. Hernádvölgyi, I. Solving the sequential ordering problem with automatically generated lower bounds. In Operations Research Proceedings 2003; Springer: Berlin/Heidelberg, Germany, 2004; pp. 355–362. [Google Scholar]
  30. Escudero, L.; Guignard, M.; Malik, K. A Lagrangian relax-and-cut approach for the sequential ordering problem with precedence relationships. Ann. Oper. Res. 1994, 50, 219–237. [Google Scholar] [CrossRef]
  31. Gambardella, L.; Dorigo, M. An ant colony system hybridized with a new local search for the sequential ordering problem. INFORMS J. Comput. 2000, 12, 237–255. [Google Scholar] [CrossRef]
  32. Karan, M.; Skorin-Kapov, N. A branch and bound algorithm for the sequential ordering problem. In Proceedings of the MIPRO, 2011 Proceedings of the 34th International Convention, Opatija, Croatia, 23–27 May 2011; pp. 452–457. [Google Scholar]
  33. Ascheuer, N.; Escudero, L.; Grötschel, M.; Stoer, M. A cutting plane approach to the sequential ordering problem (with applications to job scheduling in manufacturing). SIAM J. Optim. 1993, 3, 25–42. [Google Scholar] [CrossRef]
  34. Ascheuer, N.; Jünger, M.; Reinelt, G. A branch & cut algorithm for the asymmetric traveling salesman problem with precedence constraints. Comput. Optim. Appl. 2000, 17, 61–84. [Google Scholar]
  35. Montemanni, R.; Smith, D.H.; Gambardella, L.M. Ant colony systems for large sequential ordering problems. In Proceedings of the IEEE Swarm Intelligence Symposium (SIS), Honolulu, HI, USA, 1–5 April 2007; pp. 60–67. [Google Scholar]
  36. Fiala Timlin, M.; Pulleyblank, W. Precedence constrained routing and helicopter scheduling: Heuristic design. Interfaces 1992, 22, 100–111. [Google Scholar] [CrossRef]
  37. Dell’Amico, M.; Jamal, J.; Montemanni, R. A mixed integer linear program for a precedence-constrained minimum-cost arborescence problem. In Proceedings of the 8th International Conference on Industrial Engineering and Applications (Europe), Online, 8–11 January 2021; pp. 216–221. [Google Scholar]
  38. Chou, X.; Dell’Amico, M.; Jamal, J.; Montemanni, R. Precedence-Constrained Arborescences. Eur. J. Oper. Res. 2022, 307, 575–589. [Google Scholar] [CrossRef]
  39. Shi, W.; Su, C. The rectilinear Steiner arborescence problem is NP-complete. SIAM J. Comput. 2005, 35, 729–740. [Google Scholar] [CrossRef]
  40. Wang, I.L. Multicommodity network flows: A survey, Part I: Applications and Formulations. Int. J. Oper. Res. 2018, 15, 145–153. [Google Scholar]
  41. Hurkensa, C.A.J.; Woeginger, G.J. On the nearest neighbor rule for the traveling salesman problem. Oper. Res. Lett. 2004, 32, 1–4. [Google Scholar] [CrossRef]
  42. Dell’Amico, M.; Jamal, J.; Montemanni, R. Compact Models for the Precedence-Constrained Minimum-Cost Arborescence Problem. In Proceedings of the 2022 The 6th International Conference on Intelligent Traffic and Transportation (ICITT), Paris, France, 25–27 September 2022; IOS Press: Amsterdam, The Netherlands, 2023; pp. 112–126. [Google Scholar]
  43. Floyd, R. Algorithm 97: Shortest Path. Commun. ACM 1962, 5, 345. [Google Scholar] [CrossRef]
  44. Google. Google OR-Tools. 2023. Available online: https://developers.google.com/optimization (accessed on 20 November 2023).
  45. Montemanni, R.; Dell’Amico, M. Solving the Parallel Drone Scheduling Traveling Salesman Problem via Constraint Programming. Algorithms 2023, 16, 40. [Google Scholar] [CrossRef]
  46. IBM. IBM CPLEX Optimizer. 2023. Available online: https://www.ibm.com/products/ilog-cplex-optimization-studio/cplex-optimizer (accessed on 20 November 2023).
  47. Reinelt, G. TSPLIB–A travelling salesman problem library. ORSA J. Comput. 1991, 3, 376–384. [Google Scholar] [CrossRef]
  48. Montemanni, R.; Smith, D.H.; Rizzoli, A.E.; Gambardella, L.M. Sequential ordering problems for crane scheduling in port terminals. Int. J. Simul. Process Model. 2009, 5, 348–361. [Google Scholar] [CrossRef]
  49. Shobaki, G.; Jamal, J. An exact algorithm for the sequential ordeing problem and its application to switching energy minimization in compilers. Comput. Optim. Appl. 2015, 61, 343–372. [Google Scholar] [CrossRef]
  50. Wolpert, D.; Macready, W. No Free Lunch Theorems for Optimization. IEEE Trans. Evol. Comput. 1997, 1, 67. [Google Scholar] [CrossRef]
Figure 1. Instance of the PCMCA-WT problem and relative solution. The instance is depicted in the top of the figure, with its arc costs. The precedence relationship ( 1 , 3 ) R is depicted through a dashed arrow. The bottom part of the figure depicts the optimal arborescence of cost 8.
Figure 1. Instance of the PCMCA-WT problem and relative solution. The instance is depicted in the top of the figure, with its arc costs. The precedence relationship ( 1 , 3 ) R is depicted through a dashed arrow. The bottom part of the figure depicts the optimal arborescence of cost 8.
Algorithms 17 00012 g001
Figure 2. Example of how a precedence-violating path is detected using constraints (4)–(6).
Figure 2. Example of how a precedence-violating path is detected using constraints (4)–(6).
Algorithms 17 00012 g002
Table 1. Summary of the results achieved with each solver/model combination.
Table 1. Summary of the results achieved with each solver/model combination.
MILP SolverCP Solver
CompleteReducedCompleteReduced
ModelModelModelModel
Average optimality gap0.2060.1530.1590.122
Average solution time690.8270.8166.436.4
New best-known lower bounds2241332
New best-known upper bounds1151054
New optimal solutions0077
Table 2. Computational results for TSPLIB instances.
Table 2. Computational results for TSPLIB instances.
MILP SolverCP Solver
InstanceComplete ModelReduced ModelComplete ModelReduced Model
NameSizeBest-Known [38]LBUBGapTime [s]LBUBGapTime [s]LBUBGapTime [s]LBUBGapTime [s]
br17.1018[35, 44]39440.114-40440.091-44440.00046.62944440.00056.628
br17.1218[35, 44]41440.068-41440.068-44440.00022.60444440.00044.550
ESC0791906190619060.0000.028190619060.0000.070190619060.0000.023190619060.0000.025
ESC11132174217421740.0000.125217421740.0000.114217421740.0000.107217421740.0000.077
ESC12141138113811380.0000.035113811380.0000.030113811380.0000.034113811380.0000.037
ESC25271158115811580.0006.185115811580.0001.945115811580.0000.910115811580.0000.833
ESC47497477477470.00059.7607477470.00022.1537477470.0003.8867477470.0002.708
ESC63655656560.00024.60056560.00057.34756560.0001.51756560.0002.465
ESC78801196119611960.0002410.483119611960.000257.609119611960.000100.511119611960.00018.971
ft53.1544089408940890.0001764.235408940890.0002023.553408940890.000215.480408940890.000291.099
ft53.254[4135, 4284]411243170.047-416143340.040-410242840.042-410342840.042-
ft53.354[4623, 5457]474654250.125-479952790.091-4493600.161-450854840.178-
ft53.454[5657, 6439]592264200.078-592364200.077-533865020.179-535764200.166-
ft70.171[33,128, 33,298]32,77733,3080.016-32,82733,3080.014-32,66933,4720.024-33,10133,2980.006-
ft70.271[33,357, 34,450]33,05733,9770.027-33,08933,9160.024-32,93833,6700.022-32,89733,6700.023-
ft70.371[33,914, 42,732]34,15238,5460.114-34,42338,3510.102-33,82536,9390.084-33,81336,9320.084-
ft70.471[36,517, 40,404]36,73739,1450.062-36,85038,7710.050-33,82536,9390.084-35,66439,8430.105-
rbg048a50[261, 264]2602650.019-2592640.019-2632630.0009.4422632630.00025.294
rbg050c522252252250.000863.6622252250.00036.6732252250.0002.5752252250.0001.234
rbg109111[354, 414]3544260.169-3664070.101-3574880.268-3594010.105-
rbg150a152[447, 541]4475110.125-4615090.094-4635910.217-4615170.108-
rbg174a176[446, 580]4526010.248-4635530.163-4575710.200-4615720.194-
rbg253a255[477, 773]52312520.582-5327180.259-----5277220.270-
rbg323a325[926, 4035]98110,1110.903-97424660.605-----100918910.466-
rbg341a343[681, 3800]76493130.918-76129070.738-----78014570.465-
rbg358a360[706, 3296]95011,5280.918-75524530.692-----78811500.315-
rbg378a380[649, 2759]67210,2420.934-64821910.704-----67811260.398-
kro124p.1101[32,858, 35,231]32,65137,1200.120-32,63036,0990.096-32,50434,1000.047-32,56133,9620.041-
kro124p.2101[33,190, 37,956]32,88642,5730.228-33,00639,9310.173-32,76437,0740.116-32,79935,8600.085-
kro124p.3101[34,217, 53,988]33,81354,1830.376-34,00546,7640.273-33,56143,9100.236-33,48842,4160.210-
kro124p.4101[39,413, 55,187]39,96958,9440.322-39,33353,4560.264-38,43350,9100.245-37,67649,5900.240-
p43.144[2827, 4470]266040850.349-265639550.328-286039550.277-285139900.285-
p43.244[2826, 4275]99144500.777-270542100.357-285641600.313-287041800.313-
p43.344[2864, 5375]106750150.787-138344400.689-296644500.333-289742550.319-
p43.444[3101, 4900]299550350.405-312546050.321-309044950.313-309446200.330-
prob.100100[674, 1008]66821250.686-6777410.086-6667840.151-6677380.096-
prob.42421711711710.000396.4581711710.000230.5061711710.00079.6671711710.00034.245
ry48p.149[13,371, 13,722]13,11414,2720.081-13,20013,6700.034-13,03613,6700.046-13,06113,6700.045-
ry48p.249[13,508, 14,659]13,29914,4150.077-13,33614,3050.068-13,21614,3050.076-13,18514,3050.078-
ry48p.349[14,371, 16,326]13,88216,1930.143-13,99415,8400.117-13,76415,5460.115-13,72815,4770.113-
ry48p.449[17,339, 19,649]17,16219,7440.131-17,18019,5830.123-16,55019,8370.166-16,48319,4950.155-
Average 0.158552.557 0.167263.000 0.10337.184 0.12836.782
Table 3. Computational results for SOPLIB instances.
Table 3. Computational results for SOPLIB instances.
MILP SolverCP Solver
InstanceComplete ModelReduced ModelComplete ModelReduced Model
NameSizeBest-Known [38]LBUBGapTime [s]LBUBGapTime [s]LBUBGapTime [s]LBUBGapTime [s]
R.200.100.12002929290.00018.39429290.0006.01729290.00028.27129290.00031.403
R.200.100.15200[505, 1271]49714310.653-52510330.492-38118640.796-5899790.398-
R.200.100.30200[669, 2011]68632520.789-77417610.560-45130010.850-83818710.552-
R.200.100.60200[8070, 18,761]876017,0040.485-886116,9300.477-601831,5610.809-844016,1970.479
R.200.1000.12008878878870.0001288.0928878870.00015.6358878870.000649.9798878870.00026.0915
R.200.1000.15200[6665, 16,496]676916,3360.586-689512,6010.453-531825,1960.789-723112,8120.436-
R.200.1000.30200[9340, 30,351]993723,2260.572-10,51222,7810.539-738138,4100.808-10,12023,2490.565-
R.200.1000.60200[10,508, 23,748]11,39921,7060.475-12,04221,9930.452-666628,5220.766-10,66519,9340.465
R.300.100.13001313130.00037.35213130.00035.01213130.000205.73113130.00056.4263
R.300.100.15300[625, 12,903]66069580.905-66922590.704-----81120560.606-
R.300.100.30300[948, 3767]100867900.852-110231630.652-----115725900.553-
R.300.100.60300[824, 3005]91947320.806-94919540.514-----99118650.469-
R.300.1000.13007157157150.0003187.0497157150.00064.6837157150.000257.0747157150.00071.6789
R.300.1000.15300[7213, 112,424]7607110,3660.931-783224,0470.674-----876829,4230.702-
R.300.1000.30300[10,385, 40,457]11,17953,8350.792-12,07140,8630.705-----12,26931,6180.612-
R.300.1000.60300[9413, 30,655]10,18038,2120.734-10,27525,3230.594-----10,40821,6230.519-
R.400.100.1400663760.984-660.000995.137660.000726.057660.00097.3851
R.400.100.15400[729, 47,117]78135,0440.978-85622,7670.962-----96335910.732-
R.400.100.30400[780, 7243]91139,0220.977-101026,4380.962-----108430610.646-
R.400.100.60400[731, 5545]83733090.747-86126520.675-----96620690.533-
R.400.1000.14007807807800.000161.0217807800.000124.9907807800.000208.5257807800.00090.9555
R.400.1000.15400[7760, 501,543]835785,8780.903-908385,8780.894-----997635,1600.716-
R.400.1000.30400[10,076, 95,523]11,030127,2900.913-11,783127,2900.907-----12,33757,2720.785-
R.400.1000.60400[8103, 55,950]936065,6150.857-987736,6620.731-----995422,3760.555-
R.500.100.15003330.0002157.743330.0001881.297330.0002333.235330.000112.086
R.500.100.15500[924, 11,452]96411,4520.916-101811,4520.911-----125055080.773-
R.500.100.30500[773, 12,225]84916,9630.950-97614,2730.932-----109948410.773-
R.500.100.60500[669, 8427]84049,1050.983-84063570.868-----93127230.658-
R.500.1000.15002972972970.00097.4732972970.00085.4592972970.00085.2812972970.00077.4382
R.500.1000.15500[8420, 107,776]8949107,7760.917-9461107,7760.912-----10,62845,3560.766
R.500.1000.30500[10,431, 181,835]11,799156,3590.925-12,694156,3590.919-----12,57657,3300.781
R.500.1000.60500[7094, 33,260]8233112,4660.927-819245,6960.821-----655920,4650.680-
R.600.100.1600[1, 379]1550.982-1550.982-110.0002710.470110.0002182.18
R.600.100.15600[670, 5949]71459310.880-84540440.791-----93824430.616-
R.600.100.30600[873, 12,875]94518,9320.950-109918,9320.942-----74064670.886-
R.600.100.60600[751, 7893]83826,7320.969-77825,2140.969-----53824940.784-
R.600.1000.16003223223220.000352.2023223220.000140.6453223220.000127.3973223220.000103.378
R.600.1000.15600[10,181, 121,877]10,753121,8770.912-10,915121,8770.910-----940165,0390.855-
R.600.1000.30600[10,151, 151,010]11,352190,1450.940-12,431190,1450.935-----935648,7750.808-
R.600.1000.60600[7604, 87,770]7962256,4640.969-816275,2690.892-----690842,6520.838-
R.700.100.17002--------220.0001649.486220.000619.22
R.700.100.15700[799, 6561]81514,4780.944-97257180.830-----65527590.763-
R.700.100.30700[762, 20,281]89669600.871-98342180.767-----58825310.768-
R.700.100.60700[516, 9030]53870330.924-55518540.701-----38315980.760-
R.700.1000.1700[611, 621]6116160.008-6116160.008-6116110.000592.1076116110.000368.139
R.700.1000.15700[4636, 147,321]4375147,3210.970-513671450.281-----278763150.559-
R.700.1000.30700[4303, 50,000]447732,7420.863-482769810.309-----265861150.565-
R.700.1000.60700[2857, 15,579]294285340.655-299758420.487-----191353570.643-
Average 0.689912.416 0.577372.097 0.268797.801 0.492319.698
Table 4. Computational results for COMPILERS instances.
Table 4. Computational results for COMPILERS instances.
MILP SolverCP Solver
InstanceComplete ModelReduced ModelComplete ModelReduced Model
NameSizeBest-Known [38]LBUBGapTime [s]LBUBGapTime [s]LBUBGapTime [s]LBUBGapTime [s]
gsm.153.124126[246, 313]2573120.176-2693110.135-2783170.123-2803110.100-
gsm.444.350353[2103, 2873]229448780.530-240548560.505-----245643100.430-
gsm.462.7779[396, 488]4024780.159-4024770.157-4194740.116-4184650.101-
jpeg.1483.25278787870.00026.04187870.00018.55687870.0001.19487870.0001.071
jpeg.3184.107109[489, 684]5066560.229-5107150.287-5187180.279-5176920.253-
jpeg.3195.8587[22, 25]17250.320-17250.320-23250.080-22250.120-
jpeg.3198.9395[172, 204]1801880.043-1801880.043-1811880.037-1811880.037-
jpeg.3203.135137[600, 750]6029800.386-6187510.177-6299130.311-6267500.165-
jpeg.3740.15173333330.0001.52333330.0000.83933330.0000.15733330.0000.095
jpeg.4154.36389090900.000556.79890900.00060.92490900.0001.27290900.0001.764
jpeg.4753.54561641641640.0002753.7521641640.0001790.2691641640.00015.3421641640.00016.877
susan.248.197199[736, 1184]79219780.600-80213700.415-80513610.409-78013200.409-
susan.260.158160[564, 876]5689370.394-5739380.389-5969910.399-5988970.333-
susan.343.182184[591, 862]6177980.227-6227760.198-63610430.390-6327920.202-
typeset.10192.123125[280, 415]2744290.361-2823790.256-2933850.239-2923870.245-
typeset.10835.2628[99, 112]991110.108-1001120.107-1101110.009-1091110.018-
typeset.12395.4345[143, 146]1401460.041-1411460.034-1461460.0002181.9421461460.0002780.121
typeset.15087.23259797970.00060.50297970.00029.11897970.0000.47797970.0000.318
typeset.15577.36381251251250.000286.2101251250.00043.1641251250.0002.1161251250.0001.713
typeset.16000.6870[77, 80]66810.185-66800.175-79800.013-71800.113-
typeset.1723.25276060600.000590.57760600.00086.06860600.0004.01360600.0003.469
typeset.19972.246248[1325, 1929]142235620.601-145225090.421-151929610.487-152528040.456-
typeset.4391.240242[1093, 1412]110825950.573-113724760.541-114925110.542-115419050.394-
typeset.4597.4547[150, 155]1501540.026-1511540.019-1541540.000209.6591541540.000128.916
typeset.4724.433435[2460, 3433]----267361310.564-----267971940.628-
typeset.5797.33351131131130.000851.4901131130.00028.5041131130.0000.5471131130.0000.574
typeset.5881.246248[1305, 1700]137822580.390-139624260.425-140623850.410-139420840.331-
Average 0.206640.862 0.191257.180 0.154241.672 0.161293.492
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dell’Amico, M.; Jamal, J.; Montemanni, R. Compact Models to Solve the Precedence-Constrained Minimum-Cost Arborescence Problem with Waiting Times. Algorithms 2024, 17, 12. https://doi.org/10.3390/a17010012

AMA Style

Dell’Amico M, Jamal J, Montemanni R. Compact Models to Solve the Precedence-Constrained Minimum-Cost Arborescence Problem with Waiting Times. Algorithms. 2024; 17(1):12. https://doi.org/10.3390/a17010012

Chicago/Turabian Style

Dell’Amico, Mauro, Jafar Jamal, and Roberto Montemanni. 2024. "Compact Models to Solve the Precedence-Constrained Minimum-Cost Arborescence Problem with Waiting Times" Algorithms 17, no. 1: 12. https://doi.org/10.3390/a17010012

APA Style

Dell’Amico, M., Jamal, J., & Montemanni, R. (2024). Compact Models to Solve the Precedence-Constrained Minimum-Cost Arborescence Problem with Waiting Times. Algorithms, 17(1), 12. https://doi.org/10.3390/a17010012

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop