Next Article in Journal / Special Issue
Simple Constructive, Insertion, and Improvement Heuristics Based on the Girding Polygon for the Euclidean Traveling Salesman Problem
Previous Article in Journal
Parameterized Optimization in Uncertain Graphs—A Survey and Some Results
Previous Article in Special Issue
Linking Scheduling Criteria to Shop Floor Performance in Permutation Flowshops
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two-Machine Job-Shop Scheduling Problem to Minimize the Makespan with Uncertain Job Durations

by
Yuri N. Sotskov
1,*,
Natalja M. Matsveichuk
2 and
Vadzim D. Hatsura
3
1
United Institute of Informatics Problems, National Academy of Sciences of Belarus, Surganova Street 6, 220012 Minsk, Belarus
2
Department of Automated Production Management Systems, Belarusian State Agrarian Technical University, Nezavisimosti Avenue 99, 220023 Minsk, Belarus
3
Department of Electronic Computing Machines, Belarusian State University of Informatics and Radioelectronics, P. Brovki Street 6, 220013 Minsk, Belarus
*
Author to whom correspondence should be addressed.
Algorithms 2020, 13(1), 4; https://doi.org/10.3390/a13010004
Submission received: 30 October 2019 / Revised: 14 December 2019 / Accepted: 16 December 2019 / Published: 20 December 2019
(This article belongs to the Special Issue Exact and Heuristic Scheduling Algorithms)

Abstract

:
We study two-machine shop-scheduling problems provided that lower and upper bounds on durations of n jobs are given before scheduling. An exact value of the job duration remains unknown until completing the job. The objective is to minimize the makespan (schedule length). We address the issue of how to best execute a schedule if the job duration may take any real value from the given segment. Scheduling decisions may consist of two phases: an off-line phase and an on-line phase. Using information on the lower and upper bounds for each job duration available at the off-line phase, a scheduler can determine a minimal dominant set of schedules (DS) based on sufficient conditions for schedule domination. The DS optimally covers all possible realizations (scenarios) of the uncertain job durations in the sense that, for each possible scenario, there exists at least one schedule in the DS which is optimal. The DS enables a scheduler to quickly make an on-line scheduling decision whenever additional information on completing jobs is available. A scheduler can choose a schedule which is optimal for the most possible scenarios. We developed algorithms for testing a set of conditions for a schedule dominance. These algorithms are polynomial in the number of jobs. Their time complexity does not exceed O ( n 2 ) . Computational experiments have shown the effectiveness of the developed algorithms. If there were no more than 600 jobs, then all 1000 instances in each tested series were solved in one second at most. An instance with 10,000 jobs was solved in 0.4 s on average. The most instances from nine tested classes were optimally solved. If the maximum relative error of the job duration was not greater than 20 % , then more than 80 % of the tested instances were optimally solved. If the maximum relative error was equal to 50 % , then 45 % of the tested instances from the nine classes were optimally solved.

1. Introduction

A lot of real-life scheduling problems involve different forms of uncertainties. For dealing with uncertain scheduling problems, several approaches have been developed in the literature. A stochastic approach assumes that durations of the jobs are random variables with specific probability distributions known before scheduling. There are two types of stochastic scheduling problems [1], where one is on stochastic jobs and another is on stochastic machines. In the stochastic job problem, each job duration is assumed to be a random variable following a certain probability distribution. With an objective of minimizing the expected makespan, the flow-shop problem was considered in References [2,3,4]. In the stochastic machine problem, each job duration is a constant, while each completion time of the job is a random variable due to the machine breakdown or nonavailability. In References [5,6,7], flow-shop problems to stochastically minimize the makespan or total completion time have been considered.
If there is no information to determine a probability distribution for each random duration of the job, other approaches have to be used [8,9,10]. In the approach of seeking a robust schedule [8,11,12,13], a decision maker prefers a schedule that hedges against the worst-case scenario. A fuzzy approach [14,15,16] allows a scheduler to find best schedules with respect to fuzzy durations of the jobs. A stability approach [17,18,19,20] is based on the stability analysis of optimal schedules to possible variations of the durations. In this paper, we apply the stability approach to the two-machine job-shop scheduling problem with given segments of job durations. We have to emphasize that uncertainties of the job durations considered in this paper are due to external forces in contrast to scheduling problems with controllable durations [21,22,23], where the objective is to determine optimal durations (which are under the control of a decision maker) and to find an optimal schedule for the jobs with optimal durations.

2. Contributions and New Results

We study the two-machine job-shop scheduling problem with uncertain job durations and address the issue of how to best execute a schedule if each duration may take any value from the given segment. The main aim is to determine a minimal dominant set of schedules (DS) that would contain at least one optimal schedule for each feasible scenario of the distribution of durations of the jobs.
It is shown how an uncertain two-machine job-shop problem may be decomposed into two uncertain two-machine flow-shop problems. We prove several sufficient conditions for the existence of a small dominant set of schedules. In particular, the sufficient and necessary conditions are proven for the existence of a single pair of job permutations, which is optimal for the two-machine job-shop problem with any possible scenario. We investigated properties of the optimal pairs of job permutations for the uncertain two-machine job-shop problem.
In the stability approach, scheduling decisions may consist of two phases: an off-line phase and an on-line phase. Using information on the lower and upper bounds on each job duration available at the off-line phase, a scheduler can determine a small (or minimal) dominant set of schedules based on sufficient conditions for schedule dominance. The DS optimally covers all scenarios in the sense that, for each possible scenario, there exists at least one schedule in the DS that is optimal. The DS enables a scheduler to quickly make an on-line scheduling decision whenever additional information on completing some jobs becomes available. The stability approach enables a scheduler to choose a schedule, which is optimal for the most possible scenarios.
In this paper, we develop algorithms for testing a set of conditions for a schedule dominance. The developed algorithms are polynomial in the number of jobs. Their asymptotic complexities do not exceed O ( n 2 ) , where n is a number of the jobs. Computational experiments have shown effectiveness of the developed algorithms: if there were no more than 600 jobs, then all 1000 instances in each tested series were solved in no more than one second. For the tested series of instances with 10,000 jobs, all 1000 instances of a series were solved in 344 seconds at most (on average, 0.4 s per one instance).
The paper is organized as follows. In Section 3, we present settings of the uncertain scheduling problems. The related literature and closed results are discussed in Section 4. In Section 4.2, we describe in detail the results published for the uncertain two-machine flow-shop problem. These results are used in Section 5, where we investigate properties of the optimal job permutations used for processing a set of the given jobs. Some proofs of the claims are given in Appendix A. In Section 6, we develop algorithms for constructing optimal schedules if the proven dominance conditions hold. In Section 7, we report on the wide computational experiments for solving a lot of randomly generated instances. Tables with the obtained computational results are presented in Appendix B. The paper is concluded in Section 8, where several directions for further researches are outlined.

3. Problem Settings and Notations

Using the notation α | β | γ [24], the two-machine job-shop scheduling problem with minimizing the makespan is denoted as J 2 | n i 2 | C m a x , where α = J 2 denotes a job-shop system with two available machines, n i is the number of stages for processing a job, and γ = C m a x denotes the criterion of minimizing the makespan. In the problem J 2 | n i 2 | C m a x , the set J = { J 1 , J 2 , , J n } of the given jobs have to be processed on machines from the set M = { M 1 , M 2 } . All jobs are available for processing from the initial time t = 0 . Let O i j denote an operation of the job J i J processed on machine M j M . Each machine can process a job J i J no more than once provided that preemption of each operation O i j is not allowed. Each job J i J has its own processing order (machine route) on the machines in M .
Let J 1 , 2 denote a subset of the set J of the jobs with the same machine route ( M 1 , M 2 ), i.e., each job J i J 1 , 2 has to be processed first on machine M 1 and then on machine M 2 . Let J 2 , 1 J denote a subset of the jobs with the opposite machine route ( M 2 , M 1 ). Let J k J denote a set of the jobs, which has to be processed only on machine M k M . The partition J = J 1 J 2 J 1 , 2 J 2 , 1 holds. We denote m h = | J h | , where h { 1 ; 2 ; 1 , 2 ; 2 , 1 } .
We first assume that the duration p i j of each operation O i j is fixed before scheduling. The considered criterion C m a x is the minimization of the makespan (schedule length) as follows:
C m a x : = min s S C m a x ( s ) = min s S { max { C i ( s ) : J i J } } ,
where C i ( s ) denotes a completion time of the job J i J in the schedule s and S denotes a set of semi-active schedules existing for the problem J 2 | n i 2 | C m a x . A schedule is called semi-active if no job (operation) can be processed earlier without changing the processing order or violating some given constraints [1,25,26].
Jackson [27] proved that the problem J 2 | n i 2 | C m a x is polynomially solvable and that the optimal schedule for this problem may be determined as a pair ( π , π ) of the job permutations (calling it a Jackson’s pair of permutations) such that π = ( π 1 , 2 , π 1 , π 2 , 1 ) is a sequence of all jobs from the set J 1 J 1 , 2 J 2 , 1 processed on machine M 1 and π = ( π 2 , 1 , π 2 , π 1 , 2 ) is a sequence of all jobs from the set J 2 J 1 , 2 J 2 , 1 processed on machine M 2 . Job J j belongs to the permutation π h if J j J h .
In a Jackson’s pair ( π , π ) of the job permutations, the order for processing jobs from set J 1 (from set J 2 , respectively) may be arbitrary, while for the permutation π 1 , 2 , the following inequality holds for all indexes k and m, 1 k < m m 1 , 2 :
min { p i k 1 , p i m 2 } min { p i m 1 , p i k 2 }
(for the permutation π 2 , 1 , the following inequality holds for all indexes k and m, 1 k < m m 2 , 1 ) [28]:
min { p j k 2 , p j m 1 } min { p j m 2 p j k 1 }
The aim of this paper is to investigate the uncertain two-machine job-shop scheduling problem. Therefore, we next assume that duration p i j of each operation O i j is unknown before scheduling; namely, in the realization of a schedule, a value of p i j may be equal to any real number no less than the given lower bound l i j and no greater than the given upper bound u i j . Furthermore, it is assumed that probability distributions of random durations of the jobs are unknown before scheduling. Such a job-shop scheduling problem is denoted as J 2 | l i j p i j u i j , n i 2 | C m a x . The problem J 2 | l i j p i j u i j , n i 2 | C m a x is called an uncertain scheduling problem in contrast to the deterministic scheduling problem J 2 | n i 2 | C m a x . Let a set of all possible vectors p = ( p 1 , 1 , p 1 , 2 , , p n 1 , p n 2 ) of the job durations be determined as follows: T = { p : l i j p i j u i j , J i J , M j M } . Such a vector p = ( p 1 , 1 , p 1 , 2 , , p n 1 , p n 2 ) T of the possible durations of the jobs is called a scenario.
It should be noted that the problem J 2 | l i j p i j u i j , n i 2 | C m a x is mathematically incorrect. Indeed, in most cases, a single pair of job permutations which is optimal for all possible scenarios p T for the uncertain problem J 2 | l i j p i j u i j , n i 2 | C m a x does not exist. Therefore, in the general case, one cannot find an optimal solution for this uncertain scheduling problem.
For a fixed scenario p T , the uncertain problem J 2 | l i j p i j u i j , n i 2 | C m a x turns into the deterministic problem J 2 | n i 2 | C m a x associated with scenario p. The latter deterministic problem is an individual one and we denote it as J 2 | p , n i 2 | C m a x . For any fixed scenario p T , there exists a Jackson’s pair of the job permutations that is optimal for the individual deterministic problem J 2 | p , n i 2 | C m a x associated with scenario p.
Let S 1 , 2 denote a set of all permutations of m 1 , 2 jobs from the set J 1 , 2 , where | S 1 , 2 | = m 1 , 2 ! . Let S 2 , 1 denote a set of all permutations of m 2 , 1 jobs from the set J 2 , 1 , where | S 2 , 1 | = m 2 , 1 ! . Let S = < S 1 , 2 , S 2 , 1 > be a subset of the Cartesian product ( S 1 , 2 , π 1 , S 2 , 1 ) × ( S 2 , 1 , π 2 , S 1 , 2 ) such that each element of the set S is a pair of job permutations ( π , π ) S , where π = ( π 1 , 2 i , π 1 , π 2 , 1 j ) and π = ( π 2 , 1 j , π 2 , π 1 , 2 i ) , 1 i m 1 , 2 ! , 1 j m 2 , 1 ! . The set S determines all semi-active schedules and vice versa.
Remark 1.
As an order for processing jobs from set J 1 (from set J 2 ) may be arbitrary in the Jackson’s pair of job permutations ( π , π ) , in what follows, we fix both permutations π 1 and π 2 in the increasing order of the indexes of their jobs. Thus, both permutations π 1 and π 2 are now fixed, and so their upper indexes are omitted in each permutation from the pair ( π , π ) = ( ( π 1 , 2 i , π 1 , π 2 , 1 j ) , ( π 2 , 1 j , π 2 , π 1 , 2 i ) ) .
Due to Remark 1, the equality | S | = m 1 , 2 ! · m 2 , 1 ! holds. The following definition is used for a J-solution for the uncertain problem J 2 | l i j p i j u i j , n i 2 | C m a x .
Definition 1.
A minimal (with respect to the inclusion) set of pairs of job permutations S ( T ) S is called a J-solution for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with set J of the given jobs if, for each scenario p T , the set S ( T ) contains at least one pair ( π , π ) S of the job permutations, which is optimal for the individual deterministic problem J 2 | p , n i 2 | C m a x associated with scenario p.
From Definition 1, it follows that, for any proper subset S of the set S ( T ) S S ( T ) , there exists at least one scenario p T such that set S does not contain an optimal pair of job permutations for the individual deterministic problem J 2 | p , n i 2 | C m a x associated with scenario p , i.e., set S ( T ) is a minimal (with respect to the inclusion) set possessing the property indicated in Definition 1.
The uncertain job-shop problem J 2 | l i j p i j u i j , n i 2 | C m a x is a generalization of the uncertain flow-shop problem F 2 | l i j p i j u i j | C m a x , where all jobs from the set J have the same machine route. Two flow-shop problems are associated with the individual job-shop problem J 2 | l i j p i j u i j , n i 2 | C m a x . In one of these flow-shop problems, an optimal schedule for processing jobs J 1 , 2 has to be determined, i.e., J 2 , 1 = J 1 = J 2 = . In another flow-shop problem, an optimal schedule for processing jobs J 2 , 1 has to be determined, i.e., J 1 , 2 = J 1 = J 2 = . Thus, a solution of the problem J 2 | l i j p i j u i j , n i 2 | C m a x may be based on solutions of the two associated problems F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 and with job set J 2 , 1 .
The permutation π 1 , 2 of all jobs from set J 1 , 2 (the permutation π 2 , 1 of all jobs from set J 2 , 1 , respectively) is called a Johnson’s permutation, if the inequality in Equation (1) holds for the permutation π 1 , 2 (the inequality in Equation (2) holds for the permutation π 2 , 1 , respectively). As it is proven in Reference [28], a Johnson’s permutation is optimal for the deterministic problem F 2 | | C m a x .

4. A Literature Review and Closed Results

In this section, we address uncertain shop-scheduling problems if it is impossible to obtain probability distributions for random durations of the given jobs. In particular, we consider the uncertain two-machine flow-shop problem with the objective of minimizing the makespan. This problem is well studied and there are a lot of results published in the literature, unlike the uncertain job-shop problem.

4.1. Uncertain Shop-Scheduling Problems

The stability approach was proposed in Reference [17] and developed in Reference [18,29,30,31] for the C m a x criterion, and in References [19,32,33,34,35] for the total completion time criterion C i : = min s S J i J C i ( s ) . The stability approach combines a stability analysis of the optimal schedules, a multi-stage decision framework, and the solution concept of a minimal dominant set S ( T ) of schedules, which optimally covers all possible scenarios. The main aim of the stability approach is to construct a schedule which remains optimal for most scenarios of the set T. The minimality of the dominant set S ( T ) is useful for the two-phase scheduling described in Reference [36].
At the off-line phase, one can construct set S ( T ) , which enables a scheduler to make a quick scheduling decision at the on-line phase whenever additional local information becomes available. The knowledge of the minimal dominant set S ( T ) enables a scheduler to execute best a schedule and may end up executing a schedule optimally in many cases of the problem F 2 | l i j p i j u i j | C m a x [36]. In Reference [17], a formula for calculating the stability radius of an optimal schedule is proven, i.e., the largest value of independent variations of the job durations in a schedule such that this schedule remains optimal. In Reference [19], a stability analysis of a schedule minimizing the total completion time was exploited in the branch-and-bound method for solving the job-shop problem J m | l i j p i j u i j | C i with m machines. In Reference [29], for the two-machine flow-shop problem F 2 | l i j p i j u i j | C m a x , sufficient conditions have been identified when the transposition of two jobs minimizes the makespan.
Reference [37] addresses the total completion time objective in the flow-shop problem with uncertain durations of the jobs. A geometrical algorithm has been developed for solving the flow-shop problem F m | l i j p i j u i j , n = 2 | C i with m machines and two jobs. For this problem with two or three machines, sufficient conditions are determined such that the transposition of two jobs minimizes C i . Reference [38] is devoted to the case of separate setup times with the criterion of minimizing the makespan or total completion time. The job durations are fixed while each setup time is relaxed to be a distribution-free random variable within the given lower and upper bounds. Local and global dominance relations have been determined for the flow-shop problem with two machines.
Since, for the problem F 2 | l i j p i j u i j | C m a x there often does not exist a single permutation of n jobs J = J 1 , 2 which remains optimal for all possible scenarios, an additional criterion may be introduced for dealing with uncertain scheduling problems. In Reference [39], a robust solution minimizing the worst-case deviation from optimality was proposed to hedge against uncertainties. While the deterministic problem F 2 | | C m a x is polynomially solvable (the optimal Johnson’s permutation may be constructed for the problem F 2 | | C m a x in O ( n log n ) time), finding a job permutation minimizing the worst-case regret for the uncertain counterpart with a finite set of possible scenarios is NP hard.
In Reference [40], a binary NP hardness has been proven for finding a pair ( π k , π k ) S of identical job permutations that minimizes the worst-case absolute regret for the uncertain two-machine flow-shop problem with the criterion C m a x even for two possible scenarios. Minimizing the worst-case regret implies a time-consuming search over the set of n ! job permutations. In order to overcome this computational complexity in some cases, it is useful to consider a minimal dominant set of schedules S ( T ) instead of the whole set S. To solve the flow-shop problem F 2 | l i j p i j u i j | C m a x with job set J , one can restrict a search within the set S ( T ) .
We next describe in detail the results published for the flow-shop problem F 2 | l i j p i j u i j | C m a x since we use them for solving the job-shop problem J 2 | l i j p i j u i j , n i 2 | C m a x in Section 5, Section 6 and Section 7.

4.2. Closed Results

Since each permutation π uniquely determines a set of the earliest completion times C i ( π ) of the jobs J i J for the problem F 2 | | C m a x , one can identify the permutation π , ( ( π , π ) S ), with the semi-active schedule [1,25,26] determined by the permutation π . Thus, the set S becomes a set of n ! pairs ( π , π ) of identical permutations of n = m 1 , 2 jobs from the set J = J 1 , 2 since the order for processing jobs J 1 , 2 on both machines may be the same in the optimal schedule [28]. Therefore, the above Definition 1 is supplemented by the following remark.
Remark 2.
For the problem F 2 | l i j p i j u i j | C m a x considered in this section, it is assumed that a J-solution S ( T ) is a minimal dominant set of Johnson’s permutations of all jobs from the set J 1 , 2 , i.e., for each scenario p T , the set S ( T ) contains at least one optimal pair ( π k , π k ) of identical Johnson’s permutations π k such that the inequality in Equation (1) holds.
In Reference [36], it is shown how to delete redundant pairs of (identical) permutations from the set S for constructing a J-solution for the problem F 2 | l i j p i j u i j | C m a x with job set J = J 1 , 2 . The order of jobs J v J 1 , 2 and J w J 1 , 2 is fixed in the J-solution if there exists at least one Johnson’s permutation of the form π k = ( s 1 , J v , s 2 , J w , s 3 ) for any scenario p T . In Reference [29], the sufficient conditions are proven for fixing the order of two jobs from set J = J 1 , 2 . If one of the following conditions holds, then for each scenario p T , there exists a permutation π k = ( s 1 , J v , s 2 , J w , s 3 ) that is a Johnson’s one for the problem F 2 | p | C m a x associated with scenario p:
u v 1 l v 2   and   u w 2 l w 1 ,
u v 1 l v 2   and   u v 1 l w 1 ,
u w 2 l w 1   and   u w 2 l v 2 .
If at least one condition in Inequalities (3)–(5) holds, then there exists a J-solution S ( T ) for the problem F 2 | l i j p i j u i j | C m a x with fixed order J v J w of jobs, i.e., job J v has to be located before job J w in any permutation π i , ( π i , π i ) S ( T ) . If both conditions in Inequalities (4) and (5) do not hold, then there is no J-solution S ( T ) with fixed order J v J w in all permutations π i , ( π i , π i ) S ( T ) . If no analogous condition holds for the opposite order J w J v , then at least one permutation with job J v located before job J w or that with job J w located before job J v have to be included in any J-solution S ( T ) for the problem F 2 | l i j p i j u i j | C m a x . Theorem 1 is proven in Reference [41].
Theorem 1.
There exists a J-solution S ( T ) for the problem F 2 | l i j p i j u i j | C m a x with fixed order J v J w of the jobs J v and J w in all permutations π k , ( π k , π k ) S ( T ) if and only if at least one condition of Inequalities (4) or (5) holds.
In Reference [41], the necessary and sufficient conditions have been proven for the case when a single-element J-solution S ( T ) = { ( π k , π k ) } exists for the problem F 2 | l j m p j m u j m | C m a x . The partition J = J 0 J 1 J 2 J * of the set J = J 1 , 2 is considered, where
J 0 = { J i J : u i 1 l i 2 , u i 2 l i 1 } ,
J 1 = { J i J : u i 1 l i 2 , u i 2 > l i 1 } = { J i J \ J 0 : u i 1 l i 2 } ,
J 2 = { J i J : u i 1 > l i 2 , u i 2 l i 1 } = { J i J \ J 0 : u i 2 l i 1 } ,
J * = { J i J : u i 1 > l i 2 , u i 2 > l i 1 } .
For each job J k J 0 , inequalities u k 1 l k 2 and u k 2 l k 1 imply inequalities l k 1 = u k 1 = l k 2 = u k 2 . Since both segments of the possible durations of the job J k on machines M 1 and M 2 become a point, the durations p k 1 and p k 2 are fixed and equal for both machines M 1 and M 2 : p k 1 = p k 2 = : p k . In Reference [41], Theorems 2 and 3 have been proven.
Theorem 2.
There exists a single-element J-solution S ( T ) S , | S ( T ) | = 1 , for the problem F 2 | l i j p i j u i j | C m a x if and only if
(a) for any pair of jobs J i and J j from the set J 1 (from the set J 2 , respectively), either u i 1 l j 1 or u j 1 l i 1 (either u i 2 l j 2 or u j 2 l i 2 ),
(b) | J * | 1 ; for job J i * J * , the inequalities l i * 1 max { u i 1 : J i J 1 } , l i * 2 max { u j 2 : J j J 2 } hold; and max { l i * 1 , l i * 2 } p k for each job J k J 0 .
Theorem 2 characterizes the simplest case of the problem F 2 | l i j p i j u i j | C m a x when one permutation π k of the jobs J = J 1 , 2 dominates all other job permutations. The hardest case of this problem is characterized by the following theorem.
Theorem 3.
If max { l i k : J i J , M k M } < min { u i k : J i J , M k M } , then S ( T ) = S .
The J-solution S ( T ) may be represented in a compact form using the dominance digraph which may be constructed in O ( n 2 ) time. Let J × J denote the Cartesian product of two sets J . One can construct the following binary relation A J × J over set J = J 1 , 2 .
Definition 2.
For the two jobs J v J and J w J , the inclusion ( J v , J w ) A holds if and only if there exists a J-solution S ( T ) for the problem F 2 | l i j p i j u i j | C m a x such that job J v J is located before job J w J , v w , in all permutations π k , where ( π k , π k ) S ( T ) .
The binary relation ( J v , J w ) A is represented as follows: J v J w . Due to Theorem 1, if for the jobs J v J and J w J the relation J v J w , v w , holds, then for the jobs J v and J w , at least one of conditions in Inequalities (4) and (5) holds. To construct the binary relation A of the jobs on the set J , it is sufficient to check Inequalities (4) and (5) for each pair of jobs J v and J w . The binary relation A determines the digraph ( J , A ) with vertex set J and arc set A . It takes O ( n 2 ) time to construct the digraph ( J , A ) . In the general case, the binary relation A may be not transitive. In Reference [42], it is proven that, if the binary relation A is not transitive, then J 0 . We next consider the case with the equality J 0 = , i.e., J = J * J 1 J 2 (the case with J 0 has been considered in Reference [41]). For a pair of jobs J v J 1 and J w J 1 (for a pair of jobs J v J 2 and J w J 2 , respectively), it may happen that there exist both J-solution S ( T ) with job J v located before job J w in all permutations π k , ( π k , π k ) S ( T ) and J-solution S ( T ) with job J w located before job J v in all permutations π l , ( π l , π l ) S ( T ) .
In Reference [42], the following claim has been proven.
Theorem 4.
The digraph ( J , A ) has no circuits if and only if the set J = J * J 1 J 2 includes no pair of jobs J i J k and J j J k with k { 1 , 2 } such that l i k = u i k = l j k = u j k .
The binary relation A A J × J is defined as follows.
Definition 3.
For the jobs J v J and J w J , the inclusion ( J v , J w ) A holds if and only if J v J w and J w J v , or J v J w and J w J v with v < w .
The relation ( J v , J w ) A is represented as follows: J v J w . As it is shown in Reference [42], the relation J v J w implies that J v J w and that at least one condition in Inequalities (4) or (5) must hold. The relation J v J w implies exactly one of the relations J v J w or J w J v .
Since it is assumed that set J 0 is empty, the binary relation A is an antireflective, antisymmetric, and transitive relation, i.e., the binary relation A is a strict order. The strict order A determines the digraph G = ( J , A ) with arc set A . The digraph G = ( J , A ) has neither a circuit nor a loop. Properties of the dominance digraph G were studied in Reference [42]. The permutation π k = ( J k 1 , J k 2 , , J k n ) , ( π k , π k ) S , may be considered as a total strict order of all jobs of the set J . The total strict order determined by permutation π k is a linear extension of the partial strict order A if each inclusion ( J k v , J k w ) A implies inequality v < w . Let Π ( G ) denote a set of permutations π k S 1 , 2 defining all linear extensions of the partial strict order A . The cases when Π ( G ) = S 1 , 2 and Π ( G ) = { π k } are characterized in Theorems 2 and 3. In the latter case, the strict order A over set J can be represented as follows: J k 1 J k i J k i + 1 J k n 1 , 2 . In Reference [42], the following claims have been proven.
Theorem 5.
Let J = J * J 1 J 2 . For any scenario p T , the set Π ( G ) contains a Johnson’s permutation for the problem F 2 | p | C m a x .
Corollary 1.
If J = J * J 1 J 2 , then there exists a J-solution S ( T ) for the problem F 2 | l i j p i j u i j | C m a x such that π Π ( G ) for all pairs of job permutations, { ( π , π ) } S ( T ) .
In Reference [42], it was studied how to construct a minimal dominant set S ( T ) = { ( π , π ) } , π Π ( G ) . Two types of redundant permutations were examined, and the following claim was proven.
Lemma 1.
Let J = J * J 1 J 2 . If permutation π t Π ( G ) is redundant in the set Π ( G ) , then π t is a redundant permutation either of type 1 or type 2.
Testing whether set Π ( G ) contains a redundant permutation of type 1 takes O ( n 2 ) time, and testing whether permutation π g Π ( G ) is a redundant permutation of type 2 takes O ( n ) time. In Reference [42], it is shown how to delete all redundant permutations from the set Π ( G ) . Let Π * ( G ) denote a set of permutations remaining in the set Π ( G ) after deleting all redundant permutations of type 1 and type 2.
Theorem 6.
Assume the following condition:
max { l i , 3 k , l j , 3 k } < l i k = u i k = l j k = u j k < min { u i , 3 k , u j , 3 k } .
If set J = J * J 1 J 2 does not contain a pair of jobs J i J k and J j J k , k { 1 , 2 } , such that the above condition holds, then S ( T ) = < Π * ( G ) , Π * ( G ) > .
To test conditions of Theorem 6 takes O ( n ) time. Due to Theorem 6 and Lemma 1, if there are no jobs such that condition (6) holds, then a J-solution can be constructed via deleting redundant permutations from set Π ( G ) . Since the set Π * ( G ) is uniquely determined [42], we obtain Corollary 2.
Corollary 2
([42]). If set J = J * J 1 J 2 does not contain a pair of jobs J i and J j such that condition (6) holds, then the binary relation A determines a unique J-solution S ( T ) = < Π * ( G ) , Π * ( G ) > for the problem F 2 | l i j p i j u i j | C m a x .
The condition of Theorem 6 is sufficient for the uniqueness of a J-solution Π * ( G ) = S ( T ) for the problem F 2 | l i j p i j u i j | C m a x . Due to Theorem 1, one can construct a digraph G = ( J , A ) in O ( n 2 ) time. The digraph G = ( J , A ) determines a set S ( T ) and may be considered a condensed form of a J-solution for the problem F 2 | l i j p i j u i j | C m a x . The results presented in this section are used in Section 5 for constructing precedence digraphs for the problem J 2 | l i j p i j u i j , n i 2 | C m a x .

5. Properties of the Optimal Pairs of Job Permutations

We consider the uncertain job-shop problem J 2 | l i j p i j u i j , n i 2 | C m a x and prove sufficient conditions for determining a small dominant set of schedules for this problem. In what follows, we use Definition 4 of the dominant set D S ( T ) S along with Definition 1 of the J-solution S ( T ) S .
Definition 4.
A set of the pairs of job permutations D S ( T ) S is called a dominant set (of schedules) for the uncertain problem J 2 | l i j p i j u i j , n i 2 | C m a x if, for each scenario p T , the set D S ( T ) contains at least one optimal pair of job permutations for the individual deterministic problem J 2 | p , n i 2 | C m a x with scenario p.
Every J-solution (Definition 1) is a dominant set for the problem J 2 | l i j p i j u i j , n i 2 | C m a x . Before processing jobs of the set J (before the realization of a schedule s S ), a scheduler does not know exact values of the job durations. Nevertheless, it is needed to choose a pair of permutations of the jobs J , i.e., it is needed to determine orders of jobs for processing them on machine M 1 and machine M 2 . When all jobs will be processed on machines M (a schedule will be realized) and the job durations will take on exact values p i j * , l i j p i j * u i j , and so a factual scenario p * T will be determined. A schedule s chosen for the realization should be optimal for the obtained factual scenario p * . In the stability approach, one can use two phases of scheduling for solving an uncertain scheduling problem: the off-line phase and the on-line phase. The off-line phase of scheduling is finished before starting the realization of a schedule. At this phase, a scheduler knows only given segments of the job durations and the aim is to find a pair of job permutations ( π , π ) which is optimal for the most scenarios p T . After constructing a small dominant set of schedules D S ( T ) , a scheduler can choose a pair of job permutations in the set D S ( T ) , which dominates the most pairs of job permutations ( π , π ) S for the given scenarios T. Note that making a decision at the off-line phase may be time-consuming since the realization of a schedule is not started.
The on-line phase of scheduling can begin once the earliest job in the schedule ( π , π ) starts. At this phase, a scheduler can use additional on-line information on the job duration since, for each operation O i j , the exact value p i j * becomes known at the time of the completion of this operation. At the on-line phase, the selection of a next job for processing should be quick.
In Section 5.1, we investigate sufficient conditions for a pair of job permutations ( π , π ) such that equality D S ( T ) = { ( π , π ) } holds. In Section 5.2, the sufficient conditions allowing to construct a single optimal schedule dominating all other schedules in the set S are proven. If a single-element dominant set D S ( T ) does not exist, then one should construct two partial strict orders A 1 , 2 and A 2 , 1 on the set J 1 , 2 and on the set J 2 , 1 of jobs as it is described in Section 4.2. These orders may be constructed in the form of the two precedence digraphs allowing a scheduler to reduce a size of the dominant set D S ( T ) . Section 5.4 presents Algorithm 1 for constructing a semi-active schedule, which is optimal for the problem F 2 | l i j p i j u i j | C m a x for all possible scenarios T provided that such a schedule exists. Otherwise, Algorithm 1 constructs the precedence digraphs determining a minimal dominant set of schedules for the problem F 2 | l i j p i j u i j | C m a x .

5.1. Sufficient Conditions for an Optimal Pair of Job Permutations

In the proofs of several claims, we use a notion of the main machine, which is introduced within the proof of the following theorem.
Theorem 7.
Consider the following conditions in Inequalities (7) or (8):
J i J 1 , 2 u i 1 J i J 2 , 1 J 2 l i 2   a n d   J i J 1 , 2 l i 2 J i J 2 , 1 J 1 u i 1
J i J 2 , 1 u i 2 J i J 1 , 2 J 1 l i 1   a n d   J i J 2 , 1 l i 1 J i J 1 , 2 J 2 u i 2
If one of the above conditions holds, then any pair of job permutations ( π , π ) S is a single-element dominant set D S ( T ) = { ( π , π ) } for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with set J = J 1 J 2 J 1 , 2 J 2 , 1 of the given jobs.
Proof. 
Let the condition in Inequalities (7) hold. Then, we consider an arbitrary pair of job permutations ( π , π ) S with any fixed scenario p T and show that this pair of job permutations ( π , π ) is optimal for the individual deterministic problem J 2 | p , n i 2 | C m a x with scenario p, i.e., C m a x ( π , π ) = C m a x .
Let c 1 ( π ) ( c 2 ( π ) ) denote a completion time of all jobs J 1 J 1 , 2 J 2 , 1 (jobs J 2 J 1 , 2 J 2 , 1 ) on machine M 1 (machine M 2 ) in the schedule ( π , π ) , where π = ( π 1 , 2 , π 1 , π 2 , 1 ) and π = ( π 2 , 1 , π 2 , π 1 , 2 ) . For the problem J 2 | p , n i 2 | C m a x , the maximal completion time of the jobs in schedule ( π , π ) may be calculated as follows: C m a x ( π , π ) = max { c 1 ( π ) , c 2 ( π ) } .
Machine M 1 (machine M 2 ) is called a main machine for the schedule ( π , π ) if equality C m a x ( π , π ) = c 1 ( π ) holds (equality C m a x ( π , π ) = c 2 ( π ) holds, respectively).
For schedule ( π , π ) S , the following equality holds:
c 1 ( π ) = J i J 1 , 2 J 2 , 1 J 1 p i 1 + I 1 ; c 2 ( π ) = J i J 1 , 2 J 2 , 1 J 2 p i 2 + I 2 ,
where I 1 and I 2 denote total idle times of machine M 1 and machine M 2 in the schedule ( π , π ) , respectively. We next show that, if the condition in Inequalities (7) holds, then machine M 2 is a main machine for schedule ( π , π ) and machine M 2 has no idle time, i.e., machine M 2 is completely filled in the segment [ 0 , c 2 ( π ) ] for processing jobs from the set J 1 , 2 J 2 , 1 J 2 . At the initial time t = 0 , machine M 2 begins to process jobs from the set J 2 , 1 J 2 without idle times until the time moment t 1 = J i J 2 , 1 J 2 p i 2 .
From the first inequality in (7), we obtain the following relations:
J i J 1 , 2 p i 1 J i J 1 , 2 u i 1 J i J 2 , 1 J 2 l i 2 J i J 2 , 1 J 2 p i 2 = t 1 .
Therefore, at the time moment t 1 , machine M 2 begins to process jobs from the set J 1 , 2 without idle times and we obtain the following equality: c 2 ( π ) = J i J 1 , 2 J 2 , 1 J 2 p i 2 , where I 2 = 0 and machine M 2 has no idle time. We next show that machine M 2 is a main machine for the schedule ( π , π ) . To this end, we consider the following two possible cases.
(a) Let machine M 1 have no idle time.
By summing Inequalities (7), we obtain the following inequality:
J i J 1 , 2 J 2 , 1 J 1 u i 1 J i J 1 , 2 J 2 , 1 J 2 l i 2 .
Thus, the following relations hold:
c 1 ( π ) = J i J 1 , 2 J 2 , 1 J 1 p i 1 J i J 1 , 2 J 2 , 1 J 1 u i 1 J i J 1 , 2 J 2 , 1 J 2 l i 2 J i J 1 , 2 J 2 , 1 J 2 p i 2 = c 2 ( π ) .
Hence, machine M 2 is a main machine for the schedule ( π , π ) .
(b) Let machine M 1 have an idle time.
An idle time of machine M 1 is only possible if some job J j from set J 2 , 1 is processed on machine M 2 at the time moment t 2 when this job J j could be processed on machine M 1 .
Obviously, after the time moment J i J 2 , 1 p i 2 when machine M 2 completes all jobs from set J 2 , 1 , machine M 1 can process some jobs from set J 2 , 1 without an idle time. Therefore, the inequality t 2 + I 1 J i J 2 , 1 p i 2 holds and we obtain the following relations:
c 1 ( π ) t 2 + I 1 + J i J 2 , 1 p i 1 J i J 2 , 1 p i 2 + J i J 2 , 1 J 1 p i 1 J i J 2 , 1 p i 2 + J i J 2 , 1 J 1 u i 1
J i J 2 , 1 p i 2 + J i J 1 , 2 l i 2 J i J 2 , 1 p i 2 + J i J 1 , 2 p i 2 J i J 2 , 1 J 2 J 1 , 2 p i 2 = c 2 ( π ) .
We conclude that, in case (b), machine M 2 is a main machine for the schedule ( π , π ) . Thus, if the condition in Inequalities (7) holds, then machine M 2 is a main machine for the schedule ( π , π ) and machine M 2 has no idle time, i.e., equality C m a x ( π , π ) = c 2 ( π ) holds and machine M 2 is completely filled in the segment [ 0 , c 2 ( π ) ] with processing jobs from the set J 1 , 2 J 2 , 1 J 2 .
Thus, the pair of permutations ( π , π ) is optimal for scenario p T . Since scenario p was chosen arbitrarily in the set T, we conclude that the pair of job permutations ( π , π ) is a singleton D S ( T ) = { ( π , π ) } for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with set J = J 1 J 2 J 1 , 2 J 2 , 1 of the given jobs. As a pair of permutations ( π , π ) is an arbitrary pair of job permutations in the set S, any pair of job permutations ( π , π ) S is a singleton D S ( T ) = { ( π , π ) } for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J = J 1 J 2 J 1 , 2 J 2 , 1 .
The case when the condition in Inequalities (8) holds may be analyzed similarly via replacing machine M 1 by machine M 2 and vice versa. □
If conditions of Theorem 7 hold, then in the optimal pair of job permutations ( π , π ) existing for the problem J 2 | l i j p i j u i j , n i 2 | C m a x , the orders of jobs from sets J 1 , 2 J and J 2 , 1 J may be chosen arbitrarily. Theorem 7 implies the following two corollaries.
Corollary 3.
If the following inequality holds:
J j J 1 , 2 u i 1 J j J 2 , 1 J 2 l i 2 ,
then set < { π 1 , 2 } , S 2 , 1 > S , where π 1 , 2 is an arbitrary permutation in set S 1 , 2 , is a dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with set J = J 1 J 2 J 1 , 2 J 2 , 1 of the given jobs.
Proof. 
We consider an arbitrary vector p T of the job durations and an arbitrary permutation π 1 , 2 in the set S 1 , 2 . The set S 2 , 1 contains at least one Johnson’s permutation π 2 , 1 * for the deterministic problem F 2 | p 2 , 1 | C m a x with job set J 2 , 1 and scenario p 2 , 1 (the components of vector p 2 , 1 are equal to the corresponding components of vector p). We consider a pair of job permutations ( π , π ) = ( ( π 1 , 2 , π 1 , π 2 , 1 * ) , ( π 2 , 1 * , π 2 , π 1 , 2 ) ) < { π 1 , 2 } , S 2 , 1 > S and show that it is an optimal pair of job permutations for the problem J 2 | p , n i 2 | C m a x with job set J and scenario p. Without loss of generality, both permutations π 1 and π 2 are ordered in increasing order of the indexes of their jobs.
Similar to the proof of Theorem 7, one can show that, if the condition in Inequalities (9) holds, then machine M 2 processes jobs without idle times and equality c 2 ( π ) = J i J 1 , 2 J 2 , 1 J 2 p i 2 holds, where the value of c 2 ( π ) cannot be reduced. If machine M 1 has no idle time, we obtain equalities
C m a x ( π , π ) = max { c 1 ( π ) , c 2 ( π ) } = max { J i J 1 , 2 J 2 , 1 J 1 p i 1 , J i J 1 , 2 J 2 , 1 J 2 p i 2 } = C m a x .
On the other hand, an idle time of machine M 1 is only possible if some job J j from set J 2 , 1 is processed on machine M 2 at the time moment t 2 when job J j could be processed on machine M 1 . In such a case, the value of c 1 ( π ) is equal to the makespan C m a x ( π 2 , 1 * ) for the problem F 2 | p 2 , 1 | C m a x with job set J 2 , 1 and scenario p 2 , 1 . As the permutation π 2 , 1 * is a Johnson’s permutation, the value of C m a x ( π 2 , 1 * ) cannot be reduced and we obtain the following equalities:
C m a x ( π , π ) = max { c 1 ( π ) , c 2 ( π ) } = max { C m a x ( π 2 , 1 * ) , J i J 1 , 2 J 2 , 1 J 2 p i 2 } = C m a x .
Thus, the pair of job permutation ( π , π ) = ( ( π 1 , 2 , π 1 , π 2 , 1 * ) , ( π 2 , 1 * , π 2 , π 1 , 2 ) ) < { π 1 , 2 } , S 2 , 1 > S is optimal for the problem J 2 | p , n i 2 | C m a x with scenario p T . The optimal pair of job permutations for the problem J 2 | p , n i 2 | C m a x with scenario p T belongs to the set < { π 1 , 2 } , S 2 , 1 > . As vector p is an arbitrary vector in the set T, the set < { π 1 , 2 } , S 2 , 1 > contains an optimal pair of job permutations for all scenarios from set T. Due to Definition 4, the set < { π 1 , 2 } , S 2 , 1 > S is a dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J . □
Corollary 4.
Consider the following inequality:
J j J 2 , 1 u i 2 J j J 1 , 2 J 1 l i 1 .
If the above inequality holds, then set < S 1 , 2 , { π 2 , 1 } > , where π 2 , 1 is an arbitrary permutation in set S 2 , 1 , is a dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with set J = J 1 J 2 J 1 , 2 J 2 , 1 of the given jobs.
This claim may be proven similar to Corollary 3. If the conditions of Corollary 3 (Corollary 4) hold, then the order for processing jobs from set J 1 , 2 J (set J 2 , 1 J , respectively) in the optimal schedule ( π , π ) = ( ( π 1 , 2 , π 1 , π 2 , 1 ) , ( π 2 , 1 , π 2 , π 1 , 2 ) ) for the problem J 2 | l i j p i j u i j , n i 2 | C m a x may be arbitrary. Since the orders of jobs from the sets J 1 and J 2 are fixed in the optimal schedule (Remark 1), we need to determine only orders for processing jobs from set J 2 , 1 (set J 1 , 2 , respectively). To do this, we will consider two uncertain problems F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 J and with the machine route ( M 1 , M 2 ) and that with job set J 2 , 1 J and with the opposite machine route ( M 2 , M 1 ) .
Lemma 2.
If S 1 , 2 S 1 , 2 is a set of permutations from the dominant set for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 , then < S 1 , 2 , S 2 , 1 > S is a dominant set for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J = J 1 J 2 J 1 , 2 J 2 , 1 .
The proof of Lemma 2 and those for other statements in this section are given in Appendix A.
Lemma 3.
Let S 2 , 1 S 2 , 1 be a set of permutations from the dominant set for the problem F 2 | l i j p i j u i j | C m a x with job set J 2 , 1 , S 2 , 1 S 2 , 1 . Then, < S 1 , 2 , S 2 , 1 > is a dominant set for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J .
The proof of this claim is similar to that for Lemma 2 (see Appendix A).
Theorem 8.
Let S 1 , 2 S 1 , 2 be a set of permutations from the dominant set for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 , and let S 2 , 1 S 2 , 1 be a set of permutations from the dominant set for the problem F 2 | l i j p i j u i j | C m a x with job set J 2 , 1 . Then, < S 1 , 2 , S 2 , 1 > S is a dominant set for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J = J 1 J 2 J 1 , 2 J 2 , 1 .
Theorem 9.
Let a pair of identical permutations ( π 1 , 2 , π 1 , 2 ) determine a single-element J-solution for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 , and let a pair of identical permutations ( π 2 , 1 , π 2 , 1 ) determine a single-element J-solution for the problem F 2 | l i j p i j u i j | C m a x with job set J 2 , 1 . Then, the pairs of permutations { ( π 1 , 2 , π 1 , π 2 , 1 ) a n d ( π 1 , 2 , π 2 , π 2 , 1 ) } are a single-element dominant set D S ( T ) for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J = J 1 J 2 J 1 , 2 J 2 , 1 .
The following claim follows directly from Theorem 9.
Corollary 5.
If the conditions of Theorem 9 hold, then there exists a single pair of job permutations, which is an optimal pair of job permutations for the problem J 2 | p , n i 2 | C m a x with job set J and any scenario p T .
Theorem 9 implies also the following corollary proven in Appendix A.
Corollary 6.
If the conditions of Theorem 9 hold, then there exists a single pair of job permutations which is a J-solution for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J = J 1 J 2 J 1 , 2 J 2 , 1 .
Note that the criterion for a single-element J-solution for the problem F 2 | l i j p i j u i j | C m a x is given in Theorem 2.

5.2. Precedence Digraphs Determining a Minimal Dominant Set of Schedules

In Section 4.2, it is assumed that J 1 , 2 = J 1 , 2 1 J 1 , 2 2 J 1 , 2 * and J 2 , 1 = J 2 , 1 1 J 2 , 1 2 J 2 , 1 * , i.e., J 1 , 2 0 = J 2 , 1 0 = . Based on the results presented in Section 4.2, we can determine a binary relation A 1 , 2 for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 and a binary relation A 2 , 1 for this problem with job set J 2 , 1 . For job set J 1 , 2 , the binary relation A 1 , 2 determines the digraph G 1 , 2 = ( J 1 , 2 , A 1 , 2 ) with the vertex set J 1 , 2 and the arc set A 1 , 2 . For job set J 2 , 1 , the binary relation A 2 , 1 determines the digraph G 2 , 1 = ( J 2 , 1 , A 2 , 1 ) with the vertex set J 2 , 1 and the arc set A 2 , 1 .
Let us consider the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 and the corresponding digraph G 1 , 2 = ( J 1 , 2 , A 1 , 2 ) (the same results for the problem F 2 | l i j p i j u i j | C m a x with job set J 2 , 1 can be derived in a similar way).
Definition 5.
Two jobs, J x J 1 , 2 and J y J 1 , 2 , x y , are called conflict jobs if they are not in the relation A 1 , 2 , i.e., ( J x , J y ) A 1 , 2 and ( J y , J x ) A 1 , 2 .
Due to Definitions 2 and 3, for the conflict jobs J x J 1 , 2 and J y J 1 , 2 , x y , Inequalities (4) and (5) do not hold either for the case v = x with w = y or for the case v = y with w = x .
Definition 6.
The subset J x J 1 , 2 is called a conflict set of jobs if, for any job J y J 1 , 2 \ J x , either relation ( J x , J y ) A 1 , 2 or relation ( J y , J x ) A 1 , 2 holds for each job J x J x (provided that any proper subset of the set J x does not possess such a property).
From Definition 6, it follows that the conflict set J x is a minimal set (with respect to the inclusion). Obviously, there may exist several conflict sets in the set J 1 , 2 . (A conflict set of the jobs J x J 2 , 1 can be determined similarly.) Let the strict order A 1 , 2 for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 be represented as follows:
J 1 J 2 J k { J k + 1 , J k + 2 , , J k + r } J k + r + 1 J k + r + 2 J m 1 , 2 ,
where all jobs between braces are conflict ones and each of these jobs is in relation A 1 , 2 with any job located outside the brackets in Relation (10). In such a case, an optimal order for processing jobs from the set { J 1 , J 2 , , J k } is determined as follows: ( J 1 , J 2 , , J k ) .
Due to Theorem 5, we obtain that set Π ( G 1 , 2 ) of the permutations generated by the digraph G 1 , 2 contains an optimal Johnson’s permutation for each vector p 1 , 2 of the durations of jobs from the set J 1 , 2 . Thus, due to Definition 1, the singleton { ( π 1 , 2 , π 1 , 2 ) } , where π 1 , 2 Π ( G 1 , 2 ) , is a J-solution for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 . Analogously, the singleton { ( π 2 , 1 , π 2 , 1 ) } , where π 2 , 1 Π ( G 2 , 1 ) , is a J-solution for the problem F 2 | l i j p i j u i j | C m a x with job set J 2 , 1 . We can determine a dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J as follows: < Π ( G 1 , 2 ) , Π ( G 2 , 1 ) > S . The following theorems allow us to reduce a dominant set for the problem J 2 | l i j p i j u i j , n i 2 | C m a x . We use the following notation: L 2 = J i J 2 , 1 J 2 l i 2 .
Theorem 10.
Let the strict order A 1 , 2 over set J 1 , 2 = J 1 , 2 * J 1 , 2 1 J 1 , 2 2 be determined as follows: J 1 J k { J k + 1 , J k + 2 , , J k + r } J k + r + 1 J m 1 , 2 . Consider the following inequality:
i = 1 k + r u i 1 L 2 + i = 1 k l i 2 ,
If the above inequality holds, then set S = < { π } , Π ( G 2 , 1 ) > S with π Π ( G 1 , 2 ) is a dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J .
Proof. 
We consider an arbitrary vector p T of the job durations and an arbitrary permutation π from the set Π ( G 1 , 2 ) . The set Π ( G 2 , 1 ) contains at least one optimal Johnson’s permutation π 2 , 1 * for the problem F 2 | p 2 , 1 | C m a x with job set J 2 , 1 and vector p 2 , 1 of the job durations (components of this vector are equal to the corresponding components of the vector p).
We consider a pair of job permutations ( π , π ) = ( ( π , π 1 , π 2 , 1 * ) , ( π 2 , 1 * , π 2 , π ) ) S and show that it is an optimal pair of job permutations for the problem J 2 | p , n i 2 | C m a x with set J of the jobs and scenario p. To this end, we show that the value of C m a x ( π , π ) = max { c 1 ( π ) , c 2 ( π ) } cannot be reduced. Indeed, an idle time for machine M 1 is only possible if some job J j from the set J 2 , 1 is processed on machine M 2 at the same time when job J j could be processed on machine M 1 . In such a case, c 1 ( π ) is equal to the makespan C m a x ( π 2 , 1 * ) for the problem F 2 | p 2 , 1 | C m a x with job set J 2 , 1 and vector p 2 , 1 of the job durations. As permutation π 2 , 1 * is a Johnson’s permutation, the value of
c 1 ( π ) = max { J i J 1 , 2 J 2 , 1 J 1 p i 1 , C m a x ( π 2 , 1 * ) }
cannot be reduced. In the beginning of the permutation π , the jobs of set { J 1 , J 2 , , J k } are arranged in the Johnson’s order. Thus, if machine M 2 has an idle time while processing these jobs, this idle time cannot be reduced. From Inequality (11), it follows that machine M 2 has no idle time while processing jobs from the conflict set.
In the end of the permutation π , jobs of set { J k + r + 1 , , J m 1 , 2 } are arranged in Johnson’s order. Therefore, if machine M 2 has an idle time while processing these jobs, this idle time cannot be reduced. Thus, the value of c 2 ( π ) cannot be reduced by changing the order of jobs in the conflict set.
We obtain the qualities C m a x ( π , π ) = max { c 1 ( π ) , c 2 ( π ) } = C m a x . The pair of job permutations ( π , π ) = ( ( π , π 1 , π 2 , 1 * ) , ( π 2 , 1 * , π 2 , π ) ) is optimal for the problem J 2 | p , n i 2 | C m a x with scenario p T . Thus, set S = < { π } , Π ( G 2 , 1 ) > contains an optimal pair of job permutations for the problem J 2 | p , n i 2 | C m a x with scenario p T . As vector p is an arbitrary vector in set T, set S contains an optimal pair of job permutations for each vector from set T. Due to Definition 4, set S is a dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J . □
Theorem 11.
Let the partial strict order A 1 , 2 over set J 1 , 2 = J 1 , 2 * J 1 , 2 1 J 1 , 2 2 be determined as follows: J 1 J k { J k + 1 , J k + 2 , , J k + r } J k + r + 1 J m 1 , 2 . Consider the following inequality:
u k + s , 1 L 2 + i = 1 k + s 1 ( l i 2 u i 1 )
If the above inequality holds for all s { 1 , 2 , , r } , then the set S = < { π } , S 2 , 1 > , where π = ( J 1 , , J k 1 , J k , J k + 1 , J k + 2 , , J k + r , J k + r + 1 , , J m 1 , 2 ) Π ( G 1 , 2 ) , is a dominant set for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J .
Proof. 
We consider an arbitrary scenario p T and a pair of job permutations ( π , π ) = ( ( π , π 1 , π 2 , 1 * ) , ( π 2 , 1 * , π 2 , π ) ) S , where π 2 , 1 * S 2 , 1 is a Johnson’s permutation of the jobs from the set J 2 , 1 with vector p 2 , 1 of the job durations (components of this vector are equal to the corresponding components of vector p). We next show that this pair of job permutations ( π , π ) is optimal for the individual deterministic problem J 2 | p , n i 2 | C m a x with scenario p, i.e., C m a x ( π , π ) = C m a x .
If conditions of Theorem 11 hold, then machine M 2 processes jobs from the conflict set { J k + 1 , J k + 2 , , J k + r } without idle times. At the initial time t = 0 , machine M 1 begins to process jobs from the permutation π without idle times. Let a time moment t 1 be as follows: t 1 = i = 1 k + 1 p i 1 . At the time moment t 1 , job J k + 1 is ready for processing on machine M 2 .
On the other hand, at the time t = 0 , machine M 2 begins to process jobs from the set J 2 , 1 J 2 without idle times and then jobs from the permutation ( J 1 , J 2 , , J k + 1 ) . Let t 2 denote the first time moment when machine M 2 is ready for processing job J k + 1 . Obviously, the following inequality holds: t 2 L 2 + i = 1 k + 1 p i 2 . From the condition in Inequality (12) with s = 1 , we obtain inequality i = 1 k + 1 u i 1 L 2 + i = 1 k l i 2 .
Therefore, the following relations hold:
t 1 = i = 1 k + 1 p i 1 i = 1 k + 1 u i 1 L 2 + i = 1 k l i 2 L 2 + i = 1 k + 1 p i 2 = t 2 .
Machine M 2 processes job J k + 1 without an idle time between job J k and job J k + 1 .
Analogously, using s { 2 , 3 , , r } , one can show that machine M 2 processes jobs from the conflict set { J k + 1 , J k + 2 , , J k + r } without idle times between jobs J k + 1 and J k + 2 , between jobs J k + 2 and J k + 3 , and so on to between jobs J k + r 1 and J k + r . To end this proof, we have to show that the value of C m a x ( π , π ) = max { c 1 ( π ) , c 2 ( π ) } cannot be reduced.
An idle time for machine M 1 is only possible between some jobs from the set J 2 , 1 . However, the permutation π 2 , 1 * is a Johnson’s permutation of the jobs from the set J 2 , 1 for the vector p 2 , 1 of the job durations. Therefore, the value of c 1 ( π ) cannot be reduced. On the other hand, in the permutation π , all jobs J 1 , J 2 , , J k and all jobs J k + r + 1 , , J m 1 , 2 are arranged in Johnson’s orders. Therefore, if machine M 2 has an idle time while processing these jobs, this idle time cannot be reduced. It is clear that machine M 2 has no idle time while processing jobs from the conflict set. Thus, the value of c 2 ( π ) cannot be reduced by changing the order of jobs from the conflict set. We obtain the equalities C m a x ( π , π ) = max { c 1 ( π ) , c 2 ( π ) } = C m a x .
It is shown that the pair of job permutations ( π , π ) = ( ( π , π 1 , π 2 , 1 * ) , ( π 2 , 1 * , π 2 , π ) ) S is optimal for the problem J 2 | p , n i 2 | C m a x with vector p T of job durations. As vector p is an arbitrary one in set T, the set S contains an optimal pair of job permutations for each scenario from set T. Due to Definition 4, the set S is a dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J . □
The proof of the following theorem is given in Appendix A.
Theorem 12.
Let the partial strict order A 1 , 2 over set J 1 , 2 = J 1 , 2 * J 1 , 2 1 J 1 , 2 2 have the form J 1 J k { J k + 1 , J k + 2 , , J k + r } J k + r + 1 J m 1 , 2 . If inequalities
i = r s + 2 r + 1 l k + i , 1 j = r s + 1 r u k + j , 2
hold for all indexes s { 1 , 2 , , r } , then the set S = < { π } , S 2 , 1 > , where π = ( J 1 , , J k 1 , J k , J k + 1 , J k + 2 , , J k + r , J k + r + 1 , , J m 1 , 2 ) Π ( G 1 , 2 ) , is a dominant set of pairs of permutations for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J .
Similarly, one can prove sufficient conditions for the existence of an optimal job permutation for the problem F 2 | l i j p i j u i j | C m a x with job set J 2 , 1 , when the partial strict order A 2 , 1 on the set J 2 , 1 = J 2 , 1 * J 2 , 1 1 J 2 , 1 2 has the following form: J 1 J k { J k + 1 , J k + 2 , , J k + r } J k + r + 1 J m 2 , 1 .
To apply Theorems 11 and 12, one can construct a job permutation that satisfies the strict order A 1 , 2 . Then, one can check the conditions of Theorems 11 and 12 for the constructed permutation. If the set of jobs { J 1 , J 2 , , J k } is empty in the constructed permutation, one needs to check conditions of Theorem 12. If the set of jobs { J k + r + 1 , , J m 1 , 2 } is empty, one needs to check the conditions of Theorem 11. It is needed to construct only one permutation to check Theorem 11 and only one permutation to check Theorem 12.

5.3. Two Illustrative Examples

Example 1.
We consider the uncertain job-shop scheduling problem J 2 | l i j p i j u i j , n i 2 | C m a x with lower and upper bounds of the job durations given in Table 1.
These bounds determine the set T of possible scenarios. In Example 1, jobs J 1 , J 2 , and J 3 have the machine route ( M 1 , M 2 ) ; jobs J 6 , J 7 , and J 8 have the machine route ( M 2 , M 1 ) ; and job J 4 (job J 5 , respectively) has to be processed only on machine M 1 (on machine M 2 , respectively). Thus, J 1 , 2 = { J 1 , J 2 , J 3 } , J 2 , 1 = { J 6 , J 7 , J 8 } , J 1 = { J 4 } , J 2 = { J 5 } .
We check the conditions of Theorem 7 for a single pair of job permutations, which is optimal for all scenarios T. For the given jobs, the condition in Inequalities (7) of Theorem 7 holds due to the following relations:
J i J 1 , 2 u i 1 = u 1 , 1 + u 2 , 1 + u 3 , 1 = 7 + 9 + 9 = 25 J i J 2 , 1 J 2 l i 2 = l 6 , 2 + l 7 , 2 + l 8 , 2 + l 5 , 2 = 3 + 3 + 3 + 16 = 25 ;
J i J 1 , 2 l i 2 = l 1 , 2 + l 2 , 2 + l 3 , 2 = 6 + 5 + 5 = 16 J i J 2 , 1 J 1 u i 1 = u 6 , 1 + u 7 , 1 + u 8 , 1 + u 4 , 1 = 3 + 3 + 3 + 3 = 12 .
Due to Theorem 7, the order of jobs from the set J 1 , 2 = { J 1 , J 2 , J 3 } and the order of jobs from the set J 2 , 1 = { J 6 , J 7 , J 8 } may be arbitrary in the optimal pair of job permutations for the problem J 2 | l i j p i j u i j , n i 2 | C m a x under consideration. Thus, any pair of job permutations ( π , π ) S is a single-element dominant set D S ( T ) = { ( π , π ) } for Example 1.
Example 2.
Let us now consider the problem J 2 | l i j p i j u i j , n i 2 | C m a x with numerical input data given in Table 1 with the following two exceptions: l 5 , 2 = 2 and u 5 , 2 = 3 .
We check the condition in Inequalities (7) of Theorem 7 and obtain
J i J 1 , 2 u i 1 = u 1 , 1 + u 2 , 1 + u 3 , 1 = 7 + 9 + 9 = 25 J i J 2 , 1 J 2 l i 2 = l 6 , 2 + l 7 , 2 + l 8 , 2 + l 5 , 2 = 3 + 3 + 3 + 2 = 11 .
Thus, the condition of Inequalities (7) does not hold for Example 2. We check the condition of Inequalities (8) of Theorem 7 and obtain
J i J 2 , 1 u i 2 = u 6 , 2 + u 7 , 2 + u 8 , 2 = 4 + 4 + 4 = 12 J i J 1 , 2 J 1 l i 1 = l 1 , 1 + l 2 , 1 + l 3 , 1 + l 4 , 1 = 6 + 8 + 7 + 2 = 23 .
However, we see that the condition of Equation (8) does not hold:
J i J 2 , 1 l i 1 = l 6 , 1 + l 7 , 1 + l 8 , 1 = 1 + 1 + 1 = 3 J i J 1 , 2 J 2 u i 2 = u 1 , 2 + u 2 , 2 + u 3 , 2 + u 5 , 2 = 7 + 6 + 6 + 3 = 22 .
From Equation (14), it follows that the condition of Inequalities (9) of Corollary 3 does not hold. On the other hand, due to Equation (15), conditions of Corollary 4 hold. Thus, the order for processing jobs from set J 2 , 1 J in the optimal schedule ( π , π ) = ( ( π 1 , 2 , π 1 , π 2 , 1 ) , ( π 2 , 1 , π 2 , π 1 , 2 ) ) for the problem J 2 | l i j p i j u i j , n i 2 | C m a x may be arbitrary. One can fix permutation π 2 , 1 with the increasing order of the indexes of their jobs: π 2 , 1 = ( J 6 , J 7 , J 8 ) . Since the orders of jobs from the sets J 1 and J 2 are fixed in the optimal schedule (Remark 1), i.e., π 1 = ( J 4 ) and π 2 = ( J 5 ) , we need to determine the order for processing jobs in set J 1 , 2 . To this end, we consider the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 . We see that conditions of Theorem 2 do not hold for the jobs in set J 1 , 2 since J 1 J 1 , 2 * , J 2 J 1 , 2 2 , and J 3 J 1 , 2 2 ; however the following inequalities hold: u 2 , 2 > l 3 , 2 and u 3 , 2 > l 2 , 2 .
We next construct the binary relation A 1 , 2 over set J 1 , 2 based on Definition 3 and Theorem 1. Due to checking Inequalities (4) and (5), we conclude that the inequality in Equation (5) holds for the pair of jobs J 1 and J 2 . We obtain the relation J 1 J 2 . Analogously, we obtain the relation J 1 J 3 . For the pair of jobs J 2 and J 3 , neither Inequality (4) nor Inequality (5) hold. Therefore, the partial strict order A 1 , 2 over set J 1 , 2 has the following form: J 1 { J 2 , J 3 } . The job set { J 2 , J 3 } is a conflict set of these jobs (Definition 6).
Let us check whether the sufficient conditions given in Section 5.2 hold.
We check the conditions of Theorem 10 for the jobs from set J 1 , 2 . For k = 1 and r = 2 , we obtain the following equalities: L 2 = J i J 2 , 1 J 2 l i 2 = l 6 , 2 + l 7 , 2 + l 8 , 2 + l 5 , 2 = 3 + 3 + 3 + 2 = 11 . The condition of Theorem 10 does not hold since the following relations hold:
i = 1 k + r u i 1 = u 1 , 1 + u 2 , 1 + u 3 , 1 = 7 + 9 + 9 = 25 L 2 + i = 1 k l i 2 = L 2 + l 1 , 2 = 11 + 6 = 17 .
For checking the conditions of Theorem 11, we need to check both permutations of the jobs from set J 1 , 2 , which satisfy the partial strict order A 1 , 2 : Π ( G 1 , 2 ) = { π 1 , 2 1 , π 1 , 2 2 } , where π 1 , 2 1 = { J 1 , J 2 , J 3 } and π 1 , 2 2 = { J 1 , J 3 , J 2 } .
We consider permutation π 1 , 2 1 . As in the previous case, L 2 = 11 , k = 1 , r = 2 , and we must consider two inequalities in the condition in Equaiton (12) with s = 1 and s = 2 . For s = 1 , we obtain the following:
u 1 + 1 , 1 = u 2 , 1 = 9 L 2 + i = 1 1 + 1 1 ( l i 2 u i 1 ) = L 2 + i = 1 1 ( l i 2 u i 1 ) = 11 + ( l 1 , 2 u 1 , 1 ) = 11 + ( 6 7 ) = 10 .
However, for s = 2 , we obtain
u 1 + 2 , 1 = u 3 , 1 = 9 L 2 + i = 1 1 + 2 1 ( l i 2 u i 1 ) = L 2 + i = 1 2 ( l i 2 u i 1 )
= 11 + ( l 1 , 2 u 1 , 1 ) + ( l 2 , 2 u 2 , 1 ) = 11 + ( 6 7 ) + ( 5 9 ) = 6 .
Thus, the conditions of Theorem 11 do not hold for permutation π 1 , 2 1 .
We consider permutation π 1 , 2 2 , where J k + 1 = J 3 and J k + 2 = J 2 . Again, we must test the two inequalities in Equation (12), where either s = 1 or s = 2 . For s = 1 , we obtain
u k + 1 , 1 = u 3 , 1 = 9 L 2 + i = 1 k + 1 1 ( l i 2 u i 1 ) = L 2 + i = 1 1 ( l i 2 u i 1 ) = 11 + ( l 1 , 2 u 1 , 1 ) = 11 + ( 6 7 ) = 10 .
However, for s = 2 , we obtain
u k + 2 , 1 = u 2 , 1 = 9 L 2 + i = 1 k + 2 1 ( l i 2 u i 1 ) = L 2 + i = 1 k + 1 ( l i 2 u i 1 ) = 11 + ( l 1 , 2 u 1 , 1 ) + ( l 3 , 2 u 3 , 1 )
= 11 + ( 6 7 ) + ( 5 9 ) = 6 .
Thus, the conditions of Theorem 11 do not hold for permutation π 1 , 2 2 .
Note that we do not check the conditions of Theorem 12 since the conflict set of jobs { J 2 , J 3 } is located at the end of the partial strict order A 1 , 2 . We conclude that none of the proven sufficient conditions are satisfied for a schedule optimality. Thus, there does not exist a pair of permutations of the jobs in set J = J 1 , 2 J 2 , 1 J 1 J 2 which is optimal for any scenario p T . The J-solution S ( T ) for Example 2 consists of the following two pairs of job permutations: { ( π 1 , π 1 ) , ( π 2 , π 2 ) } = S ( T ) , where
π 1 = ( π 1 , 2 1 , π 1 , π 2 , 1 ) = ( J 1 , J 2 , J 3 , J 4 , J 6 , J 7 , J 8 ) , π 1 = ( π 2 , 1 , π 2 , π 1 , 2 1 ) = ( J 6 , J 7 , J 8 , J 5 , J 1 , J 2 , J 3 ) ,
π 2 = ( π 1 , 2 2 , π 1 , π 2 , 1 ) = ( J 1 , J 3 , J 2 , J 4 , J 6 , J 7 , J 8 ) , π 2 = ( π 2 , 1 , π 2 , π 1 , 2 2 ) = ( J 6 , J 7 , J 8 , J 5 , J 1 , J 3 , J 2 ) .
We next show that none of these two pairs of job permutations is optimal for all scenarios p T using the following two scenarios: p = ( 7 , 6 , 9 , 5 , 9 , 6 , 2 , 0 , 0 , 2 , 1 , 3 , 1 , 3 , 1 , 3 ) T and p = ( 7 , 6 , 9 , 6 , 9 , 5 , 2 , 0 , 0 , 2 , 1 , 3 , 1 , 3 , 1 , 3 ) T . For scenario p , only pair of permutations ( π 2 , π 2 ) is optimal with C max ( π 2 , π 2 ) = 30 since C max ( π 1 , π 1 ) = 31 > 30 . On the other hand, for scenario p , only the pair of permutations ( π 1 , π 1 ) is optimal with C max ( π 1 , π 1 ) = 30 since C max ( π 2 , π 2 ) = 31 > 30 .
Note that the whole set S of the semi-active schedules has the cardinality | S | = m 1 , 2 ! · m 2 , 1 ! = 3 ! · 3 ! = 6 · 6 = 36 . Thus, for solving Example 2, one needs to consider only two pairs of job permutations { ( π 1 , π 1 ) , ( π 2 , π 2 ) } = S ( T ) S instead of 36 semi-active schedules.

5.4. An Algorithm for Checking Conditions for the Existence of a Single-Element Dominant Set

We describe Algorithm 1 for checking the existence of an optimal permutation for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 if the partial strict order A 1 , 2 on the set J 1 , 2 has the following form: J 1 J k { J k + 1 , J k + 2 , , J k + r } J k + r + 1 J m 1 , 2 . Algorithm 1 considers a set of conflict jobs and checks whether the sufficient conditions given in Section 5.2 hold. For a conflict set of jobs, it is needed to construct two permutations and to check the condition in Inequality (12) for the first permutation and the condition in Inequality (13) for the second one. If at least one of these conditions holds, Algorithm 1 constructs a permutation which is optimal for the problem F 2 | l i j p i j u i j | C m a x with any scenario p T .
Obviously, testing the conditions of Theorems 11 and 12 takes O ( r ) , where the conflict set contains r jobs. The construction of the permutation of r jobs takes O ( r log r ) . Therefore, the total complexity of Algorithm 1 is O ( r log r ) .
Remark 3.
If Algorithm 1 is completed at Step 7 (STOP 1), we suggest to consider a set of conflict jobs { J k + 1 , J k + 2 , , J k + r } and construct a Johnson’s permutation for the deterministic problem F 2 | p | C m a x with job set J = { J k + 1 , J k + 2 , J k + r } , where vector p = ( p k + 1 , 1 , p k + 1 , 2 , p k + r , 1 , p k + r , 2 ) of the durations of conflict jobs { J k + 1 , J k + 2 , J k + r } is calculated for each operation O i j of the conflict job J i { J k + 1 , J k + 2 , J k + r } on the corresponding machine M j M as folows:
p i j = ( u i j + l i j ) / 2
Theorem 11 and Theorem 12 imply the following claim.
Corollary 7.
Algorithm 1 constructs a permutation π * either satisfying conditions of Theorem 11 or Theorem 12 (such permutation π * is optimal for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 and any scenario p T ) or establishes that an optimal job permutation for the problem F 2 | l i j p i j u i j | C m a x with any scenario p T does not exist.
The set of jobs J 2 , 1 for the problem F 2 | l i j p i j u i j | C m a x with job set J = J 2 , 1 can be tested similarly to the set of jobs J 1 , 2 .
Algorithm 1: Checking conditions for the existence of a single-element dominant set of schedules for the problem F 2 | l i j p i j u i j | C m a x
Input: 
Segments [ l i j , u i j ] for all jobs J i J and machines M j M ,
a partial strict order A 1 , 2 on the set J 1 , 2 = J 1 , 2 * J 1 , 2 1 J 1 , 2 2 in the form
J 1 J k { J k + 1 , J k + 2 , , J k + r } J k + r + 1 J m 1 , 2 .
Output: 
EITHER an optimal job permutation for the problem
    F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 and any scenario p T , (see STOP 0)
OR there no permutation π 1 , 2 of jobs from set J 1 , 2 , which is optimal
   for all scenarios p T , (see STOP 1).
Step 1:
Set δ s = l k + s , 2 u k + s , 1 for all s { 1 , 2 , , r }
  construct a partition of the set of conflicting jobs into two subsets X 1 and X 2 ,
  where J k + s X 1 if δ s 0 , and J k + s X 2 , otherwise.
Step 2:
Construct a permutation π 1 = ( J 1 , J 2 , , J k , π 1 , π 2 , J k + r + 1 , , J m 1 , 2 ) , where the permutation
   π 1 contains jobs from the set X 1 in the non-decreasing order of the values u k + i , 1 and the
  permutation π 2 contains jobs from the set X 2 in the non-increasing order of the values
   l k + i , 2 , renumber jobs in the permutations π 1 and π 2 based on their orders.
Step 3:
IF for the permutation π 1 conditions of Theorem 11 hold THEN GOTO step 8.
Step 4:
Set δ s = l k + s , 1 u k + s , 2 for all s { 1 , 2 , , r }
  construct a partition of the set of conflicting jobs into two subsets
   Y 1 and Y 2 , where J k + s Y 1 if δ s 0, and J k + s Y 2 , otherwise.
Step 5:
Construct a permutation π 2 = ( J 1 , J 2 , , J k , π 2 , π 1 , J k + r + 1 , , J m 1 , 2 ) , where the permutation
   π 1 contains jobs from the set Y 1 in the non-increasing order of the values u k + i , 2 , and the
  permutation π 2 contains jobs from the set Y 2 in the non-decreasing order of the
  values l k + i , 1 , renumber jobs in the permutations π 1 and π 2 based on their orders.
Step 6:
IF for the permutation π 2 conditions of Theorem 12 hold THEN GOTO step 9.
Step 7:
ELSE there is no a single dominant permutation for problem
   F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 and any scenario p T STOP 1.
Step 8:
RETURN permutation π 1 , which is a single dominant permutation
  for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 STOP 0.
Step 9:
RETURN permutation π 2 , which is a single dominant permutation
  for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 STOP 0.

6. Algorithms for Constructing a Small Dominant Set of Schedules for the Problem J 2 | l i j p i j u i j , n i 2 | C m a x

In this section, we describe Algorithm 2 for constructing a small dominant set D S ( T ) of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x . Algorithm 2 is developed for use at the off-line phase of scheduling (before processing any job from the set J ). Based on the initial data, Algorithm 2 checks the conditions of Theorem 7 for a single optimal pair of job permutations for the uncertain problem J 2 | l i j p i j u i j , n i 2 | C m a x . If the sufficient conditions of Theorem 7 do not hold, Algorithm 2 proceeds to consider the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 and the problem F 2 | l i j p i j u i j | C m a x with job set J 2 , 1 . For each of these problems, the conditions of Theorem 2 are checked. If these conditions do not hold, then strict orders of the jobs J based on Inequalities (4) and (5) are constructed. In this general case, Algorithm 2 constructs a partial strict order A 1 , 2 of the jobs from set J 1 , 2 and a partial strict order A 2 , 1 of the jobs from set J 2 , 1 . Each of these partial orders may contain one or several conflict sets of jobs. For each such conflict set of jobs, Algorithm 2 checks whether the sufficient conditions given in Section 5.2 hold. Thus, if some sufficient conditions for a schedule optimality presented in Section 4 and Section 5 are satisfied, then there exists a pair of permutations of jobs from set J which is optimal for any scenario p T . Algorithm 2 constructs such a pair of job permutations { ( π , π ) } = D S ( T ) . Otherwise, the precedence digraphs determining a minimal dominant set D S ( T ) of schedules is constructed by Algorithm 2. The more job pairs are involved in the binary relations A 1 , 2 and A 2 , 1 , the more job permutations will be deleted from set S while constructing a J-solution S ( T ) S for the problems F 2 | l i j p i j u i j | C m a x with job sets J 1 , 2 and J 2 , 1 .
Algorithm 2: Construction of a small dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x
Input: 
Lower bounds l i j and upper bounds u i j , 0 < l i j u i j , of the durations
of all operations O i j of jobs J i J processed on machines M j M = { M 1 , M 2 } .
Output: 
EITHER pair of permutations ( π , π ) = ( ( π 1 , 2 , π 1 , π 2 , 1 ) , ( π 2 , 1 , π 2 , π 1 , 2 ) ) ,
   where π is a permutation of jobs from set J 1 , 2 J 1 J 2 , 1 on machine
    M 1 , π is a permutation of jobs from set J 1 , 2 J 2 J 2 , 1 on machine M 2 ,
   such that { ( π , π ) } = D S ( T ) , (see STOP 0),
OR permutation π 2 , 1 of jobs from set J 2 , 1 on machines M 1 and M 2 and
   a partial strict order A 1 , 2 of jobs from set J 1 , 2 ,
OR permutation π 1 , 2 of jobs from set J 1 , 2 on machines M 1 and M 2 and
   a partial strict order A 2 , 1 of jobs from set J 2 , 1 ,
OR a partial strict order A 1 , 2 of jobs from set J 1 , 2 and
   a partial strict order A 2 , 1 of jobs from set J 2 , 1 , (see STOP 1).
Step 1:
Determine a partition J = J 1 J 2 J 1 , 2 J 2 , 1 of the job set J ,
  permutation π 1 of jobs from set J 1 and permutation π 2 of jobs from
  set J 2 , arrange the jobs in the increasing order of their indexes.
Step 2:
IF the first inequality in condition (7) of Theorem 7 holds THEN BEGIN
  Construct a permutation π 1 , 2 of jobs from set J 1 , 2 ,
    arrange them in the increasing order of their indexes;
  IF the second inequality in condition (7) of Theorem 7 holds
    THEN construct a permutation π 2 , 1 of jobs from set J 2 , 1 ,
    arrange them in the increasing order of their indexes GOTO Step 10 END
Step 3:
IF the first inequality in condition (8) of Theorem 7 holds THEN BEGIN
  Construct a permutation π 2 , 1 of jobs from set J 2 , 1 ,
    arrange them in the increasing order of their indexes;
  IF the second inequality in condition (8) of Theorem 7 holds THEN
    construct a permutation π 1 , 2 of jobs from set J 1 , 2 ,
    arrange the jobs in the increasing order of their indexes END
Step 4:
IF both permutations π 1 , 2 and π 2 , 1 are constructed THEN GOTO Step 10.
Step 5:
IF permutation π 1 , 2 is not constructed THEN fulfill Algorithm 3.
Step 6:
IF permutation π 2 , 1 is not constructed THEN fulfill Algorithm 4.
Step 7:
IF both permutations π 1 , 2 and π 2 , 1 are constructed THEN GOTO Step 10.
Step 8:
IF permutation π 2 , 1 is constructed THEN GOTO Step 11.
Step 9:
IF permutation π 1 , 2 is constructed THEN GOTO Step 12 ELSE GOTO Step 13.
Step 10:
RETURN pair of permutations ( π , π ) , where π is the permutation
  of jobs from set J 1 , 2 J 1 J 2 , 1 processed on machine M 1 and π is
  the permutation of jobs from set J 1 , 2 J 2 J 2 , 1 processed
  on machine M 2 such that { ( π , π ) } = D S ( T ) STOP 0.
Step 11:
RETURN the permutation π 2 , 1 of jobs from set J 2 , 1 processed on machines M 1 and M 2 ,
  the partial strict order A 1 , 2 of jobs from set J 1 , 2 GOTO Step 14.
Step 12:
RETURN the permutation π 1 , 2 of jobs from set J 1 , 2 processed on machines M 1 and M 2 ,
  the partial strict order A 2 , 1 of jobs from set J 2 , 1 GOTO Step 14.
Step 13:
RETURN the partial strict order A 1 , 2 of jobs from set J 1 , 2
  and the partial strict order A 2 , 1 of jobs from set J 2 , 1
Step 14:
STOP 1.
Algorithm 2 may be applied for solving the problem J 2 | l i j p i j u i j , n i 2 | C m a x exactly or approximately as follows. If at least one of the sufficient conditions proven in Section 5.1 hold, then Algorithm 2 constructs a pair of job permutations ( π , π ) = ( ( π 1 , 2 , π 1 , π 2 , 1 ) , ( π 2 , 1 , π 2 , π 1 , 2 ) ) , which is optimal for any scenario p T (Step 10).
It may happen that the constructed strict order on the set J 1 , 2 or on the set J 2 , 1 is not a linear strict order. If for at least one of the sets J 1 , 2 or J 2 , 1 , the constructed partial strict order is not a linear one, a heuristic solution for the problem J 2 | l i j p i j u i j , n i 2 | C m a x is constructed similar to that for the problem F 2 | l i j p i j u i j | C m a x solved by Algorithm 1 (see Section 5.4). If Algorithm 2 is completed at Steps 11-13 (STOP 1), we consider a set of conflict jobs { J k + 1 , J k + 2 , , J k + r } and construct a Jackson’s pair of job permutation for the deterministic problem J 2 | p , n i 2 | C m a x with job set J = { J k + 1 , J k + 2 , J k + r } , where the vector p = ( p k + 1 , 1 , p k + 1 , 2 , p k + r , 1 , p k + r , 2 ) of the durations of conflict jobs { J k + 1 , J k + 2 , J k + r } is calculated using the equality of Equation (16) for each operation O i j of the conflict job J i { J k + 1 , J k + 2 , J k + r } on the corresponding machine M j M (Remark 3).
Algorithm 3: Construction of a strict order A 1 , 2 on the set J 1 , 2
Input: 
Lower bounds l i j and upper bounds u i j , 0 < l i j u i j , of the durations
of all operations O i j of jobs J i J on machines M j M = { M 1 , M 2 } .
Output: 
EITHER permutation π 1 , 2 , which is optimal for the problem
   F 2 | l i j p i j u i j | C m a x with any scenario p T for the jobs J 1 , 2 ,
OR partial strict order A 1 , 2 on the set J 1 , 2 .
Step 1:
Construct a partition J 1 , 2 = J 1 , 2 1 J 1 , 2 2 J 1 , 2 * of the set J 1 , 2 of the jobs.
Step 2:
IF conditions of Theorem 2 hold THEN
Step 3:
Construct permutation π 1 , 2 = ( π 1 , 2 1 , J 1 , 2 * , π 1 , 2 2 ) , where π 1 , 2 1 is a permutation for
    processing jobs from the set J 1 , 2 1 in the non-decreasing order of the values u i 1 ,
     π 1 , 2 2 is a permutation for processing jobs from the set J 1 , 2 2
    in the non-increasing order of the values u i 2 GOTO Step 7 ELSE
Step 4:
FOR each pair of jobs J v J 1 , 2 and J w J 1 , 2 , v w , DO
    IF at least one of two conditions (4) and (5) holds THEN
    determine the relation J v J w
END FOR
Step 5:
Renumber jobs in the set J 1 , 2 such that relation v < w holds if J v J w .
Step 6:
FOR each conflict set of jobs DO
    IF condition of Theorem 10 holds THEN
    Order jobs in the conflict set in the increasing order of their indexes GOTO Step 7
    ELSE fulfill Algorithm 1
END FOR
Step 7:
IF the partial strict order A 1 , 2 is linear THEN
    construct a permutation π 1 , 2 generated by the linear order A 1 , 2
STOP.
Algorithm 4 is obtained from the above Algorithm 3 by replacing the set J 1 , 2 of jobs by the set J 2 , 1 of jobs, machine M 1 by machine M 2 , and vice versa. Obviously, testing the conditions of Theorems 11 and 12 takes O ( r ) , where conflict set contains r jobs. Construction of permutation of r jobs takes O ( r log r ) . Therefore, the total complexity of Algorithm 1 is O ( r log r ) .
Testing the conditions of Theorem 2 takes O ( m 1 , 2 log m 1 , 2 ) time. A strict order A 1 , 2 on the set J 1 , 2 is constructed by comparing no more than m 1 , 2 ( m 1 , 2 1 ) pairs of jobs in the set J 1 , 2 . Thus, it takes O ( m 1 , 2 ( m 1 , 2 1 ) ) time. The complexity of Algorithm 1 is O ( r log r ) time provided that the conflict set contains r jobs, where r m 1 , 2 . Since a strict order A 1 , 2 is constructed once in Algorithm 3, we conclude that a total complexity of Algorithm 3 (and Algorithm 4) is O ( n 2 ) time.
In Algorithm 2, testing the condition of Theorem 7 takes O ( max { m 1 , 2 , m 2 , 1 } ) time. Every Algorithm 3 or Algorithm 4 is fulfilled at most once. Therefore, the complexity of Algorithm 2 is O ( n 2 ) time.

7. Computational Experiments

We describe the conducted computational experiments and discuss the results obtained for randomly generated instances of the problem J 2 | l i j p i j u i j , n i 2 | C m a x . In the computational experiments, each tested series consisted of 1000 randomly generated instances with the same numbers n { 10 , 20 , , 100 , 200 , , 1000 , 2000 , , 10.000 } of jobs in the set J provided that a maximum relative length δ of the given segment of the possible durations of the operations O i j takes the following values: { 5 % , 10 % , 15 % , 20 % , 30 % , 40 % , and 50 % } . The lower bounds l i j and upper bounds u i j for possible values of the durations p i j of the operations O i j , p i j [ l i j , u i j ] using the value δ have been determined as follows. First, a value of the lower bound l i j is randomly chosen from the segment [ 10 , 1000 ] using a uniform distribution. Then, the upper bound u i j is calculated using the following equality:
u i j = l i j 1 + δ 100
For example, we assume that δ = 5 % . Then, for the lower bounds l i j = 50 and l i j = 500 , the upper bounds u i j = 52.5 and u i j = 525 are calculated using Reference (17). If δ = 50 % , then based on the lower bounds l i j = 50 and l i j = 500 and on Reference (17), we obtain the upper bounds u i j = 75 and u i j = 750 . Thus, rather wide ranges for the tested durations of the jobs J were considered.
In the experiments, the bounds l i j and u i j were decimal fractions with the maximum possible number of digits after the decimal point. For all tested instances of the problem J 2 | l i j p i j u i j , n i 2 | C m a x , a strict inequality l i j < u i j was guarantied for each job J i J and each machine M j M .
We used Algorithms 1 – 4 described in Section 5.4 and Section 6 for solving the problem J 2 | l i j p i j u i j , n i 2 | C m a x . These algorithms were coded in C# and tested on a PC with Intel Core i7-7700 (TM) 4 Quad, 3.6 GHz, and 32.00 GB RAM. Since Algorithms 1 – 4 are polynomial in number n jobs in set J , the calculations were carried out quickly. In the experiments, we tested 15 classes of randomly generated instances of the problem J 2 | l i j p i j u i j , n i 2 | C m a x with different ratios between numbers m 1 , m 2 , m 1 , 2 , and m 2 , 1 of the jobs in subsets J 1 , J 2 , J 1 , 2 , and J 2 , 1 of the set J . The obtained computational results are presented in Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8, Table A9, Table A10, Table A11, Table A12, Table A13, Table A14 and Table A15 for 15 classes of the solved instances. Each tested class of the instances of the problem J 2 | l i j p i j u i j , n i 2 | C m a x is characterized by the following ratio of the percentages of the number of jobs in the subsets J 1 , J 2 , J 1 , 2 , and J 2 , 1 of the set J :
m 1 n · 100 % : m 2 n · 100 % : m 1 , 2 n · 100 % : m 2 , 1 n · 100 %
Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8 and Table A9 present the computational results obtained for classes 1–9 of the tested instances characterized by the following ratios (Equation (18)):
25 % : 25 % : 25 % : 25 % (Table A1); 10 % : 10 % : 40 % : 40 % (Table A2);
10 % : 40 % : 10 % : 40 % (Table A3); 10 % : 30 % : 10 % : 50 % (Table A4);
10 % : 20 % : 10 % : 60 % (Table A5); 10 % : 10 % : 10 % : 70 % (Table A6);
5 % : 20 % : 5 % : 70 % (Table A7); 5 % : 15 % : 5 % : 75 % (Table A8);
5 % : 5 % : 5 % : 85 % (Table A9).
Note that all instances from class 1 of the instances with the ratio from Equation (18), 25 % : 25 % : 25 % : 25 % , were optimally solved by Algorithms 1–4 for all values of δ { 5 % , 10 % , 15 % , 20 % , 30 % , 40 % , and 50 % } . We also tested classes 10–15 of the hard instances of the problem J 2 | l i j p i j u i j , n i 2 | C m a x characterized by the following ratios (Equation (18)):
3 % : 2 % : 5 % : 90 % (Table A10); 2 % : 3 % : 5 % : 90 % (Table A11);
2 % : 2 % : 1 % : 95 % (Table A12); 1 % : 2 % : 2 % : 95 % (Table A13);
1 % : 1 % : 3 % : 95 % (Table A14); 1 % : 1 % : 1 % : 97 % (Table A15).
All Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8, Table A9, Table A10, Table A11, Table A12, Table A13, Table A14 and Table A15 are organized as follows. Number n of given jobs J in the instances of the problem J 2 | l i j p i j u i j , n i 2 | C m a x are presented in column 1. The values of δ (a maximum relative length of the given segment of the job durations) in percentages are presented in the first line of each table. For the fixed value of δ , the obtained computational results are presented in four columns called O p t , N C , S C , and t. The column O p t determines the percentage of instances from the series of 1000 randomly generated instances which were optimally solved using Algorithms 1–4. For each such instance, an optimal pair ( π , π ) of the job permutations was constructed in spite of the uncertain durations of the given jobs J . In other words, the equality C m a x ( π , π ) = C m a x ( π * , π * * ) holds, where ( π * , π * * ) S is a pair of job permutations which is optimal for the deterministic problem J 2 | p * , n i 2 | C m a x associated with the factual scenario p * T . The factual scenario p * T for the instance of the uncertain problem J 2 | l i j p i j u i j , n i 2 | C m a x is assumed to be unknown until completing the jobs J .
Column N C presents total number of conflict sets of the jobs in the partial strict orders A 1 , 2 on the job sets J 1 , 2 and partial strict orders A 2 , 1 on the job sets J 2 , 1 constructed by Algorithm 2. The value of N C is equal to the total number of decision points, where Algorithm 2 has to select an order for processing jobs from the corresponding conflict set. To make a correct decision for such an order means to construct a permutation of all jobs from the conflict set, which is optimal for the factual scenario (which is unknown before scheduling). In particular, if all conflict sets have received correct decisions in Algorithm 2, then the constructed pair of job permutations will be optimal for the problem J 2 | p * , n i 2 | C m a x , where p * T is the factual scenario.
Column S C presents a percentage of the correct decisions made for determining optimal orders of the conflict jobs by Algorithm 2 with Algorithms 3 and 4. Column t presents a total CPU time (in seconds) for solving all 1000 instances of the corresponding series.
Average percentages of the instances which were optimally solved ( O p t ) are presented in Figure 1 for classes 1–9 of the tested instances and in Figure 2 for classes 10–15 of the hard-tested instances.
Percentages of the average values of the correct decisions ( S C ) made for determining optimal orders of the conflict jobs for classes 1–9 are presented in Figure 3. Most instances from these nine classes were optimally solved (Table 2). If the values of δ were no greater than 20 % , i.e., δ { 5 % , 10 % , 15 % , 20 % } , then more than 80 % of the tested instances were optimally solved in spite of the data uncertainty. If the value δ is increased, the percentage of the optimally solved instances decreased. If the value δ was equal to 50 % , then 45 % of the tested instances was optimally solved.
For all series of the hard instances presented in Table A10, Table A11, Table A12, Table A13, Table A14 and Table A15 (see the third line in Table 2), only a few instances were optimally solved. If δ = 5 % , then 70 % of the tested instances was optimally solved. If value δ belongs to the set { 20 % , 30 % , 40 % , 50 % } , then only 1 % of the tested instances was optimally solved. There were no hard-tested instances optimally solved for the value of δ = 50 % .
Percentages of the average values of the correct decisions made for determining optimal orders of the conflict jobs by Algorithm 2, Algorithm 3 and Algorithm 4 for the hard classes 10–15 of the tested instances are presented in Figure 4. Note that there is a correlation between values of O p t and S C presented in Figure 1 and Figure 3 for classes 1–9 of the tested instances and those presented in Figure 2 and Figure 4 for classes 10–15 of the hard-tested instances.

8. Concluding Remarks and Future Works

The uncertain flow-shop scheduling problem F 2 | l i j p i j u i j | C m a x and its generalization the job-shop problem J 2 | l i j p i j u i j , n i 2 | C m a x attract the attention of researchers since these problems are applicable in many real-life scheduling systems. The optimal scheduling decisions for these problems allow the plant to reduce the costs of productions due to a better utilization of the available machines and other resources. In Section 5, we proved several properties of the optimal pairs ( π , π ) of job permutations (Theorems 7–12). Using these properties, we derived Algorithms 1–4 for constructing optimal pairs ( π , π ) of job permutations or a small dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x . If it is impossible to construct a single pair ( π , π ) of job permutations, which dominates all other pairs of job permutations for all possible scenarios T, then Algorithm 2 determines the partial strict order A 1 , 2 on the job set J 1 , 2 (Algorithm 3) and the partial strict order A 2 , 1 on the job set J 2 , 1 (Algorithm 4). The precedence digraphs ( J 1 , 2 , A 1 , 2 ) and ( J 2 , 1 , A 2 , 1 ) determine a minimal dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x .
From the conducted extensive computational experiments, it follows that pairs of job permutations constructed using Algorithm 2 are close to the optimal pairs of job permutations, which may be determined after completing all jobs J when factual operation durations become known. We tested 15 classes of the randomly generated instances J 2 | l i j p i j u i j , n i 2 | C m a x . Most instances from tested classes 1–9 were optimally solved at the off-line phase of scheduling. If the values of δ were no greater than 20 % , i.e., δ { 5 % , 10 % , 15 % , 20 % } , then more than 80 % of the tested instances was optimally solved in spite of the uncertainty of the input data. If δ = 50 % , then 45 % of the tested instances was optimally solved. However, less than 5 % of the instances with δ 20 % from hard classes 10–15 were optimally solved at the off-line phase of scheduling (Figure 2). There were no tested hard instances optimally solved for the value δ = 50 % .
In future research, the on-line phase of scheduling will be studied for the problem J 2 | l i j p i j u i j , n i 2 | C m a x . To this end, it will be useful to find sufficient conditions for existing a dominant pair of job permutations at the on-line phase of scheduling. The additional information on the factual value of the job duration becomes available once the processing of the job on the corresponding machine is completed. Using this additional information, a scheduler can determine a smaller dominant set DS of schedules, which is based on sufficient conditions for schedule dominance. The smaller DS enables a scheduler to quickly make an on-line scheduling decision whenever additional information on processing the job becomes available. To solve the problem J 2 | l i j p i j u i j , n i 2 | C m a x at the on-line phase, a scheduler needs to use fast (better polynomial) algorithms. The investigation of the on-line phase of scheduling for the uncertain job-shop problem is under development.
We suggest to investigate properties of the optimality box and optimality region for a pair ( π , π ) of the job permutations and to develop algorithms for constructing a pair ( π , π ) of the job permutations that have the largest optimality box (or the largest optimality region). We also suggest to apply the stability approach for solving the uncertain flow-shop and job-shop scheduling problems with | M | > 2 available machines.

Author Contributions

Methodology, Y.N.S.; software, V.D.H.; validation, Y.N.S., N.M.M. and V.D.H.; formal analysis, Y.N.S. and N.M.M.; investigation, Y.N.S. and N.M.M; writing—original draft preparation, Y.N.S. and N.M.M.; writing—review and editing, Y.N.S.; visualization, N.M.M.; supervision, Y.N.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We are thankful for useful remarks and suggestions provided by the editors and three anonymous reviewers on the earlier draft of our paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proofs of the Statements

Appendix A.1. Proof of Lemma 2

We choose an arbitrary vector p in the set T, p T , and show that set < S 1 , 2 , S 2 , 1 > contains at least one optimal pair of job permutations for the problem J 2 | p , n i 2 | C m a x with scenario p T .
Let ( π * , π * * ) = ( ( π 1 , 2 * , π 1 , π 2 , 1 * ) , ( π 2 , 1 * , π 2 , π 1 , 2 * ) ) be a Jackson’s pair of job permutations for the problem J 2 | p , n i 2 | C m a x with scenario p T , i.e., C m a x ( π * , π * * ) = C m a x . Without loss of generality, one can assume that jobs in both permutations π 1 and π 2 are ordered in increasing order of their indexes. It is clear that π 2 , 1 * S 2 , 1 . If inclusion π 1 , 2 * S 1 , 2 holds as well, then ( π * , π * * ) < S 1 , 2 , S 2 , 1 > and set < S 1 , 2 , S 2 , 1 > contains an optimal pair of job permutations for the problem J 2 | p , n i 2 | C m a x with scenario p T . We now assume that π 1 , 2 * S 1 , 2 . The set S 1 , 2 contains at least one optimal permutation for the problem F 2 | p 1 , 2 | C m a x with job set J 1 , 2 and scenario p 1 , 2 (the components of vector p 1 , 2 are equal to the corresponding components of vector p). We denote this permutation as π 1 , 2 . Remember that permutation π 1 , 2 may be not a Johnson’s permutation for the problem F 2 | p 1 , 2 | C m a x with job set J 1 , 2 and scenario p 1 , 2 . We consider a pair of job permutations ( π , π * * ) = ( ( π 1 , 2 , π 1 , π 2 , 1 * ) , ( π 2 , 1 * , π 2 , π 1 , 2 ) ) < S 1 , 2 , S 2 , 1 > and show that equality C m a x ( π , π * * ) = C m a x holds. We consider the following two possible cases.
(j) C m a x ( π , π * * ) = c 1 ( π ) .
If equality c 1 ( π ) = J i J 1 , 2 J 2 , 1 J 1 p i 1 holds, then c 1 ( π ) c 1 ( π * ) .
We now assume that inequality c 1 ( π ) > J i J 1 , 2 J 2 , 1 J 1 p i 1 holds. Then, machine M 1 has an idle time. As it is mentioned in the proof of Theorem 7, an idle time for machine M 1 is only possible if some job J j from the set J 2 , 1 is processed on machine M 2 at the time moment t 2 when job J j could be processed on machine M 1 . Thus, the value of c 1 ( π ) is equal to the makespan C m a x ( π 2 , 1 * ) for the problem F 2 | p 2 , 1 | C m a x with job set J 2 , 1 and scenario p 2 , 1 (the components of vector p 2 , 1 are equal to the corresponding components of vector p). As jobs from the set J 2 , 1 are processed as in the permutation π 2 , 1 * , which is a Johnson’s permutation, the value of c 1 ( π ) cannot be reduced and so c 1 ( π ) c 1 ( π * ) . We obtain the following relations: C m a x ( π , π * * ) = c 1 ( π ) c 1 ( π * ) max { c 1 ( π * ) , c 2 ( π * * ) } = C m a x ( π * , π * * ) = C m a x . Thus, equality C m a x ( π , π * * ) = C m a x holds.
(jj) C m a x ( π , π * * ) = c 2 ( π * * ) .
Similarly to case (j), we obtain the following equality:
c 2 ( π * * ) = max { J i J 1 , 2 J 2 , 1 J 2 p i 2 , C m a x ( π 1 , 2 ) } ,
where C m a x ( π 1 , 2 ) is the makespan for the problem F 2 | p 1 , 2 | C m a x with job set J 1 , 2 and vector p 1 , 2 of the job durations (it is assumed that π 1 , 2 is an optimal permutation for this problem). Thus, the value of c 2 ( π * * ) cannot be reduced and equality C m a x ( π , π * * ) = C m a x holds.
In both considered cases, the pair of job permutations ( π , π * * ) is an optimal schedule for the problem J 2 | p , n i 2 | C m a x with scenario p T . Therefore, an optimal pair of job permutations for the problem J 2 | p , n i 2 | C m a x with scenario p T belongs to the set < S 1 , 2 , S 2 , 1 > . As vector p is an arbitrary vector in set T, the set < S 1 , 2 , S 2 , 1 > contains an optimal pair of job permutations for each scenario from set T. Due to Definition 4, the set < S 1 , 2 , S 2 , 1 > is a dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J .

Appendix A.2. Proof of Theorem 8

We consider an arbitrary vector p T of the job durations from set T and relevant vectors p 1 , 2 and p 2 , 1 of the durations of jobs from set J 1 , 2 and set J 2 , 1 , respectively. Set S 1 , 2 contains an optimal permutation π 1 , 2 for the problem F 2 | p 1 , 2 | C m a x with job set J 1 , 2 and with vector p 1 , 2 of the job durations. Set S 2 , 1 contains an optimal permutation π 2 , 1 for the problem F 2 | p 2 , 1 | C m a x with job set J 2 , 1 and with vector p 2 , 1 of the job durations. We next show that the pair of job permutations ( π , π ) = ( ( π 1 , 2 , π 1 , π 2 , 1 ) , ( π 2 , 1 , π 2 , π 1 , 2 ) ) is an optimal pair of job permutations for the problem J 2 | p , n i 2 | C m a x with scenario p T (the jobs in permutations π 1 and π 2 are ordered in increasing order of their indexes). From the proofs of Lemmas 2 and 3, we obtain the value of C m a x ( π , π ) = max { c 1 ( π ) , c 2 ( π ) }
= max { max { J i J 1 , 2 J 2 , 1 J 1 p i 1 , C m a x ( π 2 , 1 ) } , max { J i J 1 , 2 J 2 , 1 J 2 p i 2 , C m a x ( π 1 , 2 ) } } ,
which cannot be reduced. Therefore, C m a x ( π , π ) = C m a x . An optimal pair of job permutations for the problem J 2 | p , n i 2 | C m a x with vector p T of the job durations belongs to the set < S 1 , 2 , S 2 , 1 > . As vector p is arbitrary in set T, the set < S 1 , 2 , S 2 , 1 > contains an optimal pair of job permutations for all vectors from set T. Due to Definition 4, the set < S 1 , 2 , S 2 , 1 > S is a dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J .

Appendix A.3. Proof of Theorem 9

We consider an arbitrary scenario p T . Due to Definition 1, the permutation π 1 , 2 is a Johnson’s permutation for the problem F 2 | p 1 , 2 | C m a x with job set J 1 , 2 and scenario p 1 , 2 (the components of this vector are equal to the corresponding components of vector p). Due to Definition 4, the singleton { ( π 1 , 2 , π 1 , 2 ) } is a minimal dominant set of schedules for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 .
Similarly, the singleton { ( π 2 , 1 , π 2 , 1 ) } is a minimal dominant set of schedules for the problem F 2 | l i j p i j u i j | C m a x with job set J 2 , 1 . We consider permutations π 1 and π 2 of the jobs J 1 and J 2 , respectively (due to Remark 1, the jobs in permutations π 1 and π 2 are ordered in increasing order of their indexes). Due to Theorem 8, the pair of permutations ( ( π 1 , 2 , π 1 , π 2 , 1 ) , ( π 1 , 2 , π 2 , π 2 , 1 ) ) is a single-element dominant set (DS(T)) for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J = J 1 J 2 J 1 , 2 J 2 , 1 .

Appendix A.4. Proof of Corollary 6

In the proof of Theorem 9, it is shown that the pair of job permutations ( ( π 1 , 2 , π 1 , π 2 , 1 ) , ( π 1 , 2 , π 2 , π 2 , 1 ) ) is a single-element dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J = J 1 J 2 J 1 , 2 J 2 , 1 . We next show that the pair of permutations ( ( π 1 , 2 , π 1 , π 2 , 1 ) , ( π 1 , 2 , π 2 , π 2 , 1 ) ) satisfies to Definition 1, i.e., this pair of permutations is a Jackson’s pair of job permutations for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J (the minimality condition is obvious). Indeed, due to conditions of Theorem 9, the permutation π 1 , 2 is a Johnson’s permutation for the problem F 2 | l i j p i j u i j | C m a x with job set J 1 , 2 and the permutation π 2 , 1 is a Johnson’s permutation for the problem F 2 | l i j p i j u i j | C m a x with job set J 2 , 1 . Therefore, pair ( ( π 1 , 2 , π 1 , π 2 , 1 ) , ( π 1 , 2 , π 2 , π 2 , 1 ) ) is a Jackson’s pair of permutations for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J . Due to Definition 1, the pair of job permutations ( ( π 1 , 2 , π 1 , π 2 , 1 ) , ( π 1 , 2 , π 2 , π 2 , 1 ) ) is a single-element J-solution for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J = J 1 J 2 J 1 , 2 J 2 , 1 .

Appendix A.5. Proof of Theorem 12

We consider any fixed scenario p T and a pair of job permutations ( π , π ) = ( ( π , π 1 , π 2 , 1 * ) , ( π 2 , 1 * , π 2 , π ) ) S , where π 2 , 1 * S 2 , 1 is a Johnson’s permutation of the jobs from the set J 2 , 1 with vector p 2 , 1 of the job durations (components of this vector are equal to the corresponding components of vector p). We next show that this pair of job permutations ( π , π ) is optimal for the individual problem J 2 | p , n i 2 | C m a x with scenario p, i.e., C m a x ( π , π ) = C m a x .
At time t = 0 , machine M 1 begins to process jobs from the permutation π without idle times. We denote t 1 = i = 1 k + r + 1 p i 1 . At time moment t 1 , job J k + r + 1 is ready for processing on machine M 2 . From the condition of Inequality (13) with s = 1 , it follows that, even if machine M 2 has an idle time before processing job J k + r + 1 , machine M 2 is available for processing this job at time t 1 . If in addition, the condition of Inequality (13) holds with s { 2 , 3 , , r } , then machine M 2 may also have idle times between processing jobs from the conflict set { J k + 1 , J k + 2 , , J k + r } . However, machine M 2 is available for processing job J k + r + 1 from the time moment t 1 = i = 1 k + r + 1 p i 1 .
In permutation π , jobs J k + r + 1 , , J m 1 , 2 are arranged in Johnson’s order. Therefore, if machine M 2 has an idle time while processing these jobs, this idle time cannot be reduced.
Thus, the value of c 2 ( π ) cannot be reduced by changing the order of jobs from the conflict set. Note that an idle time for machine M 1 is only possible between some jobs from the set J 2 , 1 . Since the permutation π 2 , 1 * is a Johnson’s permutation of the jobs from set J 2 , 1 with scenario p 2 , 1 , the value of c 1 ( π ) cannot be reduced. Thus, we obtain C m a x ( π , π ) = max { c 1 ( π ) , c 2 ( π ) } = C m a x and the pair of permutations ( π , π ) = ( ( π , π 1 , π 2 , 1 * ) , ( π 2 , 1 * , π 2 , π ) ) S is optimal for the problem J 2 | p , n i 2 | C m a x with scenario p T . As the vector p is an arbitrary vector in the set T, set S contains an optimal pair of job permutations for each vector from the set T. Due to Definition 4, set S is a dominant set of schedules for the problem J 2 | l i j p i j u i j , n i 2 | C m a x with job set J .

Appendix B. Tables with Computations Results

Table A1. Computational results for randomly generated instances with ratio 25 % : 25 % : 25 % : 25 % of the number of jobs in the subsets.
Table A1. Computational results for randomly generated instances with ratio 25 % : 25 % : 25 % : 25 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
2010061000100191000100351000100701000100139100010025010001003391000
401000-010041000100201000100331000100101100010013610001003331000
5010071000100310001001610001008100010050100010011410001002241000
701000-01000-010021000100310001001110001007110001001491000
801000-01000-01000-0100310001000-01003510001001221000
1001000-01000-01000-01000-01000-0100191000100841000
2001000-01000-01000-01000-01000-01000-010051000
3001000-01000-01000-01000-01000-01000-01000-0
4001000-01000-01000-01000-01000-01000-01000-0
5001000-01000-01000-01000-01000-01000-01000-0
6001000-01000-01000-01000-01000-01000-01000-0
7001000-01000-01000-01000-01000-01000-01000-0
8001000-01000-01000-01000-01000-01000-01000-0
9001000-01000-01000-01000-01000-01000-01000-0
10001000-01000-01000-01000-01000-01000-01000-0
20001000-01000-01000-01000-01000-01000-01000-0
30001000-01000-01000-01000-01000-01000-01000-0
40001000-01000-01000-01000-01000-01000-01000-0
50001000-01000-01000-01000-01000-01000-01000-0
60001000-01000-01000-01000-01000-01000-01000-0
70001000-01000-01000-01000-01000-01000-01000-0
80001000-01000-01000-01000-01000-01000-01000-0
90001000-01000-01000-01000-01000-01000-01000-0
10,0001000-01000-01000-01000-01000-01000-01000-0
Aver.1000.5410001001.0810001003.0410001004.88100010012.54100010026.04100010052.331000
Table A2. Computational results for randomly generated instances with ratio 10 % : 10 % : 40 % : 40 % of the number of jobs in the subsets.
Table A2. Computational results for randomly generated instances with ratio 10 % : 10 % : 40 % : 40 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
10100711000100153100099.624198.34099.432098.13098.148196.05095.861893.20091.971388.640
20100235100010053110001008111000100103210001001341100099.9145099.93099.5142499.650
301004601000100887100010013341000100164310001001912100099.9189399.95010018081000
401006361000100118510001001659100010020681000100235210001002162100099.9195399.950
501008241000100149610001002074100010024111000100254610001002211100010020091000
601008931000100154210001002222100010026191000100275810001002440100010021091000
701008411000100158910001002285100010027751000100293510001002477100010021061000
801009811000100157010001002342100010028961000100299510001002567100010022491000
901008781000100166010001002310100010029051000100310310001002598100010022731000
1001008261000100163310001002368100010030561000100311410001002585100010023211000
2001004111000100114510001001999100010030651000100339210001002709100010022501000
300100181100010072110001001708100010028881000100336510001002579100010021171000
4001005110001003021000100981100010024661000100326310001002469100010019661000
5001001110001002401000100813100010023071000100313810001002362100010018381000
6001000-0100881000100499100010020761000100295110001002202100010016921000
7001000-0100451000100528100010018941000100277910011002015100110015851001
8001000-0100361000100294100010017071000100265610001001866100110014851001
9001000-01000-0100318100010014421000100239210011001677100110014201001
10001000-01000-0100196100010012751000100225510011001630100110012981001
20001000-01000-0100310001004411000100145210031001137100310010441002
30001000-01000-01000-01001601000100112710061001025100510010111004
40001000-01000-01000-0100861000100103210091001005100810010001007
50001000-01000-01000-0100341000100101110014100100010012100100010010
60001000-01000-01000-0100231000100100210021100100110017100100110014
70001000-01000-01000-010081000100100010028100100010023100100010019
80001000-01000-01000-010061000100100110037100100010031100100010025
90001000-01000-01000-010031000100100010048100100010039100100010032
10,0001000-01000-01000-010041001100100010061100100010049100100010040
Aver.1002611000100529100099.9989299.92099.98148699.930.0499.93212099.868.2199.84177499.756.8299.69156099.585.61
Table A3. Computational results for randomly generated instances with ratio 10 % : 40 % : 10 % : 40 % of the number of jobs in the subsets.
Table A3. Computational results for randomly generated instances with ratio 10 % : 40 % : 10 % : 40 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
1098.710487.5097.320786.96096.928789.20094.535984.68092.554386.19088.559780.74081.571374.050
2099.746699.36099.683099.52099.3105799.24098.5124498.79095.3143396.65091146793.73084.6150189.610
301001070100099.8159799.87099.7187999.84098.2200799.10096.7205398.30091.4194595.58082.1179089.830
40100170010001002355100099.8266099.92099.7262899.89097.6246299.03092.4213596.35085.5197992.520
501002425100099.7317499.91099.9319799.97099.8304899.93097.8266499.17092.5228396.54082.2207591.420
6010032181000100380810001003702100099.9339499.97098.5288199.41091.9244296.60080.5217290.750
70100391110001004385100099.9406399.98099.9364899.97098.7295999.56091.9251796.62078214689.470
80100481710001004902100099.8437099.95099.9382999.97098.5310399.52092.2262796.99077.3225789.770
901005518100010053981000100465610001003910100097.7315499.21092.5267597.12079.5224990.840
100100619510001005697100099.9485399.9801004047100098.2320799.44091.3270696.78075.6232889.480
20010010,620100010076081000100571710001004645100098.9332099.64090.8271796.61072.7228187.900
30010013,110100010082591000100607010001004782100099.5336999.85094.3260597.81074.5211787.810
40010014,309100110086341000100611310001004650100099.8324799.94094.3246097.68073.1200286.560
50010014,935100010086581000100610210001004630100099.9313799.97095.1229797.87078.5180888.110
60010015,78010001008832100010060211000100449210001002911100097.2215398.70077.8170586.980
70010015,97110001008753100010057891000100437910001002786100097.7199698.85082.4161389.090
80010016,43910001008806100010057931000100417610001002533100098.8184699.35084148789.240
90010016,26810001008574100110056081001100400510011002379100098.8171799.30089.1136692.020
100010016,61410011008419100110054001001100380710011002279100199.6165599.76190.9130293.010
200010015,5391002100690610021003715100210024221002100140110011001135100198.4104098.461
300010013,8841004100525910041002599100410016241003100110910031001021100399.8100699.803
400010012,3021007100391110071001874100610012911006100104410051001004100499.8100199.804
500010010,4211001310029351001110014851001010011261009100100810081001000100610010001005
60001008822100171002299100161001262100141001043100131001004100101001000100910010001008
7000100742610024100185510022100114510020100102610017100100010014100100110011100100010010
8000100634610033100156910030100108410026100100710023100100010018100100010015100100010013
9000100537810042100136210038100103810033100100210030100100010024100100010019100100010017
10,000100452910054100123710048100102810042100100010038100100010029100100010024100100010020
Aver.99.94886199.537.0799.87486599.516.4399.83352099.575.6899.66282999.375.1198.91214299.144.0495.79178697.613.3286.71156992.382.89
Table A4. Computational results for randomly generated instances with ratio 10 % : 30 % : 10 % : 50 % of the number of jobs in the subsets.
Table A4. Computational results for randomly generated instances with ratio 10 % : 30 % : 10 % : 50 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
1099.419496.91097.433492.22095.547490.51094.356689.93088.374984.38079.184975.03074.193471.950
2099.976799.87099.4123299.51098.4148998.86097.5168098.39094.7176296.99085.6177091.64075.1172485.150
3099.7153299.80099.5215099.77099.2239999.67099249999.52096236798.31086.8220093.77071.7194185.270
4099.8242299.92099.7315199.90099.8329599.94098.9297299.63095.2258198.06083232692.48069.5205585.060
501003422100099.9401399.98099.9384499.97099.5345199.86095.3285798.35082.1248692.76064.2220283.470
60100442510001004681100099.9418999.98099.4375099.84094.8298198.26083.7256693.53064223883.600
7010053381000100518110001004569100099.3402799.83094.3318398.21079.6259492.14061.1228482.750
80100616910001005770100099.9491599.98099.6411299.90095.6325798.62080.5262592.42058226081.150
9010069981000100601810001004984100099.6421399.91094.6333298.38078.5268091.87052.3227078.720
10010077141000100629810001005197100099.6435899.91094.5336798.31075.8264290.61054.6229979.990
20010012,228100010079201000100595110001004748100095.7333098.71073.5266590.02045.4223375.410
30010014,375100010087351000100609610001004723100096.4328598.90069.6246487.66043.4206472.530
40010015,022100010087621000100613510001004712100096.4303698.81070.5228687.05048.1182071.430
50010015,705100010088231000100587610001004497100097.5284299.12073.5210087.38055.1168673.370
60010016,442100010087121000100575310001004252100097.6267499.10074.8194187.02062.6153075.560
70010015,910100010086701000100560910001003976100099.1251099.64076.5177686.77069.3142078.310
80010016,215100110084191001100549210011003773100199.3227199.69181.9164889.02175.9131981.731
90010016,347100110082681001100525410011003597100199.2217399.63184.8157590.35180.7124584.501
100010016,355100110081331001100506410011003369100199.7199899.85186.8142690.74184.9118987.301
200010014,67910031005955100310030951003100197210021001243100298.6105698.67297.7101797.742
300010012,64310061004207100610020361005100135410051001038100499.9100399.90410010011003
400010010,3751001210029271001110014671001010011521009100101110071001000100610010001006
50001008524100191002140100181001205100151001032100141001003100111001000100910010001008
6000100694210028100172410026100109510023100101410020100100010016100100010014100100010012
7000100546310040100139810035100105010032100100710028100100210023100100010018100100010016
8000100454310054100124010048100102810043100100210045100100010030100100010024100100010021
9000100375110069100113510062100100510055100100110048100100010038100100010032100100010026
10,000100305610086100107810077100100310068100100010060100100010048100100010040100100010033
Aver.99.96884199.8711.4399.85489699.6910.3299.74355699.609.1899.53285099.538.3697.29213898.626.5085.90177492.895.4375.28156286.254.64
Table A5. Computational results for randomly generated instances with ratio 10 % : 20 % : 10 % : 60 % of the number of jobs in the subsets.
Table A5. Computational results for randomly generated instances with ratio 10 % : 20 % : 10 % : 60 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
1098.826395.44097.549194.91094.565791.4809179188.50081.5101481.36074.9110276.59066118670.490
2099.7103499.71098.8160199.19098.8190499.32097.6198498.74089.9211395.13080.4196889.94063184579.730
3099.9213199.95099.8277899.93099305999.67098.3287899.37089.6260895.82074.3227088.46060.3206080.340
401003224100099.9374799.97099.8369899.95098.3339299.50090.9287796.63071.8246988.42052217477.600
50100437010001004704100099.6417499.90099370199.73089.4298896.32066.6256686.83047.1222875.940
60100547310001005368100099.9460899.98098.2398799.55089.3309896.51067.2264387.48042.6227974.590
70100645410001005985100099.9496899.98099.4412599.85087.5321496.02062.6266985.69038.6229172.850
801007498100099.9623599.98099.8519499.94098.8433399.70087.3337296.14061.4271685.64033.6226070.400
9099.9828199.9901006560100099.8524399.96098.9450299.76087.8338896.40061.4275785.78032236171.030
100100916910001007056100099.7550799.95099.2458699.83083.7336094.97058.2268984.27027.7228768.210
20010013,36610001008131100099.9602999.98099481499.79083.1332994.89043.9254177.80010.3217258.520
30010014,9991000100886910001006010100098.9467599.76082312794.18032.2232970.8504.6187048.660
40010015,7041000100884810001006048100099.7449099.93082.5289993.96028.3212066.0801.3171042.280
50010015,7751000100872010001005825100099.7429099.93083.2263893.63021.6188558.3000.7154135.560
60010016,336100010084201001100558210001003938100087.7242094.88118172752.5200140828.980
70010016,298100110084661001100536010011003733100188.4220394.73115157446.0010128222.001
80010016,7071001100803010011005023100199.9347999.97188.7207794.56112.9145740.1510.2120717.321
90010016,135100110079361001100480810011003265100189.5193494.52110.8136834.8010117214.681
100010016,015100110076651001100452810011003049100191.6173795.1619.6131431.2010114412.591
200010013,921100410051011004100254910041001622100398.6113898.7731.210243.523010020.203
300010011,344100910034001009100163610081001210100799.8101499.8060.410030.7050100005
40001008769100171002283100161001245100141001054100131001003100100.110000.190100008
50001006948100281001691100251001102100271001023100211001001100160100001401000011
60001005409100421001362100381001041100341001003100301001000100270100002001000017
70001004121100591001214100531001016100471001001100421001000100340100002701000023
80001003368100801001093100711001006100621001000100661001000100430100003601000030
900010026461001021001048100901001000100801001000100711001000100560100004701000039
10,0001002248100126100102410011910010021001001001000100891001000100700100006301000048
Aver.99.94885799.8216.8299.85492399.7915.3699.67356599.6513.5799.14285499.4312.3691.14212796.239.6431.17175747.908.1417.14155336.506.71
Table A6. Computational results for randomly generated instances with ratio 10 % : 10 % : 10 % : 70 % of the number of jobs in the subsets.
Table A6. Computational results for randomly generated instances with ratio 10 % : 10 % : 10 % : 70 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
1098.637196.23096.261293.63093.685392.15089.5101589.06082.1123184.73067.1130173.10054.3134863.650
2099.8143299.86099.3198899.65097.2227398.77095.5234598.04085.9222293.61067.5211984.29050.7195274.130
3099.7271399.89099.6334199.88099.2339799.76096.9320499.03085.9277794.63065.2234484.98043.5213972.980
4099.9405699.98099.7443999.93099.2401599.80097.1372299.19085298394.90058.4252582.89036.5219170.610
501005231100099.9513099.98099.7459399.93097.7395299.42080.4316593.52055259582.35028.6220067.180
601006574100099.9580499.98099.3493499.84097.1418299.26084.1328394.94049.9265680.80025.6225566.390
7099.9744499.9901006365100099.4511599.88097.8432899.49080.4326193.90048.3270680.60021233065.750
80100850510001006737100099.7541599.94096.5442299.21075.3330492.37042.7266778.29016.8225863.020
901009185100099.8733399.97099.8562399.96097.7455599.50076.4341792.98037.6269676.34013.3228861.320
1001009909100099.9730599.99099.6557199.93098.2454699.60074.4344992.49035.8269575.92011.8231461.620
20010013,80610001008387100099.8614699.97096.2473699.18063.5326188.75016.1252766.6802.7200651.400
30010015,55010001008870100099.9608499.98097.3456399.41053.6306784.8406.9221557.9700.6176543.630
40010015,85610001008573100099.9585299.98096.9430499.28048.3273781.1103.4204952.7600159637.340
50010016,1581000100857610001005760100097.9406799.48144.9247177.7001.8172743.1400.1140228.740
60010016,21610011008425100199.9541699.98198.8372499.68142.9221774.2411.6153936.0010127921.811
70010016,3381001100814210011005055100199.1343299.74140.4205971.0111.4142030.5610119716.461
80010016,5481001100790910011004744100199320699.69136.8182165.2910.4131924.4910113311.741
90010016,0001001100749410011004477100199292999.66130.5171659.5010.1126420.971010988.931
100010015,8061001100729410011004123100199.3269499.74136.1153558.3710.1117715.121010575.391
200010013,179100610044241005100220010051001431100424.3105928.524010070.703010020.203
300010099601001310027461001210014071001110011031001020.2100720.7580100007010010.107
400010074021002410018431002210011301002010010271001716.2100016.2150100001201000010
500010056161004010014501003610010421003110010081002812.2100012.2230100001901000015
60001004234100591001204100531001019100561001000100429.810009.8340100002801000023
70001003236100821001083100641001007100651001000100587.710007.7470100003801000031
800010025111001211001040100981001002100901001000100767.610007.6610100005101000043
90001002059100140100101510012410010011001101001000100966.210006.2790100006501000052
10,0001001728100174100101110015410010011001591001000100120510005970100008001000065
Aver.99.93884499.8523.7199.80494899.7520.4699.51358199.6419.7198.13283999.2016.3247.00210960.8213.3219.98173441.001110.91152931.879.07
Table A7. Computational results for randomly generated instances with ratio 5 % : 20 % : 5 % : 70 % of the number of jobs in the subsets.
Table A7. Computational results for randomly generated instances with ratio 5 % : 20 % : 5 % : 70 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
2098.8135399.04096.7199398.29093.3232597.08086.9234394.15071.4221086.38052.1209176.09034.8194165.330
4099.2410199.80098.7442699.71093.5391898.34088.8368596.93060.9300786.63033.2253973.02019218962.270
6099.1662899.86098.7571299.75093.1489398.51082.5410695.64048.8317683.44021.1258468.89010.9226760.170
8099.5861499.93098.3670199.75094.7533598.93079.5444795.21038.9327280.96014.5273868.0405.7225457.590
10099.910,16599.99098.7720299.81092.6574098.69076.5463494.80029.4333878.16010.5271166.5802.7222455.940
20010013,856100098.6850999.84087.5608397.93059.5469191.20013.4333373.6904.8254662.3300.7204151.200
30010015,201100099.3870599.92082.4618797.07044.5456687.6907.3302569.1602.4228757.2400.1183345.500
40099.915,92499.99098.4896499.82075.3588895.77032.3433884.2309.2272766.5901.9197050.2000.1159237.250
50010016,186100098858899.76071565294.80028.6398781.89112.6246864.5501.5179445.0900137927.480
60010016,531100197.5843799.70165.2539193.51126.8366079.92116.4218461.7210.6155636.0510128722.301
70010016,251100198.7828299.84163.8496792.65124.6344178.00117.1208760.2820.3145231.3410118615.681
80010016,462100198.3793799.79162.1473691.91126.5319276.94119.3180655.2610133825.2610113111.581
90010016,099100197761399.61158.8443990.70128.7288575.22123.5169454.8410.4123719.4810111810.551
100010015,750100196.6715799.52159.2418690.18129.7270874.04123155150.3510.1121117.511010807.411
200010013,055100697.8452199.51565.8216484.20576.3141683.26523.5106328.034010121.193010030.303
300010010,0381001399.5276699.821286.4140390.311194110994.591017.9100718.47801000070100006
400010075681002599.6182399.782296.3111896.692098.5102198.531815.9100115.98140100001201000010
500010056131004010014301003597.4105697.543299.5100899.502911.4100011.4230100001901000016
600010041571005910011871005499.5101499.514899.9100299.904610.4100010.4340100003301000023
700010031081008510010761007499.7100799.70661001000100608.610008.6470100003901000031
8000100258110011010010511009899.7100499.70881001000100786.410006.4650100005001000041
90001002029100140100101410013110010001001151001000100996.710006.7790100006301000052
10,00010016721001751001010100166100100010013810010001001225.910005.9990100008001000066
Aver.99.84969399.9428.6198.71504899.7526.1784.23350095.8122.9668.85270690.5120.5721.65195447.1316.486.23161230.3613.523.22141423.0711
Table A8. Computational results for randomly generated instances with ratio 5 % : 15 % : 5 % : 75 % of the number of jobs in the subsets.
Table A8. Computational results for randomly generated instances with ratio 5 % : 15 % : 5 % : 75 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
2099147599.32096.5224298.44093.4245097.18084.7251493.52068.7234086.07049.2210374.99030.8195263.730
4099.4448799.84098.2451699.53093.5424598.35085.9385096.21052.9311284.38030.1256571.85013.7214959.190
601006924100097.6594299.58091.7493698.30080.3422595.24043.5324282.02017.7262067.8205.9228258.150
8099.8911399.98098.7686399.78092.2551298.55075.3448794.25032.8335079.37012.3268666.5703.1230157.320
10099.910,50399.99098.3770799.78091.4588498.45073460693.90025.9341177.8404.4274464.7201.1229256.680
20010014,104100097.7873199.73083.2609097.14053474489.80010.2330072.4200.9248659.9000200950.020
30099.815,76799.99098.2898799.80077.5616796.25036.8456285.8802.3296966.8600214453.3600173942.500
40099.915,94899.99097.5870499.71069577594.53025.9423982.3302.1267563.1800193948.3800152234.300
50099.916,45699.99196.7856699.61163.5546093.28120.8392679.6711.3239358.7110169540.9400135726.310
60010016,327100196.6838699.59158.8514991.96117361676.9411.5210253.0910151033.7710128021.881
70010016,349100195.8815299.48152.1483190.04114.2328773.8110.6191147.8310137327.1710119015.971
80010016,329100194.8775399.33152458489.49115.7304472.1710.4173242.4410131123.7210112310.951
90010016,036100196.4743099.52250.3420688.11117.4276070.0410.4155035.6810122318.231010847.751
100010015,802100294.9701899.27250.3385987.12218.8258868.6220.2146031.6410115513.421010575.391
200010012,622100694.3422198.63665.2200682.65636.5134052.615010444.214010090.8940100003
300010093781001597.3252498.931484.7130888.301336.7108041.3912010080.791101000080100007
400010068131002899.7167799.822596.4109296.702340.8101241.5020010000160100001401000011
500010048741004899.8131899.854298.6103198.643838.1100438.3533010000270100002101000018
600010037371006910011781006299.4101099.415537.8100337.9949010000390100003101000026
700010027821009610010621008499.8100499.807735.9100035.968010000540100004101000036
800010022801001281001024100112100100010010133.2100033.290010000710100005701000052
90001001802100163100100610014299.9100199.912733.9100033.9112010000940100007301000061
10,000100153610019910010101001761001000100160311000311410100001140100009101000075
Aver.99.9962899.963397.78504499.5829.1781.00346194.9626.3940.99269166.0123.3510.56193938.5518.964.98159028.95152.37140622.1812.78
Table A9. Computational results for randomly generated instances with ratio 5 % : 5 % : 5 % : 85 % of the number of jobs in the subsets.
Table A9. Computational results for randomly generated instances with ratio 5 % : 5 % : 5 % : 85 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
2098.8189699.31094.9258597.95090.4276396.20082.2277793.30061.1253283.37041.9222072.03024.2206361.660
4099.5513899.90097.4506699.49092.9446898.30079.9392894.63046.8304081.71021.4258868.6608.7223158.090
6099.8796399.97098640999.69091.5520998.33072.5433293.44035.5330180.01011.5268166.2103.2234958.490
8099.7988999.97097.8727299.68088.7556497.90062.9454591.44025341577.3906.7271865.1600.6227355.700
10099.411,14299.95097.5787899.67086.9574997.67064.1473292.16017.4332674.7403.5270564.0700.2225755.380
20099.814,58299.99096.9877899.64073.5608195.54038.2477386.8002.9319069.2800.2240758.3300191547.780
30010015,862100095.8885699.51063.9590793.84021.9444082.2700.6280864.3200207051.6400165139.370
40010016,410100093.8868599.25051.1571591.32011.6412978.4000.1249659.9000177743.7300143830.390
50010016,610100192.8843399.13143.5528689.1418.4366574.9210217053.9810151433.9510128021.881
60010016,129100190.5810998.79137.3497587.3214.5335571.5110196649.1410140728.8610117815.111
70010016,180100190.6779598.78134.7455785.6312.6308968.3710175743.0910128922.3410111710.471
80010015,972100190.8740498.76129.3428183.4412.1272363.9010160737.7710120617.081010908.261
90010015,602100288.3686798.28224.7391880.7820.7254761.0120144830.9410114512.661010444.211
100010015,428100286.6657097.95219.3358677.4420.3227956.1620131423.902010938.512010313.011
200010011,672100883361095.2986.2171645.3470122018.037010292.826010100.9950100004
300010083201001979.4217190.51181.9118617.2817010444.2115010010.1012010000100100008
400010056481003780.9147287.02341.510506.1930010070.7027010000220100001801000015
500010041101006182.8118685.50560.310192.1649010020.2044010000350100002901000023
600010029761009183.8109185.15800.310030.607201000091010000510100004401000034
7000100228010012589103689.381100.110010.209901000088010000700100005701000053
8000100180310016488.8101088.91145010000136010000115010000930100007501000062
9000100151710020989.5100389.531840.110000.11640100001470100001170100009601000079
10,000100129010027090.5100390.5323001000020401000018201000014701000011901000098
Aver.99.87949699.9643.1390.41496995.1537.9636.44339358.4734.1719.65263449.1931.488.23188736.1924.353.70155826.70201.60138820.4316.61
Table A10. Computational results for randomly generated instances with ratio 3 % : 2 % : 5 % : 90 % of the number of jobs in the subsets.
Table A10. Computational results for randomly generated instances with ratio 3 % : 2 % : 5 % : 90 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
10099.7115110.9997095.480160.99407959140.963050.546970.89200.632170.6880127120.6310022820.5590
20099.715,0600.999809487980.993060.462440.935024.346060.83300.128120.6430023520.5740018970.4720
30099.615,6560.999708887970.98604158660.898010.243940.7920023950.5821020020.4990015720.3640
40099.516,2080.9997081.685030.97812656410.86813.539600.7561021420.5331017210.4180013760.2731
50099.716,3970.9998177.283600.973120.351400.84411.436210.7261018680.4641015550.3571012420.1951
60099.616,0830.9998171.979830.964114.248010.82110.732940.6971016880.4081013740.2721011710.1461
70099.615,8730.9997169.176230.959111.244380.80010.429560.6631015250.3441012210.1811010930.0851
80099.615,9350.9997264.671400.95028.141160.77620.526710.6272014030.2872011780.1511010660.0621
90099.615,4900.9997263.667600.94524.937350.74520.123680.5782013240.2452011100.0992010450.0431
100099.715,1210.9998260.564130.93734.434050.7192021540.5362010190.0197010910.0832010270.0262
200099.911,2770.99991030.533090.7909016240.3849011950.1638010030.00314010040.00470100005
300099.878730.99972519.219490.58521011540.13321010270.0261701000025010010.0011201000010
400099.752750.99944315.613800.38839010420.04035010070.00731010000400100002201000017
500099.936410.9997709.811410.20964010080.00857010010.00150010000590100003301000027
6000100266711029.110410.12798010010.0018201000073010000620100004801000039
7000100204511417.610240.098127010010.001118010000100010000920100007301000054
8000100163311866.210050.0671660100001470100001630100001050100008501000071
9000100140712364.510030.04821001000018601000016501000013301000010901000089
10,000100124812942.510000.025274010000231010000206010000168010000137010000112
Aver.99.77100211.0058.7445.8448030.6353.6314.1831120.4747.164.8223130.3843.320.0414950.2237.580.0513330.1728.11011980.1222.74
Table A11. Computational results for randomly generated instances with ratio 2 % : 3 % : 5 % : 90 % of the number of jobs in the subsets.
Table A11. Computational results for randomly generated instances with ratio 2 % : 3 % : 5 % : 90 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
1009911,5560.99909479140.99207558590.956046.846310.88109.334180.73100.926790.6240022060.5440
20099.314,8621.000090.186190.988051.361770.919018.546440.82000.430930.6760023240.5670018890.4680
30099.415,9241.000085.888920.984035.359010.88806.743600.78300.227490.6350020340.5080016110.3790
40099.616,3121.000077.886220.974122.456520.86112.639360.7521024120.5851017330.4221013710.2700
50099.716,3051.000172.684150.967113.251350.82911.236530.7281021220.5281015280.3461012320.1881
60099.216,2971.000166.179560.95619.549170.81410.531990.6881018500.4591013880.2801011790.1521
70099.716,2521.000158.274720.94317.644800.7932028950.6551016600.3981012330.1891010800.0741
80099.516,0751.000251.871860.93323.340880.7632026620.6242015060.3361011780.1511010610.0571
90099.415,6171.000248.667710.92322.638280.7442023980.5832013790.2752011410.1243010360.0351
100099.415,2971.000242.563810.90932.234350.7154021140.5272013520.2602010920.0843010250.0242
200098.811,3490.9991020.133320.7609016460.3929011820.1548010240.02370100006010010.0015
300098.177770.998235.419380.51221011660.14219010380.03717010020.002140100001301000010
400097.452970.995453.913910.30947010370.03635010030.00331010000260100002001000017
500099.138020.998702.211310.13564010090.0095601000050010000460100003301000027
60009926860.9961020.810470.05399010030.0038201000084010000580100004701000039
700099.619960.9981420.510110.016126010010.001112010000100010000870100006601000054
800099.916410.9991850.210040.0061780100001500100001320100001060100008601000070
900099.813870.9992360.210050.00721101000018601000018001000013401000010901000089
10,00099.912470.999293010010.001260010000231010000205010000165010000134010000111
Aver.99.2510,0881.0058.6837.9447940.605411.7131230.47474.0223010.38430.5216090.2634.320.0513330.1727.63011940.1222.58
Table A12. Computational results for randomly generated instances with ratio 2 % : 2 % : 1 % : 95 % of the numbers of jobs in the subsets.
Table A12. Computational results for randomly generated instances with ratio 2 % : 2 % : 1 % : 95 % of the numbers of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
10097.412,1240.998086.578710.982058.458670.926027.346890.84002.833770.70500.227270.6300018980.4730
20097.115,3220.998074.286120.969030.862470.88606.147740.7990030860.6740023260.5680015860.3690
30095.515,8720.99705787280.950010.857870.84300.943790.7710026640.6240019690.4920013480.2580
4009516,3540.997140.986080.93015.455180.8271039180.7431023340.5701016850.4071012320.1881
50094.116,1030.996130.981940.91411.150120.8011036060.7201020760.5181014390.3051011350.1191
60091.416,2410.99512177920.89810.746900.7871030960.6771017500.4291013220.2442010670.0631
70092.316,1380.995116.774030.8862042650.7652028370.6472016040.3771012010.1671010580.0551
80089.615,7850.993210.770810.87320.239090.7442025440.6062014640.3172011500.1301010390.0382
90087.315,3680.99226.965050.8562035950.7212023210.5692013470.2582010940.0862010190.0192
100084.515,2000.99035.262030.8473032660.6943020820.5203012540.2032010690.06520100006
200048.310,8670.952110.131410.68211015360.34910011630.1409010150.0157010020.002601000011
300027.972070.90026018470.45924011240.11021010110.01119010010.001160100001301000021
400014.149130.82549013040.23346010190.01941010030.00335010000280100002301000030
50006.933530.72280011250.11171010020.00263010010.00157010000480100003601000047
60003.423590.591118010370.036104010030.0039201000081010000650100005301000044
70001.918950.482175010170.017142010000126010010.001112010000900100007301000065
80000.815090.343210010030.0031840100001640100001460100001380100009501000078
90000.713080.241269010010.001235010000210010000187010000150010000121010000164
10,0000.511550.139330010010.001290010000257010000230010000183010000151010000123
Aver.54.1499510.8067.3218.4347090.5658.895.6530440.4552.421.8122860.3746.740.1515770.2538.680.0113150.1630.58011250.0831.42
Table A13. Computational results for randomly generated instances with ratio 1 % : 2 % : 2 % : 95 % of the number of jobs in the subsets.
Table A13. Computational results for randomly generated instances with ratio 1 % : 2 % : 2 % : 95 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
10097.311,8830.998085.480600.981058.558610.928030.447180.84603.234490.71400.526460.6210021670.5370
20096.115,2550.997072.787270.96802960580.88104.946980.7950030700.6730023010.5650018280.4510
30096.916,0300.998057.987870.950010.658690.84501.743040.7700027300.6320019560.4880015450.3530
40093.715,8290.996144.786230.93413.954620.8231039310.7451023490.5741016290.3861013890.2801
50092.816,2120.995131.281780.91511.451610.80610.235900.7211020380.5081014720.3211012410.1941
6009416,4140.996122.977660.89910.847330.7891032120.6871018060.4461013100.2371011600.1381
70091.415,9440.99411472960.88020.241710.7602027990.6431015620.3601012350.1901010770.0711
80089.715,7720.99329.370790.87020.139180.7452025100.6022014520.3112011350.1191010640.0601
90085.415,2670.99026.266180.8582035840.7202022590.5572013140.2392011060.0962010250.0242
10008415,0050.98933.361350.8413032510.6923020740.5173012630.2082010660.0622010140.0142
200048.710,9380.95313031570.68310015010.3349011450.1279010070.0077010010.00160100006
300027.672870.90028017920.44124011150.10321010210.02122010010.001160100001301000011
400013.749250.82549012900.22544010290.02839010060.00635010000280100002401000021
50007.933780.72780011020.09372010060.0066301000056010000450100003601000034
60004.424490.610117010490.0471040100009101000081010000670100005301000043
70002.118780.479160010190.019141010010.001126010000112010000890100007401000060
80000.915140.345210010040.0041850100001690100001470100001170100009601000078
90000.412980.233286010010.001236010000208010000186010000149010000121010000100
10,0000.211650.143330010000290010000303010000230010000184010000149010000123
Aver.54.0699180.8067.5818.2947200.5658.845.5030340.4554.791.9622770.3746.790.1715810.2537.470.0313080.1630.58011850.1125.53
Table A14. Computational results for randomly generated instances with ratio 1 % : 1 % : 3 % : 95 % of the number of jobs in the subsets.
Table A14. Computational results for randomly generated instances with ratio 1 % : 1 % : 3 % : 95 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
10098.611,3770.999088.578970.985061.858940.93303046950.84603.533760.70800.426760.6200022220.5450
2009815,0120.999076.887290.972030.860030.88108.447010.8000031700.6820023270.5690018560.4600
30097.415,8470.998067.887270.963014.358420.85101.542580.7670027460.6340019590.4890015580.3580
40096.716,2610.998149.684860.93916.155000.82810.439220.7451023390.5721017240.4201013300.2480
50095.916,2660.997138.182400.92312.751420.8101035400.7161020600.5151014570.3141012300.1871
6009515,9110.997129.378290.90911.347590.7911031510.6821018300.4541013250.2451011260.1121
70092.615,9220.995123.173150.89420.343170.7692028060.6432016430.3911012080.1721010930.0851
80094.415,7270.996215.870050.87920.439060.7452025080.6012014690.3192011540.1331010510.0491
90092.415,3110.995212.766280.8682036060.7212022620.5572013800.2752011020.0932010300.0292
100090.214,9790.99338.162150.8523032310.6903020110.5033012550.2033010600.0572010200.0202
200067.910,7780.970110.530680.67610015670.3629011470.1289010130.013701000060100005
300051.373980.93427017750.43724011290.11422010240.02319010000170100001301000018
400034.249030.86649013060.23448010170.01739010030.00335010000280100002301000019
500022.933120.76779011130.10271010090.0096301000056010000450100003601000030
600017.624400.662116010300.029103010010.0019201000082010000650100005301000044
700015.418380.540160010160.016141010000125010000111010000890100007301000060
80001014700.388208010070.0071850100001650100001460100001170100009501000079
90008.413390.316266010000236010000208010000185010000148010000121010000100
10,0007.811710.213331010000292010000322010000229010000200010000162010000124
Aver.62.4698560.8266.2121.5947050.5659.056.1930490.4555.632.1222650.3746.530.1815940.2538.260.0213150.1631.11011850.1125.63
Table A15. Computational results for randomly generated instances with ratio 1 % : 1 % : 1 % : 97 % of the number of jobs in the subsets.
Table A15. Computational results for randomly generated instances with ratio 1 % : 1 % : 1 % : 97 % of the number of jobs in the subsets.
δ % 5%10%15%20%30%40%50%
n O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t O p t N C S C t
10096.912,0630.997084.979010.980052.559210.916026.147130.83702.433660.7050026180.6120021800.5340
20096.115,3000.997068.186630.962022.160510.86804.146650.7910031160.6770023550.5740018370.4560
30093.615,9130.996050.487900.94206.958830.84000.842670.7650027410.6340019390.4830015540.3560
40093.416,7110.996132.385770.91913.155480.8221039350.7451023570.5761016670.4001013280.2471
50090.616,1570.994120.681800.90211.250540.8031034300.7081020480.5121014470.3091011980.1651
60089.516,0690.993113.278710.88910.147470.7881030890.6761018080.4471013260.2461011280.1131
70087.116,0380.99228.473200.8742042600.7632027100.6301016410.3891012000.1671010930.0851
80084.615,7700.99024.570350.8632037600.7332025100.6022014580.3142011480.1291010400.0381
90078.415,3600.98623.865250.8513034670.7112022330.5522013540.2612010910.0832010240.0232
100073.715,0500.98232.160450.8373031860.6863020410.5103012360.1912010730.0682010160.0162
200030.110,8110.93512030340.67011015380.35010011470.1289010120.0128010030.00370100006
300011.870950.87528017900.44125011100.09922010200.02020010000160100001301000011
40003.547110.79551012740.21561010300.02941010010.00136010000290100002401000021
50000.732740.69684011040.09474010030.0036701000059010000470100004101000032
60000.323350.573121010390.038107010010.0019501000085010000680100005501000046
70000.117470.428166010090.009147010000130010000117010000930100007601000063
8000014870.328218010050.00519401000017001000015101000012301000010001000081
9000012810.219280010000266010000217010000194010000155010000126010000105
10,000011640.141344010000305010000268010000239010000298010000156010000130
Aver.48.9799120.7869.2615.1746930.5563.324.5230290.4454.321.6322510.3748.470.1315860.2544.58013090.1631.95011790.1126.53

References

  1. Pinedo, M. Scheduling: Theory, Algorithms, and Systems; Prentice-Hall: Englewood Cliffs, NJ, USA, 2002. [Google Scholar]
  2. Elmaghraby, S.; Thoney, K.A. Two-machine flowshop problem with arbitrary processing time distributions. IIE Trans. 2000, 31, 467–477. [Google Scholar] [CrossRef]
  3. Kamburowski, J. Stochastically minimizing the makespan in two-machine flow shops without blocking. Eur. J. Oper. Res. 1999, 112, 304–309. [Google Scholar] [CrossRef]
  4. Ku, P.S.; Niu, S.C. On Johnson’s two-machine flow-shop with random processing times. Oper. Res. 1986, 34, 130–136. [Google Scholar] [CrossRef]
  5. Allahverdi, A. Stochastically minimizing total flowtime in flowshops with no waiting space. Eur. J. Oper. Res. 1999, 113, 101–112. [Google Scholar] [CrossRef]
  6. Allahverdi, A.; Mittenthal, J. Two-machine ordered flowshop scheduling under random breakdowns. Math. Comput. Model. 1999, 20, 9–17. [Google Scholar] [CrossRef]
  7. Portougal, V.; Trietsch, D. Johnson’s problem with stochastic processing times and optimal service level. Eur. J. Oper. Res. 2006, 169, 751–760. [Google Scholar] [CrossRef]
  8. Daniels, R.L.; Kouvelis, P. Robust scheduling to hedge against processing time uncertainty in single stage production. Manag. Sci. 1995, 41, 363–376. [Google Scholar] [CrossRef]
  9. Sabuncuoglu, I.; Goren, S. Hedging production schedules against uncertainty in manufacturing environment with a review of robustness and stability research. Int. J. Comput. Integr. Manuf. 2009, 22, 138–157. [Google Scholar] [CrossRef] [Green Version]
  10. Sotskov, Y.N.; Werner, F. Sequencing and Scheduling with Inaccurate Data; Nova Science Publishers: Hauppauge, NY, USA, 2014. [Google Scholar]
  11. Pereira, J. The robust (minmax regret) single machine scheduling with interval processing times and total weighted completion time objective. Comput. Oper. Res. 2016, 66, 141–152. [Google Scholar] [CrossRef]
  12. Kasperski, A.; Zielinski, P. A 2-approximation algorithm for interval data minmax regret sequencing problems with total flow time criterion. Oper. Res. Lett. 2008, 36, 343–344. [Google Scholar] [CrossRef]
  13. Wu, Z.; Yu, S.; Li, T. A meta-model-based multi-objective evolutionary approach to robust job shop scheduling. Mathematics 2019, 7, 529. [Google Scholar] [CrossRef] [Green Version]
  14. Grabot, B.; Geneste, L. Dispatching rules in scheduling: A fuzzy approach. Int. J. Prod. Res. 1994, 32, 903–915. [Google Scholar] [CrossRef]
  15. Özelkan, E.C.; Duckstein, L. Optimal fuzzy counterparts of scheduling rules. Eur. J. Oper. Res. 1999, 113, 593–609. [Google Scholar] [CrossRef]
  16. Kasperski, A.; Zielinski, P. Possibilistic minmax regret sequencing problems with fuzzy parameteres. IEEE Trans. Fuzzy Syst. 2011, 19, 1072–1082. [Google Scholar] [CrossRef]
  17. Lai, T.-C.; Sotskov, Y.N.; Sotskova, N.Y.; Werner, F. Optimal makespan scheduling with given bounds of processing times. Math. Comput. Model. 1997, 26, 67–86. [Google Scholar]
  18. Lai, T.-C.; Sotskov, Y.N. Sequencing with uncertain numerical data for makespan minimization. J. Oper. Res. Soc. 1999, 50, 230–243. [Google Scholar] [CrossRef]
  19. Lai, T.-C.; Sotskov, Y.N.; Sotskova, N.Y.; Werner, F. Mean flow time minimization with given bounds of processing times. Eur. J. Oper. Res. 2004, 159, 558–573. [Google Scholar] [CrossRef]
  20. Sotskov, Y.N.; Egorova, N.M.; Lai, T.-C. Minimizing total weighted flow time of a set of jobs with interval processing times. Math. Comput. Model. 2009, 50, 556–573. [Google Scholar] [CrossRef]
  21. Cheng, T.C.E.; Shakhlevich, N.V. Proportionate flow shop with controllable processing times. J. Sched. 1999, 27, 253–265. [Google Scholar] [CrossRef]
  22. Cheng, T.C.E.; Kovalyov, M.Y.; Shakhlevich, N.V. Scheduling with controllable release dates and processing times: Makespan minimization. Eur. J. Oper. Res. 2006, 175, 751–768. [Google Scholar] [CrossRef]
  23. Jansen, K.; Mastrolilli, M.; Solis-Oba, R. Approximation schemes for job shop scheduling problems with controllable processing times. Eur. J. Oper. Res. 2005, 167, 297–319. [Google Scholar] [CrossRef] [Green Version]
  24. Graham, R.L.; Lawler, E.L.; Lenstra, J.K.; Rinnoy Kan, A.H.G. Optimization and approximation in deterministic sequencing and scheduling. Ann. Discret. Appl. Math. 1979, 5, 287–326. [Google Scholar]
  25. Brucker, P. Scheduling Algorithms; Springer: Berlin, Germany, 1995. [Google Scholar]
  26. Tanaev, V.S.; Sotskov, Y.N.; Strusevich, V.A. Scheduling Theory: Multi-Stage Systems; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1994. [Google Scholar]
  27. Jackson, J.R. An extension of Johnson’s results on job lot scheduling. Nav. Res. Logist. Q. 1956, 3, 201–203. [Google Scholar] [CrossRef]
  28. Johnson, S.M. Optimal two and three stage production schedules with set up times included. Nav. Res. Logist. Q. 1954, 1, 61–68. [Google Scholar] [CrossRef]
  29. Allahverdi, A.; Sotskov, Y.N. Two-machine flowshop minimum-lenght scheduling problem with random and bounded processing times. Int. Trans. Oper. Res. 2003, 10, 65–76. [Google Scholar] [CrossRef]
  30. Matsveichuk, N.M.; Sotskov, Y.N. A stability approach to two-stage scheduling problems with uncertain processing times In Sequencing and Scheduling with Inaccurate Data; Sotskov, Y.N., Werner, F., Eds.; Nova Science Publishers: Hauppauge, NY, USA, 2014; pp. 377–407. [Google Scholar]
  31. Shafranski, Y.M. On the existence of globally optimal schedules for the Bellman-Johnson problem for two machines under uncertainty. Informatika 2009, 3, 100–110. (In Russian) [Google Scholar]
  32. Lai, T.-C.; Sotskov, Y.N.; Egorova, N.G.; Werner, F. The optimality box in uncertain data for minimising the sum of the weighted job completion times. Int. J. Prod. Res. 2018, 56, 6336–6362. [Google Scholar] [CrossRef]
  33. Sotskov, Y.N.; Egorova, N.G. Single machine scheduling problem with interval processing times and total completion time objective. Algorithms 2018, 11, 66. [Google Scholar] [CrossRef] [Green Version]
  34. Sotskov, Y.N.; Lai, T.-C. Minimizing total weighted flow time under uncertainty using dominance and a stability box. Comput. Oper. Res. 2012, 39, 1271–1289. [Google Scholar] [CrossRef]
  35. Sotskov, Y.N.; Egorova, N.G. Minimizing total flow time under uncertainty using optimality and stability boxes. In Sequencing and Scheduling with Inaccurate Data; Sotskov, Y.N., Werner, F., Eds.; Nova Science Publishers: Hauppauge, NY, USA, 2014; pp. 345–376. [Google Scholar]
  36. Matsveichuk, N.M.; Sotskov, Y.N.; Egorova, N.G.; Lai, T.-C. Schedule execution for two-machine flow-shop with interval processing times. Math. Comput. Model. 2009, 49, 991–1011. [Google Scholar] [CrossRef]
  37. Sotskov, Y.N.; Allahverdi, A.; Lai, T.-C. Flowshop scheduling problem to minimize total completion time with random and bounded processing times. J. Oper. Res. Soc. 2004, 55, 277–286. [Google Scholar] [CrossRef]
  38. Allahverdi, A.; Aldowaisan, T.; Sotskov, Y.N. Two-machine flowshop scheduling problem to minimize makespan or total completion time with random and bounded setup times. Int. J. Math. Math. Sci. 2003, 39, 2475–2486. [Google Scholar] [CrossRef] [Green Version]
  39. Kouvelis, P.; Yu, G. Robust Discrete Optimization and Its Application; Kluwer Academic Publishers: Boston, MA, USA, 1997. [Google Scholar]
  40. Kouvelis, P.; Daniels, R.L.; Vairaktarakis, G. Robust scheduling of a two-machine flow shop with uncertain processing times. IEEE Trans. 2011, 32, 421–432. [Google Scholar] [CrossRef]
  41. Ng, C.T.; Matsveichuk, N.M.; Sotskov, Y.N.; Cheng, T.C.E. Two-machine flow-shop minimum-length scheduling with interval processing times. Asia-Pac. J. Oper. Res. 2009, 26, 1–20. [Google Scholar] [CrossRef]
  42. Matsveichuk, N.M.; Sotskov, Y.N.; Werner, F. The dominance digraph as a solution to the two-machine flow-shop problem with interval processing times. Optimization 2011, 60, 1493–1517. [Google Scholar] [CrossRef]
Figure 1. Average percentages of the instances presented in Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8 and Table A9, which were optimally solved at the off-line phase of scheduling.
Figure 1. Average percentages of the instances presented in Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8 and Table A9, which were optimally solved at the off-line phase of scheduling.
Algorithms 13 00004 g001
Figure 2. Average percentages of the instances presented in Table A10, Table A11, Table A12, Table A13, Table A14 and Table A15, which were optimally solved at the off-line phase of scheduling.
Figure 2. Average percentages of the instances presented in Table A10, Table A11, Table A12, Table A13, Table A14 and Table A15, which were optimally solved at the off-line phase of scheduling.
Algorithms 13 00004 g002
Figure 3. Average percentages of the correct decisions made for constructing permutations of the conflict jobs for the instances presented in Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8, Table A9.
Figure 3. Average percentages of the correct decisions made for constructing permutations of the conflict jobs for the instances presented in Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8, Table A9.
Algorithms 13 00004 g003
Figure 4. Average percentages of the correct decisions made for constructing permutations of the conflict jobs for the hard instances presented in Table A10, Table A11, Table A12, Table A13, Table A14 and Table A15.
Figure 4. Average percentages of the correct decisions made for constructing permutations of the conflict jobs for the hard instances presented in Table A10, Table A11, Table A12, Table A13, Table A14 and Table A15.
Algorithms 13 00004 g004
Table 1. Input data for Example 1.
Table 1. Input data for Example 1.
J i l i 1 u i 1 l i 2 u i 2
J 1 6767
J 2 8956
J 3 7956
J 4 23--
J 5 --1620
J 6 1334
J 7 1334
J 8 1334
Table 2. Average percentage of the instances which were optimally solved.
Table 2. Average percentage of the instances which were optimally solved.
δ % 5%10%15%20%30%40%50%Average
Instances from Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8 and Table A999.9398.4888.9380.6663.9750.1844.1075.18
Instances from Table A10, Table A11, Table A12, Table A13, Table A14 and Table A1569.7824.897.962.730.200.030.0015.08

Share and Cite

MDPI and ACS Style

Sotskov, Y.N.; Matsveichuk, N.M.; Hatsura, V.D. Two-Machine Job-Shop Scheduling Problem to Minimize the Makespan with Uncertain Job Durations. Algorithms 2020, 13, 4. https://doi.org/10.3390/a13010004

AMA Style

Sotskov YN, Matsveichuk NM, Hatsura VD. Two-Machine Job-Shop Scheduling Problem to Minimize the Makespan with Uncertain Job Durations. Algorithms. 2020; 13(1):4. https://doi.org/10.3390/a13010004

Chicago/Turabian Style

Sotskov, Yuri N., Natalja M. Matsveichuk, and Vadzim D. Hatsura. 2020. "Two-Machine Job-Shop Scheduling Problem to Minimize the Makespan with Uncertain Job Durations" Algorithms 13, no. 1: 4. https://doi.org/10.3390/a13010004

APA Style

Sotskov, Y. N., Matsveichuk, N. M., & Hatsura, V. D. (2020). Two-Machine Job-Shop Scheduling Problem to Minimize the Makespan with Uncertain Job Durations. Algorithms, 13(1), 4. https://doi.org/10.3390/a13010004

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