Next Article in Journal
Optimizing OCR Performance for Programming Videos: The Role of Image Super-Resolution and Large Language Models
Previous Article in Journal
Geometric Approximation of Point Interactions in Three-Dimensional Domains
Previous Article in Special Issue
The Expected Competitive Ratio on a Kind of Stochastic-Online Flowtime Scheduling with Machine Subject to an Uncertain Breakdown
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Approaches to the Formulation of Scheduling Problems

by
José Manuel García-Sánchez
* and
Plácido Moreno
Department of Industrial Management I, University of Seville, E41092 Seville, Spain
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(7), 1035; https://doi.org/10.3390/math12071035
Submission received: 4 March 2024 / Revised: 21 March 2024 / Accepted: 27 March 2024 / Published: 30 March 2024
(This article belongs to the Special Issue Mathematical Models and Methods of Scheduling Theory)

Abstract

:
This paper presents two novel formulations for scheduling problems, namely order-position hybrid formulation (OPH) and order-disjunctive hybrid formulation (ODH), which extend and combine parts of existing formulation strategies. The first strategy (OPH) is based on sequence position and linear ordering formulations, adding relationships between constraints that allow relaxing some decision variables. The second approach (ODH) is based on linear ordering and disjunctive formulations. In this work, we prove ODH to be the most efficient formulation known so far. The experiments have been carried out with a large set of problems, which consider single machines and identical parallel machines. Computational results show that OPH is better than the rest of the existing formulations for the case of weighted completion objectives, while ODH turns out to be the best approach for most scenarios studied.

1. Introduction

In this paper, two new approaches to formulate machine scheduling problems are proposed. These new formulations extend and mix parts of previous models. We focus on two scenarios, namely a single machine and identical parallel machines, in order to be able to study a wide range of problems. The set of problems to be considered is 1||LMAX, 1||ΣUi, 1||ΣwiCi, 1||ΣTi, 1||ΣwiTi, Pm||ΣwiCi, and Pm||ΣwiTi.
From a general point of view, the main formulations for single-machine scheduling have been the following:
  • DC: disjunctive constraints formulation [1];
  • LO: linear ordering formulation [2,3];
  • TI: time-indexed formulation [4];
  • SP: sequence position formulation [5,6].
An experimental study of the previous formulations on a single machine was carried out [7] by considering objectives such as 1||ΣwjCj, 1||LMAX, 1||ΣwjTj, and 1||ΣUj, as well as the same problems with release dates. Another experimental study of these formulations on the single-machine total tardiness problem was performed in [8], including a hybrid formulation between the LO formulation and the SP formulation, with an adaptation of the traveling salesman problem formulation with Desrochers and Laporte anti-loop constraints [9]. Although these last two approaches are more recent, neither of them achieves any relevant improvements over the other formulations. The DC and SP formulations were also compared empirically in other single-machine scenarios [10].
However, as far as we are concerned, there has not been a relevant contribution to the formulation of scheduling problems for years. Therefore, the main novelty of this paper is to develop new strategies for modeling scheduling problems, aimed at both single-machine and parallel-machine scenarios.
Considering the results from previous studies and the results presented in this work, it can be concluded that the TI formulation is the best formulation for short job sizes, such as the uniform distribution from the interval [1, 10]. However, with the increase in job size, it becomes a very inefficient formulation. For instance, the time to build the model in a single-machine problem with no more than 100 jobs can take several hours if the job size is in the uniform distribution from the interval [1, 100].
When the problem objective does not imply making the job completion time explicit, SP is the best formulation for any job size [8], but the objective of the problem can be modeled using the completion time of each position, which happens in problems like 1||ΣCj, 1||ΣTj, 1||LMAX, or 1||ΣUj. However, weighted objectives in linear programming cannot be modeled without obtaining the completion time for each job. For these cases, the SP formulation loses efficiency, although its behavior depends on the specific weighted objective. Regarding weighted completion times, its behavior becomes very inefficient, whereas its results improve for weighted tardiness purposes.
LO is a good formulation for problems with few jobs [7]. The set of θ(n3) transitivity constraints, where n refers to the number of jobs, significantly worsens its performance when the number of jobs increases. Finally, DC is the formulation with the worst behavior, even for problems with the computational complexity of class P [8].
The next section shows the modeling of these formulations for single-machine problems. In Section 3, we present the new approaches to formulate single-machine problems. In Section 4, we adapt these formulations for a parallel-machine scenario. Finally, in Section 5, an empirical analysis of the formulations is developed.

2. Overview of Previous Formulations

In this section, a summary is presented for each of the four to-date best approaches for weighted objectives, namely the DC, LO, TI, and SP formulations.
Let N = {1..n} be the set of jobs to be sequenced.
j N ; let pj be the processing time, wj be the weight of the job, and dj be the due date of the job. Also, R = Σ j = 1 n p j .

2.1. Disjunctive Constraints Formulation (DC)

The variables of the formulation are defined as follows:
j , k N / j k : α j k = 1 if job j precedes job k; 0 otherwise;
j N : C j = completion time of job j; L M A X = maximum lateness;
j N : T j = tardiness time of job j; j N : U j = 1 if job j is late; 0 otherwise.
Common constraints in all DC single-machine models are (1)–(3), i.e., constraints (1)–(3) are included for all objectives:
j , k N / j k : α j k + α k j = 1 ,
j N : C j p j ,
j , k N / j k : C j + p k C k + R 1 α j k ,
For   1 | | L M A X : j N : L M A X C j d j ,
M i n   L M A X ,
For   1 | | Uj :   j N : C j d j Σ j = 1 n p j d j U j ,
M i n   Σ j = 1 n U j ,
For   1 | | Σ w j C :   M i n   Σ j = 1 n w j C j ,
For   1 | | Σ T j :     j N : T j C j d j ,
M i n   Σ j = 1 n T j ,
For   1 | | Σ w j T j : c o n s t r a i n t s   ( 9 ) ,
M i n   Σ j = 1 n w j T j ,

2.2. Linear Ordering Formulation (LO)

It uses the same variables that were defined for the DC formulation. Common constraints in all LO single-machine models are the following:
Precedence   constraints :   c o n s t r a i n t s   ( 1 ) ,
Calculation of lower bound for the time completion:
j N : C j Σ k N k j p k α k j + p j ,
Transitivity constraints:
j , k , i N / j k i : α j k + α k i + α i j 2 ,
All objectives for scheduling problems are modeled in the same way as that for the DC formulation (Section 2). The only difference is the replacement of common constraints (2) and (3) by constraints (12) and (13).

2.3. Time-Indexed Formulation (TI)

For job processing, let t = 1, 2, …T be the set of time periods, where T = Σ j = 1 n p j .
The variables of the formulation are defined as follows:
j N , t 1 . . T p j + 1 : x j t = 1 if job j starts processing at time t; 0 otherwise;
j N : C j = completion time of job j.
Common constraints in all TI single machine models are:
Process   each   job   once :     j N : Σ t = 1 T p j + 1 x j t = 1 ,
No more than one job is processed at each time:
t = 1 T : Σ j = 1 n Σ s = max ( 1 , t p j + 1 ) t x j s 1 ,
Objective functions for each type of problem are:
For   1 | | L M A X : j N : C j = Σ t = 1 T p j + 1 ( t + p j 1 ) x j t ,
e x p r e s s i o n s   4   a n d   ( 5 ) ,
For   1 | | U j : M i n   Σ j = 1 n Σ t = d j p j + 1 T p j + 1 x j t ,
For   1 | | Σ w j C :   M i n   Σ j = 1 n Σ t = 1 T p j + 1 w j t + p j 1 x j t ,
For   1 | | Σ T j :         M i n   Σ j = 1 n Σ t = 1 T p j + 1 max 0 , t + p j 1 d j x j t ,
For   1 | | Σ W j T j : M i n   Σ j = 1 n Σ t = 1 T p j + 1 w j max 0 , t + p j 1 d j x j t ,

2.4. Sequence Position Formulation (SP)

Let S = {1…n} be the set of positions in the sequence. Since n is the number of jobs, n is also the number of positions in the sequence.
The variables of the SP formulation are defined as follows:
j N , s S : β j s = 1 if job j is processed at position s; 0 otherwise;
s S : C s P = completion time of position s;
s S : U s = 1 if job in position s is late, 0 otherwise;
j N : C j = completion time of job j;
s S : T s P = tardiness time of position s;
j N : T j = tardiness time of job j.
Common constraints in all SP single-machine models are:
Each   job   at   a   position :   j N : Σ s = 1 n β j s = 1 ,
In   each   position   a   job :   s S : Σ j = 1 n β j s = 1 .
Calculation of the lower bound for the completion time of each position:
s S : C s P Σ j = 1 n p j Σ k = 1 s β j k ,
For   1 | | L MAX : s S : L M A X C s P Σ j = 1 n d j β j s ,
a n d   e x p r e s s i o n   ( 5 ) .
For   1 | | Uj :   s S : C s P Σ j = 1 n d j β j s + R U s ,
M i n   Σ s = 1 n U s .
In the previous cases, where no weighted objectives were imposed, there was no need to work with variables considering the job. In other words, the model could include only variables relative to the position. However, in this case, we need to calculate the lower bound for the completion time of each job from the completion time of each position:
For   1 | | Σ w j C jj : j N , s S : C j C s P R 1 β j s ,
a n d   e x p r e s s i o n   ( 8 ) .
For   1 | | Σ T j : s S : T s P C s P Σ j = 1 n d j β j s ,
M i n   Σ s = 1 n T s P .
For   1 | | Σ w j T j : e x p r e s s i o n s   ( 9 ) ,   ( 11 ) ,   and   ( 26 ) .

3. Novel Formulations—Single-Machine Scenario

In this section, the formulations of the novel approaches are presented for the case of single-machine problems. The parallel-machine scenario will be developed in Section 4.

3.1. Order-Position Hybrid Formulation (OPH)

This novel approach takes elements from both the LO formulation and SP formulation, so that the transitive constraints can be replaced. It uses the following variables:
j , k N / j k : α j k = 1 if job j precedes job k; 0 otherwise;
j N , s S : β j s = 1 if job j is processed at position s; 0 otherwise;
j N : C j = completion time of job j;
j N : T j = tardiness time of job j.
Common constraints in all DC single-machine models are expressions (1) and (12) from the LO formulation and (22) from the SP formulation:
j , k N / j k : α j k + α k j = 1 , j N : C j Σ k N k j p k α k j + p j , s S : Σ j = 1 n β j s = 1 .
To establish a connection between formulations and prevent precedents from forming a loop, there is a need to incorporate the following relationship between both formulations into the model:
j N : Σ k N k j α k j + 1 = Σ s = 1 n s β j s
Constraint (30) ensures that the position in the sequence occupied by job j is equal to the number of jobs preceding job j, plus one. Constraint (30) also forces each job to have a position, thus making it unnecessary to include constraint (21).
Proof. 
Because of constraint (22), there must be a job in each position. If a job were to be in more than one position, some job j would not be able to be placed in any position; hence, Σ k = 1 n β j k = 0 . By constraint (30), for that job j: Σ k N k j α k j = 1 , which is impossible. Therefore, each job has a unique integer position in the sequence, which corresponds to the integer Σ k = 1 n k β j k ,   j N . □
In addition, the set of constraint (30) allows us to relax variables α j k .
Proof. 
Let the established order be:
k = 1 Job o1; k = 2 Job o2; …; k = i Job oi; …; k = n Job on
There is a bijective mapping between {o1, o2, …, on} and the set of jobs {1…n}.
Position(o1) = 1:
o 1 : b y     c o n s t r a i n t s   30 Σ o k N o k o 1 α o k o 1 + 1 = 1   Σ o k N o k o 1 α o k o 1 = 0   k / k > 1 : α o k o 1 = 0 b y     c o n s t r a i n t s   1 k / k > 1 : α o 1 o k = 1
Position(o2) = 2:
b y   c o n s t r a i n t s   30 Σ o k N o k o 2 α o k o 2 = 1 ;   b y   e x p r e s s i o n s   31 . a α o 1 o 2 = 1 k / o k o 2 & o k o 1 ( k > 2 ) : α o k o 2 = 0 b y   c o n s t r a i n t s   ( 1 ) k / k > 2 : α o 2 o k = 1
Position(oi) = i:
b y   c o n s t r a i n t s   30 Σ o k N o k o i α o k o i = i 1 ;   b y   e x p r e s s i o n s   31 . a α o 1 o i = 1 b y   e x p r e s s i o n s   31 . b α o 2 o i = 1 b y   e x p r e s s i o n s   31 . i 1 α o i 1 o i = 1   k / k > i : α o k o i = 0 b y   c o n s t r a i n t s   ( 1 ) k / k > i : α o i o k = 1
Position(on) = n:
b y   c o n s t r a i n t s   30 Σ o k N o k o n α o k o n = n 1 ;   b y   e x p r e s s i o n s   31 . a α o 1 o n = 1 b y   e x p r e s s i o n s   31 . b α o 2 o n = 1 b y   e x p r e s s i o n s   31 . i α o i o n = 1 b y   e x p r e s s i o n s   31 . n 1 α o n 1 o n = 1  
Therefore, all variables α j k are assigned a value of 1 or 0, even if they are defined as continuous. □
For all objectives, the constraints and expressions to be included correspond to those proposed for the DC or LO formulations.

3.2. Order-Disjunctive Hybrid Formulation (ODH)

The LO formulation problem shows a high number of O(n3) of transitive constraints (constraint (13)), which prevents precedence loops amongst jobs. We propose to replace transitive constraints with the following disjunctive constraint (32), which only employs n2 constraints while guaranteeing the proper sequence of the jobs:
j , k N / j k : C j + p k α j k C k + R α k j ,
which are equivalent to constraint (3).
Therefore, common constraints of the ODH formulation will be expressions (1) and (12). The other expressions correspond to those of the DC or LO formulation.

4. Formulations for Parallel-Machine Scenario

In this scenario, we need to incorporate set M={1…m} of parallel machines as data. Precedence relationships for LO, OPH, and ODH formulations must only happen amongst jobs processed on the same machine. The machine where each job is processed needs to be considered a decision variable for the LO formulation, but not for the OPH formulation, where it can be treated as an auxiliary calculation of the position in the sequence of each machine, as in the SP formulation. The DC formulation is not going to be studied in this section because of its poor behavior. Finally, we are going to focus on the problems with objectives Pm||ΣwjCj and Pm||ΣwjTj. Other objectives for the study of identical parallel machines were considered in [11,12,13].

4.1. Linear Ordering Formulation (LO)

The variables of the formulation are defined as follows:
j , k N / j k : α j k = 1 if job j precedes job k; 0 otherwise;
j N , i M : ω j i = 1 if job j is processed in machine i; 0 otherwise;
j , k N / j < k : δ j k = 1 if jobs j and k are processed in the same machine; 0 otherwise;
j N : C j = completion time of job j;
j N : T j = Tardiness time of job j.
Common constraints in all LO parallel-machine models are:
Precedence   constraints : j , k N / j k : α j k + α k j = δ j k ,
Activation   of   δ j k : i M , j , k N / j < k : ω j i + ω k i 1 + δ j k ,
A   machine   of   each   job : j N : Σ i = 1 m ω j i = 1 .
In addition, constraints (12) and (13) need to be incorporated. The objective functions for every problem are the same as those of the single-machine scenario.

4.2. Time-Indexed Formulation (TI)

For job processing, let t = 1, 2, …T be the set of time periods, where T = Σ j = 1 n p j .
The variables of the formulation are defined as follows:
i M , j N , t 1 . . T p j + 1 : x j t i = 1 ;
If job j starts processing at time t in machine i; 0 otherwise;
j N : C j = completion time of job j.
Common constraints in all TI parallel-machine models are:
Process   each   job   once :       j N : Σ i = 1 m Σ t = 1 T p j + 1 x j t i = 1 .
No more than one job is processed at each time:
i M , t = 1 T : Σ j = 1 n Σ s = max ( 1 , t p j + 1 ) t x j s i 1 .
Objective functions for each type of problem are:
For   1 | | Σ w j C : M i n   Σ i = 1 m Σ j = 1 n Σ t = 1 T p j + 1 w j t + p j 1 x j t i ,
For   1 | | Σ w j T j : M i n   Σ i = 1 m Σ j = 1 n Σ t = 1 T p j + 1 w j max 0 , t + p j 1 d j x j t i .

4.3. Sequence Position Formulation (SP)

The variables of the SP formulation are defined as follows:
j N , s S , i M : β j s i = 1 if job j is processed at position s in machine i; 0 otherwise;
s S : C s P = completion time of position s;
j N : C j = completion time of job j;
j N : T j = tardiness time of job j.
Common constraints in all SP parallel-machine models are:
Each   job   at   a   position :   j N : Σ s = 1 n Σ i = 1 m β j s i = 1 .
No more than one job in each position:
s S , i M , : Σ j = 1 n β j s i 1 .
Calculation of lower bound for the time completion of each position:
s S , i M : C s i P Σ j = 1 n p j Σ k = 1 s β j k i .
Calculation of lower bound for the time completion of each job:
s S , i M : C s i P Σ j = 1 n p j Σ k = 1 s β j k i .
The objective functions for each type of problem are the same as those defined in the LO formulation.

4.4. Order-Position Hybrid Formulation (OPH)

All variables from the LO formulation and variables β j s i from the SP formulation are used in the OPH formulation.
Expressions (33) and (34) from the LO formulation and (40) and (41) from the SP formulation need to be incorporated. The relationship between the LO and SP formulations is defined as:
s S , i M , : Σ j = 1 n β j s i 1 .
As in expression (30), variables α j k can be considered continuous. Regarding variables δ j k , their value is set to 0 whenever any variable ω j i or ω j k is equal to 0. Therefore, the next constraints must be added:
i M , j , k N / j < k : ω j i ω k i δ j k 1   ,
i M , j , k N / j < k : ω k i ω j i δ j k 1 .
However, variables ω j i are not decision variables in this model but are calculated from variables β j s i by using the following expression:
j N , i M : ω j i = Σ j = 1 n β j s i
The objective functions for each type of problem are the same as those defined in the LO formulation.

4.5. Order-Disjunctive Hybrid Formulation (ODH)

The ODH formulation with parallel machines uses variables from the LO formulation, along with its constraints, (33)–(35) and (12). In addition, constraint (3) must be added to guarantee proper sequences on each machine. Constraint (32), which was defined for the ODH formulation in a single machine, cannot be applied in the case of parallel machines, since two jobs on different machines do not have precedence relationships.
The objective functions for each type of problem are the same as those defined in the LO formulation.

5. Computational Results

Experiments have been performed by following the experimental analysis proposed by [14,15]. Therefore, the following parameters have been considered:
  • Processing time pj: generated from two uniform distributions: [1, 100] and [1, 10];
  • Weight wj: generated from the uniform distribution [1, 10];
  • Due date dj: an integer generated from the uniform distribution [P(L − R/2),P(L + R/2)], where P is the sum of the processing times of all jobs. Parameters L and R correspond to relative measures of the location and range of the distribution, respectively, and take three values each: {0.2, 0.4, 0.6} and {0.6, 1, 1.4};
  • Number of jobs (n): six sizes have been considered {20, 40, 60, 80, 100, and 200}, although no problem includes the six sizes.
From each combination of parameters, four instances are generated. Experiments were run on an Intel(R) Core(TM) i7-10700K CPU @ 3.80 GH with 16 Gb RAM. The optimization library GUROBI Optimizer 10.0.1 was used.
All instances for all problems will be executed during a number of seconds equal to the number of jobs times the number of machines in the corresponding problem, which is a significant difference compared to the analysis of previous scheduling models, where the execution time is usually one hour.

5.1. Linear Programming Relaxation

The linear programming relaxation has been analyzed for 1||ΣwjCj and 1||ΣwjTj. This relaxation of the integer models produces optimal solutions for the formulations TI, LO, OPH, and ODH in the case of problem 1||ΣwjCj. However, the SP formulation obtains solutions with an integrality gap of 97.5%.
Regarding problem 1||ΣwjTj, Table 1 shows the lower bound gap. In this case, the TI formulation presents the best bounds. ODH also shows a better gap than the rest of the formulations. LO and OPH have the same lower bounds.

5.2. Results for 1||ΣwjCj

Regarding the computational results of the integer resolution of the problem, Table 2 shows the average resolution times for the four instances of each size and pi interval. The behaviors of the LO, TI, OPH, and ODH models can be observed. However, results from the sequencing-position model (SP) were not included because, with just 20 jobs and after an hour, no instance had finished, and the best solution found so far was not optimal.
Since the building time of a model is equivalent in all scheduling problems that are to be analyzed, the analysis of the TI formulation for the interval pj ∈ [1, 100] is discarded. For any scheduling problem, preliminary tests for the interval pj ∈ [1, 100] conclude that the TI models do not yield good feasible solutions for most of the cases.
As can be seen in Table 2, OPH and ODH are clearly the fastest formulations, equivalent to the resolution time of the LP relaxation. Although TI is slightly better than LO for pj ∈ [1, 10], the growth of the TI model may make its resolution unfeasible for the pj ∈ [1, 100] distribution interval.
To evaluate the growth of the models, Table 3 shows how to calculate the size of the integer formulations for all four alternatives. In the TI formulation, T corresponds to the number of time periods. Table 4 displays the size for N = {20, 40, 60, 100} and shows that TI contains the smallest number of constraints for pj∈ [1, 10], whereas the OPH formulation contains the smallest number of constraints for pj∈ [1, 100]. LO and ODH contain the smallest number of variables, but ODH considerably reduces the number of constraints.

5.3. Results for 1||LMAX

It was proved that the best formulation to solve the 1||LMAX problem was the SP formulation [9], where the maximum execution time was one hour. In our study, the maximum time is set to n seconds. Table 5 shows the number of instances where the branch and bound execution finished (A), the number of instances where the optimal solution was reached (B), and the number of instances where each model obtained optimal or best feasible solutions compared to other formulations (C). Please note that we have considered the values of n = {40, 60, 80, 100}.
Each row of Table 5 comprises 36 instances, i.e., 4 instances for each combination of [R, L]. Therefore, 288 problems have been solved. The ODH formulation achieves the best solutions, and pj ∈ [1, 100] reaches optimal values for all instances. It can also be seen that the SP formulation is the one that most often completes the resolution, thus guaranteeing the optimal solution. Although the TI formulation shows very good behavior, it significantly decreases convergence with 100 jobs. LO never provides optimal solutions when it does not finish. Finally, OPH presents the worst performance here.

5.4. Results for 1||ΣUj

Although ODH yields very good results for problem 1||ΣUj, the SP model provides the best overall performance (Table 6). However, if the results are analyzed based on parameter R (Table 7), ODH turns out to be the best model with R equal to 1 and 1.4, whereas its performance decreases with R equal to 0.6. At that value, SP works very well. In other words, SP is the best formulation for problems where delivery dates are very similar, but if we increase the range of delivery date values, ODH becomes the best formulation.
Table 6 also shows that TI loses convergence as the number of jobs increases, while the ODH formulation presents an overall better behavior for a greater interval of pj.

5.5. Results for 1||ΣTj

Henceforth, formulations are studied for the NP problems. Results are analyzed for n = {40, 60, 80} and are shown in Table 8. Since 4 instances for each combination [L, R] are considered, there is a total of 216 instances. The TI formulation presents the best overall results, provided that the formulation is feasible. If all ranges are considered, ODH presents the best results by a considerable margin. Although the LO formulation comes in second place, it loses efficiency when the number of jobs increases, due to its constraints to avoid loops in the precedents of order O(n3).
As in the previous models, SP continues to have a direct formulation of the objective function with the use of positions. It is not necessary to calculate the delay of the job, but the delay of the position. However, SP formulation does not show good results, and there are no important differences regarding the length of the process times.

5.6. Results for 1||ΣwjTj

The problem of minimizing the total weighted tardiness in single-machine scheduling is a well-known, strongly NP-hard problem [16]. Therefore, the resolution of instances until reaching an optimal solution becomes harder with size. We have decided to maintain the same instances as in the previous problem 1||ΣTj.
In this problem, SP needs to calculate the delay of each job from the position delay, so this formulation will lose efficiency. Results are displayed in Table 9 and follow a similar trend as problem 1||ΣTj.
The TI formulation yields the best results. Regarding the rest of the formulations, ODH is the one presenting the best convergence. LO and OPH improve the results of the SP formulation, with LO being slightly better. There are no notable differences with respect to the length of the process times.

5.7. Results for Pm||ΣwjCj

To minimize completion times with parallel machines, we have considered the number of jobs to be n = {20, 40, 60} and the number of machines m = {2, 3, 5}. Four instances for each combination of the processing time interval, {pj ∈ [1, 100], pj ∈ [1, 10]}, n, and m, have been generated, leading to a total of 72 instances.
The TI formulation reaches the optimum in all 36 problems with pj ∈ [1, 10]. If any range of pj is considered, ODH outperforms the rest of the formulations, displaying the best performance in 45 problems. The OPH and LO formulations show similar results, with OPH being slightly better. Like the single-machine problem, SP exhibits very poor convergence. Since there were no differences regarding the length of the process times, results of both intervals have been integrated in Table 10. Therefore, each row contains eight problems, except for the TI formulation, which has only executed four problems.

5.8. Results for Pm||ΣwjTj

In the case of problem Pm||ΣwjTj, 648 problems are tested. As in the Pm||ΣwjCj problem, there were no differences regarding the interval of the process times, so the results of both intervals have been integrated in Table 11. Therefore, each row contains 72 problems, except for the TI formulation, which has only executed 36 problems in each row.
Like the previous problem, TI achieves optimal solutions for all problems. Although the results of the SP and ODH formulations are roughly similar, ODH provides a greater number of optimal solutions. Finally, OPH improves the LO formulation.

5.9. Summary and Applications

Scheduling models are used in almost all production systems and environments. Their formulations are needed in many applications as a tool for the experimental study of heuristic allocation techniques and even as a way to achieve exact optimal solutions in environments with a restrained number of jobs.
As a summary of results, the total number of instances where each formulation obtained the optimal or the best feasible solution in each problem is shown in Figure 1. In other words, Figure 1 depicts the summation of the C column from Table 5, Table 6, Table 7, Table 8, Table 9, Table 10 and Table 11. It can be seen that ODH outperforms the other formulations in most of the problems. Regarding the Pm||ΣwjTj problem, the TI formulation achieves a great performance only in scenarios with a small number of jobs.

6. Conclusions

In this paper, two novel modeling approaches for scheduling problems have been presented by extending already existing modeling strategies. The strategy for the analysis has been to use short execution times: a number of seconds equal to the number of jobs times the number of machines. Afterward, we collected the number of completions reached and the number of optimal solutions, and instead of the average error, we have computed the number of times that each formulation achieved the best result amongst the five approaches.
The first formulation, OPH, based on the combination of the linear order and sequence position formulations, shows very good results in weighted completion time objectives. The second formulation, ODH, based on the combination of the linear order formulation and the disjunctive formulation, yields the best results in our experimental study. Therefore, ODH should be considered the main formulation for further work and other study scenarios.
If only processing times belonging to the interval [1, 10] are assessed, the time index formulation produces the best results. However, when the processing time is extended to interval [1, 100], the model drastically increases its size. For instance, it takes more than 15 min to generate a problem with 40 jobs and one machine. In addition, its results are very poor compared to the processing time.
Future lines of research involve extending the two novel formulations to flowshop environments. Concerning limitations of the research methodology, mathematical models can be inefficient in complex scenarios with a large number of jobs. However, as a tool for studying heuristic techniques, ODH achieves optimal solutions while showing a better behavior than previous formulations.

Author Contributions

Conceptualization, J.M.G.-S.; methodology, J.M.G.-S.; software, J.M.G.-S.; validation, J.M.G.-S. and P.M.; formal analysis, J.M.G.-S.; investigation, J.M.G.-S.; resources, J.M.G.-S.; writing—original draft preparation, J.M.G.-S.; writing—review and editing, P.M.; visualization, J.M.G.-S. and P.M.; supervision, J.M.G.-S. and P.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Dataset available on request from the authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Balas, E. On the facial structure of scheduling polyhedra. Math. Program. 1985, 24, 179–218. [Google Scholar] [CrossRef]
  2. Dyer, M.E.; Wolsey, L.A. Formulating the single machine sequencing problem with release dates as a mixed integer program. Discret. Appl. Math. 1990, 26, 255–270. [Google Scholar] [CrossRef]
  3. Nemhauser, G.L.; Savelsbergh, M.W.P. A cutting plane algorithm for the single machine scheduling problem with release times. In Combinatorial Optimization: New Frontiers in Theory and Practice; Akgul, M., Hamacher, H.W., Tufekci, S., Eds.; Springer: Berlin/Heidelberg, Germany, 1992; pp. 63–84. [Google Scholar] [CrossRef]
  4. Sousa, J.P.; Wolsey, L.A. A time-indexed formulation of non-preemptive single-machine scheduling problems. Math. Program. 1992, 54, 353–367. [Google Scholar] [CrossRef]
  5. Wagner, H.M. An integer programming model for machine scheduling. Nav. Res. Logist. 1959, 6, 131–140. [Google Scholar] [CrossRef]
  6. Lasserre, J.B.; Queyranne, M. Generic scheduling polyhedral and a new mixed-integer formulation for single-machine scheduling. In Proceedings of the 2nd IPCO (Integer Programming and Combinatorial Optimization) Conference, Pittsburgh, PA, USA, 25–27 May 1992; pp. 136–149. [Google Scholar]
  7. Keha, A.B.; Khowala, K.; Fowler, J.W. Mixed integer programming formulations for single machine scheduling problems. Comput. Ind. Eng. 2009, 56, 357–367. [Google Scholar] [CrossRef]
  8. Baker, K.R.; Keller, B. Solving the single-machine sequencing problem using integer programming. Comput. Ind. Eng. 2010, 59, 730–735. [Google Scholar] [CrossRef]
  9. Desrochers, M.; Laporte, G. Improvements and extensions to the Miller-Tucker-Zem-lin subtour elimination constraints. Oper. Res. Lett. 1991, 10, 27–36. [Google Scholar] [CrossRef]
  10. Chao, M.-T.; Lin, B.M.T. Scheduling of Software Test to Minimize the Total Completion Time. Mathematics 2023, 11, 4705. [Google Scholar] [CrossRef]
  11. Ouazene, Y.; Nguyen, N.-Q.; Yalaoui, F. Workload Balancing on Identical Parallel Machines: Theoretical and Computational Analysis. Appl. Sci. 2021, 11, 3677. [Google Scholar] [CrossRef]
  12. Kramer, A.; Dell’Amico, M.; Iori, M. Enhanced arc-flow formulations to minimize weighted completion time on identical parallel machines. Eur. J. Oper. Res. 2019, 275, 67–79. [Google Scholar] [CrossRef]
  13. Kramer, A.; Iori, M.; Lacomme, P. Mathematical formulations for scheduling jobs on identical parallel machines with family setup times and total weighted completion time minimization. Eur. J. Oper. Res. 2021, 289, 825–840. [Google Scholar] [CrossRef]
  14. Potts, C.N.; Van Wassenhove, L.N. An algorithm for single machine sequencing with deadlines to minimize total weighted completion time. Eur. J. Oper. Res. 1983, 12, 379–389. [Google Scholar] [CrossRef]
  15. Hariri, A.M.A.; Potts, C.N. An algorithm for single machine sequencing with release dates to minimize total weighted completion time. Discret. Appl. Math. 1983, 5, 99–109. [Google Scholar] [CrossRef]
  16. Lenstra, J.K.; Rinnooy Kan, A.H.G.; Brucker, P. Complexity of machine scheduling problems. Ann. Discret. Math. 1977, 1, 343–362. [Google Scholar] [CrossRef]
Figure 1. Total results for each formulation.
Figure 1. Total results for each formulation.
Mathematics 12 01035 g001
Table 1. Lower Bound GAP (%).
Table 1. Lower Bound GAP (%).
LOTISPOPHODH
GAP (%)55.91.565.455.918.5
Table 2. Summary of computational results on problem 1||ΣwjCj.
Table 2. Summary of computational results on problem 1||ΣwjCj.
NLOTIOPHODH
pj ∈ [1, 100]200.023.210.020.02
400.1817.620.040.11
600.7839.620.090.21
1004.31115.640.260.49
20040.6611.040.51
pj ∈ [1, 10]200.020.100.030.02
400.210.410.050.12
600.771.200.100.31
1003.654.230.230.63
20040.3318.271.013.55
1 To generate the TI models, in the case of pj ∈ [1, 100], it took 8 s for 20 jobs, 33 s for 40 jobs, 937 s for 60 jobs, and 11.185 s for 100 jobs. Therefore, the creation of the TI model for 200 jobs was ruled out.
Table 3. Formulation size expressions for 1||ΣwjCj.
Table 3. Formulation size expressions for 1||ΣwjCj.
Jobs LOSPOPHODHTI
nConstraintsn3  5n2/2 + 3n/2n2+ 3nn2/2 + 7n/2n2T + n
Variablesn2n2+ 2n2n2n2T(n − 1) + n
Integer variablesn2nn2n2n2nT(n − 1) + n
Table 4. Formulation size for 1||ΣwjCj.
Table 4. Formulation size for 1||ΣwjCj.
TI
Jobs LOSPOPHODHpj ∈ [1, 10]pj ∈ [1, 100]
20Constraints7030460270400136922
Variables400440800400232018,580
Integer variables380400400380232018,580
40Constraints60,060172094016002601874
Variables1600168032001600880073,360
Integer variables1560160016001560880073,360
60Constraints207,0903780201036003912820
Variables360037207200360019,860165,580
Integer variables354036003600354019,860165,580
80Constraints975,15010,300535010,0006615087
Variables10,00010,20020,00010,00056,133498,733
Integer variables990010,00010,000990056,133498,733
Table 5. Computational results for 1||LMAX problem.
Table 5. Computational results for 1||LMAX problem.
LOSPTIOPHODH
nABCABCABCABCABC
pj ∈ [1, 100]40212121363636---223636233636
60363636293636---233131253636
80242424273434---121919213636
100171717262929---133203636
pj ∈ [1, 10]40363636363636363636203636243636
60353535363636363636202828213636
80262628303634353535101615223634
100202020273030252525166223434
Total 217217217247273273132132132109175175178286286
Greatest number of optimal solutions in each scenario (A, B and C) has been highlighted in bold.
Table 6. Computational results for 1||ΣUj problem.
Table 6. Computational results for 1||ΣUj problem.
LOSPTIOPHODH
nABCABCABCABCABC
pj ∈ [1, 100]40122424313636---172828233636
60161818283535---131818233333
80111111263636---577213131
100888202930---233182525
pj ∈ [1, 10]40152525313636363636203131243636
60121313303636282929111515203434
80999303232121515777192828
100888243233788333202829
Total 9111611622027227483888878112112168251252
Greatest number of optimal solutions in each scenario (A, B and C) has been highlighted in bold.
Table 7. Computational results for 1||ΣUj problem with respect to parameter R.
Table 7. Computational results for 1||ΣUj problem with respect to parameter R.
LOSPTIOPHODH
RABCABCABCABCABC
0.60131389949621242411232345960
1384848548484272929152929689696
1.4535555779494353535526060969696
Total9111611622027227483888878112112168251252
Greatest number of optimal solutions in each scenario (A, B and C), as well as for each value of parameter R, has been highlighted in bold.
Table 8. Computational results for 1||ΣTj problem.
Table 8. Computational results for 1||ΣTj problem.
LOSPTIOPHODH
nABCABCABCABCABC
pj ∈ [1, 100]40102727102525---52626223333
607212271417---42325172628
80617175715---0151792324
pj ∈ [1, 10]40727277363636363652424213131
608202110191932343661920152728
80619196151530363601010112727
Total 4413113345116127981061082011712295167171
Greatest number of optimal solutions in each scenario (A, B and C) has been highlighted in bold.
Table 9. Computational results for 1||ΣwjTj problem.
Table 9. Computational results for 1||ΣwjTj problem.
LOSPTIOPHODH
nABCABCABCABCABC
pj ∈ [1, 100]4019293091313---52627233435
609212371013---51617183029
80413195710---066152425
pj ∈ [1, 10]401526266141436363672626223232
609212110121235363671515182626
8051313566333336099152527
Total 611231324262681041051082498100111171174
Greatest number of optimal solutions in each scenario (A, B and C) has been highlighted in bold.
Table 10. Computational results for Pm||ΣwjCj problem.
Table 10. Computational results for Pm||ΣwjCj problem.
LOSPTIOPHODH
mnABCABCABCABCABC
220888000444888888
40000000444000004
60003000444001002
320055000444055088
40000000444004004
60000000444005003
520024000444001046
40002000444002004
60000000444002006
Total 815220003636368132882045
Greatest number of optimal solutions in each scenario (A, B and C) has been highlighted in bold.
Table 11. Computational results for Pm||ΣwjTj problem.
Table 11. Computational results for Pm||ΣwjTj problem.
LOSPTIOPHODH
mnABCABCABCABCABC
22038454782222363636253840465860
4022916183636362292413
602221111223636364450020
32015212561319363636101922233245
4010321203636362014206
6000109925363636327003
5201423276131836363692025203441
4020700173636362072013
6000063173236361132124
Total 74931304978178320324324588613297129225
Greatest number of optimal solutions in each scenario (A, B and C) has been highlighted in bold.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

García-Sánchez, J.M.; Moreno, P. Novel Approaches to the Formulation of Scheduling Problems. Mathematics 2024, 12, 1035. https://doi.org/10.3390/math12071035

AMA Style

García-Sánchez JM, Moreno P. Novel Approaches to the Formulation of Scheduling Problems. Mathematics. 2024; 12(7):1035. https://doi.org/10.3390/math12071035

Chicago/Turabian Style

García-Sánchez, José Manuel, and Plácido Moreno. 2024. "Novel Approaches to the Formulation of Scheduling Problems" Mathematics 12, no. 7: 1035. https://doi.org/10.3390/math12071035

APA Style

García-Sánchez, J. M., & Moreno, P. (2024). Novel Approaches to the Formulation of Scheduling Problems. Mathematics, 12(7), 1035. https://doi.org/10.3390/math12071035

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