Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2015 / Article
Special Issue

Recent Advances on Modeling, Control, and Optimization for Complex Engineering Systems

View this Special Issue

Research Article | Open Access

Volume 2015 |Article ID 387082 |

Jinwei Gu, Manzhan Gu, Xingsheng Gu, "A Mutualism Quantum Genetic Algorithm to Optimize the Flow Shop Scheduling with Pickup and Delivery Considerations", Mathematical Problems in Engineering, vol. 2015, Article ID 387082, 17 pages, 2015.

A Mutualism Quantum Genetic Algorithm to Optimize the Flow Shop Scheduling with Pickup and Delivery Considerations

Academic Editor: Kang Li
Received05 Jun 2014
Revised03 Nov 2014
Accepted11 Nov 2014
Published19 Mar 2015


A mutualism quantum genetic algorithm (MQGA) is proposed for an integrated supply chain scheduling with the materials pickup, flow shop scheduling, and the finished products delivery. The objective is to minimize the makespan, that is, the arrival time of the last finished product to the customer. In MQGA, a new symbiosis strategy named mutualism is proposed to adjust the size of each population dynamically by regarding the mutual influence relation of the two subpopulations. A hybrid Q-bit coding method and a local speeding-up method are designed to increase the diversity of genes, and a checking routine is carried out to ensure the feasibility of each solution; that is, the total physical space of each delivery batch could not exceed the capacity of the vehicle. Compared with the modified genetic algorithm (MGA) and the quantum-inspired genetic algorithm (QGA), the effectiveness and efficiency of the MQGA are validated by numerical experiments.

1. Introduction

The coordination of logistics activities in a supply chain has received a lot of attention recently. From the manufacturer’s point of view, the important problem in the supply chain scheduling is the coordination of the three stages including material supply, production scheduling, and product delivery. Traditionally, research on scheduling generally focuses on the models with various machine setting, job characteristics, and performance measures [1], and the transportation arrangements of the materials and the finished products are normally ignored in these models.

In the past decades, many researchers studied the machine scheduling problems with transportation under consideration. Some of the results consider the transportation of the jobs between machines in the flow shop model, which incorporates transport times when the jobs are transferred from one machine to another. Maggu et al. [24] firstly considered the scheduling problem with transportation between machines. In addition, Langston [5] designed algorithms for the problem of planning and coordinating movement within a deterministic flow shop system and analyzed the worst-case performances of the algorithms. Haouari and Ladhari [6] presented an effective branch and bound algorithm for the model proposed by Langston [5]. Hurink and Knust [7] considered the flow shop problems with transportation and a single robot with the objective of minimizing the makespan. They derived complexity results for the problems. Hurink and Knust [8] developed a tabu search algorithm for the problem with time windows. Hurink and Knust [9] extended some of the results in [8] to the job-shop scheduling models. For the model with a single robot, Lee and Strusevich [10] studied the two machine flow shop and open shop problems. They presented a best possible approximation algorithm for each of the two problems with some constraints. Naderi et al. [11] proposed simulated annealing algorithms for flow shop problems with the objective of minimizing the total weighted tardiness and makespan. For more new results on the scheduling problems with transportation between machines, please see [1215], among others.

Another type of scheduling model with transportation focuses on the delivery of finished jobs to customers. Lee and Chen [16] considered some scheduling models that incorporate the delivery decisions of the finished jobs. They researched the computational complexity of some problems and proposed polynomial or pseudopolynomial algorithms for them. Chang and Lee [17] extended one of the models in [16] to the environment in which each job has an individual amount of space. Li et al. [18] studied a single-machine scheduling problem with routing decisions of a vehicle that serves customers at different addresses. In addition, Chen and Vairaktarakis [19] researched eight problems with production and distribution under consideration and either designed a polynomial time algorithm or proved the NP-hardness for each of the problems. Geismar et al. [20] extended the Chen and Vairaktarakis model to the problem with a short shelf life and no inventory in the process. In addition, Soukhal et al. [21] investigated the scheduling models with constraints on both transportation and buffer capacities to minimize the makespan. They proved new complexity results for special cases of the problems considered. We refer to, for example, [2225], for more recent results on the model with delivery operations for finished products.

This paper considers the scheduling model that integrates the pickup of materials, flow shop scheduling, and the delivery of finished jobs. In this model, the material warehouse, the factory, and the customer are located at three different places. There are two vehicles (namely, conveyor and truck) each with a limited capacity. One vehicle (conveyor) travels between the factory and the warehouse for material transportation, and the other vehicle (truck) travels between the customer and the factory for finished products delivery. This model applies to many situations in supply chain business activities. For example, a shoes manufacturer purchases materials from India and arranges the production in China. Finally, the finished products are delivered to the USA.

To the best of our knowledge, research on this model includes Hall and Potts [26], Li and Ou [27], and Wang and Cheng [28]. Hall and Potts [26] provided algorithms for models in which jobs are produced on machines and formed in batches for delivery. Li and Ou [27] studied the problem with a vehicle traveling between a warehouse and the factory. They proved the problem is NP-hard in strong sense and developed polynomial time and effective heuristic algorithms, respectively, for special case and the general setting. Wang and Cheng [28] studied the complexity of the model in which the warehouse, the factory, and the customer locations are different. They proved the problem is strongly NP-hard and proposed a heuristic with a tight bound of 2.

The model considered in this paper is different from that in [26]. Each job in our model has its own transportation time between the warehouse and the factory and between the factory and the customer. Our model also differs from Li and Ou’s model in which the locations of the warehouse and the customer are the same. Furthermore, this paper considers flow shop during the production scheduling, while the model in [28] considers processing jobs in a single machine.

Quantum genetic algorithms (QGA) are heuristic search techniques inspired from the principles of survival of the fittest in natural genetic evolution and quantum theory. They are known to be efficient in a large search space, without explicitly requiring additional information (such as convexity or derivative information) about the objective. In addition, the -bit encoding method has the superposition and probability expression characteristics, and each individual can express more states to increase the diversity of populations. For these reasons, in the last few years, they have been applied to many combinatorial problems, including scheduling and vehicle routing applications that are partially related to our problem.

In order to solve the integrated supply chain scheduling model with the materials pickup, flow shop scheduling, and the finished products delivery, a technique called mutualism quantum genetic algorithm (MQGA) is developed, which differs from traditional QGA method in two aspects. The first aspect is the encoding method. With the encoding method in [29], a chromosome including two segments, the job sequence and the job-to-batch assignment, is designed, and this encoding method could equally and uniquely represent all possible solutions. Furthermore, in order to improve the encoding performance, a local speeding-up method is introduced to perform the deep search operation. After that, a checking routine is carried out to ensure the feasibility of each solution, that is, ensure that the total physical space of each delivery batch could not exceed the capacity of the vehicle. The second aspect is as follows: based on the relationship between two mutually symbiotic species, a population growth model (namely, mutualism) is introduced to improve the performance of the algorithm. In this model, a symbiosis strategy is developed to dynamically adjust the sizes of two subpopulations. This strategy helps to increase the diversity of genes and avoid premature convergence.

The remainder of this paper is organized as follows. In Section 2, the problem including the assumptions and notations is introduced. In Section 3, the two species growth model and the corresponding evolutionary strategy are proposed. The mechanism of MQGA is introduced in Section 4, and the experimental results are discussed in Section 5. Section 6 contains some concluding remarks.

2. Problem Description

In this paper, it is assumed that the material warehouse, the factory, and the customer are located at different addresses, and there is a set of jobs, , processed in a flow shop scheduling, which consists of machines . All unprocessed jobs are initially located at the supplier’s warehouse and need to be transported to the factory for processing by a vehicle (conveyor). The conveyor is located at the warehouse at the beginning and is available for transportation of jobs from warehouse to the factory and comes back. Each pickup journey takes a constant amount of time, and in a trip the conveyor loads a limited number of jobs due to the limitation of its space capacity. In addition, the finished jobs need to be delivered to a customer by another vehicle (truck), which is initially located at the factory. Because of the constraint of its space capacity, the truck transports limited finished jobs to the customer and returns to the factory. Each trip also takes a constant amount of time. The objective is to minimize the makespan, that is, the arrival time of the last finished job to the customer. Figure 1 gives a description of the integrated supply chain problem including the pickup of materials, flow shop scheduling, and the delivery of finished products.

2.1. Assumptions

Before the introduction of the notations, a number of assumptions are given as follows:(i)all the facilities including conveyor, machines, and truck are available from time zero;(ii)there is no idle time between any consecutive two pickup journeys;(iii)the time of loading and unloading jobs is included in the pickup time and delivery time;(iv)the time of transporting jobs between machines is negligible;(v)the storage or buffer capacities between successive machines are unlimited;(vi)there is no priority among jobs;(vii)machine failure is not considered.

2.2. Notations

The following notations are used in the problem model:: physical space of job ;: physical space capacity of the conveyor (truck);: travel time from the warehouse to the factory (from the factory to the warehouse);: travel time from the factory to the customer (from the customer to the factory);: the arrival time of to the factory;: the arrival time of to the customer;: processing time of on machine ;: the earliest start time of the th position job on machine ;: departure time of the th batch from the warehouse to the workshop;: departure time of the th batch from the workshop to the customer;: equal to 1 if job is scheduled on the th position on machine ; 0, otherwise;: the latest completion time of the jobs in the th delivery batch;: equal to 1 if job is in the th pickup batch; 0, otherwise;: equal to 1 if job is in the th delivery batch; 0, otherwise;: very large positive constant.

2.3. The Optimization Model

Tang and Gong [30] proposed a mixed integer programming (MIP) model for a single-machine batch scheduling with pickup and delivery transportations. Based on it, this paper studies a different integrated supply chain problem, which consists of the pickup of materials, flow shop scheduling, and the delivery of finished products. For this problem, a MIP model is constructed below. This model is to decide the following policy: how and when the materials are assigned to an appropriate delivery batch, what the scheduling of n jobs on a series of machines is, how to dispatch the finished product as soon as possible to satisfy the customer’s demand:Objective function (1) represents the arrival time of the last finished product to the customer. Constraints (2) ensure that job should be placed in one and only one position on each machine. Constraints (3) guarantee that there is one and only one job processed on the th position on each machine. Constraints (4) define that the earliest start time of the job on the th position should not be earlier than the completion time of the th position job on the same machine. Constraints (5) ensure that, for the same job, the earliest start time on machine should not be earlier than its completion time on machine . Constraints (6) define that the start processing time of a job should not be earlier than its release time. Constraints (7) indicate that the ready time for delivering the th job should not be earlier than its completion time on the last machine. Constraints (8) and (9) define that the number of jobs in each batch should satisfy the space capacity of the conveyor or the truck. Constraints (10) and (11) determine each job to be assigned to one and only one batch. Constraints (12) and (13) indicate that if no job is assigned to batch , then no job will be assigned to batch . Constraints (14) and (15) define the departure time . Constraints (16) and (17) define that the departure time of the th batch should not be earlier than the returning time of the th batch. Constrains (18) and (19) specify the nonnegativity of and . Equations (20)–(25) define the variables.

3. The Mutualism Strategy for Population Growth

As the long evolutionary process of nature, the relationship between living beings is complicated. The phenomenon that two species live together is generally referred to as the symbiosis. This paper studies a kind of symbiosis called mutualism, which refers to two species living together and depending on each other over a long period of time. The nutrition of one species is the food source of the other one.

Mutualism brings two advantages to our algorithm. Firstly, there is no need to design the fitness function, and the fitness value of an individual is obtained by the cooperation between two populations, which could reduce the dependence on the domain knowledge. Secondly, different from the traditional framework in which the evolution of populations depends on the fitness values of individuals, mutualism incorporates the cooperative behavior between an individual and its surroundings, so as to postpone the premature convergence and improve the convergence speed.

3.1. The Mutualism Population Growth Model

A differential equation of mutualism population growth model is introduced in this subsection. Assume two species A and B live in the same environment and they do not take each other as food, but the existence of one species can promote the population growth of the other one. For example, algae and fungi are two species living in lichen. Algae provide nutrient to fungi through photosynthesis, and fungi offer algae water and inorganic substance. If algae and fungi are separated in lichen, both of them will die.

Let and , respectively, be the numbers of individuals of species A and B at the th generation. Define the maximal individual number of species A (B) by , and the cooperation degree is . In this way, the differential equation of the cooperation system can be written as follows:where denotes the accrual rate, that is, birth rate subtracting death rate. In time interval , define as the modified number of species A. So . Once the population number reaches the maximal value, species A will stop growing, and the accrual rate will become zero. denotes the block function on the population size of species A. Suppose the total amount of food is 1; then is the amount of food consumed by species A. When species A and B live together in the same environment, they cooperate with each other and bring direct advantages for population growth, and denotes the influence degree of species B on the population size of species A. We have similar explanation for (27).

3.2. The Mutualism Strategy

In the standard genetic algorithm, the fitness function will be given as an input, and the fitness is easy to reach the peak value. However, in the real environment the adaptability of a species is dynamic and is affected by its surroundings. Therefore, in MQGA, the fitness value of an individual is replaced by its cooperation degree with individuals in other species, and the details on the computation of the cooperation degree are given below.

3.2.1. Expert Evaluation Method

Given two individuals and , respectively, from populations A and B, a method is proposed to evaluate the cooperation degree of individual . Let and be the objective values of and . Denote, respectively, by , , and the minimum value, the maximum value, and the average objective value of all individuals in A and B. By cooperating with individual , define as the cooperative degree of individual .

Case 1. Consider , . In this case, each of the two individuals has a high survival rate and the “matching” ability to cooperate. They can offer nutrition to each other, and the collaborative effect is the best. So, the cooperative degree of individual is set to be “4.”

Case 2. Consider , . In this case, compared to individual , individual has a smaller objective value, which indicates that individual is more capable of growing and developing. If the two individuals work together, individual cannot obtain enough nutrition from individual , and the cooperation is a “mismatch.” Therefore, the cooperative degree of individual is set to be “2.”

Case 3. Consider , . Known from Case 2, in this case individual is more capable of growing and developing, and individual will absorb more nutrition from individual . That will eventually impede the development and growth of individual . Therefore, the cooperative degree of individual is defined to be “0.”

Case 4. Consider , . In this case, the two individuals have the similar low levels of growing, and a long period of time is needed for the development of the individuals. But they have the “matching” abilities and can cooperate with each other. Therefore, the cooperative degree of individual is set to be “1.” The four cases are summarized in Table 1.


, 4
, 2
, 0
, 1

3.2.2. Cooperative Degree of Individual

Given two populations A and B, denote by the cooperative degree of individual in A. Firstly choose some individuals from B to form a delegation , and can be chosen by ways including random selection, greedy selection, or tournament. Then let individual cooperate with each individual in with the expert evaluation method. The fitness of is calculated as follows:where denotes the cooperative degree of individual by cooperating with individual .

3.2.3. Cooperative Degree of a Population

The cooperative degree of a population stands for the cooperative level of the population, and it can be calculated as follows:where Pop1 (Pop2) is the delegation of population A (B). If , A is more cooperative and obtains more advantage or nutrition from B. Conversely, if , B is more cooperative than A. The two parameters and show the mutual influence degree on population sizes.

3.2.4. Adjustment of Population Sizes

In our mutualism population growth model, assumewhere are constants. Based on (26) and (27), at each generation , the population sizes of A and B will be modified by and , respectively, which are given as follows: The newly added individuals are generated randomly, and the deleted individuals are the ones with the smallest cooperative degree values.

If the cooperative degree of a population is high, its population size will increase. Thus the growth of a population size mostly depends on its ability to coordinate. If the ability is strong, the size will become bigger and bigger till the maximum value. Otherwise, the size will become small and eventually extinct. In order to survive, each population will try its best to gain good genes and cooperate with others. In this way, the quality of solution is improved and the performance of algorithm is enhanced.

3.3. Framework of the Mutualism Strategy for Population Growth

A description of the mutualism strategy for population growth is given in Figure 2.

4. The Mutualism Quantum Genetic Algorithm

The genetic algorithm (GA) proposed by Holland has achieved a great success for solving complex combinational problem. To improve GA, Han and Kim [31] recently introduced a quantum-inspired genetic algorithm, which attracted many researchers’ attention. However, the experiment results are still not satisfactory for large-scale problems. In this section, a method named mutualism quantum genetic algorithm (MQGA) is introduced to solve the integrated supply chain problem including the pickup of materials, the flow shop scheduling, and the delivery of finished jobs. The symbiosis evolutionary strategy is embedded in the quantum genetic algorithm (QGA).

4.1. Hybrid -Bit Coding

The encoding in [29] is adopted in this paper, as it has obtained the success on the coordinated scheduling problem with batching machine scheduling and two-stage transportation. Their encoding divided a chromosome into two segments, the job sequence and the job-to-batch assignment, which can represent all possible solutions. Besides, the -bit representation is combined in our algorithm for its advantage of linear superposition of solutions and increase of diversity of genes. The details of the encoding are illustrated as follows.

Firstly, based on the concept and principles of quantum computing, the smallest unit of information stored in a two-state quantum computer is called a -bit. A -bit may be in the “1” state, “0” state, or any superposition of the two. A -bit can be represented as , where and are complex numbers that specify the probability amplitudes of the corresponding states. gives the probability that the -bit will be found in the “0” state and gives the probability that the -bit will be found in the “1” state. Normalization of the state to unity guarantees . Therefore, in MQGA, a population in -bit representation is initially randomly generated, and an individual is defined as a string of -bits , where , .

Secondly, the individual is converted into binary representation in the following way. For each -bit, let be a random number generated in . If , the -bit is converted to be 1; otherwise it is 0.

Thirdly, the -bit representation needs to be changed into two segments, the job sequence and the job-to-batch assignment, which include the comprehensive information of production and delivery. Thus we give a converting rule to form three parts, flow shop part (the code size is ), material pickup part (the code size is ), and product delivery part (the code size is ), in which we set the last -bit of part 2 and part 3 as “1” to ensure the last job belongs to the last batch.

For example, given a 3-job, 3-machine problem, let 15 -bits represent a code, and every three -bits form one group. So a code contains five groups, which are then separated into 3 parts. The first three groups are defined as Part 1, the fourth group is Part 2, and the fifth group is Part 3.

For Part 1, suppose the binary representation is , which is already converted from -bit representation; then get . This job permutation is regarded as a random key representation. If values of two elements in the representation are different, let the smaller one denote the job with smaller index; otherwise, let the first one denote the job with smaller index. So, the above random key representation is corresponding to job permutation .

For Part 2 and Part 3, suppose the binary representations are and , which represent the job-to-batch assignment of the conveyor and the truck, respectively. The final code is the combination of Part 1, Part 2, and Part 3. An example is given in Figure 3.

The final code for another problem with 8 jobs and 2 machines is presented in (32) which is the individual structure of the MQGA, where the first line represents a permutation of 8 jobs and the second (third) line denotes the job-to-batch assignment in the conveyor (truck):

The rule is that the jobs in the positions after each “1” to the next “1” belong to a new batch. Therefore, in the pickup transportation, job 5 is in the first batch, job 8 and job 2 belong to the second batch, job 1, job 3, job 6, and job 4 are in the third batch, and job 7 is in the last batch. The third line is the delivery batch of truck. Similarly, in the delivery transportation, job 5, job 8, and job 2 are in the same batch, job 1 and job 3 belong to the next batch, and job 6, job 4, and job 7 are in the last batch. Figure 4 illustrates the job sequence and the job-to-batch assignment.

Assume the processing times of the 8 jobs on the two machines are given in Table 2. The pickup time and returning time of the conveyor are 4 time units and 2 time units, respectively. The delivery time and returning time of truck are 8 time units and 4 time units, respectively. Figure 5 gives the visual assignment of jobs to be transported and produced.

Processing time


4.2. Further Processing Mechanisms Solutions
4.2.1. The Speeding-Up Method of Inserting Neighborhood

When applying local search techniques to search better solutions, the neighborhood structure is very important, since it will directly affect the results. For the flow shop problem, the results in [3234] show that the inserting operation is better than exchanging when searching in the local neighborhood. Therefore, this paper employs the inserting operation during the process of search in the neighborhood. The inserting neighborhood of a permutation is a set of solutions, that is, the permutations generated by deleting a job from the original position and inserting the job into a new position.

Denote by the job in the th position in the permutation , and define as the subpermutation after removing job . Assuming job is inserted in the th position of , let be the new permutation and the position where job is inserted to achieve the minimal makespan; that is,

Take the two-machine flow shop problem as an example. Let be the makespan for a given permutation , and let and , respectively, be the earliest possible completion times of job on machines and in the forward channel. In addition, define , , respectively, denoting the earliest possible completion times of job on machines and in the backward channel. The iterative formulas are given as follows:

For the flow shop problem with more than two machines, the details of obtaining the local optimal solution are given in Procedure 1.

While () do
 Obtain the subpermutation by removing job in the original permutation .
Calculate , , and .
For and
  Calculate .
  Calculate .
 Denote by the position where job is inserted to get the minimal objective value.
Output the local optimal individual.

4.2.2. A Modified Mechanism of the Solution

Because the total physical space of all jobs in a batch may exceed the capacity of the vehicle, our encoding could not ensure a feasible solution, and a modified mechanism is needed to update the solution, which is an important step, not only in encoding process, but also in the steps of crossover and mutation, as an unfeasible solution may affect the search process in local neighborhood.

Let TVol be the capacity of the vehicle and the total space of all jobs in batch . For simplicity, denote by the jobs in batch , where is the number of jobs in . The modified mechanism is described in Procedure 2.

For batch
 If ,
 Find the job in such that
   and .
 Let and . Update the index of jobs in .

4.3. Other Operations

The operations introduced below are performed at each generation of the quantum genetic algorithm (QGA). The fitness of each individual can be calculated by the mutualism strategy in Section 3.2. In order to ensure the individual with good genes can be chosen in the offspring, the roulette selection is applied here, which could maintain the diversity of genes. In the crossover operation, two decoded individuals are selected as parents to operate single-point crossover. See Figure 6 for the illustration of this operation, in which the corresponding genes of -bit parents are swapped to produce two new offspring.

In the mutation operation, a NOT Gate is used as the mutation operator. Firstly, select individuals with mutation probability ; then randomly generate a mutation position. Quantum rotation gate is used to maintain diversity of genes, which could change the probability amplitude of quantum states. See more details in our previous work [35]. Besides, a catastrophe operator is used to avoid premature convergence. If the solution does not change in two consecutive generations, it could be regarded to be trapped in local optimal solution, and the local best solution will be reserved while others will be replaced by solutions generated randomly.

4.4. The Main Procedure of MQGA

In this subsection, the details of the MQGA are introduced in Algorithm 1.

Step  1. Set parameters: crossover probability , mutation probability , iterative
   generation GN, two subpopulations sizes , .
Step  2. Initialize subpopulations in hybrid -bit representation.
While () do
Step 3. If (subpopulation is not empty)
     Apply mutualism strategy, rescale the size of subpopulations according to population growth model.
Step  4. Perform selection, crossover, mutation operations in each subpopulation;
Step  5. If (catastrophe condition is satisfied)
     Perform catastrophe operation.
     Apply quantum rotation operation.
Step  6. Calculate objective value and fitness value of each individual in the current generation,
   and update the best solution if possible.
Output the global best result.

5. Experiment

In this paper, all algorithms are programmed with MATLAB language, and all the computations are conducted on a Pentium PC 1.66 GHZ with 512 MB memory. In order to evaluate the performance of MQGA, an extensive set of instances with different characteristics are generated based on the flow shop benchmark problems, including FT benchmark and ABZ benchmark. In order to test the relative large-scale problem, three test instances are designed, where there are 50 jobs and 5 machines, 50 jobs and 10 machines, and 100 jobs and 5 machines, seen in the Appendix. Table 3 shows the transporting times of conveyor and truck in different problems. Table 4 shows the physical spaces of the 100 jobs. These test instances are named N(FT06), N(FT10), N(FT20), N(ABZ3), N(ABZ6), N(ABZ7), N(ABZ8), N(ABZ9), N(TA1), N(TA2), and N(TA3).



Note: 1 and 2 are the capacities of conveyor and truck, respectively. Tr1 and Tr2 are the pickup and delivery times. 1 and 2 are the returning times of conveyor and truck.

Jobs 1–20
2 6 1 0.3 3 1.3 2 3.3 1.2 4 2.8 3.2 0.7 0.1 0.9 2 6 1 0.3 3

Jobs 21–40
1.3 2 3.3 1.2 4 2.8 3.2 0.7 0.1 0.9 2 6 1 0.3 3 1.3 2 3.3 1.2 4

Jobs 41–60
2.8 3.2 0.7 0.1 0.9 3 4 3 6 1 2 6 1 0.3 3 1.3 2 3.3 1.2 4

Jobs 61–80
2.8 3.2 0.7 0.1 0.9 2 6 1 0.3 3 1.3 2 3.3 1.2 4 2.8 3.2 0.7 0.1 0.9

Jobs 81–100
2 6 1 0.3 3 1.3 2 3.3 1.2 4 2.8 3.2 0.7 0.1 0.9 3 4 3 6 1

5.1. Parameter Settings

Two parameters including crossover rate and mutation rate are discussed in this section. Assume the two rates are set with the following levels:

In this way, there are a total of 25 combinations. For each combination, MGA [13] is used to test the 11 testing instances in Table 3. After the experiment, the data is transformed to be the relative percentage deviation (RPD), which is calculated as follows:where is the corresponding objective value of the solution achieved by proposed algorithm and is the minimum objective value under multiple simulation of the same test problem. In addition, for a combination of and values, the analysis of variance (AOV) model is employed to eliminate the different RPD values caused by the block (each test problem is regarded as a block). The main experiment results of the AOV are shown in Table 5.

SourceSum of squaresDfMean square value

A: 0.008340.00213.3334
B: 0.007240.00184.3732
Interactions AB0.0438160.00297.3679

According to the theory of statistical analysis, the greater the value is, the more significantly the corresponding factors affect. Table 5 shows the interaction between the two parameters has a significant effect on the performance of the algorithm. Based on the experiment, we find the best value of is 0.8 and the best value of is 0.2, since they could find the most times of optimal objective value.

5.2. Simulation Work
5.2.1. The Detailed Simulation Result of N(TA1) Problem

In order to validate the performance of MQGA, MGA [13] and QGA [36] are compared. Firstly problem N(TA1) is chosen as the representative instance to analyze the performance of algorithms. Three combinations of iterative generations and population sizes are 500 × 10 (i.e., the maximum number of iterations is 500, population size is 10, and subpopulation size is 5), 1000 × 10, and 2000 × 20. For each combination, the results are presented in Table 6, and more results are given in Table 7, where “BV” and “WV,” respectively, denote the best and worst objective values in the 10 run times. “AV” is the average objective value.

RUN500 × 101000 × 102000 × 20




500 × 10315531923165.4326733823325.7319932293216
1000 × 10313831663150.2321833513296.2314532223196.3
2000 × 20312831673147.3321033453294.6313132263186.3

Known from the results in Table 6, for each calculation, MQGA always obtains a better solution than that obtained by MGA and QGA. When the iteration is 500, the range of objective values of MQGA is (3155, 3192) and (3267, 3382) and (3199, 3229) for MGA and QGA. When the iteration is 2000, the ranges of objective values are, respectively, (3128, 3167), (3210, 3345), and (3131, 3226). Therefore, the search ability of MQGA increases faster than that of MGA and QGA as the iteration increases, and MGA and QGA are more likely to be trapped in local optimal solutions.

The convergence curves of MQGA, MGA, and QGA obtained for each combination are depicted in Figure 7, which indicates the convergence speed of MQGA is faster than that of MGA and QGA, especially for the combination with larger iteration.

The solutions distributions of MQGA, QGA, and MGA are given in Figure 8, which shows that there is no overlap in the solution space of each algorithm. Moreover, the change in the solutions achieved by MQGA is smaller than that of MGA and QGA, which implies that the quality of solution obtained by MQGA is better than that of MGA and QGA.

5.2.2. The Comparison Results of MQGA with MGA and QGA

The maximum number of iterations is 1000, and the size of population is 100 (in MQGA, there are two subpopulations, and the size of each subpopulation is 50). To be fair, each algorithm runs 10 times for each instance. The experiment results are shown in Table 8, in which each item is an average value. The “Improved%” column represents the percentage difference between the average objective values obtained by the current algorithm with MGA and is calculated by the following formula:where represents the objective value achieved by the current algorithm and represents the objective value achieved by MGA. Boldface and italic indicate the best results for each problem.


(6 * 6)

(10 * 10)

(20 * 5)

(10 * 10)

(10 * 10)

(15 * 20)

(15 * 20)

(15 * 20)

(50 * 5)

(50 * 10)

(100 * 5)

Note: represents the converged generation and represents the average computation time.

In the performance of solution, from Table 8, MQGA is superior to MGA and QGA for all the test instances, and QGA always outperforms MGA. The improvement becomes more obvious as the scale of the problem increases. For example, for medium-scale problem N(FT20), the optimal ability of MQGA improves 0.32%. For the large-scale problem N(TA3), the MQGA improves 14.76%. The reason is that when the number of jobs increases the transportation of vehicle will become the bottleneck for MGA and QGA. Therefore, MGA and QGA are not applicable to the relative large-scale problem instances. In addition, in the performance of convergence, MQGA can obtain a better solution than MGA and QGA with less iterative generations. For the relative large-scale problems like N(TA1), N(TA2), and N(TA3), MGA and QGA converge faster to the local optimal solution. However, MQGA could jump out of the local optimum with more iterative generations. Finally, in the performance of time, MQGA needs more computation time than MGA and QGA but obtains an improved solution.

6. Conclusions and Future Research

This paper studies an integrated scheduling with the materials pickup, flow shop scheduling, and the finished products delivery. The objective is to find a coordinated schedule to minimize the arrival time of the last completed product to the customer. In order to solve the problem, a biologically inspired quantum genetic algorithm is proposed with a new mutualism strategy. The experiment results demonstrate that MQGA can find a satisfactory solution with an acceptable amount of computation time.

Future research could address problems with multiple customers or multiple transport vehicles or different shop environments, including flexible scheduling and job-shop. Problems with other performance measures, including minimum mean tardiness, and multimeasures should also be studied.


See Table 9.


MachinesProcessing times

73 87 13 11 41 43 93 69 80 13 24 72 38 81 83 88 26 6 89 67 70 30 89 30 68 21 78 46 99 10 17 23 83 47 86 18 67 46 4 14 4 20 88 30 84 38 93 76 30 30

26 37 23 93 49 12 39 17 46 20 32 44 92 73 93 33 10 43 2 62 62 82 29 29 94 20 42 80 94 33 8 41 63 4 71 30 14 32 30 30 27 98 39 84 63 12 38 43 49 13

48 4 92 92 72 43 3 98 93 17 79 11 16 89 81 92 43 61 39 28 94 87 23 1 33 91 67 91 4 60 38 23 90 93 13 63 23 34 47 98 91 11 46 30 77 3 14 47 80 43

26 67 4 14 93 34 21 20 6 18 73 23 16 77 28 24 13 77 36 16 32 46 21 81 28 70 89 34 96 62 46 60 19 97 13 7 44 7 73 13 66 70 97 33 97 64 73 28 4 87

77 94 9 37 29 79 33 73 63 86 23 39 76 24 38 3 91 29 22 27 39 31 46 18 93 38 83 38 97 10 79 93 2 87 17 18 10 30 8 26 14 21 13 10 83 46 42 18 36 2


MachinesProcessing times

46 32 79 43 97 10 44 24 83 73 66 49 93 61 19 47 84 13 11 19 98 2 83 44 7 73 19 69 12 73 83 23 33 16 88 8 26 42 38 63 7 2 44 38 24 76 83 61 32 90

61 87 31 23 73 93 28 90 94 39 64 2 16 33 33 40 81 26 83 4 4 10 63 96 33 71 66 94 7 13 11 99 37 30 36 69 22 36 67 63 96 74 4 42 40 30 93 36 23 87

31 38 83 33 71 38 36 64 43 48 69 96 33 82 33 64 11 61 36 33 87 88 10 32 38 23 24 90 7 11 49 2 76 17 32 39 9 83 69 67 28 88 23 91 71 3 26 41 96

31 24 21 37 69 31 30 31 21 19 63 91 11 6 31 63 36 39 37 47 36 63 39 4 10 12 62 43 49 34 87 29 2 18 73 39 77 69 13 78 68 37 22 41 92 67 24 87 91 31

37 16 42 47 94 14 94 34 72 36 88 31 41 71 94 99 11 97 44 77 69 91 38 23 87 7 66 34 86 49 3 48 44 93 37 82 31 39 78 33 36 3 38 10 98 6 44 62 24 94

79 93 68 73 37 44 34 39 76 62 74 28 78 43 98 83 91 27 6 82 60 44 43 76 99 66 11 33 32 8 40 62 23 24 30 1 73 27 16 91 33 11 99 2 60 90 36 62 13 3

83 87 38 38 86 67 23 19 97 78 66 67 7 23 67 8 77 71 83 29 49 3 94 76 93 48 4 37 82 37 61 6 97 3 27 93 46 92 46 32 8 11 7 34 72 37 83 22 87 63

22 29 99 23 98 33 80 82 33 68 47 74 26 61 93 33 11 42 72 14 8 98 90 36 73 69 26 24 33 98 86 30 92 94 66 47 3 41 41 47 89 28 39 80 47 37 74 38 39 3

27 92 73 94 18 41 37 38 36 20 2 39 91 81 33 14 88 22 36 63 79 23 66 3 13 31 2 81 12 40 39 32 16 87 78 41 43 94 1 93 22 93 62 33 30 34 27 30 34 77

24 47 39 66 41 46 24 23 68 30 93 22 64 81 94 97 34 82 11 91 23 32 26 22 12 23 34 87 39 2 38 84 62 10 11 93 37 81 10 40 62 49 90 34 11 81 31 21 39 27


MachinesProcessing times






Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This work is supported by the National Natural Science Foundation of China (Grant nos. 11201282 and 61304209), Humanity and Social Science Youth Foundation of Ministry of Education of China (Grant no. 10YJCZH032), Innovation Program of Shanghai Municipal Education Commission (14YZ127), and the Fund of Scheme for Training Young Teachers in Colleges and Universities in Shanghai (ZZCD12006).


  1. M. L. Pinedo, Scheduling: Theory, Algorithms, and Systems, Prentice-Hall, Upper Saddle River, NJ, USA, 2nd edition, 2002. View at: MathSciNet
  2. P. L. Maggu and G. Das, “On 2n sequencing problem with transportation times of jobs,” Pure and Applied Mathematica Sciences, vol. 12, pp. 1–6, 1980. View at: Google Scholar
  3. P. L. Maggu, G. Das, and R. Kumar, “On equivalent-job for job-block in 2×n sequencing problem with transportation-times,” Journal of the Operations Research Society of Japan, vol. 24, no. 2, pp. 136–146, 1981. View at: Google Scholar | MathSciNet
  4. P. L. Maggu, M. L. Singhal, N. Mohammad, and S. K. Yadav, “On n-job, 2-machine flow-shop scheduling problem with arbitrary time lags and transportation times of jobs,” Journal of the OR Society of Japan, vol. 23, pp. 219–227, 1982. View at: Google Scholar
  5. M. A. Langston, “Interstage transportation planning in the deterministic flow-shop environment,” Operations Research, vol. 35, no. 4, pp. 556–564, 1987. View at: Publisher Site | Google Scholar
  6. M. Haouari and T. Ladhari, “Minimising maximum lateness in a two-machine flowshop,” Journal of the Operational Research Society, vol. 51, no. 9, pp. 1100–1106, 2000. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  7. J. Hurink and S. Knust, “Makespan minimization for flow-shop problems with transportation times and a single robot,” Discrete Applied Mathematics, vol. 112, no. 1–3, pp. 199–216, 2001. View at: Publisher Site | Google Scholar | MathSciNet
  8. J. Hurink and S. Knust, “A tabu search algorithm for scheduling a single robot in a job-shop environment,” Discrete Applied Mathematics, vol. 119, no. 1-2, pp. 181–203, 2002. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  9. J. Hurink and S. Knust, “Tabu search algorithms for job-shop problems with a single transport robot,” European Journal of Operational Research, vol. 162, no. 1, pp. 99–111, 2005. View at: Publisher Site | Google Scholar
  10. C.-Y. Lee and V. A. Strusevich, “Two-machine shop scheduling with an uncapacitated interstage transporter,” IIE Transactions, vol. 37, no. 8, pp. 725–736, 2005. View at: Publisher Site | Google Scholar
  11. B. Naderi, R. Tavakkoli-Moghaddam, and M. Khalili, “Electromagnetism-like mechanism and simulated annealing algorithms for flowshop scheduling problems minimizing the total weighted tardiness and makespan,” Knowledge-Based Systems, vol. 23, no. 2, pp. 77–85, 2010. View at: Publisher Site | Google Scholar
  12. M. Khalili and R. Tavakkoli-Moghaddam, “A multi-objective electromagnetism algorithm for a bi-objective flowshop scheduling problem,” Journal of Manufacturing Systems, vol. 31, no. 2, pp. 232–239, 2012. View at: Publisher Site | Google Scholar
  13. J. Behnamian, S. M. T. Fatemi Ghomi, F. Jolai, and O. Amirtaheri, “Minimizing makespan on a three-machine flowshop batch scheduling problem with transportation using genetic algorithm,” Applied Soft Computing Journal, vol. 12, no. 2, pp. 768–777, 2012. View at: Publisher Site | Google Scholar
  14. L. Tang, J. Guan, and G. Hu, “Steelmaking and refining coordinated scheduling problem with waiting time and transportation consideration,” Computers & Industrial Engineering, vol. 58, no. 2, pp. 239–248, 2010. View at: Publisher Site | Google Scholar
  15. A. Elmi and S. Topaloglu, “A scheduling problem in blocking hybrid flow shop robotic cells with multiple robots,” Computers and Operations Research, vol. 40, no. 10, pp. 2543–2555, 2013. View at: Publisher Site | Google Scholar | MathSciNet
  16. C.-Y. Lee and Z.-L. Chen, “Machine scheduling with transportation considerations,” Journal of Scheduling, vol. 4, no. 1, pp. 3–24, 2001. View at: Google Scholar | MathSciNet
  17. Y.-C. Chang and C.-Y. Lee, “Machine scheduling with job delivery coordination,” European Journal of Operational Research, vol. 158, no. 2, pp. 470–487, 2004. View at: Publisher Site | Google Scholar | MathSciNet
  18. C.-L. Li, G. Vairaktarakis, and C.-Y. Lee, “Machine scheduling with deliveries to multiple customer locations,” European Journal of Operational Research, vol. 164, no. 1, pp. 39–51, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  19. Z.-L. Chen and G. L. Vairaktarakis, “Integrated scheduling of production and distribution operations,” Management Science, vol. 51, no. 4, pp. 614–628, 2005. View at: Publisher Site | Google Scholar
  20. H. N. Geismar, G. Laporte, L. Lei, and C. Sriskandarajah, “The integrated production and transportation scheduling problem for a product with a short lifespan,” INFORMS Journal on Computing, vol. 20, no. 1, pp. 21–33, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  21. A. Soukhal, A. Oulamara, and P. Martineau, “Complexity of flow shop scheduling problems with transportation constraints,” European Journal of Operational Research, vol. 161, no. 1, pp. 32–41, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  22. M. M. Mazdeh and M. Rostami, “A branch-and-bound algorithm for two-machine flow-shop scheduling problems with batch delivery costs,” International Journal of Systems Science: Operations and Logistics, vol. 1, no. 2, pp. 94–104, 2014. View at: Google Scholar
  23. P. J. Kalczynski and J. Kamburowski, “An empirical analysis of heuristics for solving the two-machine flow shop problem with job release times,” Computers & Operations Research, vol. 39, no. 11, pp. 2659–2665, 2012. View at: Publisher Site | Google Scholar | MathSciNet
  24. M. Rasti-Barzoki, S. R. Hejazi, and M. Mahdavi Mazdeh, “A branch and bound algorithm to minimize the total weighted number of tardy jobs and delivery costs,” Applied Mathematical Modelling, vol. 37, no. 7, pp. 4924–4937, 2013. View at: Publisher Site | Google Scholar | MathSciNet
  25. C. Low, C.-M. Chang, R.-K. Li, and C.-L. Huang, “Coordination of production scheduling and delivery problems with heterogeneous fleet,” International Journal of Production Economics, vol. 153, pp. 139–148, 2014. View at: Publisher Site | Google Scholar
  26. N. G. Hall and C. N. Potts, “Supply chain scheduling: batching and delivery,” Operations Research, vol. 51, no. 4, pp. 566–584, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  27. C.-L. Li and J. Ou, “Machine scheduling with pickup and delivery,” Naval Research Logistics, vol. 52, no. 7, pp. 617–630, 2005. View at: Publisher Site | Google Scholar | MathSciNet
  28. X. Wang and T. C. Cheng, “Production scheduling with supply and delivery considerations to minimize the makespan,” European Journal of Operational Research, vol. 194, no. 3, pp. 743–752, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  29. C.-H. Liu, “Using genetic algorithms for the coordinated scheduling problem of a batching machine and two-stage transportation,” Applied Mathematics and Computation, vol. 217, no. 24, pp. 10095–10104, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  30. L. Tang and H. Gong, “A hybrid two-stage transportation and batch scheduling problem,” Applied Mathematical Modelling, vol. 32, no. 12, pp. 2467–2479, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  31. K.-H. Han and J.-H. Kim, “Quantum-inspired evolutionary algorithm for a class of combinatorial optimization,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 6, pp. 580–593, 2002. View at: Publisher Site | Google Scholar
  32. G. I. Zobolas, C. D. Tarantilis, and G. Ioannou, “Minimizing makespan in permutation flow shop scheduling problems using a hybrid metaheuristic algorithm,” Computers & Operations Research, vol. 36, no. 4, pp. 1249–1267, 2009. View at: Publisher Site | Google Scholar
  33. J. Grabowski and J. Pempera, “Some local search algorithms for no-wait flow-shop problem with makespan criterion,” Computers and Operations Research, vol. 32, no. 8, pp. 2197–2212, 2005. View at: Publisher Site | Google Scholar
  34. G. L. Deng and X. S. Gu, “A hybrid discrete differential evolution algorithm for the no-idle permutation flow shop scheduling problem with makespan griterion,” Computers and Operations Research, vol. 39, pp. 2132–2160, 2012. View at: Google Scholar
  35. J. Gu, X. Gu, and M. Gu, “A novel parallel quantum genetic algorithm for stochastic job shop scheduling,” Journal of Mathematical Analysis and Applications, vol. 355, no. 1, pp. 63–81, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  36. L. Wang, H. Wu, F. Tang, and A. Z. Zheng, “A hybrid quantum-inspired genetic algorithm for flow shop scheduling,” in Advances in Intelligent Computing, vol. 3654 of Lecture Notes in Computer Science, pp. 636–644, Springer, 2005. View at: Google Scholar

Copyright © 2015 Jinwei Gu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.