Abstract

In this paper, we investigate a static stochastic single machine JIT scheduling problem in which the jobs’ processing times are stochastically independent and follow geometric distributions whose mean is provided, due dates are geometrically distributed with a common mean, and both the unit penalty of earliness/tardiness and the fixed penalty of earliness/tardiness are deterministic and different. The objective is to minimize the expected total penalties for quadratic earliness, quadratic tardiness, and early and tardy jobs. We prove that the optimal schedule to minimize this problem is V-shaped with respect to the ratio of mean processing time to unit tardiness penalty under the specific condition. Also, we show a special case and two theorems related to this JIT scheduling problem under specific situations where the optimal solutions exist. Finally, based on the V-shaped characteristic, a dynamic programming algorithm is designed to achieve an optimal V-shaped schedule in pseudopolynomial time.

1. Introduction

In just-in-time (JIT) scheduling, the decision maker intends to find a rational scheduling scheme so that all jobs will not be completed either too early or too late. In general, it is assumed that there are job due dates in the JIT scheduling model, and both earliness and tardiness are penalized. This paper tackles a scheduling model involving earliness and tardiness penalties in both quadratic and fixed form, that is, these penalties contain not only the variable charges hinging upon the quadratic earliness/tardiness of the jobs but also the fixed charges generated once the job is not punctually finished no matter how much the earliness/tardiness is. Hence, in order to avoid these penalties as much as possible, it is of great importance for decision maker to allocate the job due dates.

The single machine environment plays a crucial role in scheduling theory. The performance of complex systems such as parallel machines generally depends on the schedule quality of a single bottleneck machine. Also, some results and insights obtained from the research of the single machine problems can be regarded to be valuable for solving the scheduling problems on more complex systems. The scheduling models with a single processor seem to appear frequently in practice. For the deterministic version of the models, researchers have studied many scheduling problems and obtained corresponding optimization results. Unfortunately, a large quantity of scheduling problems we encounter in practical are filled with many uncertain factors, so extending the model to the field of stochastic scheduling is an attempt to make the scheduling theory more practical. However, the research progress related to scheduling problem in the stochastic model has been slow. Hence, the scheduling problems under stochastic environment are still a significant and challenging research direction.

The deterministic models related to single machine JIT scheduling problem have become an extensive research topic discussed since the later 1970s. Sidney [1] was one of the pioneers to study this class of scheduling problem. Kanet [2] proposed a constructive algorithm to find an optimal solution in polynomial time for the scheduling problem, where earliness and tardiness penalties are unit weights. Lann and Mosheiov [3], firstly, investigated the problem with the number of early and tardy jobs related to several cost structures by some heuristic algorithms. Later, there were many scheduling problems with earliness and tardiness penalties which have been studied in [46] by different algorithms, whilst this type of scheduling problem with controllable processing times was discussed in [7, 8], and the notion of job preemption in this type of scheduling problem was also examined in [9, 10]. The problems of minimizing the quadratic earliness and tardiness have also been considered. In [1113], researchers, under different assumptions, considered the problem in which all quadratic earliness and tardiness penalties are equal, and this problem with symmetric weighted quadratic penalty function was also studied in [14, 15], while the asymmetric weighted version of this problem was investigated in [1618]. The relevant scheduling problems with quadratic cost function have also been studied in [19, 20]. For the problems with more than two types of penalties, Lee et al. [21] developed a dynamic programming algorithm to solve two instances of the problem which includes a combination of three types of penalties, namely, earliness, tardiness, and tardy jobs. Shabtay [22] considered a case with the objective function , where , , and are nondecreasing functions. Koulamas [23] presented an algorithm for the problem. Baker and Scudder [24] presented an excellent survey for the initial research situations on both linear E-T and quadratic E-T scheduling problems, whilst a brilliant review of scheduling problems with multiple types of penalties was also provided in [25]. In [26, 27], the concept of batch in JIT scheduling problems with different objectives was investigated by several authors.

Although there are most of JIT problems in which processing times and due dates are certain, jobs usually have random processing times and due dates in many practical situations. Surprisingly, there is little literature studied on scheduling problems in stochastic environments. Soroush and Fredendall [28] analyzed the single machine scheduling problem with normally distributed processing times and deterministic due dates to minimize the total expected earliness and tardiness penalties. Baker [29] studied the same case and obtained optimal solutions by designing a branch and bound algorithm. However, this problem was also considered in [30] when both processing times and due dates follow exponential distributions, and they derived the optimal V-shaped schedule by giving a dynamic programming algorithm. Liu and Liu [31] dealt with this type of problem in which release times and processing times are random variables. Cai and Zhou [32] investigated a problem with the three types of penalties, namely, earliness, tardiness, and flow time under normally distributed processing times and identically distributed due dates. Soroush [33] firstly explored the problem whose goal is to minimize the total expected weighted number of early and tardy jobs, where processing times are arbitrarily distributed, due dates are certain and jobs are penalized by fixed earliness weights and tardiness weights. The majority of the corresponding problems with E-T penalties have also been studied in [3436] under distinct constraint conditions and assumptions. The stochastic version of problems with a relevant quadratic cost function has also been studied by several authors. Mittenthal and Raghavachari [37] analyzed a problem with quadratic earliness and tardiness penalties and stochastic machine breakdowns, whilst processing times are deterministic and due dates are common with certainty, and their objective is to minimize the expected value of the weighted sum of symmetric quadratic earliness and tardiness penalties of all jobs. Cai and Zhou [38] considered a similar case with uncertain uptimes and downtimes, but where both processing times and due dates are exponentially distributed. The problem of stochastically minimizing maximum lateness was analyzed in [39], where processing times and due dates follow exponential distribution. In [40], Soroush and Alqallaf investigated the stochastic problem of finding an optimal sequence that minimizes the total expected costs with weight quadratic tardiness penalties on a single machine, where idle time is not allowed, whilst the extending problem with initial idle time was examined in [41].

In this paper, we, under stochastic scenario, address a static single machine JIT scheduling problem with the objective of minimizing the expected total penalties for quadratic earliness, quadratic tardiness, and early and tardy jobs. In this scheduling problem, it is assumed that processing times follow geometric distributions with distinct parameters and due dates follow geometric distributions with a common parameter. Geometric distribution is a commonly used discrete distribution, where there is only one parameter which represents the probability of success of the event, and the value range of geometric random variable is all positive integers. It is also one of the simplest forms in processing time distributions. Moreover, the jobs are penalized by variable weights depending on quadratic earliness/tardiness and fixed earliness/tardiness weights. These weights are distinct and job-dependent. We show that the optimal solution of this problem has V-shaped characteristic with respect to the ratio of mean processing time to unit tardiness penalty, that is, the schedule will first arrange jobs in nonincreasing order of the ratio of mean processing time to unit tardiness penalty and then arrange jobs in nondecreasing order of the ratio of mean processing time to unit tardiness penalty. To the best of the author’s knowledge, no other research, under geometrically distributed processing times and due dates, has studied this objective function which simultaneously incorporated the four types of penalty.

This paper is structured as follows. Section 2 defines and formulates a static single machine stochastic JIT scheduling problem. In Section 3, we obtain the V-shaped property of the optimal schedule for this problem and show a special case and two theorems related to this problem under specific situations where the optimal solutions exist. The designed dynamic programming is presented to obtain an optimal V-shaped schedule in Section 4. Section 5 is the conclusions.

2. Problem Definition and Formulation

Consider stochastic JIT scheduling problem that a set of jobs are to be scheduled on a single machine. Assume that all jobs prepare for processing at time zero, and job preemption is not allowed that, as soon as a job starts to process on a single machine, it can not be disrupted until the process is finished. In addition, it is assumed that the machine which is always available can handle at most one job at a time, no machine breakdown occurs, and the machine inserted idle time are not allowed. In this problem, the processing time of each job is assumed to be independent random variable which follows geometric distribution with parameter . Each job has a due date which has geometric distribution with common parameter , and is independent of .s

The objective is to determine a schedule that minimizes the expected total penalties for quadratic earliness, quadratic tardiness, and early and tardy jobs which can be formulated as follows:where determines the order of processing the jobs is a sequence denotes the completion time of job under denotes the due date of job denotes the earliness of job , where denotes the tardiness of job , where denotes the early indictor variable of job , where  = 1 if , otherwise, denotes the tardy indictor variable of job , where if , otherwise, is unit earliness penalty of job is unit tardiness penalty of job is fixed earliness penalty of job is fixed tardiness penalty of job E represents the expectation of a random variable

3. The V-Shaped Characteristic of the Optimal Schedule

In this section, our goal is to obtain the structure characteristic of the optimal sequence that minimizes , and the objective in (1) can be equivalently expressed intowhere is an indicative function, which takes value 1 when event occurs; otherwise, it takes 0, and , , are defined by

In order to tackle this scheduling problem, we first present a lemma below.

Lemma 1. If the processing times are stochastically independent and geometrically distributed random variables with the mean , the due dates follow geometric distributions with a common mean , independent of . Then, given an arbitrary schedule , let denote the schedule that interchanges order of job and job in . Hence,

Proof. Due to the geometrically distributed of and , we have , , , and . Therefore, we have the expected completion time of job under :where is the set of jobs which is to be processed before job in . Since and are geometrically distributed, which are stochastically independent, we haveandwhere , which impliesSubsequently, by simple calculation, we haveand then, by (11) and (13), we can defineHence, according to (14a) and (14b), the objective in (3) and (4) can, respectively, be rewritten byandFollowing expectation and variance of , we haveand by (9), (17), and being independent of , the objective in (5) can be further written aswhere , , are similarly defined byThus, due to (15), (16), and (18) and , (2) can be rewritten byLet denote the set of jobs scheduled before job in (or before job in ). Therefore, by (12) and (14a), under and , we haveand according to (12), (14b), and similar process with (21),Substituting (21) into (15), (6) can be obtained, and substituting (22) into (16), (7) can be obtained.
Similarly, by (19a) and (19b),andSubstituting (23), (24), and (25) into (18), (8) can be obtained.
According to Lemma 1, when both the processing times and the due dates are geometrically distributed, we have the main result on the expected penalties of quadratic earliness, quadratic tardiness, and early and tardy jobs below.

Theorem 1. Assume that the processing times are stochastically independent and geometrically distributed random variables with means , individually, and the due dates follow geometric distributions with a common mean , independent of . Definewhere , , and are defined unless or . If , , and satisfy the following condition,and the jobs are firstly arranged such that in nonincreasing order of , then they are arranged such that, in nondecreasing order of , for , the schedule minimizes . In other words, the optimal solutions given by the schedule are V-shaped with respect to the order of .

Proof. By combining (6), (7), and (8) with (26a) and (26b), we obtainThen, by extracting common factor, Now, given a sequence , then we denote the sequence which interchanges the order of job and job in and the sequence which interchanges the order of job and job in . Definefor all , and then, we obtainandIf and according to the assumption of Theorem 1, the jobs are firstly handled such thatThen, (31) impliesIf and combining (27) and (30), we deduce thatwhere . Hence, we obtain . Due to the assumption of Theorem 1, the jobs are secondly handled such thatThen, (32) impliesThus, the result of Theorem 1 is proved that an optimal sequence to minimize is V-shaped with respect to .
For distinct , , , and , a large quantity of cases are covered by Theorem 1. When the four weights are some constants, we discuss a special case below.
A special case: assume that , , and are some constants, then condition (27) is rewritten byIf , , and , the above condition evidently holds. Therefore, , , , and have the same order. Hence, by Theorem 1, the optimal schedule is V-shaped with respect to , , , or .
In addition, we, under different assumptions, show two theorems related to this problem where the optimal solutions exist.

Theorem 2. Assume that , and , , and are some constants.(i)If , , and , that is, , , , and have same orders, then is minimized by the schedule in nondecreasing order of , , , or (ii)If , ,and , that is, and have the same orders, and have the same orders, and and have opposite orders, then is minimized by the schedule in nondecreasing order of or , or, equivalently, in nonincreasing order of or (iii)If , and , that is, , , and have same orders and and have opposite orders, then is minimized by the schedule in nondecreasing order of , or, equivalently, in nonincreasing order of , , or

Theorem 3. Assume that , and , , and are some constants.(i)If , , and , that is, , , , and have the same orders, then is minimized by the schedule in nonincreasing order of , , , or (ii)If , , and , that is, and have the same orders, and have the same orders, and and have opposite orders, then is minimized by the schedule in nonincreasing order of or , or, equivalently, in nondecreasing order of or (iii)If , , and , that is, , , and have the same orders and and have opposite orders, then is minimized by the schedule in nonincreasing order of , , or , or, equivalently, in nondecreasing order of

Proof. Let , , , and ; according to the assumption that and the conditions of Part (i), we have and . It follows from (29) that . This proves that if and only if , , , or .
According to the similar proof procedure of Part (i), we can prove Part (ii), Part (iii), Part (iv), Part (v), and Part (vi).

4. Dynamic Programming Algorithm

In this section, in order to obtain the optimal solution of this stochastic problem among all V-shaped schedules, we propose dynamic programming algorithm based on Theorem 1 for the objective function (20). In general, it is necessary for the jobs to be numbered in nondecreasing order of , i.e., . Then, think about a set of the jobs for . Hence, job can be scheduled as the first or the last in this schedule for all the jobs in Ai if a schedule has a V-shaped characteristic. Suppose that is denoted as the optimal V-shaped schedule for the objective function (20) and represents a set of the jobs is scheduled before those jobs in under . Let . Then, provided , , and , is defined as the contribution of all the jobs in to the objective function (20). Afterwards, it is obvious that the contribution of the jobs if job is the first and last job scheduled among all jobs in is, respectively,andwhere , and .

On the basis of the optimal principle of dynamic programming, and and are expressed in (39a) and (39b), respectively, hence the jobs are handled by the optimal V-shaped schedule such that

Thus, we can summarize the dynamic programming algorithm designed to achieve the optimal V-shaped schedule for this stochastic JIT problem.

Since this is a dynamic programming algorithm, we can obtain according to its standard procedure, and the detailed process of the backward tracking is omitted. The algorithm is used to calculate all , at most O steps, for each . Thus, Algorithm 1 can derive the optimal schedule in O pseudopolynomial time for the stochastic single machine JIT problem.

Step 1: for , we can calculate by (39a), (39b), and (40) for all possible values , , and .
 Step 2: define
 Step 3: structure the optimal V-shaped schedule that obtains due to the backward tracking process

5. Conclusions

In this paper, we study a static stochastic single machine JIT scheduling problem in which the processing times are stochastically independent and geometrically distributed with distinct parameters, the due dates are geometrically distributed with a common parameter, and both the unit penalty of earliness/tardiness and the fixed penalty of earliness/tardiness are certain and distinct, where the JIT scheduling concept means that all the jobs are scheduled to complete as close to their due dates as possible. The objective is to find a schedule of jobs that minimizes the expected total penalties for quadratic earliness, quadratic tardiness, and early and tardy jobs. We show that an optimal schedule minimizing this problem has a V-shaped characteristic with respect to by mathematical derivation under the specific condition (27). That is, the schedule will first arrange jobs in nonincreasing order of and then arrange jobs in nondecreasing order of . Nevertheless, there is a special case under three constants with specific value range, namely, , , and , and then, the optimal solution obtained is V-shaped with respect to , , , or . Furthermore, we also give two theorems which contain a total of six situations related to this problem under different assumptions, and then, the optimal solution obtained is either monotonous nonincreasing or monotonous nondecreasing with respect to the ratio of mean processing time to different penalties. Finally, based on the V-shaped characteristic, we develop a dynamic programming algorithm which is designed to obtain an optimal V-shaped schedule in pseudopolynomial time.

As a direction for the future research, it may be an interesting and significant challenge to consider this problem with other distributed processing times and due dates. In addition, this problem would be worth extending towards more complex systems such as parallel processors or dedicated processors.

Data Availability

No data were used to support the findings of the study.

Conflicts of Interest

The author declares that there are no conflicts of interest reported in this paper.