Next Article in Journal
Statistical Indicators of the Scientific Publications Importance: A Stochastic Model and Critical Look
Next Article in Special Issue
Comparison of Three Computational Approaches for Tree Crop Irrigation Decision Support
Previous Article in Journal
Optimal One-Point Iterative Function Free from Derivatives for Multiple Roots
Previous Article in Special Issue
Innovative Platform for Designing Hybrid Collaborative & Context-Aware Data Mining Scenarios
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Efficient Hybrid Genetic Approach for Solving the Two-Stage Supply Chain Network Design Problem with Fixed Costs

Department of Mathematics and Computer Science, North University Center of Baia Mare, Technical University of Cluj-Napoca, 430083 Baia Mare, Romania
*
Authors to whom correspondence should be addressed.
Mathematics 2020, 8(5), 712; https://doi.org/10.3390/math8050712
Submission received: 31 March 2020 / Revised: 23 April 2020 / Accepted: 24 April 2020 / Published: 3 May 2020
(This article belongs to the Special Issue Computational Intelligence)

Abstract

:
This paper deals with a complex optimization problem, more specifically the two-stage transportation problem with fixed costs. In our investigated transportation problem, we are modeling a distribution network in a two-stage supply chain. The considered two-stage supply chain includes manufacturers, distribution centers, and customers, and its principal feature is that in addition to the variable transportation costs, we have fixed costs for the opening of the distribution centers, as well as associated with the routes. In this paper, we describe a different approach for solving the problem, which is an effective hybrid genetic algorithm. Our proposed hybrid genetic algorithm is constructed to fit the challenges of the investigated supply chain network design problem, and it is achieved by incorporating a linear programming optimization problem within the framework of a genetic algorithm. Our achieved computational results are compared with the existing solution approaches on a set of 150 benchmark instances from the literature and on a set of 50 new randomly generated instances of larger sizes. The outputs proved that we have developed a very competitive approach as compared to the methods that one can find in the literature.

1. Introduction

Supply chains (SCs) are considered to be worldwide networks in which the actors are: suppliers, manufacturer plants, distribution centers (DCs), retailers, and customers, and their principal objective is the fulfillment of the customer needs. In order to obtain an effective management of SC systems, researchers have emphasized the transportation system design, as it plays a significant role within the SC. A supply chain is represented usually as a multi-level structure, while its optimal design has been acknowledged to be an NP-hard problem. For more information on supply chain network design, we refer to Govindan et al. [1], Klibi et al. [2], Melo et al. [3], Wang [4], Dotoli [5,6], etc.
The goal of a transportation model is to minimize the total cost in transporting goods from a set of sources to a set of destinations, fulfilling the request of the destinations and using the capacities of the sources. In practical applications, fixed costs are associated with the arcs connecting sources to destinations, supplementary to the variable transportation costs, which are proportional to the amount of goods distributed along the arcs, resulting in a problem known in the literature as the fixed cost transportation problem (FCTP). Obviously, FCTP is a generalization of the classical transportation problem, and it was introduced by Balinski [7]. Guisewite and Pardalos [8] showed that the FCTP is N P -hard. For more information on the FCTP including a review of exact and heuristic approaches developed for solving the problem, we refer to Buson et al. [9].
This paper focuses on a variation of the fixed cost transportation problem in a supply chain network, namely the two-stage supply chain network design problem in which we consider two kinds of fixed costs: ones for opening the DCs and the others associated with the routes between manufacturers and DCs and between DCs and customers. In the form considered in our paper, the two-stage supply chain network design problem with fixed costs (TSSCNDP-FC) was defined by Hong et al. [10]. The same authors proposed an integer linear programming model of the problem, as well as a solution approach based on ant colony optimization tested on a set of 150 instances split into three classes: small, medium, and large sized instances. Recently, Sabo et al. [11] described a valid model of the problem and pointed out some inaccuracies regarding the paper published by Hong et al. [10].
Some other two-stage transportation problems with fixed costs considered in the literature and related to the investigated problem are:
  • The two-stage transportation problem with fixed costs associated with the routes: Raj and Rajendran [12] proposed two scenarios of the two-stage transportation problem: the first one, called Scenario 1, takes into consideration fixed costs associated with the routes in addition to unit transportation costs and boundless capacities of the DCs, while the second one, called Scenario 2, considers the opening costs of the DCs in addition to unit transportation costs. The same authors developed a genetic algorithm (GA) with a particular coding scheme applicable for two-stage transportation problems, and also, they provided a set of 20 benchmark instances. Another GA dealing with the two-stage transportation problem with fixed charge associated with the routes from plants to customers through DCs was proposed by Jawahar and Balaji [13]. Pop et al. [14] proposed a hybrid method that combines a steady-state GA with a powerful local search procedure. Cosma et al. [15] described an efficient multi-start iterated local search (ILS) procedure for the total transportation cost minimization of the two-stage transportation problem, which begins with a feasible solution of the problem, makes use of a local search procedure with the goal of increasing the exploration, a perturbation mechanism, and a neighborhood operator with the scope of diversifying the search.
  • The two-stage transportation problem with fixed costs for opening the distribution centers (DCs): This two-stage transportation problem was introduced by Gen et al. [16]. The present literature regarding the two-stage transportation problem with fixed costs for opening the DCs is rather limited. This optimization problem has also been investigated by Raj and Rajendran [12], who called it Scenario 2. Calvete et al. [17] proposed a hybrid evolutionary algorithm whose principal characteristic is the employment of a chromosome encoding that offers information about the DCs used within the transportation system. Cosma et al. [18] described an effective heuristic algorithm that reduces the solution search space to a subspace with a reasonable size, without losing optimal or sub-optimal solutions by means of a perturbation mechanism that allows the reconsideration of the feasible solutions that are discarded and that might lead to such solutions. Lately, Cosma et al. [19] proposed a matheuristic approach for solving the two-stage transportation problem with fixed costs associated with the routes by incorporating a linear programming optimization problem within the framework of a genetic algorithm.
  • A particular case is where there exists only one plant manufacturer, and this version was considered by Molla et al. [20]. They proposed an integer linear programming mathematical model of the problem, and in addition, they described two solution approaches for solving it: a spanning tree-based genetic algorithm with a Prüfer number representation and an artificial immune algorithm. Some remarks regarding the mathematical model of the problem were published by El-Sherbiny [21]. Pintea et al. [22] proposed some hybrid algorithms, and Pintea and Pop [23] described an efficient hybrid approach combining the nearest neighbor search heuristic with a local search procedure for solving this particular two-stage transportation problem with fixed costs. Pop et al. [24] developed an innovative hybrid heuristic method achieved by combining a genetic algorithm based on a hash table coding of the individuals with a powerful local search procedure. Recently, Cosma et al. [25] described an effective hybrid heuristic approach that builds an initial feasible solution, then uses a local search procedure whose goal is to increase the exploration and a neighborhood structure for diversifying the search.
  • Another two-stage transportation problem takes into consideration its effect on the environment by reducing the greenhouse gas emissions and was introduced by Santibanez-Gonzales [26] for dealing with a practical application from the public sector. For this version of the problem, Pintea et al. [27] described a set of hybrid heuristic methods, and Pop et al. [28] proposed an effective reverse distribution system for solving it.
As we can observe, the investigated supply chain network design problem generalizes the previously mentioned transportation problems by considering simultaneously two types of fixed costs: ones associated with the transportation routes and the others for opening the DCs in addition to the variable transportation costs, which are proportional to the amount of goods distributed along the arcs.
We organize the remainder of the paper as follows: in Section 2, we define the investigated two-stage supply chain problem with fixed costs and present a set of notations that will be used throughout the paper, and in Section 3, we describe a mixed integer linear formulation of the problem. The novel hybrid method, which incorporates a linear programming problem within the framework of a genetic algorithm, is presented in Section 4, and the comprehensive computational experiments with their outcomes are showcased and analyzed in Section 5. Finally, we conclude our work and discuss our plans for future work in Section 6.

2. Definition of the Two-Stage Supply Chain Network Design Problem with Fixed Costs for Opening the Distribution Centers and Transportation Routes

In order to define and model the two-stage supply chain network design problem with fixed costs for opening the DCs and transportation routes, we consider a tripartite directed graph G = ( V , A ) that consists of a set of vertices V = V 1 V 2 V 3 and a set of arcs A = A 1 A 2 defined as follows:
A 1 = { ( i , j ) | i V 1 a n d j V 2 } and A 2 = { ( j , k ) | j V 2 and k V 3 }
The entire set of nodes V is partitioned into three mutually exclusive sets corresponding to the set of manufacturers denoted by V 1 with | V 1 | = m , the set of distribution centers denoted by V 2 with | V 2 | = d , and the set of customers denoted by V 3 with | V 3 | = r .
In addition, we suppose that:
  • Every manufacturer i V 1 has S i units of supply; every distribution center j V 2 has a given capacity S C j ; each customer k V 3 has a demand D k ;
  • Every manufacturer may transport to any of the q distribution centers at a transportation cost c i j per unit from manufacturer i V 1 to DC j V 2 ;
  • Every DC may transport to any of the r customers at a transportation cost c j k per unit from DC j V 2 to customer k V 3 ;
  • In order to open any of the DCs, we have to pay a given fixed cost denoted by f j , and there exist fixed transportation costs from each manufacturer to each distribution center, denoted by f i j , where i V 1 and j V 2 , and from each DC to each customer, denoted by f j k , where j V 2 and k V 3 .
The aim of the two-stage supply chain network design problem with fixed costs associated with the transportation routes and for opening the DCs is to select the DCs and the routes to be opened and the corresponding transported quantities on these routes, such that the demands of the customers are satisfied, all transportation restrictions are fulfilled, and the total transportation costs are minimized.
Figure 1 illustrates the two-stage supply chain network design problem with fixed costs that we investigated.

3. A Valid Mathematical Model of the Two-Stage Supply Chain Network Design Problem with Fixed Costs

In this section, we present a valid mathematical formulation based on mixed integer programming of the investigated two-stage supply chain network design problem with fixed costs associated with the transportation routes and for opening the DCs.
We introduce the following decision variables:
  • Linear variables:
    -
    x if , specifying the number of units shipped from plant i to the DC j;
    -
    x j k , specifying the number of units shipped from DC j to the customer k;
  • Binary variables:
    -
    y if , specifying if there are units transported from manufacturer plant i to the DC j ( y i j = 1 , if x i j > 0 , and y i j = 0 , otherwise);
    -
    y j k , specifying if there are units transported from DC j to the customer k ( y j k = 1 , if x j k > 0 , and y j k = 0 , otherwise);
    -
    z j , specifying if the DC j is open ( z j = 1 , if the DC j is open, and z j = 0 , otherwise).
Then, the two-stage supply chain network design problem with fixed costs can be formulated as the following mixed integer problem, described by Sabo et al. [11]:
Z = min i = 1 m j = 1 d ( c if x if + f if y if ) + j = 1 d k = 1 r ( c j k x j k + f j k y j k ) + j = 1 d f j z j
s . t . j = 1 d x i j S i , i V 1
j = 1 d x j k = D k , k V 3
i = 1 m x i j = k = 1 r x j k , j V 2
k = 1 r x j k S C j · z j , j V 2
x i j 0 , i V 1 , j V 2
x j k 0 , j V 2 , k V 3
y i j { 0 , 1 } , i V 1 , j V 2
y j k { 0 , 1 } , j V 2 , k V 3
z j { 0 , 1 } , j V 2
Our objective is to minimize the total transportation cost including the unit transportation costs and the fixed costs (for opening DC’s and associated with the routes). Constraint (2) guarantees that the capacity of the manufacturers is not surpassed. Constraint (3) guarantees that the customers demands are fulfilled. Constraint (4) is the flow conservation conditions and assures that the units collected by a given distribution center from manufacturers are equal to the units transported from that distribution center to the customers. Constraint (5) guarantees that the storage capacities of the distribution centers are not surpassed. Finally, the last constraint sets the ranges of the decision variables.
Hong et al. [10] described an illustrative example consisting of 2 manufacturing plants, 4 DCs, and 6 customers, whose characteristics are presented in Figure 2.
We solved this example using our proposed mixed integer programming formulation of the problem with CPLEX Version 12.7.0, and in Figure 3, we present the obtained optimal solution.
In Figure 3, we point out the capacities of the manufacturers, the storage capacities of the DCs, the demands of the customers, the transportation costs, and with red color, the fixed costs for opening the DCs and fixed costs associated with the selected transportation routes.
In order to obtain the optimal solution for this example, using our model by means of CPLEX only took 0.05 s and 30 iterations, in contrast to the model proposed by Hong et al. [10], which used the LINGO solver in order to find the optimal solution of a cost of 449,050 within 93 iterations, and their proposed ant colony approach, which provided an suboptimal solution cost of 476,138 within 0.24 s.

4. Description of the Novel Solution Approach

To solve TSSCNDP-FC, we propose a genetic algorithm, hybridized with a linear programming procedure. Genetic algorithms (GAs) were introduced by Holland [29] and are search metaheuristic techniques inspired by the Darwinian evolutionary theory based on the “survival of the fittest” concept. GAs have the capability to deliver “good” sub-optimal solutions within reasonable computational running times, making them very attractive for solving optimization problems characterized by a large feasible solution space. A genetic algorithm begins with a collection of feasible solutions, called the initial population, which are represented by chromosomes. Solutions from the current population are selected and employed to compose a new population. This is motivated by the belief that the newly created population will be better than the old one in terms of the quality of the solutions. The solutions that are picked to generate the offspring are selected according to their fitness; the more appropriate they are, the more opportunities they have to reproduce. This is repeated until some criteria (for example, the number of populations, improvement of the best solution, etc.) are fulfilled.
Our genetic algorithm builds different breeds of chromosomes that evolve separately from random populations, until evolution stagnates. Then, the breeds are merged together, hoping that the newly formed hybrid chromosomes will be better.
The operating principle of our hybrid genetic algorithm is shown in Figure 4, and the description of its blocks are described within this section.
The proposed GA makes use of an effective representation in which the genes of the chromosome represent the estimates of the flows corresponding to the optimal solution of the problem. Therefore, the chromosome contains two parts: the first one is an m × d matrix associated with the links from manufacturers to distribution centers, and the second one is a d × r matrix associated with the links from distribution centers to customers. We denote by x ˜ i j the gene corresponding to the link between manufacturer i and DC j and x ˜ j k the gene corresponding to the link between DC j and customer k. Initially, the chromosomes are generated randomly, with the only natural conditions that the estimates must not exceed the capacities of the manufactures, respectively the demands of the customers, i.e., x ˜ i j [ 0 , S i ] , for all i V 1 and j V 2 , and x ˜ j k [ 0 , D k ] , for all j V 2 and k V 3 .
It is unlikely that such a random chromosome would represent a correct estimate of a feasible solution of the two-stage supply chain network design problem with fixed costs for opening the DCs and associated with the routes. However, each chromosome has an associated feasible solution that can be effectively determined by solving the following linear programming problem, which is a simplified variant of the mathematical model of the TSSCNDP-FC:
min i = 1 m j = 1 d c ˜ i j x i j + j = 1 d k = 1 r c ˜ j k x j k s . t . ( 2 ) ( 7 )
where:
c ˜ i j = { c i j + f i j x ˜ i j + f j x ˜ j , if x ˜ i j 0 c if + f if , if x ˜ if = 0 and x ˜ j 0 c if + f if + f i , if x ˜ if = 0 and x ˜ j = 0
c ˜ i j = { c i j + f i j x ˜ i j + f j x ˜ j , if x ˜ i j 0 c if + f if , if x ˜ if = 0 and x ˜ j 0 c if + f if + f i , if x ˜ if = 0 and x ˜ j = 0
x ˜ j = i = 1 m x ˜ i j + k = 1 r x ˜ i j
The simplified variant of the mathematical model of the TSSCNDP-FC is actually a minimum cost flow problem for which there are well-known algorithms that solve it optimally in an efficient manner. We used the network simplex algorithm for solving this linear programming problem. The cost of the TSSCNDP-FC solution associated with the considered chromosome can be obtained using Relation (1) and the flows x i j and x j k determined by solving the simplified variant of the problem.
In order to boost the chances of discovering the optimal solution of the investigated supply chain problem with fixed costs, we developed a chromosome enhancement procedure shown in Algorithm 1, that processes all the chromosomes created throughout our genetic algorithm.
Algorithm 1: Procedure Chromosome enhancement
Mathematics 08 00712 i001
Each iteration of the loop in the chromosome enhancement procedure involves solving a simplified model of the TSSCNDP-FC, based on chromosome c (Step 3). The flows x i j and x j k determined in Step 3 are used in Step 4 to calculate the cost of the TSSCNDP-FC solution associated with chromosome c. If the solution has been improved, then the chromosome genes are updated in Step 8, using the streams determined in step 3. However, there is no guarantee that each iteration of the algorithm will enhance the TSSCNDP-FC solution associated with chromosome c. The algorithm stops when the solution worsens or a duplicate chromosome is reached. Finally, the genes of c are replaced with those of the last chromosome saved in Step 6. A chromosome c 2 is considered a duplicate of c 1 if the two chromosomes have the same corresponding TSSCNDP-FC solution.
For the newly created random chromosomes, the condition in Step 5 is changed to z ˜ z . This produces cleaner chromosomes in the initial population. Experiments showed that for the other chromosomes, it was better that the decision remained unchanged, because thus, there would be more diversity among the offspring.
In Figure 5, we illustrate the representation of an estimated chromosome associated with the example provided by Hong et al. [10], whose characteristics were presented in the previous section, where the entries of the matrix must not overcome the capacities of the manufacturers and the demands of the customers, and the representation of the chromosome after applying the chromosome enhancement procedure:
We remark that the chromosome C 1 correspond to an unfeasible solution, while the enhanced chromosome P 1 corresponds to a feasible solution of the problem.
For initializing the algorithm, an initial population of 2 × ( m × d + d × r ) enhanced chromosomes is generated by the random generator block; see Figure 4.
The admission block chooses the chromosomes that will form the current population (generation) from a pool of chromosomes that can be either from the current population or the offspring resulting from applying genetic operators (crossover and mutation). The size of each generation is kept constant at ( m × d + d × r ) / 2 chromosomes. The admission block applies the following rules, which were adjusted based on computational experiments:
  • The first two-thirds of the current population will be completed with the best chromosomes in the pool. At least half of these chromosomes must be newborn, i.e., not being part of the current population in the previous stages of evolution.
  • The other chromosomes in the current population are randomly chosen from the pool.
Selection is the phase of a GA in which individual chromosomes are selected from a population for later breeding. Our selection block uses the tournament selection strategy for choosing the two chromosomes that will be selected to undergo crossover in order to form an offspring. The number of participants for each tournament was randomly settled between two and 10.
Crossover is a genetic operator that combines the genetic information of two parents in order to achieve new offspring. Our crossover block groups the chromosomes supplied by the selection block two-by-two. Then, with each pair, it forms an offspring. Each gene of the offspring is taken with the same probability either from the first parent or from the second one. The way our proposed crossover operator works is illustrated in Figure 6.
The chromosomes C 1 and C 2 in Figure 6 are two random chromosomes associated with the illustrative example presented in the previous section. By applying the chromosome enhancement procedure, they are transformed into two new chromosomes denoted by P 1 and P 2 , which are feasible solutions of the TSSCNDP-FC. The offspring O 1 is obtained by applying the crossover operator to parents P 1 and P 2 , and finally, the enhanced offspring O 1 e is obtained by applying the enhancement procedure to O 1 . We can observe that the enhanced offspring O 1 e represents a better solution than the two feasible solutions corresponding to the parents P 1 and P 2 . In Figure 6, we represent with red color the genes inherited by the offspring O 1 from each of the parents P 1 and P 2 .
Mutation is a genetic operator whose main scope is to maintain the diversity of the chromosomes between consecutive generations of the GA. It is similar to biological mutation and usually happens with low probability. Our mutation operation was applied with a probability of 1% to each newborn chromosome. To accomplish this operation, a client k and a set of minimum 1 and maximum d DCs are randomly chosen. Then, all genes corresponding to the links between the k client and the chosen DCs are changed as follows: x ˜ j k r a n d o m v a l u e [ 1 , D k ] . Then, a DC j and a set of a minimum of one and a maximum of m manufacturers are randomly chosen. Then, all genes corresponding to the links between the d DCs and the chosen manufacturers are changed as follows: x ˜ i j r a n d o m v a l u e [ 1 , S i ] . The way the proposed mutation operator works is depicted in Figure 7.
The chromosome O 2 in Figure 7 represents the result of applying the mutation operator to the offspring O 1 presented in Figure 5. In this example, the mutation operator chooses the sixth client and a group of three DCS (1, 3, and 4), after which, it chooses the fourth DC and only one manufacturer (2). The enhanced offspring O 2 e is obtained by applying the enhancement procedure to O 2 . This actually represents the optimal solution of the considered TSSCNDP-FC.
The genes of the chromosomes involved in the crossover and mutation examples are presented in Figure 8. The red entries in the representations of P 1 and P 2 are the genes inherited by the offspring O 1 and in O 2 are the genes that are changed by applying the mutation operator.
The evolution of a chromosome population ends when the best chromosome has not been improved in the last 25 generations. If a breed pool is available, then it is merged with the current population, using hybrid selection. The hybrid selection is different from the selection operation by the fact that it processes chromosomes from two different pools: the current population and the breed pool. One of the parent chromosomes is selected from the current population, and the other one is selected from the breed pool. The merging completes when the desired number of newborn chromosomes has been added to the offspring pool or when the desired number of crossover operations has been performed.
When the merging operation finishes, the breed pool is replaced with chromosomes taken from the current population, the existing breed pool, and the offspring. This operation is performed by the hybrid admission block, which applies the same rules as the admission block.

5. Discussion

This section is dedicated to the achieved computational results with the aim of assessing the effectiveness of our developed approach for solving the two-stage supply chain network design problem with fixed costs associated with the transportation routes and for opening the DCs.
We performed our computational experiments for solving the two-stage supply chain problem with fixed costs associated with the routes and for opening the distribution centers on a set of 200 instances randomly generated with varying characteristics in the same way as Hong et al. [10] generated. Since the test instances used by Hong et al. [10] could not be obtained in the literature, we generated new instances similar to those in Hong et al. [10]: the first 150 test instances classified into three problem classes: smaller, which consisted of 2 manufacturing plants, 5 DCs, and 10 customers, medium, which consisted of 4 manufacturing plants, 8 DCs, and 15 customers, and large, which consisted of 6 manufacturing plants, 10 DCs, and 20 customers. In addition, we generated randomly 50 test instances of the larger size consisting of 8 manufacturing plants, 12 distribution centers, and 25 customers, respectively 10 manufacturing plants, 15 distribution centers, and 30 customers. All the instances used in our computational experiments are available in [30].
We coded our algorithm in Java 8, and for each instance, we carried out 10 independent trials, on a PC with Intel Core i5-4590 3.3GHz, 4GB RAM, and the Windows 10 Education 64 bit operating system. To solve the proposed mixed integer programming formulation of the problem, we used CPLEX Version 12.7.0.
Table 1, Table 2, Table 3, Table 4 and Table 5 summarize the computational experiments performed for solving the considered instances using CPLEX 12.7.0 and the proposed hybrid solution approach.
In Table 1, we present the results obtained by CPLEX and the proposed hybrid genetic algorithm for solving small instances of the problem. The experimental study implied running each instance ten times. The first column displays the number of the instance; the next two columns contain the value of the optimal solution Z o p t and the necessary computing time spent in solving the instances provided by CPLEX. The following six columns contain the results obtained by our hybrid genetic algorithm: the minimum and the maximum of the objective function achieved in the ten runs of each instance ( Z m i n and Z m a x ), the percentage gap defined as 100 × ( Z a v g . Z m i n ) / Z m i n , the minimum, maximum, and average computing times ( T m i n , T m a x , and T a v g . ) necessary for solving the instances. Z a v g . is the average of the objective function achieved in the ten runs of each instance. The last column provides the improvement gap of the time necessary to deliver the optimal solution by the hybrid genetic algorithm in comparison to CPLEX. The improvement time gap is calculated as follows: 100 × ( T C P L E X T a v g . ) / T C P L E X . The instances are ordered based on these improvements.
Analyzing the results presented in Table 1, we can remark that for all the considered small sized instances, CPLEX delivered the optimal solutions within a computational time ranging from 0.047 to 0.172 s. Our hybrid solution approach delivered as well the optimal solutions in all ten runs of each instance and the average computational time spent in solving the instances ranged from 0.000 to 0.009 s. The average improvement gap of the time spent to deliver the optimal solution by the hybrid genetic algorithm in comparison to CPLEX was greater than 95%.
In Table 2, we present the results obtained by CPLEX and the proposed hybrid genetic algorithm for solving medium instances of the problem. The columns of Table 2 are similar to those of Table 1, and we ordered the instances based on the improvement of the time necessary to deliver the optimal solution by the hybrid genetic algorithm in comparison to CPLEX.
Analyzing the results presented in Table 2, we can remark that for all the considered medium sized instances, CPLEX delivered the optimal solutions within a computational time ranging from 1.39 to 16.63 s. Our hybrid solution approach delivered as well the optimal solutions in all ten runs of each instance, i.e., Z m i n = Z m a x = Z o p t , and the average computational time spent in solving the instances ranged from 0.09 to 3.07 s. The average improvement gap of the time spent to deliver the optimal solution by the hybrid genetic algorithm in comparison to CPLEX was greater than 87%.
In Table 3, we present the results obtained by CPLEX and the proposed hybrid genetic algorithm for solving large instances of the problem. The first column displays the number of the instance, and the next two columns contain the value of the optimal solution Z C P L E X achieved by CPLEX and the necessary computing time spent in solving the instances provided by CPLEX when available, otherwise the solution determined by CPLEX when interrupting the run after 3600 s of computing time. The following six columns contain the results obtained by our hybrid genetic algorithm: the minimum and the maximum of the objective function achieved in the ten runs of each instance ( Z m i n and Z m a x ), the percentage gap defined as 100 × ( Z m a x Z m i n ) / Z m i n , the minimum, maximum, and average computing times necessary for solving the instances. The last column provides the improvements of the solutions delivered by the hybrid genetic algorithm in comparison to the solutions achieved by CPLEX. The instances are ordered based on these improvements. The improvement gap is calculated as follows: 100 × ( Z C P L E X Z m i n ) / Z C P L E X . For the instances presented in Table 3, the running time of the algorithm was limited to 200 s.
Analyzing the computational results reported in Table 3, we can remark that for the first 43 instances, the best solution provided by our hybrid genetic algorithm improved the solution delivered by CPLEX within 3600 s, and the average computational time spent in solving the instances ranged from 11.63 to 113.26 s. The average improvement gap of the solutions delivered by the hybrid genetic algorithm in comparison to the solutions achieved by CPLEX was 0.114%. In the case of the last seven instances, CPLEX delivered the optimal solutions within a computational time ranging from 2071.50 to 3595.30 s. Our hybrid solution approach delivered as well the optimal solutions in all ten runs of each instance, i.e., Z m i n = Z m a x = Z o p t , but with much less computational time effort, and the average computational time spent in solving the instances ranged from 1.76 to 31.05 s. The improvement time gap ranged from 98.52% to 99.95%, and the average was greater than 99%. For 29 out of 50 instances, our proposed approach did not achieve the same solution in all ten runs, but we could observe that the average percentage gap between the minimum and the maximum of the objective function achieved ranged from 0.01% and 0.25%, a fact that proved the stability of our proposed solution approach.
In Table 4, we report the results obtained by CPLEX and the developed hybrid genetic algorithm for solving larger instances of the problem, namely with 8 manufacturing plants, 12 distribution centers, and 25 customers. The columns of Table 4 are similar to those of Table 3, and the instances are ordered based on the achieved improvement gaps. For the instances presented in Table 4, the running time of the algorithm was limited to 500 s.
Analyzing the computational results reported in Table 4, we can remark that for the first 18 instances, the best solution provided by our hybrid genetic algorithm improved the solution delivered by CPLEX within 3600 s, and the average computational time spent in solving the corresponding instances ranged from 24.52 to 255.88 s. In the case of the last seven instances, the solutions delivered by CPLEX within 3600 s were the same as the solutions achieved by our hybrid solution approach in all ten runs of each instance, but with much less computational time effort, and the average computational time spent in solving the instances ranged from 19.81 to 260.02 s. For 12 out of 25 instances, our proposed approach did not achieve the same solution in all ten runs, but we could observe that the average percentage gap between the minimum and the maximum of the objective function achieved ranged from 0.003% and 0.11%, a fact that showed the stability of our proposed solution approach.
In Table 5, we report the results obtained by CPLEX and the proposed hybrid genetic algorithm for solving larger instances of the problem, namely with 10 manufacturing plants, 15 distribution centers, and 30 customers. The columns of Table 5 are similar to those of Table 3, and the instances are ordered based on the achieved improvement gaps. For the instances presented in Table 5, the running time of the algorithm was limited to 800 s.
Analyzing the computational results reported in Table 5, we can remark that for all 25 instances, the best solution provided by our hybrid genetic algorithm improved the solution delivered by CPLEX within 3600 s, and the average computational time spent in solving the corresponding instances ranged from 22.9 to 484.6 s. In the case of these instances, our proposed approach did not achieve the same solution in all ten runs, but we can observe that the average percentage gap between the minimum and the maximum of the objective function achieved ranged from 0.001% and 0.3%, a fact that proved the stability of our proposed solution approach.
Overall, the results achieved by our hybrid genetic algorithm can be summarized as follows:
  • In the case of small and medium sized instances, our algorithm delivered the optimal solution in all ten runs of each instance, but with much less computational time effort in comparison to CPLEX.
  • As regards the 50 large sized instances, for seven of them, our hybrid algorithm delivered the optimal solution in all ten runs of each instance, but with much less computational time effort in comparison to CPLEX, and for the remaining instances, the best solution achieved by our algorithm improved the solution provided by CPLEX within 3600 s, but with much less computational time effort in comparison to CPLEX.
  • In the case of the proposed 50 larger instances, the best solution provided by our hybrid genetic algorithm improved the solution delivered by CPLEX within 3600 s, and the average computational times spent in solving the corresponding instances were lower compared to CPLEX.
  • We can remark that our developed solution outperformed in terms of the quality of the solutions and of computational times the ACO-based heuristic approach proposed by Hong et al. [10], which according to the previously mentioned authors, provided sub-optimal solutions with a gap of about 10 % on average from the optimal solutions.

6. Conclusions

In this paper, we investigated the two-stage supply chain network design problem in which we considered two types of fixed costs: ones for opening the DCs and the others associated with the routes. This optimization problem was introduced recently by Hong et al. [10], and it generalized the previously considered two-stage transportation problems.
The main goal of our paper was to describe an effective hybrid genetic algorithm for solving the two-stage transportation network design problem with fixed costs. Our hybrid algorithm was achieved by incorporating a linear programming optimization problem within the framework of a genetic algorithm. The method we proposed had certain important features, and here are some of these: the employment of an efficient representation in which the chromosome was generated in two stages: from an estimation of the flows to a feasible solution of the problem, the building of different breeds of chromosomes that evolved separately from random populations, until evolution stagnates, then the breeds were merged together, hoping that the newly formed hybrid chromosomes would be better.
We evaluated the effectiveness of the proposed solution approach on a set of 200 instances classified into four classes: small, medium, large, and larger sized instances. The first 150 instances were generated in the same way as Hong et al. [10] suggested, and the others were larger sized instances. The computational results that we obtained proved the efficiency of our developed hybrid genetic algorithm in yielding high quality solutions within reasonable running times, besides its superiority as compared to CPLEX and the other existing solution approaches from the literature.
It is our intention to continue our research and further improve the developed hybrid genetic algorithm through a combination with the local search methods. We intend to carry out a parallel implementation of the proposed algorithm, in order to fully benefit from the processing power of the new multicore processors and to generalize the TSSCNDP-FC model so that customers can be supplied directly by manufacturers and so that the distribution centers can also supply other distribution centers. We also plan to assess the generality and scalability of the solution approach that we suggested by testing it on even larger instances.

Author Contributions

Conceptualization, P.C.P.; methodology, O.C.; software, O.C.; validation, O.C. and C.S.; formal analysis, P.C.P.; investigation, O.C. and C.S.; resources, P.C.P and C.S.; data curation, O.C.; writing, original draft preparation, O.C.; writing, review and editing, O.C. and P.C.P.; visualization, O.C. and P.C.P.; All authors read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This work was partially supported by the project 21 PFEin the frame of the program PDI-PFE-CDI2018 and by the post-doctoral program “Entrepreneurial competences and excellence research in doctoral and postdoctoral programs—ANTREDOC”, a project co-funded by the European Social Fund.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Govindan, K.; Fattahi, M.; Keyvanshokooh, E. Supply chain network design under uncertainty: A comprehensive review and future research directions. Eur. J. Oper. Res. 2017, 263, 108–141. [Google Scholar] [CrossRef]
  2. Klibi, W.; Martel, A.; Guitouni, A. The design of robust value-creating supply chain networks: A critical review. Eur. J. Oper. Res. 2010, 203, 283–293. [Google Scholar] [CrossRef]
  3. Melo, M.T.; Nickel, S.; Da Gamma, F.S. Facility location and supply chain management—A review. Eur. J. Oper. Res. 2009, 196, 401–412. [Google Scholar] [CrossRef]
  4. Wang, H.S. A two-phase ant colony algorithm for multi-echelon defective supply chain netwok design. Eur. J. Oper. Res. 2009, 192, 243–252. [Google Scholar] [CrossRef]
  5. Dotoli, M.; Fanti, M.P.; Meloni, C.; Zhou, M.C. A Multi-Level Approach for Network Design of Integrated Supply Chains. Int. J. Prod. Res. 2005, 43, 4267–4287. [Google Scholar] [CrossRef]
  6. Dotoli, M.; Fanti, M.P.; Meloni, C.; Zhou, M.C. A decision support system for the supply chain configuration. IEEE Int. Conf. Syst. Man Cybern. 2003, 3, 2667–2672. [Google Scholar]
  7. Balinski, M.I. Fixed-cost transportation problems. Nav. Res. Logist. 1961, 8, 41–54. [Google Scholar] [CrossRef]
  8. Guisewite, G.; Pardalos, P. Minimum concave-cost network flow problems: Applications, complexity, and algorithms. Ann. Oper. Res. 1990, 25, 75–99. [Google Scholar] [CrossRef]
  9. Buson, E.; Roberti, R.; Toth, P. A reduced-cost iterated local search heuristic for the fixed-charge transportation problem. Oper. Res. 2014, 62, 1095–1106. [Google Scholar] [CrossRef] [Green Version]
  10. Hong, J.; Diabat, A.; Panicker, V.V.; Rajagopalan, S. A two-stage supply chain problem with fixed costs: An ant colony optimization approach. Int. J. Prod. Econ. 2018, 204, 214–226. [Google Scholar] [CrossRef]
  11. Sabo, C.; Horvat-Marc, A.; Pop, P.C. Comments on “A two-stage supply chain problem with fixed costs: An ant colony optimization approach” by Hong et al. International Journal of Production Economics (2018). Creat. Math. Inform. 2019, 28, 183–189. [Google Scholar]
  12. Raj, K.A.A.D.; Rajendran, C. A genetic algorithm for solving the fixed-charge transportation model: Two-stage problem. Comput. Oper. Res. 2012, 39, 2016–2032. [Google Scholar]
  13. Jawahar, N.; Balaji, A.N. A genetic algorithm for the two-stage supply chain distribution problem associated with a fixed charge. Eur. J. Oper. Res. 2009, 194, 496–537. [Google Scholar] [CrossRef]
  14. Pop, P.C.; Sabo, C.; Biesinger, B.; Hu, B.; Raidl, G. Solving the Two-Stage Fixed-Charge Transportation Problem with a Hybrid Genetic Algorithm. Carpathian J. Math. 2017, 33, 365–371. [Google Scholar]
  15. Cosma, O.; Pop, P.C.; Pop Sitar, C. An efficient iterated local search heuristic algorithm for the two-stage fixed-charge transportation problem. Carpathian J. Math. 2019, 35, 153–164. [Google Scholar]
  16. Gen, M.; Altiparmak, F.; Lin, L. A genetic algorithm for two-stage transportation problem using priority based encoding. OR Spectr. 2006, 28, 337–354. [Google Scholar] [CrossRef]
  17. Calvete, H.; Gale, C.; Iranzo, J. An improved evolutionary algorithm for the two-stage transportation problem with fixed charge at depots. OR Spectr. 2016, 38, 189–206. [Google Scholar] [CrossRef]
  18. Cosma, O.; Dănciulescu, D.; Pop, P.C. On the two-stage transportation problem with fixed charge for opening the distribution centers. IEEE Access 2019, 79, 113684–113698. [Google Scholar] [CrossRef]
  19. Cosma, O.; Pop, P.C.; Dănciulescu, D. A novel matheuristic approach for a two-stage transportation problem with fixed costs associated with the routes. Comput. Oper. Res. 2020, 118, 104906. [Google Scholar] [CrossRef]
  20. Molla-Alizadeh-Zavardehi, S.; Hajiaghaei-Kesteli, M.; Tavakkoli-Moghaddam, R. Solving a capacitated fixed-cost transportation problem by artificial immune and genetic algorithms with a Prüfer number representation. Expert Syst. Appl. 2011, 38, 10462–10474. [Google Scholar] [CrossRef]
  21. El-Sherbiny, M.M. Comments on “Solving a capacitated fixed-cost transportation problem by artificial immune and genetic algorithms with a Prüfer number representation" by Molla-Alizadeh-Zavardehi, S. et al. Expert Syst. Appl. Expert Syst. Appl. 2012, 39, 11321–11322. [Google Scholar] [CrossRef]
  22. Pintea, C.-M.; Pop Sitar, C.; Hajdu-Macelaru, M.; Pop, P.C. A Hybrid Classical Approach to a Fixed-Charge Transportation Problem. In Proceedings of HAIS 2012, Part I Lecture Notes in Computer Science; Corchado, E., Ed.; Springer: Salamanca, Spain, 2012; Volume 7208, pp. 557–566. [Google Scholar]
  23. Pintea, C.M.; Pop, P.C. An improved hybrid algorithm for capacitated fixed-charge transportation problem. Log. J. IJPL 2015, 23, 369–378. [Google Scholar] [CrossRef]
  24. Pop, P.C.; Matei, O.; Pop Sitar, C.; Zelina, I. A hybrid based genetic algorithm for solving a capacitated fixed-charge transportation problem. Carpathian J. Math. 2016, 32, 225–232. [Google Scholar]
  25. Cosma, O.; Pop, P.C.; Matei, O.; Zelina, I. A hybrid iterated local search for solving a particular two-stage fixed-charge transportation problem. In Proceedings of HAIS 2018, Lecture Notes in Computer Science; Springer: Oviedo, Spain, 2018; Volume 10870, pp. 684–693. [Google Scholar]
  26. Santibanez-Gonzalez, E.D.R.; Robson Mateus, G.; Pacca Luna, H. Solving a public sector sustainable supply chain problem: A Genetic Algorithm approach. In Proceedings of the International Conference on Artificial Intelligence (ICAI), Las Vegas, NV, USA, 18–21 July 2011; pp. 507–512. [Google Scholar]
  27. Pintea, C.-M.; Pop, P.C.; Hajdu-Măcelaru, M. Classical Hybrid Approaches on a Transportation Problem with Gas Emissions Constraints. Adv. Intell. Soft Comput. 2013, 188, 449–458. [Google Scholar]
  28. Pop, P.C.; Pintea, C.-M.; Pop Sitar, C.; Hajdu-Macelaru, M. An efficient Reverse Distribution System for solving sustainable supply chain network design problem. J. Appl. Log. 2015, 13, 105–113. [Google Scholar] [CrossRef]
  29. Holland, J.H. Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology. In Control and Artificial Intelligence; MIT Press: Cambridge, MA, USA, 1992. [Google Scholar]
  30. Two Stage Transportation Problem Instances. Available online: https://sites.google.com/view/tstp-instances/ (accessed on 26 April 2020).
Figure 1. Illustration of the two-stage supply chain network design problem with fixed costs.
Figure 1. Illustration of the two-stage supply chain network design problem with fixed costs.
Mathematics 08 00712 g001
Figure 2. The characteristics of the example described by Hong et al. [10].
Figure 2. The characteristics of the example described by Hong et al. [10].
Mathematics 08 00712 g002
Figure 3. An illustrative example of the of the two-stage supply chain network design problem with fixed costs and the obtained optimal solution.
Figure 3. An illustrative example of the of the two-stage supply chain network design problem with fixed costs and the obtained optimal solution.
Mathematics 08 00712 g003
Figure 4. The operation principle of our proposed hybrid genetic algorithm for solving the two-stage supply chain network design problem with fixed costs.
Figure 4. The operation principle of our proposed hybrid genetic algorithm for solving the two-stage supply chain network design problem with fixed costs.
Mathematics 08 00712 g004
Figure 5. The representation of the estimated chromosome and the enhanced chromosome.
Figure 5. The representation of the estimated chromosome and the enhanced chromosome.
Mathematics 08 00712 g005
Figure 6. The crossover operator.
Figure 6. The crossover operator.
Mathematics 08 00712 g006
Figure 7. The mutation operator.
Figure 7. The mutation operator.
Mathematics 08 00712 g007
Figure 8. The genes of the chromosomes involved in the mutation and crossover examples.
Figure 8. The genes of the chromosomes involved in the mutation and crossover examples.
Mathematics 08 00712 g008
Table 1. Computational results achieved in the case of small sized instances.
Table 1. Computational results achieved in the case of small sized instances.
InstanceCPLEXOur Proposed Hybrid Genetic AlgorithmImprovement
Z opt T cplex Z min Z max % gap T min T max T avg . Time (%)
1.150,7870.047150,787150,7870.000.0000.0000.000100.0
2.156,6540.047156,654156,6540.000.0000.0000.000100.0
3.139,8110.047139,811139,8110.000.0000.0000.000100.0
4.128,8590.047128,859128,8590.000.0000.0000.000100.0
5.116,6370.093116,637116,6370.000.0000.0150.00298.4
6.87,6940.07887,69487,6940.000.0000.0160.00297.9
7.105,4200.063105,420105,4200.000.0000.0150.00297.6
8.120,0770.063120,077120,0770.000.0000.0150.00297.6
9.117,5900.063117,590117,5900.000.0000.0160.00297.5
10.131,2330.063131,233131,2330.000.0000.0160.00297.5
11.121,4110.062121,411121,4110.000.0000.0160.00297.4
12.149,5250.109149,525149,5250.000.0000.0160.00397.1
13.134,5910.156134,591134,5910.000.0000.0160.00597.0
14.105,0470.047105,047105,0470.000.0000.0150.00296.8
15.127,1840.047127,184127,1840.000.0000.0150.00296.8
16.122,6300.094122,630122,6300.000.0000.0160.00396.7
17.101,7220.094101,722101,7220.000.0000.0160.00396.7
18.132,5930.047132,593132,5930.000.0000.0160.00296.6
19.113,9760.047113,976113,9760.000.0000.0160.00296.6
20.129,4360.172129,436129,4360.000.0000.0160.00696.4
21.146,0780.125146,078146,0780.000.0000.0160.00596.2
22.119,6270.079119,627119,6270.000.0000.0160.00396.1
23.95,8450.07895,84595,8450.000.0000.0160.00396.0
24.104,1900.078104,190104,1900.000.0000.0160.00395.9
25.127,8690.140127,869127,8690.000.0000.0160.00695.6
26.100,4510.171100,451100,4510.000.0000.0160.00895.4
27.128,2530.094128,253128,2530.000.0000.0160.00595.1
28.150,3790.125150,379150,3790.000.0000.0160.00695.0
29.108,9420.093108,942108,9420.000.0000.0160.00594.9
30.137,9550.062137,955137,9550.000.0000.0160.00394.8
31.170,1190.109170,119170,1190.000.0000.0160.00694.3
32.126,6030.079126,603126,6030.000.0000.0160.00594.2
33.134,6530.078134,653134,6530.000.0000.0160.00594.1
34.142,5590.078142,559142,5590.000.0000.0160.00594.0
35.126,6110.125126,611126,6110.000.0000.0320.00893.8
36.116,8690.047116,869116,8690.000.0000.0150.00393.6
37.108,2010.094108,201108,2010.000.0000.0160.00693.4
38.129,0050.047129,005129,0050.000.0000.0160.00393.2
39.127,9970.110127,997127,9970.000.0000.0160.00892.9
40.104,8910.109104,891104,8910.000.0000.0470.00892.8
41.111,1040.063111,104111,1040.000.0000.0160.00592.7
42.108,0010.062108,001108,0010.000.0000.0160.00592.6
43.143,7000.063143,700143,7000.000.0000.0160.00592.5
44.138,5100.078138,510138,5100.000.0000.0160.00692.1
45.135,1010.079135,101135,1010.000.0000.0160.00692.0
46.129,8540.078129,854129,8540.000.0000.0160.00691.9
47.116,6700.094116,670116,6700.000.0000.0160.00891.8
48.141,4670.047141,467141,4670.000.0000.0160.00590.2
49.141,3060.094141,306141,3060.000.0000.0310.00990.1
50.125,5500.062125,550125,5500.000.0000.0160.00690.0
Table 2. Computational results achieved in the case of medium sized instances.
Table 2. Computational results achieved in the case of medium sized instances.
InstanceCPLEXOur Proposed Hybrid Genetic AlgorithmImprovement
Z opt T cplex Z min Z max % gap T min T max T avg . Time (%)
1.544,06211.86544,062544,0620.000.050.420.1398.90
2.552,7933.44552,793552,7930.000.050.130.0997.43
3.581,8252.25581,825581,8250.000.050.110.0995.93
4.558,0473.66558,047558,0470.000.050.620.1795.43
5.524,3466.53524,346524,3460.000.060.510.3095.33
6.536,3443.00536,344536,3440.000.050.650.1595.08
7.496,75311.83496,753496,7530.000.081.930.6194.85
8.632,8675.94632,867632,8670.000.061.290.3494.20
9.569,3514.38569,351569,3510.000.050.760.3093.07
10.584,3862.47584,386584,3860.000.050.610.1992.14
11.479,2043.27479,204479,2040.000.050.810.2692.11
12.573,4881.39573,488573,4880.000.060.420.1191.80
13.496,39912.06496,399496,3990.000.252.311.0391.44
14.507,9624.20507,962507,9620.000.081.940.3890.99
15.451,7123.55451,712451,7120.000.090.830.3390.80
16.439,67710.66439,677439,6770.000.053.890.9990.68
17.527,75616.63527,756527,7560.000.134.881.6590.09
18.512,0293.88512,029512,0290.000.111.110.3990.04
19.549,3111.63549,311549,3110.000.050.470.1689.91
20.495,9126.81495,912495,9120.000.081.940.6989.83
21.569,7623.67569,762569,7620.000.091.330.3889.78
22.546,12711.19546,127546,1270.000.063.521.1989.40
23.650,9472.94650,947650,9470.000.060.820.3189.31
24.625,2852.48625,285625,2850.000.060.540.2789.08
25.487,6652.22487,665487,6650.000.060.620.2688.29
26.465,1591.61465,159465,1590.000.060.640.1988.25
27.531,8913.72531,891531,8910.000.091.730.4488.22
28.525,3361.86525,336525,3360.000.060.510.2288.21
29.523,0617.56523,061523,0610.000.062.190.9088.05
30.549,1133.06549,113549,1130.000.081.100.3887.45
31.520,6655.02520,665520,6650.000.081.230.6387.43
32.516,2356.03516,235516,2350.000.421.250.7787.16
33.539,9782.02539,978539,9780.000.051.020.2985.55
34.476,7252.25476,725476,7250.000.060.560.3485.10
35.533,0053.33533,005533,0050.000.081.210.5085.06
36.572,4112.25572,411572,4110.000.051.120.3484.81
37.511,8751.69511,875511,8750.000.060.820.2982.95
38.484,0942.55484,094484,0940.000.061.550.4482.67
39.463,1492.83463,149463,1490.000.082.200.5381.21
40.610,3282.36610,328610,3280.000.060.930.4580.93
41.488,71515.66488,715488,7150.000.117.863.0780.38
42.524,4442.58524,444524,4440.000.031.390.5280.02
43.546,9014.02546,901546,9010.000.451.550.8179.73
44.546,7892.02546,789546,7890.000.092.030.4279.16
45.541,5213.99541,521541,5210.000.082.560.8478.90
46.458,5574.28458,557458,5570.000.082.900.9278.44
47.576,0662.70576,066576,0660.000.082.250.6077.68
48.566,7295.08566,729566,7290.000.162.761.1477.55
49.498,8913.19498,891498,8910.000.082.130.7377.12
50.541,0122.36541,012541,0120.000.092.320.5875.47
Table 3. Computational results achieved in the case of large sized instances.
Table 3. Computational results achieved in the case of large sized instances.
InstanceCPLEXOur Proposed Hybrid Genetic AlgorithmImprovement
Z CPLEX T CPLEX Z min Z max % gap T min T max T avg . Gap (%)
1.1,498,462>36001,492,9591,498,4620.047.58179.7270.540.367
2.1,526,157>36001,521,5521,525,8640.2510.5671.1434.710.302
3.1,229,762>36001,226,0591,228,9000.160.7353.5812.970.301
4.1,359,194>36001,355,7591,358,5080.122.01125.8650.130.253
5.1,506,789>36001,503,2531,506,4260.080.69171.2954.240.235
6.1,360,560>36001,357,4061,360,8710.210.73195.6276.330.232
7.1,373,481>36001,370,3051,370,3050.000.33118.0632.240.231
8.1,507,877>36001,505,0231,507,8770.1415.42190.9886.360.189
9.1,313,779>36001,311,3411,313,7790.073.95187.4856.390.186
10.1,404,192>36001,401,6471,401,6470.002.59114.9931.510.181
11.1,439,014>36001,436,5111,436,5110.000.5929.3612.940.174
12.1,241,920>36001,239,8411,242,5540.024.4548.0925.270.167
13.1,182,028>36001,180,0551,182,3590.085.36185.8167.030.167
14.1,400,016>36001,397,7491,399,6950.090.41163.1748.320.162
15.1,452,499>36001,450,2241,453,1310.040.58199.5479.300.157
16.1,210,260>36001,208,8681,210,2600.090.42150.6933.190.115
17.1,172,705>36001,171,3921,173,0700.0611.33192.2992.860.112
18.1,431,188>36001,429,6511,429,6510.002.1640.6111.630.107
19.1,391,275>36001,389,9981,389,9980.008.58175.24113.260.092
20.1,354,598>36001,353,3771,355,2220.031.24192.8869.230.090
21.1,401,338>36001,400,0901,400,0900.004.8097.9734.900.089
22.1,364,919>36001,363,7491,366,0090.133.25167.7647.420.086
23.1,425,663>36001,424,5981,425,6630.040.5971.5620.010.075
24.1,164,525>36001,163,6891,164,5250.040.41171.1945.030.072
25.1,485,811>36001,484,7451,485,5950.021.1791.4337.720.072
26.1,187,478>36001,186,6401,186,6400.005.47155.6659.620.071
27.1,493,211>36001,492,2361,493,5110.0128.94191.48106.280.065
28.1,348,937>36001,348,1401,348,9370.020.54197.8868.950.059
29.1,306,728>36001,305,9671,305,9670.005.26167.7248.310.058
30.1,362,966>36001,362,2301,362,2300.006.8143.3116.370.054
31.1,432,828>36001,432,0601,432,0600.0014.73163.4558.880.054
32.1,294,145>36001,293,4701,294,1450.017.53160.4167.000.052
33.1,367,135>36001,366,4941,367,1350.038.97185.4358.010.047
34.1,345,140>36001,344,7221,345,2870.031.63182.4375.990.031
35.1,444,708>36001,444,2831,444,7080.012.39125.7533.110.029
36.1,321,196>36001,320,8201,320,8200.006.27186.8568.680.028
37.1,488,631>36001,488,2591,488,6310.0218.45196.4199.940.025
38.1,360,169>36001,359,8441,360,5780.027.05174.5681.020.024
39.1,459,280>36001,458,9801,458,9800.004.5762.3228.490.021
40.1,320,505>36001,320,2571,320,8880.015.3481.0952.970.019
41.1,515,575>36001,515,3011,515,3010.000.3929.8712.100.018
42.1,456,015>36001,455,7821,456,0150.012.82167.4745.340.016
43.1,558,119>36001,557,9801,557,9800.000.3964.5517.190.009
44.1,299,2992071.501,299,2991,299,2990.003.8381.3924.670.000
45.1,247,1632096.721,247,1631,247,1630.001.32120.0131.050.000
46.1,427,5022686.281,427,5021,427,5020.002.2194.7925.360.000
47.1,393,4123595.301,393,4121,393,4120.001.8512.646.860.000
48.1,191,0853507.411,191,0851,191,0850.000.275.601.760.000
49.1,391,9083344.051,391,9081,391,9080.000.3016.084.500.000
50.1,285,9003348.831,285,9001,285,9000.000.6114.584.940.000
Table 4. Computational results achieved in the case of larger sized instances (8-12-25).
Table 4. Computational results achieved in the case of larger sized instances (8-12-25).
InstanceCPLEXOur Proposed Hybrid Genetic AlgorithmImprovement
Z CPLEX T CPLEX Z min Z max % gap T min T max T avg . Gap (%)
1.1,630,535>36001,622,9811,622,9810.002.94214.0564.900.46
2.1,585,431>36001,581,5791,585,2610.1128.37460.47187.290.24
3.1,564,349>36001,560,9781,564,1960.0219.81247.37129.790.22
4.1,688,457>36001,686,0191,688,4570.0216.48444.18255.880.14
5.1,795,971>36001,793,4741,795,9710.0811.07420.04210.970.14
6.1,693,447>36001,691,2031,694,0660.06136.12457.29254.000.13
7.1,468,196>36001,466,2641,466,2640.006.2149.0724.520.13
8.1,670,941>36001,668,9431,669,4560.0036.49294.80122.200.12
9.1,543,092>36001,541,2801,541,2800.001.91314.22113.650.12
10.1,850,860>36001,848,7321,848,7320.0025.34295.75159.960.11
11.1,710,500>36001,709,0781,710,5000.071.74327.8447.630.08
12.1,616,695>36001,615,4171,615,4170.0013.73494.86122.640.08
13.1,660,980>36001,659,8831,660,9800.0322.75433.65147.030.07
14.1,631,768>36001,631,2441,631,2440.0014.3399.5038.020.03
15.1,658,811>36001,658,5121,658,8110.0041.30431.92119.950.02
16.1,583,688>36001,583,4491,583,6880.016.33289.25119.270.02
17.1,570,940>36001,570,7291,570,9400.014.28429.81196.770.01
18.1,778,786>36001,778,6471779,6420.0120.73430.86166.330.01
19.1,671,547>36001,671,5471,671,5470.008.60412.75135.360.00
20.1,553,874>36001,553,8741,553,8740.0016.71404.55140.020.00
21.1,618,396>36001,618,3961,618,3960.002.03120.9425.570.00
22.1,305,978>36001,305,9781,305,9780.0033.56463.87260.020.00
23.1,682,622>36001,682,6221,682,6220.0016.56117.5659.010.00
24.1,738,929>36001,738,9291,738,9290.008.2743.0319.810.00
25.1,539,218>36001,539,2181,539,2180.0073.29398.56215.220.00
Table 5. Computational results achieved in the case of larger sized instances (10-15-30).
Table 5. Computational results achieved in the case of larger sized instances (10-15-30).
InstanceCPLEXOur Proposed Hybrid Genetic AlgorithmImprovement
Z CPLEX T CPLEX Z min Z max % gap T min T max T avg . Gap (%)
1.2,043,349>36002,031,6522,036,0910.0561.1742.1311.50.572
2.2,054,648>36002,046,8192,052,2300.09153.3800.4479.50.381
3.2,081,647>36002,073,8632,076,0820.0762.5781.5389.50.374
4.1,928,963>36001,923,0431,926,3610.0867.5782.5463.10.307
5.1,983,427>36001,977,5141,986,5230.30116.4707.7418.80.298
6.1,904,794>36001,899,3911,905,8470.18186.9797.8471.10.284
7.2,083,134>36002,077,2552,077,6770.0121.6769.1374.00.282
8.2,054,765>36002,049,0412,050,4520.02106.7763.722.90.279
9.1,879,761>36001,874,8751,879,5500.1153.8756.7381.20.260
10.2,115,385>36002,110,3772,112,0990.0478.4601.4408.80.237
11.2,160,852>36002,156,1132,160,8120.0828.2732.3244.90.219
12.2,151,335>36002,147,1652,151,6240.09162.8696.6400.60.194
13.1,893,790>36001,890,6821,894,1710.1241.7699.4286.00.164
14.2,116,378>36002,113,2122,118,7290.1859.4788.1359.10.150
15.2,211,552>36002,208,6642,211,9250.1263.2754.5300.90.131
16.2,234,741>36002,232,1172,234,7410.0745.5685.2338.90.117
17.2,073,949>36002,071,5232,073,5720.0543.3738.5347.00.117
18.2,243,637>36002,241,3362,243,6370.0514.4767.2426.00.103
19.2,176,878>36002,174,8432,177,4320.04116.7666.9345.80.093
20.2,112,460>36002,110,6532,113,9820.07170.9648.1383.60.086
21.1,783,878>36001,782,4391,782,6930.0014.3343.7130.60.081
22.1,910,476>36001,909,2061,913,7530.0516.8647.0279.00.066
23.1,778,394>36001,777,2351,778,8340.05148.9695.6484.60.065
24.2,101,223>36002,099,9762,103,0720.0624.7764.0415.60.059
25.2,203,959>36002,202,6762,204,8920.06112.8709.2334.30.058

Share and Cite

MDPI and ACS Style

Cosma, O.; Pop, P.C.; Sabo, C. An Efficient Hybrid Genetic Approach for Solving the Two-Stage Supply Chain Network Design Problem with Fixed Costs. Mathematics 2020, 8, 712. https://doi.org/10.3390/math8050712

AMA Style

Cosma O, Pop PC, Sabo C. An Efficient Hybrid Genetic Approach for Solving the Two-Stage Supply Chain Network Design Problem with Fixed Costs. Mathematics. 2020; 8(5):712. https://doi.org/10.3390/math8050712

Chicago/Turabian Style

Cosma, Ovidiu, Petrică C. Pop, and Cosmin Sabo. 2020. "An Efficient Hybrid Genetic Approach for Solving the Two-Stage Supply Chain Network Design Problem with Fixed Costs" Mathematics 8, no. 5: 712. https://doi.org/10.3390/math8050712

APA Style

Cosma, O., Pop, P. C., & Sabo, C. (2020). An Efficient Hybrid Genetic Approach for Solving the Two-Stage Supply Chain Network Design Problem with Fixed Costs. Mathematics, 8(5), 712. https://doi.org/10.3390/math8050712

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