[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
The Strategic Weight Manipulation Model in Uncertain Environment: A Robust Risk Optimization Approach
Next Article in Special Issue
Joint Optimization of Order Allocation and Rack Selection in the “Parts-to-Picker” Picking System Considering Multiple Stations Workload Balance
Previous Article in Journal
Analysis of Urban Electric Vehicle Adoption Based on Operating Costs in Urban Transportation Network
Previous Article in Special Issue
A New Lagrangian Problem Crossover—A Systematic Review and Meta-Analysis of Crossover Standards
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two Due-Date Assignment Scheduling with Location-Dependent Weights and a Deteriorating Maintenance Activity

School of Science, Shenyang Aerospace University, Shenyang 110136, China
*
Author to whom correspondence should be addressed.
Systems 2023, 11(3), 150; https://doi.org/10.3390/systems11030150
Submission received: 5 February 2023 / Revised: 8 March 2023 / Accepted: 9 March 2023 / Published: 15 March 2023

Abstract

:
This paper investigates single-machine scheduling with a deteriorating maintenance activity, where the processing time of a job depends on whether it is handled before or after the maintenance activity. Under common and slack due date assignments, the aim is to find the optimal job schedule, position of the maintenance activity, and optimal value of the common due date (flow-allowance) so that the linear weighted sum of earliness, tardiness and common due date (flow-allowance) value is minimized, where the weights are location-dependent (position-dependent) weights. Through a series of optimal properties, a polynomial time algorithm is proposed and it is then proven that the problem is polynomially solvable.

1. Introduction

Traditional scheduling problems consider that the machine is always available, while actually, if a maintenance (rate-modifying) activity is carried out on the machine, it becomes unavailable, and it reverts to its original condition after the maintenance activity is finished (see Lee and Leon [1]; Ma et al. [2]; Strusevich and Rustogi [3]). In 2006, Mosheiov and Oron [4] considered single-machine scheduling with a rate-modifying activity (denoted r m a ). Under the common due-date (denoted C O N ) assignment, the goal is to minimize the total earliness, tardiness and due-date cost. They proved that the problem can be solved in polynomial time. In 2010, Wang and Wang [5] addressed single-machine scheduling with r m a and slack due-date (denoted S L K ) assignment. They showed that the non-regular objective minimization is polynomially solvable. In 2012, Yin et al. [6] studied single-machine batch delivery scheduling with the r m a and C O N assignment, with the goal of minimizing the sum of earliness, tardiness, due-date, holding and delivery cost. They proved that some special cases of the problem can be solved in polynomial time. In 2014, Bai et al. [7] studied the single-machine problem with deteriorating jobs and r m a . Under the S L K allocation, the objective is to minimize the weighted sum of earliness, tardiness and the common flow allowance cost. They showed that the problem is polynomially solvable. In 2018, Cheng et al. [8] discussed single-machine batch problems with variable maintenance activities. For the minimization of makespan (total completion time), they demonstrated that the problem can be solved in polynomial time in a special case. In 2019, Detti et al. [9] addressed the single-machine problem with a flexible maintenance activity. In 2021, Wang et al. [10] considered the single-machine problem with r m a and a common due-window. For the generalized earliness/tardiness penalties, they showed that the problem is polynomially solvable. In 2022, Zhao et al. [11] studied the single-machine resource allocation problem with an aging effect as well as r m a and S L K assignment. In 2023, Sun et al. [12] delved into single-machine scheduling with worsening effects. With the maintenance activity, they showed that some regular objective minimizations are polynomially solvable. For total weighted completion time minimization, they suggested several heuristic and branch-and-bound algorithms. In addition, in 2014, Ji et al. [13] and Fan and Zhao [14] considered due-date assignment scheduling with a deteriorating maintenance activity (denoted d m a ), i.e., the maintenance duration (time) is a linear deteriorating function of the starting time of this maintenance activity. In 2015, Mor and Mosheiov [15] studied single-machine scheduling with d m a . Under the due-window assignment, they proved that some non-regular minimizations can be solved in polynomial time. In 2017, Li and Chen [16] investigated the single-machine problems with job-rejection and d m a , and they showed that the makespan, total completion time and due-date assignment minimizations are polynomially solvable. Zhu et al. [17] explored single-machine scheduling with a general (including worsening and resource-dependent) maintenance activity. They showed that the problem can be solved in polynomial time for a number of regular and non-regular objectives. Wang et al. [18] addressed identical parallel machine scheduling with d m a . They illustrated that the total completion time minimization can be solved in polynomial time. More recent papers which have studied scheduling problems with r m a and d m a include He et al. [19], Jia et al. [20], liu et al. [21] and Zou et al. [22].
At the same time, in modern manufacturing, more and more industries are adopting due-date assignment systems (see Gordon et al. [23,24], Qian and Zhan [25], Lu et al. [26], Wang [27]) and location-dependent (position-dependent) weights (Brucker [28], Liu et al. [29], Wang et al. [30], Wang et al. [31]). Recently, Jiang et al. [32] investigated proportionate flowshop problems with location-dependent weights. Under the C O N and S L K assignments, they showed that the non-regular objective minimization is polynomially solvable (i.e., time is O ( n 2 log n ) , where n is the number of jobs). Liu et al. [33] studied single-machine scheduling with resource allocation, worsening jobs and location-dependent weights. Under the C O N and S L K assignments, they provided a bi-criteria analysis for the scheduling cost (including the weighted sum of the absolute value in lateness and common due date (flow-allowance) cost) and the resource consumption cost. They proved that three versions of these two costs can be solved in polynomial time. In 2021, Lv and Wang [34] considered the same problems as Jiang et al. [32], and they demonstrated that both problems can be addressed by using a lower-order algorithm (i.e., time is O ( n log n )). Wang et al. [35] and Wang et al. [36] investigated single-machine scheduling with setup times associated with past sequences and weights related to location. Under the C O N and S L K assignments, Wang et al. [35] demonstrated that the weighted sum of earliness, tardiness and due-date cost minimization is polynomially solvable. Wang et al. [36] demonstrated that the weighted sum of lateness, number of early and delayed jobs, and due-date cost minimization is also polynomially solvable.
In real production processes, the location-dependent (position-dependent) weights can be found in services and logistics systems (e.g., in Didi taxi dispatching, Sun et al. [37]). Hence, the work of C O N ( S L K ) assignment and d m a will be continued by considering the location-dependent weights, i.e., we investigate single-machine due-date assignment scheduling with d m a and location-dependent weights. Under the C O N and S L K assignments, the aim of the article is to find the job schedule, location of the maintenance activity, and optimal value of common due-date (flow-allowance) so as to minimize the linear weighted sum of the earliness, tardiness and due date assignment costs, where the weights are related to the location. We demonstrate that the problem can be solved in polynomial time, i.e., time complexity is O ( n 4 ) . For a comparison with other similar papers (see Table 1; the related symbols are given later), this article extends the results of Mosheiova and Oron [4], Wang and Wang [5], Brucker [28], and Liu et al. [29], by scrutinizing a more general scheduling model. The rest of this article is structured as below. In Section 2, we introduce the problem. In Section 3 and Section 4, we explore the main details of common due-date and slack due-date assignments. In Section 5, an example and computational experiments are reported to verify the effectiveness of the algorithms. Section 6 presents the conclusion of the problem.

2. Problem Description

In this article, the problem can be described as: n independent jobs { T ˇ 1 , T ˇ 2 , , T ˇ n } can be usable at time 0 for handling on a single-machine and they are not preemptive. The machine stops working when a maintenance activity is being carried out, if the job T ˇ i is handled before the maintenance activity, the corresponding normal processing time (i.e., basic processing time that is not affected by other factors) is defined as p i , and its actual processing time is b i = p i . If the job T ˇ i is processed after the maintenance activity, the corresponding actual processing time is defined as b i = ε i p i , where ε i is the modified rate of the job T ˇ i , and it satisfies 0 < ε i 1 . Meanwhile, the machine only has a deteriorating maintenance activity, and its maintenance time is: φ ( t ) = t 0 + α S t , where t 0 is the basic maintenance time, α > 0 is the rate of deterioration, and S t is the starting time of d m a . We consider two different types of due-dates that include C O N and S L K . As for the C O N model, all the jobs have the same due date, i.e., d i = d o p t , in which d o p t is a decision variable. For the S L K model, the due date d i of job T ˇ i is equivalent to the sum of the actual processing time b i and the common flow-allowance q o p t , i.e., d i = b i + q o p t , and the common flow-allowance q o p t is a decision variable. For a known schedule σ , the completion (resp. starting) time of the job T ˇ i is C i (resp. S i ), the earliness (resp. tardiness) cost of the job T ˇ i is E i ˜ = max { 0 , d i C i } (resp. L i ˜ = max { 0 , C i d i } ). Let [ i ] be the job being arranged at the ith position in the schedule, then our goal is: (1) to determine the optimal job schedule; (2) to determine where to carry out the maintenance activity of the machine; (3) to determine the value of the common due date and the common flow-allowance; (4) to minimize the objective function: Z ˙ = i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i d o p t ) and Z ¨ = i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i q o p t ) , where μ i , ν i , ω i are weights of the ith location in a schedule, i.e., the weights are the location-dependent weights rather than the weights connected with the job T ˇ i . With three-field notation, the scheduling can be expressed as:
1 d m a , C O N , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i d o p t ) ;
1 d m a , S L K , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i q o p t ) .

3. Results of CON

Assuming that the maintenance activity occurs exactly before the job T ˇ [ i ] , its starting time is equivalent to the completion time of the job T ˇ [ i 1 ] , i.e., S t = C [ i 1 ] . For the sake of simplicity, let j be the location of the maintenance activity.
Lemma 1.
If C [ i ] d o p t , then C [ i + 1 ] d o p t , if C [ i ] d o p t , then C [ i 1 ] d o p t .
Proof. 
For a given schedule σ , if C [ i ] d o p t , we can get C [ i ] d o p t C [ i + 1 ] b [ i + 1 ] d o p t C [ i + 1 ] d o p t . Similarly, if C [ i ] d o p t , we can also get C [ i ] d o p t C [ i 1 ] + b [ i ] d o p t C [ i 1 ] d o p t . □
From Lemma 1, it can be seen that the jobs before the job T ˇ [ i ] are the early jobs, and the jobs after the job T ˇ [ i ] are the late jobs.
Lemma 2.
For a known schedule σ = ( T ˇ [ 1 ] , T ˇ [ 2 ] , , T ˇ [ n ] ) , the optimal value of the common due date d o p t is equal to the completion time of the hth job, i.e., d o p t = C [ h ] , where h satisfies both
i = 1 h μ i i = h + 1 n ν i + i = 1 n ω i 0 a n d i = 1 h 1 μ i i = h n ν i + i = 1 n ω i 0 .
Proof. 
See Appendix A. □
Remark 1.
For a given schedule, if h satisfies both the above inequalities, the optimal common due date can be determined by Lemma 2. But h may not meet the above both inequalities, so we need to set d o p t = 0 .
Lemma 3
(Hardy et al. [38]). For i = 1 n ξ i η i , if the sequence ξ 1 , , ξ n is in non-decreasing order and the sequence η 1 , , η n is in non-increasing order, or vice versa, it is the smallest.
Next, we consider the problem 1 d m a , C O N , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i d o p t ) . Under Lemmas 1–3, the following cases should be considered.
Case 1. If j h , the schedule is shown in Figure 1.
Where d o p t = C h = i = 1 j 1 p i + φ t + i = j h ε i p i and φ t = t 0 + α i = 1 j 1 p i . Then the objective function is expressed as
Z ˙ j = i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i d o p t ) = i = 1 h μ i C h C i + i = h + 1 n υ i C i C h + m = 1 n ω m · C h = t 0 m = 1 j 1 μ m + α i = 1 j 1 p i m = 1 j 1 μ m + i = j h ε i p i m = 1 j 1 μ m + i = 1 j 1 p i m = 1 i 1 μ m + i = j h ε i p i m = j i 1 μ m + i = h + 1 n ε i p i m = i m ν m + i = 1 j 1 p i m = 1 n ω m + t 0 m = 1 n ω m + α i = 1 j 1 p i m = 1 n ω m + i = j h ε i p i m = 1 n ω m = i = 1 j 1 p i α m = 1 j 1 μ m + m = 1 i 1 μ m + m = 1 n ω m + α m = 1 n ω m + i = j h ε i p i m = 1 i 1 μ m + m = 1 n ω m + i = h + 1 n ε i p i m = i n υ m + t 0 m = 1 j 1 μ m + m = 1 n ω m = Z ˙ 1 j + f 1 t 0
where f 1 t 0 = t 0 m = 1 j 1 μ m + m = 1 n ω m is a constant (if j is given),
Z ˙ 1 j = i = 1 j 1 p i α m = 1 j 1 μ m + m = 1 i 1 μ m + m = 1 n ω m + α m = 1 n ω m + i = j h ε i p i m = 1 i 1 μ m + m = 1 n ω m + i = h + 1 n ε i p i m = i n υ m
is only related to the location j of the maintenance activity. Therefore, for a known j, minimizing Z ˙ j is equivalent to minimizing Z ˙ 1 j . Let
λ i l = p i α m = 1 j 1 μ m + m = 1 l 1 μ m + m = 1 n ω m + α m = 1 n ω m 1 l j 1 ε i p i m = 1 l 1 μ m + m = 1 n ω m j l h ε i p i m = l n υ m h + 1 l n
Z ˙ 1 j can be minimized by solving the next assignment problem:
M i n Z ˙ 1 j = i = 1 n l = 1 n λ i l β i l
s . t i = 1 n β i l = 1 l = 1 , , n
l = 1 n β i l = 1 i = 1 , , n
β i l = 0 o r 1 1 i , l n
where β i l is a 0 or 1 variable, if the job T ˇ i is at location l, β i l = 1 , otherwise β i l = 0 .
Case 2. If j > h , the schedule is shown Figure 2.
Where d o p t = C h = i = 1 h p i , and φ t = t 0 + α i = 1 j 1 p i . Then the objective function is expressed as
Z ˙ j = i = 1 n ( μ i E [ i ] ˜ + υ i L [ i ] ˜ + ω i d o p t ) = i = 1 h μ i C h C i + i = h + 1 j 1 υ i C i C h + i = j n υ i C i C h + i = 1 n ω i i = 1 h p i = i = 1 h p i m = 1 i 1 μ m + i = h + 1 j 1 p i m = i j 1 υ m + i = h + 1 j 1 p i m = j n υ m + t 0 m = j n υ m + i = 1 j 1 p i α m = j n υ m + i = j n ε i p i m = i n υ m + i = 1 h p i m = 1 n ω m = i = 1 h p i m = 1 i 1 μ m + α m = j n υ m + m = 1 n ω m + i = h + 1 j 1 p i m = i n υ m + α m = j n υ m + i = j n ε i p i m = i n υ m + t 0 m = j n υ m = Z ˙ 2 j + f 2 t 0
where f 2 t 0 = t 0 m = j n υ m is only related to t 0 and it is a constant,
Z ˙ 2 j = i = 1 h p i m = 1 i 1 μ m + α m = j n υ m + m = 1 n ω m i = h + 1 j 1 p i m = i n υ m + α m = j n υ m + i = j n ε i p i m = i n υ m
is only related to the location j of the maintenance activity. Similarly to Case 1 (i.e., j h ), let
δ i l = p i m = 1 l 1 μ m + α m = j n υ m + m = 1 n ω m 1 l h p i m = l n υ m + α m = j n υ m h + 1 l j 1 ε i p i m = l n υ m j l n
then, the problem can be translated into the assignment problem below:
M i n Z ˙ 2 j = i = 1 n l = 1 n δ i l γ i l
s . t i = 1 n γ i l = 1 l = 1 , , n
l = 1 n γ i l = 1 i = 1 , , n
γ i l = 0 o r 1 1 i , l n
where γ i l is a 0 or 1 variable, when the job T ˇ i is at position l, γ i l = 1 , otherwise γ i l = 0 .
Case 3. If j = n + 1 , the schedule is shown Figure 3.
where d o p t = C h = i = 1 h p i and φ t = t 0 + α i = 1 n p i . This means that there is no maintenance activity on the machine. In this case, the objective function is
Z ˙ n + 1 = i = 1 n ( μ i E [ i ] ˜ + υ i L [ i ] ˜ + ω i d o p t ) = i = 1 h μ i C h C i + i = h + 1 n υ i C i C h + i = 1 n ω i i = 1 h p i = i = 1 h p i m = 1 i 1 μ m + i = h + 1 n p i m = i n υ m + i = 1 h p i m = 1 n ω m = i = 1 h p i m = 1 i 1 μ m + m = 1 n ω m + i = h + 1 n p i m = i n υ m = i = 1 n p i π i
where
π i = m = 1 i 1 μ m + m = 1 n ω m 1 i h m = i n υ m h + 1 i n
The optimal schedule σ as well as the minimum value of Z ˙ can be acquired by matching the minimum weight π i with the maximum processing time p i by using H L P rule (see Lemma 3).
Case 4. if j = 1 , the schedule is shown Figure 4.
where d o p t = C h = i = 1 h ε i p i + φ t , and φ t = t 0 . This means that the maintenance activity occurs before the jobs are processed. In this case, the objective function is
Z ˙ 1 = i = 1 n ( μ i E [ i ] ˜ + υ i L [ i ] ˜ + ω i d o p t ) = i = 1 h μ i C h C i + i = h + 1 n υ i C i C h + i = 1 h ε i p i m = 1 n ω m + t 0 m = 1 n ω m = i = 1 h ε i p i m = 1 i 1 μ m + i = h + 1 n ε i p i m = i n υ m + i = 1 h ε i p i m = 1 n ω m + t 0 m = 1 n ω m = i = 1 h ε i p i m = 1 i 1 μ m + m = 1 n ω m + i = h + 1 n ε i p i m = i n υ m + t 0 m = 1 n ω m = i = 1 n ε i p i π i + f 3 t 0
where f 3 t 0 = t 0 m = 1 n ω m is only related to t 0 and it is a constant. Similarly, the optimal schedule of this case can be obtained through Lemma 3.
From the above cases, the optimal solution of the problem 1 d m a , C O N , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i d o p t ) can be obtained by using the following algorithm.
Theorem 1.
The optimal solution of 1 d m a , C O N , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i d o p t ) can be obtained in O n 4 time by using Algorithm 1.
Proof. 
The correctness of Algorithm 1 follows from the above analysis. For each j, the complexity of the assignment problem is O n 3 , and it takes n 1 times. So the time complexity of Algorithm 1 is O n 4 . □
Algorithm 1: Solution of 1 d m a , C O N , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i d o p t )
Initialization: Let Z ˙ = , σ ˙ = 0 , d o p t = 0 and j = 0 .
Step 1: Calculate h from Lemma 2.
Step 2: For j = 1 n + 1
   If j h , then
    obtain the minimum value Z ˙ j and the schedule σ ˙ by using (3)–(6);
     If Z ˙ j < Z ˙ , then
     let Z ˙ = Z ˙ j , j = j , d o p t = d o p t and σ ˙ = σ ˙ ;
   If j > h , then
    obtain the minimum value Z ˙ j and the schedule σ ˙ by using (7)–(10);
     If Z ˙ j < Z ˙ , then
     let Z ˙ = Z ˙ j , j = j , d o p t = d o p t and σ ˙ = σ ˙ ;
   If j = n + 1 , then
    acquire the optimal value of Z ˙ = Z ˙ ( [ n + 1 ] ) and the schedule σ ˙ by
    the H L P rule;
Calculate d o p t ;
   If j = 1 , then
    obtain the optimal value of Z ˙ = Z ˙ ( [ 1 ] ) and the schedule σ ˙ by the H L P rule;
    Calculate d o p t .
Step 3: Choose the minimum value Z ˙ = min { Z ˙ [ j ] , j = 1 , 2 , , n + 1 } , and obtain the
  corresponding schedule σ ˙ , d o p t and j .

4. Results of SLK

Lemma 4.
If C [ i ] d [ i ] , then C [ i + 1 ] d [ i + 1 ] , if C [ i ] d [ i ] , then C [ i 1 ] d [ i 1 ] .
Proof. 
For a given schedule σ , under the S L K model, we have d [ i ] = b [ i ] + q o p t . If C [ i ] d [ i ] , we can get C i d i C i 1 + b i b i + q o p t C i 1 + b i q o p t C i 1 + b i + b i + 1 q o p t + b i + 1 C i + b i + 1 d i + 1 , i.e., C i + 1 d i + 1 . Similarly, if C [ i ] d [ i ] , we can also get C i d i C i 1 + b i b i + q o p t C i 1 q o p t C i 1 q o p t + b i 1 , i.e., C i 1 d i 1 . □
Lemma 5.
For a given schedule σ = ( T ˇ [ 1 ] , T ˇ [ 2 ] , , T ˇ [ n ] ) , the optimal value of the common flow-allowance q o p t is decided by the start time of the hth job, i.e., q o p t = S [ h ] , where h satisfies both
i = 1 h μ i i = h + 1 n ν i + i = 1 n ω i 0 a n d i = 1 h 1 μ i i = h n ν i + i = 1 n ω i 0 .
Proof. 
Similar to the proof of Lemma 2. □
Remark 2.
Similarly, if h does not meet the above both inequalities, we need to set q o p t = 0 .
Next, we investigate the problem 1 d m a , S L K , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i q o p t ) . Under Lemmas 3–5, the following cases will be considered.
Case 5. If j h , the schedule is identical to Case 1, where q o p t = S h = i = 1 j 1 p i + φ t + i = j h 1 ε i p i and φ t = t 0 + α i = 1 j 1 p i . Then the objective function is expressed as
Z ¨ j = i = 1 n ( μ i E i ˜ + ν i L i ˜ + ω i q o p t ) = i = 1 h μ i d i C i + i = h + 1 n υ i C i d i + i = 1 n ω i q o p t = t 0 m = 1 j 1 μ m + i = 1 j 1 p i α m = 1 j 1 μ m + i = j h 1 ε i p i m = 1 j 1 μ m + i = 1 j 1 p i m = 1 i μ m + i = j h 1 ε i p i m = j i μ m + i = h n ε i p i m = i + 1 n ν m + i = 1 j 1 p i m = 1 n ω m + t 0 m = 1 n ω m + i = 1 j 1 p i α m = 1 n ω m + i = j h 1 ε i p i m = 1 n ω m = i = 1 j 1 p i α m = 1 j 1 μ m + m = 1 i 1 μ m + m = 1 n ω m + α m = 1 n ω m + i = j h 1 ε i p i m = 1 i μ m + m = 1 n ω m + i = h n ε i p i m = i + 1 n υ m + t 0 m = 1 j 1 μ m + m = 1 n ω m = Z ¨ 3 j + f 1 t 0
where f 1 t 0 is the same as it in Case 1, and
Z ¨ 3 j = i = 1 j 1 p i m = 1 i 1 μ m + α m = 1 j 1 μ m + m = 1 n ω m + α m = 1 n ω m + i = j h 1 ε i p i m = 1 i ω m + m = 1 n ω m + i = h n ε i p i m = i + 1 n ν m
is only related to the position of j. Therefore, for a known j, minimizing Z ¨ j is identical to minimizing Z ¨ 3 j , and Z ¨ 3 j can be minimized by resolving the assignment problem. Let ρ i l represent the weight of job T ˇ i at location l, i.e.,
ρ i l = p i m = 1 l 1 μ m + α m = 1 j 1 μ m + m = 1 n ω m + α m = 1 n ω m 1 l j 1 ε i p i m = 1 l μ m + m = 1 n ω m j l h 1 ε i p i m = l + 1 n ν m h l n
The problem can be translated into the assignment problem below:
M i n Z ¨ 3 j = i = 1 n l = 1 n ρ i l χ i l
s . t i = 1 n χ i l = 1 l = 1 , , n
l = 1 n χ i l = 1 i = 1 , , n
χ i l = 0 o r 1 1 i , l n
where χ i l is a 0 or 1 variable, if the job T ˇ i is at position l, χ i l = 1 , otherwise χ i l = 0 .
Case 6. If j > h , the schedule is the same as Case 2, where q o p t = S h = i = 1 h 1 p i , and φ t = t 0 + α i = 1 j 1 p i . Then the objective function is stated as
Z ¨ j = i = 1 n ( μ i E i ˜ + ν i L i ˜ + ω i q o p t ) = i = 1 h μ i d i C i + i = h + 1 j 1 ν i C i d i + i = j n ν i C i d i + i = 1 n ω i q o p t = i = 1 h 1 p i m = 1 i μ m + i = h j 1 p i m = i + 1 j 1 ν m + i = h j 1 p i m = j n ν m + t 0 m = j n ν m + i = 1 j 1 p i α m = j n ν m + i = j n ε i p i m = i + 1 n ν m + i = 1 h 1 p i m = 1 n ω m = i = 1 h 1 p i m = 1 i μ m + α m = j n ν m + m = 1 n ω m + i = h j 1 p i m = i + 1 n ν m + α m = j n ν m + i = j n ε i p i m = i + 1 n ν m + t 0 m = j n ν m = Z ¨ 4 j + f 2 t 0
where f 2 t 0 same as it in Case 2, and
Z ¨ 4 j = i = 1 h 1 p i m = 1 i μ m + α m = j n ν m + m = 1 n ω m + i = h j 1 p i m = i + 1 n ν m + α m = j n ν m + i = j n ε i p i m = i + 1 n ν m
is only related to the location of the maintenance activity. Similarly to Case 5 (i.e., j h ), let
Ω i l = p i m = 1 l μ m + α m = j n ν m + m = 1 n ω m 1 l h 1 p i m = l + 1 n ν m + α m = j n ν m h l j 1 ε i p i m = l + 1 n ν m j l n
The problem can be translated into the assignment problem below:
M i n Z ¨ 4 j = i = 1 n l = 1 n Ω i l ς i l
s . t i = 1 n ς i l = 1 l = 1 , , n
l = 1 n ς i l = 1 i = 1 , , n
ς i l = 0 o r 1 1 i , l n
where ς i l is a 0 or 1 variable, when the job T ˇ i is at position l, ς i l = 1 , otherwise ς i l = 0 .
Case 7. If j = n + 1 , the schedule is equal to Case 3, where q o p t = S h = i = 1 h 1 p i , and φ t = t 0 + α i = 1 n p i . In this case, the objective function is
Z ¨ n + 1 = i = 1 n ( μ i E i ˜ + ν i L i ˜ + ω i q o p t ) = i = 1 h μ i C h C i + i = h + 1 n υ i C i C h + i = 1 n ω i i = 1 h 1 p i = i = 1 h 1 p i m = 1 i μ m + i = h n p i m = i + 1 n υ m + i = 1 h 1 p i m = 1 n ω m = i = 1 h 1 p i m = 1 i μ m + m = 1 n ω m + i = h n p i m = i + 1 n υ m = i = 1 n p i τ i
where
τ i = m = 1 i μ m + m = 1 n ω m 1 i h 1 m = i + 1 n υ m h i n
The optimal schedule σ and the minimum value of the objective function Z ¨ can be found by matching the minimum weight τ i with the maximum processing time p i by using Lemma 3.
Case 8. If j = 1 , the schedule is identical with Case 4, where q o p t = S h = i = 1 h 1 ε i p i + φ t and φ t = t 0 . In this case, the objective function is
Z ¨ 1 = i = 1 n ( μ i E i ˜ + ν i L i ˜ + ω i q o p t ) = i = 1 h μ i d i C i + i = h + 1 n ν i C i d i + i = 1 n ω i q o p t = i = 1 h 1 ε i p i m = 1 i μ m + i = h n ε i p i m = i + 1 n ν m + i = 1 h 1 ε i p i m = 1 n ω m + t 0 m = 1 n ω m = i = 1 h 1 ε i p i m = 1 i μ m + m = 1 n ω m + i = h n ε i p i m = i + 1 n ν m + t 0 m = 1 n ω m = i = 1 n ε i p i τ i + f 3 t 0
where f 3 t 0 is the same as it in Case 4. Similar to Case 7, the minimum weight τ i with the maximum ε i p i can minimize the objective function Z ¨ by using Lemma 3.
From above cases, the optimal solution of 1 d m a , S L K , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i q o p t ) can be obtained by the next algorithm.
Theorem 2.
The optimal solution of 1 d m a , S L K , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i q o p t ) can be obtained in O n 4 time by using Algorithm 2.
Proof. 
Similar to the proof of Theorem 1. □
Algorithm 2: Solution of 1 d m a , S L K , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i q o p t ) .
Initialization: Let Z ¨ = , σ ¨ = 0 , q o p t = 0 and j = 0 .
Step 1: Calculate h from Lemma 5.
Step 2: For j = 1 n + 1
   If j h , then
    obtain the minimum value Z ¨ j and the schedule σ ¨ by using (14)–(17);
     If Z ¨ j < Z ¨ , then
     let Z ¨ = Z ¨ j , j = j , q o p t = q o p t and σ ¨ = σ ¨ ;
   If j > h , then
    obtain the minimum value Z ¨ j and the schedule σ ¨ by using (18)–(21);
     If Z ¨ j < Z ¨ , then
     let Z ¨ = Z ¨ j , j = j , q o p t = q o p t and σ ¨ = σ ¨ ;
   If j = n + 1 , then
    acquire the optimal value of Z ¨ = Z ¨ ( [ n + 1 ] ) and the schedule σ ¨ by the H L P rule;
    Calculate q o p t ;
   If j = 1 , then
    obtain the optimal value of Z ¨ = Z ¨ ( [ 1 ] ) and the schedule σ ¨ by the H L P rule;
Calculate q o p t .
Step 3: Choose the minimum value Z ¨ = min { Z ¨ [ j ] , j = 1 , 2 , , n + 1 } , and obtain the
  corresponding schedule σ ¨ , q o p t and j .

5. An Example and Computational Experiments

5.1. An Example

Let n = 7 , t 0 = 2 , α = 0.5 , μ i = ω i = ν i , respectively are 7, 4, 3, 1, 5, 13, 2, and the remaining data are given in Table 2.
Solution: according to Lemmas 2 and 5, we can get h = 4 and μ 1 = 7 , μ 2 = 4 , μ 3 = 3 , ω 4 = 1 , ν 5 = 5 , ν 6 = 13 , ν 7 = 2 .
First, we consider the case of the common due date d o p t .
Case 9. j = 2 : from Equation (5), we can get
λ i l = p i α μ 1 + n ω m + α n ω m l = 1 ε i p i m = 1 l 1 μ m + n ω m l = 2 , 3 , 4 ε i p i m = l n υ m l = 5 , 6 , 7
and it is easy to see that
λ 11 = p 1 α μ 1 + n ω 4 + α n ω 4 = 9 × 0.5 × 7 + 7 × 1 + 0.5 × 7 × 1 = 126 ;
λ 12 = ε 1 p 1 μ 1 + n ω 4 = 0.7 × 9 × 7 + 7 × 1 = 88.2 ;
λ 13 = ε 1 p 1 μ 1 + μ 2 + n ω 4 = 0.7 × 9 × 7 + 4 + 7 × 1 = 113.4 ;
λ 14 = ε 1 p 1 μ 1 + μ 2 + μ 3 + n ω 4 = 0.7 × 9 × 7 + 4 + 3 + 7 × 1 = 132.4 ;
λ 15 = ε 1 p 1 ν 5 + ν 6 + ν 7 = 0.7 × 9 × 5 + 13 + 2 = 126 ;
λ 16 = ε 1 p 1 ν 6 + ν 7 = 0.7 × 9 × 13 + 2 = 94.5 ;
λ 17 = ε 1 p 1 ν 7 = 0.7 × 9 × 2 = 12.6 .
Similarly, the rest of the values of λ i l are given in Table 3,
We can obtain the optimal schedule σ ˙ = [ T ˇ 2 , T ˇ 7 , T ˇ 1 , T ˇ 6 , T ˇ 4 , T ˇ 3 , T ˇ 5 ] by solving the assignment problem (6), and Z ˙ 1 ( [ 2 ] ) = 88.2 + 19.8 + 84 + 75 + 63 + 64.8 + 64 = 458.8 , φ ( t ) = t 0 + α p 2 = 2 + 0.5 × 11 = 7.5 , d o p t = p 2 + φ ( t ) + ε 7 p 7 + ε 1 p 1 + ε 6 p 6 = 11 + 7.5 + 0.4 × 8 + 0.7 × 9 + 0.3 × 12 = 31.6 , f 1 ( t 0 ) = t 0 ( μ 7 + n ω 4 ) = 2 × ( 7 + 7 × 1 ) = 28 , so Z ˙ ( [ 2 ] ) = Z ˙ 1 ( [ 2 ] ) + f 1 ( t 0 ) = 458.8 + 28 = 486.8 .
j = 3 : Similarly, we can easily find the values of λ i l from Equation (5), and solve the assignment problem (6) to obtain the optimal schedule σ ˙ = [ T ˇ 1 , T ˇ 7 , T ˇ 2 , T ˇ 6 , T ˇ 4 , T ˇ 3 , T ˇ 5 ] , and Z ˙ 1 ( [ 3 ] ) = 562.6 , φ ( t ) = t 0 + α ( p 1 + p 7 ) = 2 + 0.5 × ( 9 + 8 ) = 10.5 , d o p t = p 1 + p 7 + φ ( t ) + ε 2 p 2 + ε 6 p 6 = 9 + 8 + 10.5 + 0.9 × 11 + 0.3 × 12 = 41 , f 1 ( t 0 ) = t 0 ( μ 1 + μ 2 + n ω 4 ) = 2 × ( 7 + 4 + 7 × 1 ) = 36 , so Z ˙ ( [ 3 ] ) = Z ˙ 1 ( [ 3 ] ) + f 1 ( t 0 ) = 562.6 + 36 = 598.6 .
j = 4 : Similarly, we can easily find the values of λ i l from Equation (5), and solve the assignment problem (6) to obtain the optimal schedule σ ˙ = [ T ˇ 1 , T ˇ 7 , T ˇ 3 , T ˇ 6 , T ˇ 4 , T ˇ 5 , T ˇ 2 ] , and Z ˙ 1 ( [ 4 ] ) = 754.3 , φ ( t ) = t 0 + α ( p 1 + p 7 + p 3 ) = 2 + 0.5 × ( 9 + 8 + 6 ) = 13.5 , d o p t = p 1 + p 7 + p 3 + φ ( t ) + ε 6 p 6 = 9 + 8 + 6 + 13.5 + 0.3 × 12 = 40.1 , f 1 ( t 0 ) = t 0 ( μ 1 + μ 2 + μ 3 +   n ω 4 ) = 2 × ( 7 + 4 + 3 + 7 × 1 ) = 42 , so Z ˙ ( [ 4 ] ) = Z ˙ 1 ( [ 4 ] ) + f 1 ( t 0 ) = 754.3 + 42 = 796.3 .
Case 10.  j = 5 : from Equation (9), the values of δ i l are known, similarly, we can obtain the optimal schedule σ ˙ = [ T ˇ 2 , T ˇ 7 , T ˇ 4 , T ˇ 1 , T ˇ 5 , T ˇ 6 , T ˇ 3 ] by solving the assignment problem (10), and Z ˙ 2 ( [ 5 ] ) = 929.8 , φ ( t ) = t 0 + α ( p 2 + p 7 + p 4 + p 1 ) = 2 + 0.5 × ( 11 + 8 + 10 + 9 ) = 21 , d o p t = p 2 + p 7 + p 4 + p 1 = 11 + 8 + 10 + 9 = 38 , f 2 ( t 0 ) = t 0 ( ν 5 + ν 6 + ν 7 ) = 2 × ( 5 + 13 + 2 ) = 40 , so Z ˙ ( [ 5 ] ) = Z ˙ 2 ( [ 5 ] ) + f 2 ( t 0 ) = 929.8 + 40 = 969.8 .
j = 6 : Similarly, we can obtain the values of δ i l from Equation (9) and tackle the assignment problem (10) to get the optimal schedule σ ˙ = [ T ˇ 3 , T ˇ 1 , T ˇ 4 , T ˇ 2 , T ˇ 6 , T ˇ 7 , T ˇ 5 ] , and Z ˙ 2 ( [ 6 ] ) = 1047.2 , φ ( t ) = t 0 + α ( p 3 + p 1 + p 4 + p 2 + p 6 ) = 2 + 0.5 × ( 6 + 9 + 10 + 11 + 12 ) = 26 , d o p t = p 3 + p 1 + p 4 + p 2 = 6 + 9 + 10 + 11 = 36 , f 2 ( t 0 ) = t 0 ( ν 6 + ν 7 ) = 2 × ( 13 + 2 ) = 30 , so Z ˙ ( [ 6 ] ) = Z ˙ 2 ( [ 6 ] ) + f 2 ( t 0 ) = 1047.2 + 30 = 1077.2 .
j = 7 : Similarly, we can get the values of δ i l from Equation (9) and tackle the assignment problem (10) to acquire the optimal schedule σ ˙ = [ T ˇ 3 , T ˇ 2 , T ˇ 4 , T ˇ 6 , T ˇ 7 , T ˇ 1 , T ˇ 5 ] , and Z ˙ 2 ( [ 7 ] ) = 898 , φ ( t ) = t 0 + α ( p 3 + p 2 + p 4 + p 6 + p 7 + p 1 ) = 2 + 0.5 × ( 6 + 11 + 10 + 12 + 8 + 9 ) = 30 , d o p t = p 3 + p 2 + p 4 + p 6 = 6 + 11 + 10 + 12 = 39 , f 2 ( t 0 ) = t 0 ν 7 = 2 × 2 = 4 , so Z ˙ ( [ 7 ] ) = Z ˙ 2 ( [ 7 ] ) + f 2 ( t 0 ) = 898 + 4 = 902 .
Case 11.  j = 8 : From (12), we can obtain π 1 = n ω 4 = 7 × 1 = 7 , π 2 = μ 1 + n ω 4 = 7 + 7 × 1 = 14 , π 3 = μ 1 + μ 2 + n ω 4 = 7 + 4 + 7 × 1 = 18 , π 4 = μ 1 + μ 2 + μ 3 + n ω 4 = 7 + 4 + 3 + 7 × 1 = 21 , π 5 = ν 5 + ν 6 + ν 7 = 5 + 13 + 2 = 20 , π 6 = ν 6 + ν 7 = 13 + 2 = 15 , π 7 = ν 7 = 2 , according to the H L P rule the optimal schedule is σ ˙ = [ T ˇ 6 , T ˇ 2 , T ˇ 1 , T ˇ 3 , T ˇ 7 , T ˇ 4 , T ˇ 5 ] , and φ ( t ) = t 0 + α ( p 6 + p 2 + p 1 + p 3 + p 7 + p 4 + p 5 ) = 2 + 0.5 × ( 12 + 11 + 9 + 6 + 8 + 10 + 15 ) = 37.5 , d o p t = p 6 + p 2 + p 1 + p 3 = 12 + 11 + 9 + 6 = 38 , Z ˙ ( [ 8 ] ) = p 6 π 1 + p 2 π 2 + p 1 π 3 + p 3 π 4 + p 7 π 5 + p 4 π 6 + p 5 π 7 = 12 × 7 + 11 × 14 + 9 × 18 + 6 × 21 + 8 × 20 + 10 × 15 + 15 × 2 = 866 .
Case 12.  j = 1 : Similar to Case 11, we can know π 1 = 7 , π 2 = 14 , π 3 = 18 , π 4 = 21 , π 5 = 20 , π 6 = 15 , π 6 = 2 , meanwhile ε 1 p 1 = 0.7 × 9 = 6.3 , ε 2 p 2 = 0.9 × 11 = 9.9 , ε 3 p 3 = 0.8 × 6 = 4.8 , ε 4 p 4 = 0.5 × 10 = 5 , ε 5 p 5 = 0.2 × 15 = 3 , ε 6 p 6 = 0.3 × 12 = 3.6 , ε 7 p 7 = 0.4 × 8 = 3.2 , using the H L P rule the optimal schedule is σ ˙ = [ T ˇ 1 , T ˇ 4 , T ˇ 6 , T ˇ 5 , T ˇ 7 , T ˇ 3 , T ˇ 2 ] , and f 3 ( t 0 ) = t 0 n ω 4 = 2 × 7 × 1 = 14 , φ ( t ) = t 0 = 2 , d o p t = ε 1 p 1 + ε 4 p 4 + ε 6 p 6 + ε 5 p 5 + φ ( t ) = 6.3 + 5 + 3.6 + 3 + 2 = 19.9 , so Z ˙ ( [ 1 ] ) = ε 1 p 1 π 1 + ε 4 p 4 π 2 + ε 6 p 6 π 3 + ε 5 p 5 π 4 + ε 7   p 7 π 5 + ε 3 p 3 π 6 + ε 2 p 2 π 7 + f 3 ( t 0 ) = 6.3 × 7 + 5 × 14 + 3.6 × 18 + 3 × 21 + 3.2 × 20 + 4.8 × 15 + 9.9 × 2 + 14 = 397.7 + 14 = 411.7 .
Next, we consider the case of the common flow allowance q o p t .
Case 13.  j = 2 : Similar to Case 9, the values of ρ i l can been gathered from (16), and we can obtain the optimal schedule σ ¨ = [ T ˇ 6 , T ˇ 7 , T ˇ 1 , T ˇ 5 , T ˇ 3 , T ˇ 2 , T ˇ 4 ] by solving the assignment problem (17), and Z ¨ 3 ( [ 2 ] ) = 363.4 , φ ( t ) = t 0 + α p 6 = 2 + 0.5 × 12 = 8 , q o p t = p 6 + φ ( t ) + ε 7 p 7 + ε 1 p 1 = 12 + 8 + 0.4 × 8 + 0.7 × 9 = 29.5 , f 1 ( t 0 ) = 28 , so Z ¨ ( [ 2 ] ) = Z ¨ 3 ( [ 2 ] ) + f 1 ( t 0 ) = 363.4 + 28 = 391.4 .
j = 3 : Similarly, we can easily find the values of ρ i l from (16), and solve the assignment problem (17) so that we can obtain the optimal schedule σ ¨ = [ T ˇ 1 , T ˇ 7 , T ˇ 2 , T ˇ 6 , T ˇ 3 , T ˇ 5 , T ˇ 4 ] , and Z ¨ 3 ( [ 3 ] ) = 473 , φ ( t ) = t 0 + α ( p 1 + p 7 ) = 2 + 0.5 × ( 9 + 8 ) = 10.5 , q o p t = p 1 + p 7 + φ ( t ) + ε 2 p 2 = 9 + 8 + 10.5 + 0.9 × 11 = 37.4 , f 1 ( t 0 ) = 36 , so Z ¨ ( [ 3 ] ) = Z ¨ 3 ( [ 3 ] ) + f 1 ( t 0 ) = 473 + 36 = 509 .
j = 4 : Similarly, we can obtain the values of ρ i l from (16), and tackle the assignment problem (17) to acquire the optimal schedule σ ¨ = [ T ˇ 1 , T ˇ 7 , T ˇ 3 , T ˇ 6 , T ˇ 4 , T ˇ 5 , T ˇ 2 ] , and Z ¨ 3 ( [ 4 ] ) = 648.5 , φ ( t ) = t 0 + α ( p 1 + p 7 + p 3 ) = 2 + 0.5 × ( 9 + 8 + 6 ) = 13.5 , q o p t = p 1 + p 7 + p 3 + φ ( t ) = 9 + 8 + 6 + 13.5 = 36.5 , f 1 ( t 0 ) = 42 , so Z ¨ ( [ 4 ] ) = Z ¨ 3 ( [ 4 ] ) + f 1 ( t 0 ) = 648.5 + 42 = 690.5 .
Case 14.  j = 5 : from (20), the values of Ω i l are given, similarly, we can get the optimal schedule σ ¨ = [ T ˇ 2 , T ˇ 7 , T ˇ 3 , T ˇ 1 , T ˇ 5 , T ˇ 6 , T ˇ 4 ] by solving the assignment problem (21), and Z ¨ 4 ( [ 5 ] ) = 970.2 , φ ( t ) = t 0 + α ( p 2 + p 7 + p 3 + p 1 ) = 2 + 0.5 × ( 11 + 8 + 6 + 9 ) = 19 , q o p t = p 2 + p 7 + p 3 = 11 + 8 + 6 = 25 , f 2 ( t 0 ) = 40 , so Z ¨ ( [ 5 ] ) = Z ¨ 4 ( [ 5 ] ) + f 2 ( t 0 ) = 970.2 + 40 = 1010.2 .
j = 6 : Similarly, we can obtain the values of Ω i l from (20), and tackle the assignment problem (21) to know the optimal schedule σ ¨ = [ T ˇ 2 , T ˇ 1 , T ˇ 3 , T ˇ 5 , T ˇ 6 , T ˇ 7 , T ˇ 4 ] , and Z ¨ 4 ( [ 6 ] ) = 1088 , φ ( t ) = t 0 + α ( p 2 + p 1 + p 3 + p 5 + p 6 ) = 2 + 0.5 × ( 11 + 9 + 6 + 15 + 12 ) = 28.5 , q o p t = p 2 + p 1 + p 3 = 11 + 9 + 6 = 26 , f 2 ( t 0 ) = 30 , so Z ¨ ( [ 6 ] ) = Z ¨ 4 ( [ 6 ] ) + f 2 ( t 0 ) = 1088 + 30 = 1118 .
j = 7 : Similarly, we can obtain the values of Ω i l from (20), and tackle the assignment problem (21) to acquire the optimal schedule σ ¨ = [ T ˇ 2 , T ˇ 1 , T ˇ 3 , T ˇ 5 , T ˇ 7 , T ˇ 6 , T ˇ 4 ] , and Z ¨ 4 ( [ 7 ] ) = 832 , φ ( t ) = t 0 + α ( p 2 + p 1 + p 3 + p 5 + p 7 + p 6 ) = 2 + 0.5 × ( 11 + 9 + 6 + 15 + 8 + 12 ) = 32.5 , q o p t = p 2 + p 1 + p 3 = 11 + 9 + 6 = 26 , f 2 ( t 0 ) = 4 , so Z ¨ ( [ 7 ] ) = Z ¨ 4 ( [ 7 ] ) + f 2 ( t 0 ) = 832 + 4 = 836 .
Case 15. j = 8 : From (23), we can get τ 1 = μ 1 + n ω 4 = 7 + 7 × 1 = 14 , τ 2 = μ 1 + μ 2 + n ω 4 = 7 + 4 + 7 × 1 = 18 , τ 3 = μ 1 + μ 2 + μ 3 + n ω 4 = 7 + 4 + 3 + 7 × 1 = 21 , τ 4 = ν 5 + ν 6 + ν 7 = 5 + 13 + 2 = 20 , τ 5 = ν 6 + ν 7 = 13 + 2 = 15 , τ 6 = ν 7 = 2 , τ 7 = 0 , according to the H L P rule the optimal schedule is σ ¨ = [ T ˇ 2 , T ˇ 1 , T ˇ 3 , T ˇ 7 , T ˇ 4 , T ˇ 6 , T ˇ 5 ] , and φ ( t ) = t 0 + α ( p 2 + p 1 + p 3 + p 7 + p 4 + p 6 + p 5 ) = 2 + 0.5 × ( 11 + 9 + 6 + 8 + 10 + 12 + 15 ) = 37.5 , q o p t = p 2 + p 1 + p 3 = 11 + 9 + 6 = 26 , Z ¨ ( [ 8 ] ) = p 2 τ 1 + p 1 τ 2 + p 3 τ 3 + p 7 τ 4 + p 4 τ 5 + p 6 τ 6   + p 5 τ 7 = 11 × 14 + 9 × 18 + 6 × 21 + 8 × 20 + 10 × 15 + 12 × 2 + 15 × 0 = 776 .
Case 16.  j = 1 : Similar to Cases 12 and 15, we can obtain τ 1 = 14 , τ 2 = 18 , τ 3 = 21 , τ 4 = 20 , τ 5 = 15 , τ 6 = 2 , τ 7 = 0 and ε 1 p 1 = 6.3 , ε 2 p 2 = 9.9 , ε 3 p 3 = 4.8 , ε 4 p 4 = 5 , ε 5 p 5 = 3 , ε 6 p 6 = 3.6 , ε 7 p 7 = 3.2 , and according to the H L P rule the optimal schedule is σ ¨ = [ T ˇ 4 , T ˇ 6 , T ˇ 5 , T ˇ 7 , T ˇ 3 , T ˇ 1 , T ˇ 2 ] , and f 3 ( t 0 ) = 14 , φ ( t ) = t 0 = 2 , q o p t = ε 4 p 4 + ε 6 p 6 + ε 5 p 5 + φ ( t )   = 5 + 3.6 + 3 + 2 = 13.6 , so Z ¨ ( [ 1 ] ) = ε 4 p 4 τ 1 + ε 6 p 6 τ 2 + ε 5 p 5 τ 3 + ε 7 p 7 τ 4 + ε 3 p 3 τ 5   + ε 1 p 1 τ 6 + ε 2 p 2 τ 7 + f 3 ( t 0 ) = 5 × 14 + 3.6 × 18 + 3 × 21 + 3.2 × 20 + 4.8 × 15 + 6.3 × 2 + 9.9 × 0 + 14 = 346.4 + 14 = 360.4 .
After the above calculation, the results are shown in Table 4.
We can see that the optimal solution to the question about C O N is j = 1 , σ ˙ = [ T ˇ 1 , T ˇ 4 , T ˇ 6 , T ˇ 5 , T ˇ 7 , T ˇ 3 , T ˇ 2 ] , d o p t = 19.9 and the optimal solution to the question about S L K is j = 1 , σ ¨ = [ T ˇ 4 , T ˇ 6 , T ˇ 5 , T ˇ 7 , T ˇ 3 , T ˇ 1 , T ˇ 2 ] , q o p t = 13.6 .

5.2. Computational Experiments

To test the validity of Algorithms 1 and 2, the examples are generated randomly. Microsoft visual C++ 2022 was applied to code Algorithms 1 and 2. For every problem size, 20 cases were created and coded on a PC with a 3.10 GHz CPU and 8.00 GB RAM. The features of the examples are listed below:
(1)
n = 30 , 40 , 50 , 60 , 70 , 80 , 90 , 100 , 110 , 120 , 130 , 140 , 150 , 160 , 170 , 180 , 190 , 200 , t 0 = 5 and α = 0.1 ;
(2)
p i ( i = 1 , 2 , , n ) is evenly distributed over [1, 100];
(3)
ε i ( i = 1 , 2 , , n ) is evenly distributed over [0.5, 0.95];
(4)
μ i , ν i and ω i ( i = 1 , 2 , , n ) are evenly distributed over [1, 50].
The computational tests for Algorithms 1 and 2 are given as follows. The minimum (min), average (mean) and maximum (max) CPU times (milliseconds (ms)) are shown in Table 5. From Table 5, we can see that Algorithms 1 and 2 are effective and their CPU times increase moderately as n increases from 30 to 200, and the maximum CPU time is 830,682.20 ms for n = 200 .

6. Conclusions

In this paper, we studied the single-machine due-date assignment problem with a deteriorating maintenance activity. Under the common and slack due-date assignments, the purpose of the problem is to find the optimal job schedule, the position of the maintenance activity, the optimal value of the common due date or the optimal value of the common flow-allowance so that the linear weighted sum of earliness, tardiness and due-date assignment cost is minimized. The problem is proved to be solved in polynomial time. Through computational complexity analysis and experimentation, we demonstrated that the proposed algorithm (approach) performs very well. In the future, other non-regular objectives with deteriorating maintenance activity can be studied, or multi-machine (flexible job shop, see Zhang et al. [39] and Song et al. [40]) due-date assignment problems with deteriorating maintenance activity can be addressed.

Author Contributions

Conceptualization, W.W., D.-Y.L. and J.-B.W.; methodology, W.W. and D.-Y.L.; investigation, J.-B.W.; writing, review and editing, W.W. and D.-Y.L.; supervision, J.-B.W.; writing-original draft preparation, W.W.. All authors have read and agreed to the published version of the manuscript.

Funding

This Work was supported by LiaoNing Revitalization Talents Program (grant no. XLYC2002017).

Data Availability Statement

The data used to support this paper are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof. 
When j = n + 1 , we first consider a schedule σ , and the common due date d o p t satisfies C k < d o p t < C k + 1 , Z ˙ is the corresponding value of objective function, let x ˙ = d o p t C k , y ˙ = C k + 1 d o p t ( x ˙ , y ˙ 0 ), Z ˙ and Z ˙ are the corresponding objective function values of d o p t = C k and d o p t = C k + 1 , respectively, therefore
Z ˙ = Z ˙ + x ˙ i = 1 n ν i ω i i = 1 k μ i ν i
Z ˙ = Z ˙ y ˙ i = 1 n ν i ω i i = 1 k μ i ν i .
Obviously, if i = 1 n ν i ω i i = 1 k μ i + ν i , Z ˙ Z ˙ , otherwise Z ˙ Z ˙ , this means that d o p t is equal to the completion time of some job.
We assume that the value of d o p t coincides with the completion time of the job at the hth position, i.e., d o p t = C [ h ] , then for θ 0 , the common due date shifted θ unit to the left, and it is obtained
Z ˙ C h + θ , σ Z ˙ C h , σ 0 ,
we have
i = 1 h μ i i = h + 1 n ν i + i = 1 n ω i 0 .
If the common due date shifted θ unit to the right, we have
Z ˙ C h θ , σ Z ˙ C h , σ 0 ,
and
i = 1 h 1 μ i i = h n ν i + i = 1 n ω i 0
Thus, h satisfies
i = 1 h μ i i = h + 1 n ν i + i = 1 n ω i 0 a n d i = 1 h 1 μ i i = h n ν i + i = 1 n ω i 0 .
Similarly, the above result can be obtained when j = 1 , , n , i.e., for any given schedule, d o p t = C h . □

References

  1. Lee, C.Y.; Leon, V.J. Machine scheduling with a rate-modifying activity. Eur. J. Oper. Res. 2001, 128, 119–128. [Google Scholar] [CrossRef]
  2. Ma, Y.; Chu, C.; Zuo, C. A survey of scheduling with deterministic machine availability constraints. Comput. Ind. Eng. 2010, 58, 199–211. [Google Scholar] [CrossRef]
  3. Strusevich, V.A.; Rustogi, K. Scheduling with Time-Changing Effects and Rate-Modifying Activities; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  4. Mosheiova, G.; Oron, D. Due-date assignment and maintenance activity scheduling problem. Math. Comput. Model. 2006, 44, 1053–1057. [Google Scholar] [CrossRef]
  5. Wang, X.-Y.; Wang, M.-Z. Single machine common flow allowance scheduling with a rate-modifying activity. Comput. Ind. Eng. 2010, 59, 898–902. [Google Scholar] [CrossRef]
  6. Yin, Y.; Cheng, T.C.E.; Xu, D.; Wu, C.-C. Common due date assignment and scheduling with a rate-modifying activity to minimize the due date, earliness, tardiness, holding, and batch delivery cost. Comput. Ind. Eng. 2012, 63, 223–234. [Google Scholar] [CrossRef]
  7. Bai, J.; Li, Z.-R.; Wang, J.-J.; Huang, X. Single machine common flow allowance scheduling with deteriorating jobs and a rate-modifying activity. Appl. Math. Model. 2014, 38, 5431–5438. [Google Scholar] [CrossRef]
  8. Cheng, M.; Xiao, S.; Luo, R.; Lian, Z. Single-machine scheduling problems with a batch-dependent aging effect and variable maintenance activities. Int. J. Prod. Res. 2018, 56, 7051–7063. [Google Scholar] [CrossRef]
  9. Detti, P.; Nicosia, G.; Pacifici, A.; De Lara, G.Z.M. Robust single machine scheduling with a flexible maintenance activity. Comput. Oper. Res. 2019, 107, 19–31. [Google Scholar] [CrossRef]
  10. Wang, J.-B.; Hu, Y.; Zhang, B. Common due-window assignment for single-machine scheduling with generalized earliness/tardiness penalties and a rate-modifying activity. Eng. Optim. 2021, 53, 496–512. [Google Scholar] [CrossRef]
  11. Zhao, X.; Xu, J.; Wang, J.-B.; Li, L. Bicriteria common flow allowance scheduling with aging effect, convex resource allocation, and a rate-modifying activity on a single machine. Asia-Pac. J. Oper. Res. 2022, 39, 2150046. [Google Scholar] [CrossRef]
  12. Sun, X.; Liu, T.; Geng, X.-N.; Hu, Y.; Xu, J.-X. Optimization of scheduling problems with deterioration effects and an optional maintenance activity. J. Sched. 2023. [Google Scholar] [CrossRef]
  13. Ji, P.; Li, G.; Huo, Y.; Wang, J.-B. Single-machine common flow allowance scheduling with job-dependent aging effects and a deteriorating maintenance activity. Optim. Lett. 2014, 8, 1389–1400. [Google Scholar] [CrossRef]
  14. Fan, Y.P.; Zhao, C.-L. Single machine scheduling with multiple common due date assignment and aging effect under a deteriorating maintenance activity consideration. J. Appl. Math. Comput. 2014, 46, 51–66. [Google Scholar] [CrossRef]
  15. Mor, B.; Mosheiov, G. Scheduling a deteriorating maintenance activity and due-window assignment. Comput. Oper. Res. 2015, 57, 33–40. [Google Scholar] [CrossRef]
  16. Li, S.-S.; Chen, R.X. Scheduling with rejection and a deteriorating maintenance activity on a single machine. Asia-Pac. J. Oper. Res. 2017, 34, 1750010. [Google Scholar] [CrossRef]
  17. Zhu, H.; Li, M.; Zhou, Z.J. Machine scheduling with deteriorating and resource-dependent maintenance activity. Comput. Ind. Eng. 2015, 88, 479–486. [Google Scholar] [CrossRef]
  18. Wang, J.-J.; Wang, J.-B.; Liu, F. Parallel machines scheduling with a deteriorating maintenance activity. J. Oper. Res. Soc. 2011, 62, 1898–1902. [Google Scholar] [CrossRef]
  19. He, H.; Hu, Y.; Liu, W.-W. Scheduling with deterioration effects and maintenance activities under parallel processors. Eng. Optim. 2021, 53, 2070–2087. [Google Scholar] [CrossRef]
  20. Jia, X.; Lv, D.-Y.; Hu, Y.; Wang, J.-B.; Wang, Z.; Wang, E. Slack due-window assignment scheduling problem with deterioration effects and a deteriorating maintenance activity. Asia-Pac. J. Oper. Res. 2022, 39, 2250005. [Google Scholar] [CrossRef]
  21. Liu, W.; Wang, X.; Li, L.; Zhao, P. A maintenance activity scheduling with time-and-position dependent deteriorating effects. Math. Biosci. Eng. 2022, 19, 11756–11767. [Google Scholar] [CrossRef]
  22. Zou, J.; Sui, Y.-K.; Gao, J.; Zhang, X.-Z. Parallel machines scheduling with deteriorating maintenance activities and job rejection. Asia-Pac. J. Oper. Res. 2023, 40, 2240013. [Google Scholar] [CrossRef]
  23. Gordon, V.S.; Proth, J.-M.; Chu, C.B. A survey of the state-of-the-art of common due date assignment and scheduling research. Eur. J. Oper. Res. 2002, 139, 1–25. [Google Scholar] [CrossRef]
  24. Gordon, V.S.; Proth, J.-M.; Chu, C.B. Due date assignment and scheduling: SLK, TWK and other due date assignment models. Prod. Plan. Control 2002, 13, 117–132. [Google Scholar] [CrossRef]
  25. Qian, J.; Zhan, Y. The due date assignment scheduling problem with delivery times and truncated sum-of-processing-times-based learning effect. Mathematics 2021, 9, 3085. [Google Scholar] [CrossRef]
  26. Lu, Y.-Y.; Wang, T.-T.; Wang, R.-Q.; Li, Y. A note on due-date assignment scheduling with job-dependent learning effects and convex resource allocation. Eng. Optim. 2021, 53, 1273–1281. [Google Scholar] [CrossRef]
  27. Wang, W. Single-machine due-date assignment scheduling with generalized earliness/tardiness penalties including proportional setup times. J. Appl. Math. Comput. 2022, 68, 1013–1031. [Google Scholar] [CrossRef]
  28. Brucker, P. Scheduling Algorithms; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  29. Liu, W.; Hu, X.; Wang, X. Single machine scheduling with slack due dates assignment. Eng. Optim. 2017, 49, 709–717. [Google Scholar] [CrossRef]
  30. Wang, S.-H.; Lv, D.-Y.; Wang, J.-B. Research on position-dependent weights scheduling with delivery times and truncated sum-of-processing-times-based learning effect. J. Ind. Manag. Optim. 2023, 19, 2824–2837. [Google Scholar] [CrossRef]
  31. Wang, Y.-C.; Wang, S.-H.; Wang, J.-B. Resource allocation scheduling with position-dependent weights and generalized earliness-tardiness cost. Mathematics 2023, 11, 222. [Google Scholar] [CrossRef]
  32. Jiang, C.; Zou, D.; Bai, D.; Wang, J.-B. Proportionate flowshop scheduling with position-dependent weights. Eng. Optim. 2020, 52, 37–52. [Google Scholar] [CrossRef]
  33. Liu, W.; Yao, Y.; Jiang, C. Single-machine resource allocation scheduling with due-date assignment, deterioration effect and position-dependent weights. Eng. Optim. 2020, 52, 701–714. [Google Scholar] [CrossRef]
  34. Lv, D.Y.; Wang, J.-B. Study on proportionate flowshop scheduling with due-date assignment and position-dependent weights. Optim. Lett. 2021, 15, 2311–2319. [Google Scholar] [CrossRef]
  35. Wang, L.-Y.; Huang, X.; Liu, W.-W.; Wu, Y.; Wang, J.-B. Scheduling with position-dependent weights, due-date assignment and past-sequence-dependent setup times. Rairo-Oper. Res. 2021, 55, S2747–S2758. [Google Scholar] [CrossRef]
  36. Wang, X.; Liu, W.; Li, L.; Zhao, P.; Zhang, R. Due date assignment scheduling with positional-dependent weights and proportional setup times. Math. Biosci. Eng. 2022, 19, 5104–5119. [Google Scholar] [CrossRef]
  37. Sun, X.; Geng, X.-N.; Liu, T. Due-window assignment scheduling in the proportionate flow shop setting. Ann. Oper. Res. 2020, 292, 113–131. [Google Scholar] [CrossRef]
  38. Hardy, G.H.; Littlewood, J.E.; Polya, G. Inequalities; Cambridge University Press: Cambridge, UK, 1976. [Google Scholar]
  39. Zhang, H.; Qin, C.; Zhang, W.; Xu, Z.; Xu, G.; Gao, Z. Energy-saving scheduling for flexible job shop problem with AGV transportation considering emergencies. Systems 2023, 11, 103. [Google Scholar] [CrossRef]
  40. Song, L.; Xu, Z.; Wang, C.; Su, J. A new decision method of flexible job shop rescheduling based on WOA-SVM. Systems 2023, 11, 59. [Google Scholar] [CrossRef]
Figure 1. If j h .
Figure 1. If j h .
Systems 11 00150 g001
Figure 2. If j > h .
Figure 2. If j > h .
Systems 11 00150 g002
Figure 3. If j = n + 1 .
Figure 3. If j = n + 1 .
Systems 11 00150 g003
Figure 4. If j = 1 .
Figure 4. If j = 1 .
Systems 11 00150 g004
Table 1. Summary of C O N and S L K results.
Table 1. Summary of C O N and S L K results.
ProblemComplexityRef.
1 C O N , b i = p i i = 1 n ( μ i E [ i ] ˜ + μ i L [ i ] ˜ + ω 0 d o p t ) O ( n log n ) Brucker [28]
1 S L K , b i = p i i = 1 n ( μ i E [ i ] ˜ + μ i L [ i ] ˜ + ω 0 q o p t ) O ( n log n ) Liu et al. [29]
1 S p s d , C O N , b i = p i i = 1 n ( μ i E [ i ] ˜ + μ i L [ i ] ˜ + ω 0 d o p t ) O ( n log n ) Wang et al. [36]
1 S p s d , S L K , b i = p i i = 1 n ( μ i E [ i ] ˜ + μ i L [ i ] ˜ + ω 0 q o p t ) O ( n log n ) Wang et al. [36]
1 r m a , C O N , b i = ( p i , ε i p i ) i = 1 n ( μ 0 E [ i ] ˜ + ν 0 L [ i ] ˜ + ω 0 d o p t ) O ( n 4 ) Mosheiova and Oron [4]
1 r m a , S L K , b i = ( p i , ε i p i ) i = 1 n ( μ 0 E [ i ] ˜ + ν 0 L [ i ] ˜ + ω 0 q o p t ) O ( n 4 ) Wang and Wang [5]
1 d m a , C O N , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i d o p t ) O ( n 4 ) Theorem 1
1 d m a , S L K , b i = ( p i , ε i p i ) i = 1 n ( μ i E [ i ] ˜ + ν i L [ i ] ˜ + ω i q o p t ) O ( n 4 ) Theorem 2
Table 2. The remaining data.
Table 2. The remaining data.
T ˇ i T ˇ 1 T ˇ 2 T ˇ 3 T ˇ 4 T ˇ 5 T ˇ 6 T ˇ 7
p i 91161015128
ε i 0.7 0.9 0.8 0.5 0.2 0.3 0.4
Table 3. The values of λ i l .
Table 3. The values of λ i l .
i l 1234567
1126 88 . 2 113.4 132.3 126 94.5 12.6
2154 138.6 178.2 207.9 198 148.5 19 . 8
3 84 67.2 86.4 100.8 9672 9.6
41407090105100 75 10
52104254 63 60456
6168 50.4 64 . 8 75.6 7254 7.2
7112 44.8 57.6 67.2 64 48 6.4
Table 4. The results of the example.
Table 4. The results of the example.
d opt q opt Z ˙ ( [ j ] ) Z ¨ ( [ j ] ) σ ˙ σ ¨
j = 1 19.9 13.6 411.7 360.4 [ T ˇ 1 , T ˇ 4 , T ˇ 6 , T ˇ 5 , T ˇ 7 , T ˇ 3 , T ˇ 2 ] [ T ˇ 4 , T ˇ 6 , T ˇ 5 , T ˇ 7 , T ˇ 3 , T ˇ 1 , T ˇ 2 ]
j = 2 31.6 29.5 486.8 391.4 [ T ˇ 2 , T ˇ 7 , T ˇ 1 , T ˇ 6 , T ˇ 4 , T ˇ 3 , T ˇ 5 ] [ T ˇ 6 , T ˇ 7 , T ˇ 1 , T ˇ 5 , T ˇ 3 , T ˇ 2 , T ˇ 4 ]
j = 3 41 37.4 598.6 509 [ T ˇ 1 , T ˇ 7 , T ˇ 2 , T ˇ 6 , T ˇ 4 , T ˇ 3 , T ˇ 5 ] [ T ˇ 1 , T ˇ 7 , T ˇ 2 , T ˇ 6 , T ˇ 3 , T ˇ 5 , T ˇ 4 ]
j = 4 40.1 36.5 796.3 690.5 [ T ˇ 1 , T ˇ 7 , T ˇ 3 , T ˇ 6 , T ˇ 4 , T ˇ 5 , T ˇ 2 ] [ T ˇ 1 , T ˇ 7 , T ˇ 3 , T ˇ 6 , T ˇ 4 , T ˇ 5 , T ˇ 2 ]
j = 5 3825 969.8 1010.2 [ T ˇ 2 , T ˇ 7 , T ˇ 4 , T ˇ 1 , T ˇ 5 , T ˇ 6 , T ˇ 3 ] [ T ˇ 2 , T ˇ 7 , T ˇ 3 , T ˇ 1 , T ˇ 5 , T ˇ 6 , T ˇ 4 ]
j = 6 3626 1077.2 1118 [ T ˇ 3 , T ˇ 1 , T ˇ 4 , T ˇ 2 , T ˇ 6 , T ˇ 7 , T ˇ 5 ] [ T ˇ 2 , T ˇ 1 , T ˇ 3 , T ˇ 5 , T ˇ 6 , T ˇ 7 , T ˇ 4 ]
j = 7 3926902836 [ T ˇ 3 , T ˇ 2 , T ˇ 4 , T ˇ 6 , T ˇ 7 , T ˇ 1 , T ˇ 5 ] [ T ˇ 2 , T ˇ 1 , T ˇ 3 , T ˇ 5 , T ˇ 7 , T ˇ 6 , T ˇ 4 ]
j = 8 3826866776 [ T ˇ 6 , T ˇ 2 , T ˇ 1 , T ˇ 3 , T ˇ 7 , T ˇ 4 , T ˇ 5 ] [ T ˇ 2 , T ˇ 1 , T ˇ 3 , T ˇ 7 , T ˇ 4 , T ˇ 6 , T ˇ 5 ]
Table 5. CPU times of algorithms.
Table 5. CPU times of algorithms.
Algorithm 1 (ms)Algorithm 2 (ms)
Jobs (  n )MinMeanMaxMinMeanMax
3096.15104.12116.47105.57123.42131.36
40403.13447.04468.66456.02490.37533.89
501100.621133.921207.001800.971869.211969.54
602385.212415.262480.622406.823211.603621.18
704764.214824.514962.154427.605595.506072.50
808988.259047.699174.259546.2610,260.9711,854.23
9015,607.8915,735.9415,993.6518,732.1219,081.0520,657.20
10025,808.6426,019.1226,248.5225,948.2131,214.8034,985.24
11040,791.2340,888.9941,023.3344,216.2547,840.1550,154.56
12062,351.4062,705.7962,994.4969,263.5772,753.3176,016.59
13091,992.4092,142.2692,376.4895,431.29104,475.77112,460.39
140131,937.56132,467.57134,070.77153,157.28159,192.78165,218.85
150185,206.31186,061.38187,231.45185,069.15197,686.43204,623.32
160255,436.23256,559.22259,944.61259,783.23266,256.81279,325.82
170342,530.15342,989.64344,039.89382,981.28389,495.60391,893.45
180453,007.52455,489.67457,957.65494,021.26496,965.56501,507.23
190594,550.76597,809.51602,660.32625,462.31639,839.76661,165.27
200763,412.32780,313.79791,677.65812,131.25822,461.92830,682.20
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

Wu, W.; Lv, D.-Y.; Wang, J.-B. Two Due-Date Assignment Scheduling with Location-Dependent Weights and a Deteriorating Maintenance Activity. Systems 2023, 11, 150. https://doi.org/10.3390/systems11030150

AMA Style

Wu W, Lv D-Y, Wang J-B. Two Due-Date Assignment Scheduling with Location-Dependent Weights and a Deteriorating Maintenance Activity. Systems. 2023; 11(3):150. https://doi.org/10.3390/systems11030150

Chicago/Turabian Style

Wu, Wei, Dan-Yang Lv, and Ji-Bo Wang. 2023. "Two Due-Date Assignment Scheduling with Location-Dependent Weights and a Deteriorating Maintenance Activity" Systems 11, no. 3: 150. https://doi.org/10.3390/systems11030150

APA Style

Wu, W., Lv, D. -Y., & Wang, J. -B. (2023). Two Due-Date Assignment Scheduling with Location-Dependent Weights and a Deteriorating Maintenance Activity. Systems, 11(3), 150. https://doi.org/10.3390/systems11030150

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