Asymptotic Analysis of SPTA-Based Algorithms for No-Wait Flow Shop Scheduling Problem with Release Dates

We address the scheduling problem for a no-wait flow shop to optimize total completion time with release dates. With the tool of asymptotic analysis, we prove that the objective values of two SPTA-based algorithms converge to the optimal value for sufficiently large-sized problems. To further enhance the performance of the SPTA-based algorithms, an improvement scheme based on local search is provided for moderate scale problems. New lower bound is presented for evaluating the asymptotic optimality of the algorithms. Numerical simulations demonstrate the effectiveness of the proposed algorithms.


Introduction
The no-wait flow shop is a very typical scheduling setting in steel production, food production, service industry, and so forth [1]. In this scheduling problem, ( ≥ 1) jobs have to be processed on a set of ( ≥ 2) machines following the same route. There is no intermediate storage between any two adjacent machines. The finished job has to remain on the upstream machine, until the downstream machine becomes available. The jobs to be processed arrive to the system over time. The goal is to achieve a job sequence that minimizes the sum of total completion times. For formally stating the scheduling problems, the standard three-field notation [2] is employed in the succeeding content.
Röck [3] reported the strong NP-hardness for problem F2|no-wait|Σ . As the F2|no-wait|Σ problem is a special case of the F |no-wait, |Σ problem, it implies that obtaining the optimal solution in polynomial time for the latter is impossible. Wang et al. [4] considered the F |no-wait|Σ (1 − − ) problem, where is the weight of job and is the discounted rate. The authors developed efficient polynomial time algorithms for finding the optimal schedules of the problem. Su and Lee [5] investigated a no-wait and separate setup two-machine flow shop system with a single server (i.e., F2 S1|no-wait, |Σ ). They proposed optimal solutions for several restricted cases and some properties for general case as well as established the heuristic and branch and bound algorithms to tackle the problem. Huang et al. [6] addressed a no-wait flow shop problem with two-stage multiprocessor (i.e., F2(P1, P2)|no-wait, |Σ , where is the setup time of job ). They implement integer programming model and an ant colony optimization to test, analyze, and compare simulated data. Framinan et al. [7] proposed a new constructive heuristic based on an analogy with the two-machine problem in order to select the candidate to be appended in the partial schedule. They conducted computational tests to show that the proposed heuristic outperforms existing ones. Zhao and Tang [8] extend the artificial immune system approach by proposing a new methodology termed as Psycho-Clonal algorithm to deal with problem F |nowait |Σ . Qian et al. [9] presented a hybrid algorithm based on differential evolution (HDE) to solve the F |no-wait, , |Σ , where denotes the sequence dependent setup time. The superiority of HDE in terms of searching quality, robustness, and efficiency is demonstrated by simulations.
For the unlimited buffer version of this problem (i.e., F | |Σ ) Liu et al. [10] reported the asymptotic optimality of two heuristic algorithms based on the Shortest Processing Time among Available job (SPTA) rule, namely, the SPTA-F and SPTA-A rules. In this paper, we try to prove that these 2 The Scientific World Journal rules are also asymptotically optimal for problem F |nowait, |Σ in sense of probability limit. An improvement scheme is introduced for the heuristics to enhance the quality of the original solutions. For numerically evaluating the experimental results, new lower bound is provided for the problem. At the end of the paper, computational results demonstrate the convergence of the heuristics and the performance of the improvement scheme.
The remainder of the paper is organized as follows. The formulated expression of the no-wait flow shop is given in Section 2. The asymptotic analysis on the SPTA-based algorithms is introduced in Section 3. The improvement scheme and new lower bound are presented in Section 4. Some computational results are provided in Section 5 and this paper closes with the conclusion in Section 6.

Problem Statement
Generally, a no-wait flow shop involves machines in series and jobs to be executed. Each job requires to be sequentially processed on each of the machines without preemption. The processing time of job , = 1, 2, . . . , , on machine , = 1, 2, . . . , , is ( , ). It is supposed that all the processing times are independently and identically distributed (i.i.d.) random variables defined on interval (0, 1]. The jobs arrive to the system with a release date , which is the earliest time when the job is available. A regular assumption on the release dates is that they are nonnegative. The permutation schedule is adopted; that is, the jobs pass the whole system maintaining the sequence that they are processed on the first machine. Let The objective value of algorithm is denoted as (ALG) and the optimal value is denoted as (OPT). The objective is to find a schedule that optimizes the total completion time (TCT) for all jobs on the final machine.

Asymptotic Analysis on the SPTA-Based Algorithms
In this section, we implement performance analysis on the SPTA-F and SPTA-A rules for the no-wait flow-shop TCT problem. The mathematical proof reveals the asymptotic optimality of these SPTA-based algorithms.

Asymptotic Optimality of the SPTA-F Algorithm.
The SPTA-F rule schedules the job with the minimum value of (1, ), = 1, 2, . . . , , among the available ones whenever the first machine is idle. Without regard to the extra time (caused by blocking) when jobs remain on the machines, the total time on the first machine can be treated as a lower bound (LB1) for the problem, where the asymptotic optimality of LB1 for the associated single machine problem has been provided by Kaminsky and Simchi-Levi [11]: where (LB1) denotes the objective value of LB1. Therefore, the following results are presented.
Proof. Denote the completion time of job , = 1, 2, . . . , , on the first machine in LB1 as (LB1). Therefore, we have Note that The Scientific World Journal 3 Under the assumption of Theorem 1, for every , = 1, 2, . . . , , with probability one, Xia et al. [12] proved that where = 2, 3, . . . , . With inequalities (4) and (5), dividing on both sizes of (3) and taking limit, we have For the objective values, we have Noting limit (6), dividing 2 on the both sides of (7), and taking limit, we can get the result.

Asymptotic Optimality of the SPTA-A Algorithm.
The SPTA-A rule schedules the job with the minimum value of = ∑ =1 ( , ), = 1, 2, . . . , , among the available ones whenever the first machine is idle. Similarly, the associated lower bound (LB2) is provided for the problem: where (LB2) denotes the objective value of LB2. For LB1 and LB2, Xia et al. [12] presented the following result.

Lemma 2.
Under the assumption of Theorem 1, for a given instance of the flow shop problem and its corresponding optimal solutions of LB1 and LB2, with probability one, one has where * (LB1) and * (LB2) denote the optimal values of LB1 and LB2, respectively. With this lemma, one can obtain the asymptotic optimality for the SPTA-A rule.
where (SPTA-A) is the objective value obtained by the SPTA-A rule.
Proof. For a given problem instance, we have Referring to the proof of Theorem 1, we obtain Noting Lemma 2 and limit (12), dividing 2 on the both sides of (7), and taking limit, we can get the result.

Improvement Scheme and Lower Bound
This section provides an improvement based on the properties of job sequencing to enhance the performance of the SPTA-based algorithms. For numerical verifying the asymptotic optimality of the algorithms, a new lower bound is presented.

Improvement Scheme for SPTA-Based Algorithm.
Before giving the improvement scheme, several properties for problem F2|no-wait, |Σ are introduced first.
Because sequence { , } is optimal, Hence, which leads to a contradiction.
Because sequence { , } is optimal, Hence, which leads to a contradiction.
Because sequence { , } is optimal, Hence, which leads to a contradiction.
Because sequence { , } is optimal, Hence, which leads to a contradiction.
On the basis of these properties, we design an improvement scheme to promote the solution obtained by SPTAbased rules.

Improvement Scheme
Step 1. Generate the initial sequence 0 with SPTA-based rule and calculate the objective value 0 . Reindex the jobs from 1 to .
Step 2. Divide the machines into − 1 groups. For each machine group = { − 1, }, = 1, 2, . . . , , execute the following adjustment substeps for sequence 0 . Step 2.2. If the objective value obtained in the previous substep is smaller than 0 , retain the current sequence; otherwise, hold the original sequence unchangeably.
Step 3. Stop the proceeding until all the jobs in each machine group are checked, and select the sequence with the minimum objective value as the final solution.
Obviously, without regard to the no-wait constraint, the improvement scheme can also solve the F | |Σ problem.
The Scientific World Journal 5

New Lower Bound.
For the F | |Σ problem, Bai and Ren [13] presented an asymptotically optimal lower bound, LB-BR, as a substitute for the optimal schedule: where (LB-BR) denotes the objective value obtained by LB-BR. Obviously, this lower bound can solve problem F |nowait, |Σ directly. But we find that LB-BR is not a real lower bound, and sometimes it may be larger than the optimal solution. Consider the following instance.

Example 4.
A two-machine flow shop scheduling problem involves two jobs (machines 1 and 2, jobs 1 and 2) without release date. The processing times of the jobs are listed as follows: where is an arbitrary small positive number. The LB-BR in the SPTA-A sequence { 1, 2} has the objective value of (LB-BR) = 3 + 2 , but the optimal sequence { 2, 1} has the optimal solution of (OPT) = 2 + 3 . Clearly, (LB-BR) > (OPT). To avoid such cases, a new lower bound, LB * , is introduced to deal with the flow shop total completion time problem, including no-wait constraint: where Calculating the value of LB * for sequence { 1, 2}, we have (LB * ) = 2 + 3 = (OPT) . (34)

Numerical Simulation
This section designed a series of simulation experiments to show the convergence of the SPTA-based rules and the effectiveness of the improvement scheme for different size problems. In the computational testing, the following parameters are varied: (1) machines: 3, 5, and 10 machines for testing the SPTAbased rules; 10, 15, and 20 machines for testing improvement scheme; 6 The Scientific World Journal  (2) jobs: 100, 500, 1000, and 1500 jobs for testing the SPTA-based rules; 10, 20, 50, 80, and 100 jobs for testing improvement scheme; (3) processing times: randomly generated from a discrete uniform distribution on [12,13] and a discrete normal distribution with expectation 5 and variance 25; (4) release dates: 1 = 0, +1 = + , = 2, . . . , , where is a random variable generated independently from a discrete uniform distribution on [12,13].
Ten different random trials are performed for each combination of these parameters, respectively. The averages are reported in the tables. Tables  1 and 2 listed the mean relative ratio ( (SPTA-X) − (LB * ))/ (LB * ) for the SPTA-F/A rule, where (SPTA-X) denotes the SPTA-F and SPTA-A rules. Tables 3 and 4 report the ratio (LB4)/ (LB3) for the associated lower bound obtained by the SPTA-F/A rule. The data in Tables 1 and 2 showed that the ratios approach zero as the number of jobs increases, indicating that convergence trend is independent of processing time distributions and suggesting that the SPTA-based algorithms are asymptotically optimal. For five machines and uniform distribution, for example, the ratios of the SPTA-F rule dropped from 0.05581 to 0.00895, when the number of jobs increases from 100 to 1500. The ratios in Tables 3 and 4 revealed that one lower bound dominating the other is impossible. With Lemma 2, we can conclude that the two lower bounds are asymptotically equal when the job number is sufficiently large for the same instances.

Performance of Improvement Scheme.
As the experiments in Section 5.1 demonstrate nonsignificant difference for the SPTA-F and SPTA-A rules, we select one of them for testing the performance of the improvement scheme. Table 5 listed the mean relative ratio ( (SPTA-F) − (IS))/ (IS) for the improvement scheme, where (IS) denotes the objective value obtained by the improvement scheme. The data in Table 5 indicate that the improvement scheme enhances the performance of the SPTA-F rule effectively for moderatesized problems. As the problem scale continues to enlarge, the improvement is weakened and the running time lengthens. Therefore, obtaining the near-optimal solution with the SPTA-based rules directly for large-scale problems is more practical.

Conclusions
In this paper, we first proved that the SPTA-based algorithms are asymptotically optimal as the problem size tends to infinity for problem F |no-wait, |Σ . And then, an improvement scheme based on job swap is presented to boost the performance of the SPTA-based algorithms for moderate scale problems. To evaluate the convergence trend of the The Scientific World Journal 7 algorithms, new lower bound is provided. Computational results show the effectiveness of the proposed algorithms