Abstract 1 Introduction 2 Preliminaries 3 Parameterization by the Number of Jobs 4 Parameterization by the Number of Time Steps 5 Conclusion References

Repairing Schedules by Removing Waiting Times: A Parameterized Complexity Analysis

Niels Grüttemeier ORCID Fraunhofer IOSB-INA, Lemgo, Germany Klaus Heeger ORCID Fraunhofer IOSB-INA, Lemgo, Germany
Abstract

We consider the problem of repairing production schedules in a job-shop setting by reducing pre-planned waiting times. Herein, a schedule of all jobs is given. To compensate unforeseen disturbances, this schedule contains waiting times between the execution of two consecutive tasks of a job. Further, we assume that the schedule temporarily overloads some machines, e.g. due to reduced machine capacities because of worker sickness or (partially) broken machines. We study the problem of removing as few waiting times as possible in order to eliminate the machine overloads. After formalizing this problem, we perform an extensive analysis of its parameterized complexity with respect to several natural parameters, resulting in a detailed picture of the problem’s complexity.

Keywords and phrases:
Job shop, parallel machines, reactive scheduling
Funding:
Niels Grüttemeier: Supported by the project Datenfabrik.NRW, a project by KI.NRW, funded by the Ministry for Economics, Innovation, Digitalization and Energy of the State of North Rhine-Westphalia (MWIDE).
Klaus Heeger: Supported by the German Federal Ministry for the Environment, Nature Conservation, Nuclear Safety and Consumer Protection (BMUV) under the project Smart-E-Factory.
Copyright and License:
[Uncaptioned image] © Niels Grüttemeier and Klaus Heeger; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Parameterized complexity and exact algorithms
; Theory of computation Problems, reductions and completeness ; Mathematics of computing Combinatorics
Editors:
Pat Morin and Eunjin Oh

1 Introduction

Scheduling is a crucial task in the industrial production of goods, helping to achieve key goals such as low production costs or on-time deliveries. Therefore, scheduling is an established and very important area in the intersection of theoretical computer science, operations research, and combinatorial optimization. In classical scheduling problems, one is usually given a set of jobs (and their characteristics) to be processed in the future, and the task is to compute an “optimal” (with respect to some specified criterium) schedule of the jobs. However, many real-world scenarios are not adequatly covered by this approach, for instance when a large number of unexpected obstacles (e.g., machines breaking down, workers calling in sick, arrival of urgent customer orders, etc.) may compromise the initial plan. Recomputing a new schedule from scratch after each such obstacle is hardly an option in practice, as reassignments of jobs and operations to different resources require cancellations of expensive set-ups and lead to shop floor nervousness [21].

Instead of rescheduling, there are also other ways to handle unexpected obstacles in production planning. One such approach is scheduling under uncertainties, containing several different models of uncertainties in scheduling. One of the most common such models is stochastic scheduling, where the values of job parameters (most commonly of the processing time) are independent random variables instead of numbers, see e.g. [19, 8]. Another model is multi-scenario scheduling, where a finite number of different scenarios needs to be considered (see [18] for a survey). Another approach is online scheduling, where the machine breakdowns appear one after another and one has to update the schedule immediately, see e.g. [2, 10]. Another approach is to repair a schedule. In this approach, one starts with a given schedule and applies small changes to it, adapting the schedule to the new situation. This is done by heuristically applying generic operations on the given schedule [20, 21] or by applying data-driven techniques from the field of machine learning [22, 16, 14]. Relatively new approaches combine scheduling under uncertainties with repairing policies [13]. For a very recent survey of scheduling problems in general, we refer to Agnetis et al. [1].

While the vast majority of research on repairing schedules focusses on heuristic approaches and their experimental evaluation, this work aims to lay a foundation for the theoretical study of these problems. In order to do so, we first formalize our scheduling problem as a precise mathematical problem which we call Waiting Time Removal (WTR) using notations from stringology. The idea of WTR is that in a case where unexpected resource breakdowns occur, updated information about resource availability and job durations is perfectly known and one has to adapt the old solution to the new situation. WTR is motivated by real-world production planning scenarios in German mid-sized companies. In this work, we study the computational complexity of WTR, focussing on the parameterized complexity.

t 1 2 3 4 5 6 7 8
A 3 2 2 2 2 2 0 1
B 1 1 2 0 2 2 1 1
1 A x x B
2 A B
3 A x A B
4 A x x A B
5 A A x A B

t 1 2 3 4 5 6 7 8 A 3 2 2 2 2 2 0 1 B 1 1 2 0 2 2 1 1 1 A x x B 2 A B 3 A A B 4 A A B 5 A A x A B

Figure 1: An example of an instance of WTR (left) together with a repaired schedule (right). There are time steps t=1,,8 and jobs 1,,5 that are processed on two machines A and B. The upper part shows the machine capacities at the time steps, the lower part shows the jobs aligned at their starting time steps. At time step t=5, the capacity constraints are violated, since three jobs are processed on machine A. Removing the grey marked symbols x in Job 3 and Job 4 converts the schedule into a schedule satisfying all capacity constraints.

An example of an instance of WTR together with a repaired schedule is illustrated in Figure 1. Intuitively, our scheduling problem WTR models the following scenario: there are multiple types of different machines (e.g. machine types A and B). Similar to a job shop, each job has to be processed on machines of different types in a job-specific order for a specific amount of time. However, in contrast to the standard job shop model, there may be waiting times planned between the steps, serving as a “buffer” to cover unexpected delays. To model this, each job is represented by a string enlisting the required machines (repeating a machine if it needs to be processed for multiple time units) and the planned waiting times. As an example, the string AAxB corresponds to a job that first needs to be processed on machine A for two time units, then has one time unit waiting time (encoded by the letter “x”), and finally needs to be processed on machine B for one time unit. At each time step, a different number of machines of each type is available (e.g. at time t=1, three machines of type A might be available, meaning that up to three different jobs can be processed on machines of type A); this number will be called the capacity of the machine at that time step. The capacities model e.g. the breakdown of machines or workers required to operate the machine calling in sick. Additionally, each job is assigned with a starting time. The goal is to remove a small number of waiting times from the jobs such that the resulting schedule satisfies the capacity constraints posed on the machine types in the all steps. The idea behind keeping the number of waiting time removals small is that maintaining many waiting times increases the robustness against future disruptions.

Related Work.

To the best of our knowledge, our introduced problem WTR has not yet been studied.

One framework that is conceptually related to the idea of repairing given schedules is the approach of parameterized local search: The idea is to apply a limited number of changes on a given schedule to improve a target function value. Balzereit et al. [3] study parameterized local search for a single machine scheduling problem. While their approach mostly serves as a subroutine in a hill-climbing setting, it can also be used to repair schedules.

Another framework related to repairing schedules is reoptimization [17, 4]: The idea of reoptimization is to compute approximate solutions of perturbed instances from given optimal solutions of an original instance. There are constant factor approximations for reoptimization settings of multi-machine makespan minimization under the constraint that some jobs are not allowed to be scheduled parallel [17]. Boria and Della Croce [4] provide constant factor approximizations for multiple min-sum scheduling problems with release dates.

Our Results.

We initiate the theoretical study of WTR. As WTR is NP-hard even if the number of allowed waiting time removals is unbounded (Theorem 3), the problem is inapproximable, motivating our parameterized complexity analysis. Considering the most natural parameterizations (by the number n of jobs, the number k of deletions, the time horizont T, and the alphabet size |Σ| and combinations thereof), we derive a complete picture of the parameterized complexity of WTR. An overview of the results is shown in Table 1.

Notably, while most of our hardness results hold already for alphabets of constant size, the paraNP-hardness of WTR parameterized by T (Theorem 6) is contrasted with the fixed-parameter tractability with respect to T+|Σ| (Theorem 7). A further notable result is that parameterizing with the number n of jobs is not sufficient to ensure fixed-parameter tractability (Theorem 3), but only results in an XP-algorithm (Theorem 4).

Table 1: Our results. Each cell corresponds to parameterization by the sum of the two corresponding values. The diagonal cells correspond to parameterization by n, k, |Σ|, or T alone.
n XP(Thm. 4)W[1]-h(Thm. 3)
k FPT(Cor. 5)No Poly Ker.(Cor. 10) XP(Obs. 2)W[2]-h(Cor. 9)
|Σ| XP(Thm. 4)W[1]-h(Thm. 3) XP(Obs. 2)W[2]-h(Cor. 9) Para NP-h(Thm. 3)
T FPT(Obs. 1)Poly Kernel(Obs. 1) XP(Obs. 2)W[1]-h(Thm. 11) FPT(Thm. 7)No Poly Ker.(Thm. 8) Para NP-h(Thm. 6)
n k |Σ| T

2 Preliminaries

We use standard notations from parameterized complexity [5].

For integers ab, we let [a,b]:={iaib}. We let  denote the concatenation of strings, and we let |w| denote the length of a string w. Furthermore, for i, we let w[i] be the ith symbol of w or the empty word if i[1,|w|].

Let Σ be a finite alphabet, let xΣ, and let w(Σ{x}) be a string with m occurrences of the symbol x. For R[1,m], we let wR denote the string after removing the ith occurrences of x for each iR from w. As an example, consider axbbxcx{1,3}=abbxc, where we remove the first and the third occurrence of x while keeping the second occurrence of x. Throughout this work, we call such set R an x-removal set.

Problem Definition.

Recall that we aim to study a problem where the goal is to adjust a given plan by removing waiting times. In the following, we formalize this idea.

Let Σ be a finite alphabet, let xΣ, and let T be a number of time steps. A T-step capacity bound for Σ is a mapping c:Σ×[1,T]0. Intuitively, c(σ,t) corresponds to the capacity of symbol σ at time step t. A job (w,s) is a tuple consisting of a string w(Σ{x}) and a starting time s that satisfies s+|w|1T. Intuitively, the inequality guarantees that every job is completed by time step T. Given a set of jobs 𝒥:={(wi,si)i[1,n]}, a symbol σΣ and a time step t[1,T], we define the load of symbol σ at time step t as

λσt:=|{(w,s)Jw[ts+1]=σ}|.

A set 𝒥 of jobs satisfies the capacity constraints c if λσtc(σ,t) for every combination of σΣ and t[1,T]. For a job J and some t[s,s+|J|1], we say that the (ts+1)th letter of J is processed at time t.

Given a set of jobs 𝒥:={(wi,si)i[1,n]}, we let #xi denote the number of occurrences of the symbol x in the string wi for i[1,n]. The problem WTR is defined as follows.

Waiting Time Removal (WTR)

Input:

A number of time steps T, a set 𝒥:={(wi,si)i[1,n]} of jobs, a T-step capacity bound c, and an integer k.

Question:

Does there exist a family of sets R1,,Rn with Ri[1,#xi] for each i[1,n] and i=1n|Ri|k such that the job set 𝒥:={(wiRi,si)i[1,n]} satisfies the capacity constraints c?

For notational convenience, we define for a job J=(w,s) its length |J|:=|w|, its starting time s(J):=s, and JR:=(wR,s) for R[1,#xi].

Note that Tn essentially bounds the size of an instance (T,𝒥,c,t) of WTR: Each wi has length at most T and thus, |𝒥| can be encoded by encoding at most Tn symbols. Moreover, the load of a symbol can never exceed n and thus, c can be encoded using 𝒪(T2nlog(n)) bits. Therefore, we can observe the following.

Observation 1.

WTR admits a polynomial kernel for T+n.

Furthermore, note that WTR can be solved by the following simple brut-force algorithm: Enumerate all size-k subsets of occurrences of the symbol x in the input instance, and check whether removing precisely these occurrences provides a modified set 𝒥 that satisfies the capacity constraints c. This simple brute-force strategy implies the following.

Observation 2.

WTR can be solved in |I|k+𝒪(1) time, where |I| denotes the input size.

Graph Theory.

We use standard notation from graph theory, see e.g. [12]. More precisely, all appearing graphs are undirected unless stated otherwise. A graph G=(V,E) is k-partite if its vertex set can be partitioned into k independent sets, i.e., V=V1Vk with ViVj= for each ij and for each set Vi, there is no edge whose endpoints are both contained in Vi. For a graph G=(V,E) and V1,V2V, we denote by E[V1,V2] the edges with one endpoint in V1 and the other endpoint in V2. For a vertex vV, we denote by δ(v) the set of edges incident to v.

3 Parameterization by the Number of Jobs

We first study the parameterization by the number n of jobs. While the total number of jobs naturally appears to be a large parameter in real-world applications, we start by showing that – somewhat surprisingly – WTR is unlikely to be fixed-parameter tractable for the parameter n.

Theorem 3.

For an instance of WTR, deciding the existence of a deletion family (of arbitrary size) is NP-hard and W[1]-hard parameterized by the number of jobs, even if |Σ|=2.

Proof.

We reduce from Multicolored Independent Set which is NP-hard [11] and W[1]-hard parameterized by solution size  [15]. In Multicolored Independent Set, one is given an -partite graph G=(V1V,E) with |V1|=|V2|==|V|, and the task is to find a multicolored independent set, i.e., a set of pairwise non-adjacent vertices containing exactly one vertex from each Vi. Let G=(V1V,E) be an instance of Multicolored Independent Set. We fix an arbitrary ordering e1,,em of the edges of G as well as of the vertices of each color class Vi={vi1,,vi|V1|}. For each color class Vi and each edge epE, we define the |V1|-letter substring wip to be

wip:={0j110|V1|jif epδ(vij) for some j[1,|V1|]0|V1|otherwise.

Using these substrings, we create, for each i[1,], one job Ji=(wi,1) where

wi:=x|V1|1p=1mwip.

We define the (m|V1|+|V1|1)-capacity bound c: For every p[1,m], letter 1 has capcity one at time p|V1|. All other capacities are infinite. See Figure 2 for an example.

Figure 2: The upper part shows an example for the reduction for =3 and the graph containing edges {v11,v21}, {v13,v22}, {v22,v31}, {v23,v32}, {v12,v33}, and {v13,v31}. Red columns correspond to time steps where the capacity of letter 1 equals one while all other time steps have unboundend capacity. The vertical lines separate the different “blocks” of |V1| letters corresponding to an edge. The lower part depicts a solution corresponding to the independent set {v11,v22,v33}.

The reduction clearly runs in polynomial time and the number of jobs equals .

We continue by showing correctness. First note that for any set Ri[1,|V1|1], job JiRi contributes to the load of letter 1 at time p|V1| if and only if vi|Ri|+1 is adjacent to ep. Given a multicolored independent set I={v1j1,,vj}, we construct a deletion set ={R1,,Rn} as follows: For each job Ji, we let Ri=[1,ji1]. To see that this {(wiRi,1)i[1,n]} satisfies the capacity constraints, note that only letter 1 has a finite capacity (at times t=p|V1| for p[1,m]), so it suffices to show that the load of letter 1 does not exceed its capacity. Assume towards a contradiction that there is some time t where the load of letter 1 exceeds its capacitiy, i.e., t=p|V1| for some p[1,m] and there are two jobs JiRi and JiRi both contributing to the load of 1 at time t. This implies that ep is incident to both viji and viji, a contradiction to I being an independent set.

We conclude with the reverse direction, so let ={R1,,Rn} be a deletion family satisfying the capacity constraints. Let ji:=|Ri|[1,|V1|1]. We claim that I={v1j1,,vj} is an independent set. Assume towards a contradiction that ep={viji,viji} for some p[m] and i,i[1,]. Then at time p|V1|, both JiRi and JiRi contribute to the load of letter 1, a contradiction to the feasibility of the schedule.

While WTR is presumably not fixed-parameter tractable for n, we show that – on the positive side – WTR can be solved in polynomial time for constant values of n by providing the following XP algorithm.

Theorem 4.

WTR can be solved in O(nT2n+1) time.

Proof.

We give a dynamic program solving the problem. The dynamic programming table τ has entries of the form τ[(t,p1,,pn)] for every suitable tuple (t,p1,,pn), where we call (t,p1,,pn) with t[0,T+1] and pi[0,|Ji|+1] for each i[1,n] suitable if, for every i[n], it holds that pi=0 if and only if t<s(Ji). Each entry stores the minimum size of a deletion set such that the load of any symbol for any time tt does not exceed its capacity and at time t (where each letter has capacity 0 for t{0,T+1}), the pith letter of Ji is processed. Herein, pi=0 encodes that no letter of Ji is processed until time t and pi=|Ji|+1 encodes that all letters from Ji have already been processed until time t. The table is computed as follows: We initialize the table by setting

τ[(0,p1,,pn)]={0if p1==pn=0otherwise.

For the update step, we introduce some additional notation. We say that a suitable tuple (t1,p1,,pn) is compatible with a suitable tuple (t,p1,,pn) if, for every i[n],

  • pipi,

  • if pi{0,|Ji|+1}, then pi<pi and wi[pi+1]=wi[pi+2]==wi[pi1]=x, and

  • for each σΣ, we have |{isi[pi]=σ}|λσt.

Intuitively, this means that if the pith letter of job Ji is processed at time t1, then the pith letter can be processed at time t (after deleting all x in between; the second item ensures that only x are in between the pith and pith letter). We update the table as follows:

τ[(t, p1,,pn)]=
min(t1,p1,,pn) compatible with (t,p1,,pn)τ[(t1,p1,,pn)]+i:pipi(pipi1).

An optimal solution can be read from τ[(T+1,|J1|+1,,|Jn|+1)] using traceback. More formally, we have a yes-instance of WTR if and only if τ[(T+1,|J1|+1,,|Jn|+1)]k.

The dynamic programming table has O(TΠi=1n(|Ji|+2))=O(Tn+1) entries, each of which can be computed in O(nTn) time. Thus, the dynamic program runs in O(nT2n+1) time.

It remains to show correctness. First, we show by induction on t that if τ[(t,p1,,pn)]=q<, then there is a deletion family of size q such that no load exceeds its capacity at any time tt and at time t, the pith letter of Ji is scheduled. This clearly holds for t=0. For t1, let a:=(t1,p1,,pn) such that τ[(t,p1,,pn)]=τ[a]+i:pipi(pipi1). This entry τ[a] corresponds to a deletion family . By induction, for jobs {JiRi}, no load is exceeded for any time t<t. By the second condition, the pi+1,,pi1th letters of Ji are all x, say the rith to rith occurrences of x. Setting ={Ri[ri,ri]i[n]} for each i[n] results in a deletion family of size τ[a]+i:pipi(pipi1)=q. By the third condition, no load is exceeded at time t. By induction, no load is exceeded at any time tt1 for ; consequently, also for no load is exceeded at any time tt1 for . Thus, is a deletion family of size q with the desired properties.

Next, we show that if there is a deletion family  of size q such that at time t, for each i[1,n] the pith letter of Ji is scheduled, and no load is exceeded for any tt, then τ[(t,p1,,pn)]q. We again do so by induction on t. For t=0, the statement is obvious, so consider t>0. The deletion family  is also a deletion family where at time t1, the pith letter of Ji is scheduled for some p1,,pn. By the second condition, the pi+1,,pi1th letters of Ji are all x, say the rith to rith occurrences of x. Setting Ri:=Ri[ri,ri] for i[1,n] therefore results in a deletion set ={R1,,Rn} of size qi:pipi(pipi1). By induction, we have i=1n|Ri|τ[(t1,p1,,pn)]. Thus, we have

q=i=1n|Ri|+i:pipi(pipi1) τ[(t1,p1,,pn)]+i:pipi(pipi1)
τ[(t,p1,,pn)].

The correctness follows.

Adding a simple pruning rule to the dynamic program from Theorem 4 yields fixed-parameter tractability by the combined parameter n+k:

Corollary 5.

WTR parameterized by n+k can be solved in O(Tn2k+1) time.

Proof.

For every job Ji=(wi,si), the deletion set Ri has size at most k. Thus, at each time t, at most k different letters from wi can be scheduled, implying that at most k different values for pi need to be considered. Consequently, the number of DP entries to consider is O(Tnk), and each of these entries can be computed in O(nnk) time. Therefore, WTR can be solved in O(Tn2k+1) time.

4 Parameterization by the Number of Time Steps

We next study parameterization by the number T of time steps. We start by showing that WTR cannot be solved in polynomial time for constant values of T unless P=NP. Motivated by this strong hardness result, we then study parameterizations by the combined parameters of T+|Σ| and T+k.

Theorem 6.

WTR is NP-hard even if T is constant.

Proof.

We provide a reduction from Independent Set. In Independent Set, one is given a graph G=(V,E) together with an integer . The question is whether there exists a set IV of pairwise non-adjacend vertices with |I|. Independent Set is well-known to be NP-hard even if G has maximum degree 3 [7].

Let (G=(V,E),) be an instance of Independent Set, where G has maximum degree 3. We describe how to construnct an equivalent instance of WTR in polynomial time. To this end, we assume that G has a proper edge coloring with four colors: Due to Vizing’s Theorem [23], we can partition the edge set E into four sets E1, E2, E3, and E4 such that for each i[1,4], no two edges in Ei share one endpoint. Since the proof of Vizing’s Theorem provides an algorithm computing the partition in polynomial time, we may use the sets E1,,E4 for our construction.

To construct an instance (T,𝒥,c,k) of WTR, we first set k:= and we set T:=10. Furthermore, we define the alphabet Σ:={0}E. Next, given some vV and some i[1,4], we let wvi:=00 if v has no incident edge in Ei, and we let wvi:=0e otherwise, where e is the unique incident edge in Ei. We then define wv:=xwv1wv2wv3wv40 and set 𝒥:={(wv,1)vV}. Note that for every vV, we have |wv|=10.

We complete the construction by defining the capacity c. We set c(0,10):=|V|. Furthermore, for every i{2,4,6,8} and eE we define c(e,i):=1. All remaining capacities are set to .

Before we prove the correctness, we provide some intuition. Since c(0,10)=|V| and every wv ends with the symbol 0 at time step 10, we need to remove the unique occurrence of x in  strings, which then corresponds to the choice of vertices in a solution of Independent Set. The capacities c(e,i):=1 for i{2,4,6,8} then guarantee that no pair of chosen vertices are incident with the same edge.

We conclude the proof by showing that (G,) is a yes-instance of Independent Set if and only if (T,𝒥,c,k) is a yes-instance of WTR.

Let (G,) be a yes-instance of Independent Set. Then, there is a set IV containing  vertices that are pairwise non-adjacent. We construct a solution of the WTR instance as follows: for each vI, we remove the unique occurrence of x in wv and let w~v denote the resulting string. Note that w~v[t]=wv[t+1] for t[1,9]. Note that we removed |I|==k occurrences of x. We show that the modified job set satisfies the capacity constraints c: First, observe that λ010=|V||I|=V=c(0,10). Second, let i{2,4,6,8} and assume towards a contradiction that c(e,i)2 for some symbol eΣ. Then, there exist w~v and w~u with w~v[i]=w~u[i]=e. Then, we have wv[i+1]=wu[i+1]=e and thus, by the construction of 𝒥, the vertices v and u are both incident with the same edge e in G. This contradicts the fact that I consists of pairwise non-adjacent vertices. Consequently, we have c(e,i)1 for every i{2,4,6,8} and eE. Since all other capacities are , we conclude that the modified job set satisfies the capacity constraints c and thus, (T,𝒥,c,k) is a yes-instance of WTR.

Conversely, let (T,𝒥,c,k) be a yes-instance of WTR. Then, there is a subset of jobs 𝒥𝒥 with |𝒥|k= such that removing the unique occurrence of x in the corresponding strings results in a plan that satisfies the capacity constraints c. Observe that |𝒥|= since otherwise λ010=|V||𝒥|>|V|=c(0,10). Let wv1,,wv be the strings of the jobs in 𝒥. Then, I:={v1,,v} is a set of vertices of G with |I|=. It reamins to show that the vertices in I are pairwise non-adjacent. Assume towards a contradiction that two vertices v1 and v2 of I are connected by an edge e. Recall that the edge set E is partitioned into the sets E1,,E4. Without loss of generality we assume eE1. Then, by the construction of the jobs we have wv1[3]=wv2[3]=e. Removing the single occurrence of x in wv1 and wv2 provides strings w~v1 and w~v2 with w~v1[2]=w~v2[2]=e. Since c(e,2)=1, this contradicts the fact that the modified job set satisfies the capacity constraints c. Consequently, all vertices in I are pairwise non-adjacent and thus, (G,) is a yes-instance of Independent Set.

Note that the alphabet size of the instance constructed in the proof of Theorem 6 is roughly the size of the input graph of the Independent Set instance. Recall that in our application the symbols in Σ correspond to distinct machine types on which the jobs need to be processed. In a real-world production setting it is reasonable that there are only few distinct machines types. Motivated by this observation, we study the combined parameter T+|Σ| and show that WTR is FPT for this parameterization.

Theorem 7.

WTR parameterized by T+|Σ| is FPT.

Proof.

Let (T,𝒥,c,k) be an instance of WTR. We show fixed-parameter tractability by reformulating WTR as an integer linear program (ILP) where the number of variables is bounded by a function only depending on T and |Σ|. It is well-known that ILP parameterized by the number of variables is in FPT [9]. To this end, we introduce the notion of job types. Recall that a job (w,s)𝒥 consists of a string w(Σ{x}) and a starting time s[1,T]. We say that two jobs have the same type τ=(wτ,sτ) if they have the same string wτ and the same starting time sτ. We denote the number of jobs of type τ by nτ. For each job type τ, we consider possible x-removal sets R[1,#wτ]. We provide the ILP formulation by specifying the variables, the constraints, and the target function.

Variables.

For each combination of a job type τ and a possible removal set R[1,#wτ], we introduce an integer variable Y(τ,R). The intuitive idea is that the value of Y(τ,R) equals the number of jobs with type τ whose internal x-deletions correspond to the (possibly empty) removal set R. More precisely, a solution set 𝒥 contains exactly Y(τ,R) jobs of the form (wτR,sτ).

Constraints.

For each type τ, we add the following constraint, ensuring that all jobs of type τ are scheduled:

R[1,#wτ]Y(τ,R)=nτ.

For each combination of σΣ and t[1,T], we introduce the following constraint:

Type τ=(wτ,sτ)R[1,#wτ](wτR)[tsτ+1]=σY(τ,R)c(σ,t).

The left-hand side corresponds to the load of symbol σ at time step t. Thus, these constraints guarantee that the solution satisfies the capacity constraints c: Intuitively, each string wτ has been modified by a (possibly empty) x-removal set R. For a given type τ, the inner sum counts all modified strings of type τ that have symbol σ at time step t. Since we sum up these values over all possible types, the left-hand side corresponds to the load of symbol σ at time step t.

Target Function.

The target function has the form

Type τR[1,#wτ]|R|Y(τ,R).

Recall that intuitively each value of Y(τ,R) equals the number of jobs with type τ whose internal x-deletions correspond to the (possibly empty) removal set R. Thus, for each combination of a type τ and a removal set R, we delete exactly Y(τ,R)|R| symbols. Since we sum up over all τ and R, the target function models the total number of x-deletions.

Number of Variables.

Note that the number of strings in (Σ{x}) in an instance with maximum time step T is bounded by (|Σ|+1)T, and that we have at most T starting positions. Consequently, there are at most T(|Σ|+1)T possible job types τ=(wτ,sτ). Furthermore, since |wτ|T for all τ, there are at most 2T possible x-removal sets R[1,#wτ]. Thus, the number of variables Y(τ,R) is bounded by T(|Σ|+1)T2T. By a classical result by Lenstra [9], an ILP can be solved FPT-time with respect to the number of variables. This implies that the above ILP can be solved in FPT-time with respect to T+|Σ|.

We complement the FPT result for T+|Σ| by showing that WTR is unlikely to admit a problem kernel of polynomial size for this parameterization.

Theorem 8.

WTR parameterized by T+|Σ| does not admit a polynomial kernel unless NPcoNP/poly.

Proof.

We prove the statement by providing a polynomial parameter transformation from Set Cover. In Set Cover, one is given a universe of elements U, a familiy 2U of subsets of U, and an integer . The question is whether there exists a subfamily  such that || and FF=U. Set Cover parameterized by |U| does not admit a polynomial kernel unless NPcoNP/poly [6].

Let (U,,) be an instance of Set Cover and fix an arbitrary ordering U:={u1,,u|U|} of the elements in the universe. We describe how to construct an equivalent instance (T,𝒥,c,k) of WTR for the alphabet Σ={0,1}. We set k:= and we set T:=2|U|+1. Note that T is polynomial in |U|. Next, we define the job set 𝒥. Given some F and some uU, we let wuF:=01 if uF and wuF:=00 if uF. We define the job set 𝒥:={(wF,1)F}, where

wF:=xi=1|U|wuiF.

Note that  wF[2i+1]=1 if and only if uiF and we have wF[2i]=0 for every i[1,|U|].

We finish the construction by defining the (2|U|+1)-capacity bound c: For every i[1,|U|], we set c(1,2i+1):=|{FuF}|1. All remaining c(σ,t) are set to .

We next show that (U,,) is a yes-instance of Set Cover if and only if (T,𝒥,c,) is a yes instance of WTR.

Let (U,,) be a yes-instance of Set Cover. Then, there exists a subfamily  with || such that for each i[1,|U|], there is one F with uiF. We construct a solution of the WTR instance as follows: for each F, we remove the unique occurrence of x in wF, and let w~F denote the resulting string. Note that w~F[t]=wF[t+1] for every t[1,T1]. We show that the modified job set satisfies the capacity constraints c. Since all other capacities are , it suffices to show that for each i[1,|U|] the load of symbol 1 at time step 2i+1 is not larger than c(1,2i+1). Since w~F[2i+1]=wF[2i+2]=0 for all F, we have λ12i+1=|{FuiF}|. Since each ui is contained in at least one F, we conclude that λ12i+1|{FuiF}|1=c(1,2i+1). Therefore, (T,𝒥,c,) is a yes instance of WTR.

Conversely, let (T,𝒥,c,k) be a yes-instance of WTR. Then, there is a subset 𝒥𝒥 of jobs with |𝒥| such that removing the unique occurrence of x in the corresponding strings results in a plan that satisfies the capacity constraints c. Let wF1,,wF|𝒥| be the strings of the jobs in 𝒥. Then, :={F1,,F|𝒥|} is a subfamily of  of size at most . Assume towards a contradiction that there is some uiU with uiFF. Then, for each F with uiF, the occurrence of x in wF has not been removed. Consequently, the load of symbol 1 at time step 2i+1 is |{FuiF}|>c(1,2i+1). This contradicts the fact that removing the x in the strings of 𝒥 results in a job set that satisfies the capacity constraints c. Therefore, (U,,) is a yes-instance of Set Cover.

The following two observations can be made by a taking closer look at the instance constructed in the proof of Theorem 8: First, the parameter k equals the solution size of the Set Cover instance. Since Set Cover is known to be W[2]-hard when parameterized by solution size [5], this implies the following.

Corollary 9.

WTR parameterized by k is W[2]-hard even if |Σ|=2.

Second, the parameter n+k is upper-bounded by 2||, where  is the set familiy in the Set Cover instance. Since Set Cover does not admit a polynomial kernel for || [6], this implies the following.

Corollary 10.

WTR does not admit a polynomial kernel for n+k unless NPcoNP/poly even if |Σ|=2.

Next, consider the parameterization by T+k. Recall that WTR is XP for k by Observation 2 and thus, WTR is XP for T+k. We complement this result by showing W[1]-hardness.

Theorem 11.

WTR parameterized by T+k is W[1]-hard.

Proof.

We reduce from Multicolored Clique parameterized by solution size . In Multicolored Clique, one is given an -partite graph G=(V1V,E) with |V1|=|V2|==|V|, and the task is to find a multicolored clique, i.e., a set of  pairwise adjacent vertices containing exactly one vertex from each Vi. Multicolored Clique is W[1]-hard parameterized by , even if there is some r so that each vertex has exactly r neighbors in every other color class [5]. Let G=(V1V,E) be an instance of Multicolored Clique where each vertex v has exactly r neighbors in every other color class. We let f:[1,(2)]([1,]2) be an arbitrary bijection to fix an ordering of the set ([1,]2):={X[1,]|X|=2} containing all two-element subsets of [1,]. We use the alphabet Σ=V{0,1}{zii[]}. We add the following jobs:

  • for each vVi, we add one job Jv=(wv,1) where wv starts with xzi. Afterwards, for each q[1,(2)], we append the following four letters:

    • if if(q), then we append 0000,

    • if f(q)={i,i} for some i>i, then we append 0v00, and

    • if f(q)={i,i} for some i<i, then we append 000v.

  • for each edge e={u,v}E with uVi and vVi for some i<i, we add one job Je=(xu0v,4(f({i,i})1)+2).

The capacities are as follows:

  • for each i[1,], letter zi has capacity 1 at time 1 and capacity |V1|1 at time 2,

  • for each vV, letter v has capacity r at every time step, and

  • all other capacities are set to be .

Finally, we set k:=+(2). Note that T=4(2)+2. As the reduction can clearly be computed in polynomial time, it remains to show correctness.

Before we show the correctness, we provide some intuition. The capacities of zi at time 1 ensures that at most one deletion occurs in jobs {JvvVi}, while the capacities of zi at time 2 ensure that at least one deletion occurs in {JvvVi}. Therefore, exactly one deletion occurs in {JvvVi} for every i[1,]. These vertices with a deletion correspond to a multicolored clique: Since there exists some unique viVi such that the x in Jvi is deleted, this implies that all occurrences of vi in Jvi move one time step earlier. Thereby, for each q[1,(2)] with if(q), letter vi from Jvi “collides” with letter vi in Je the r edges e incident to vi in E[Vi,Vi] where i=f(q){i}. Because the capacity of vi is r at each time step, this implies that a deletion has to occur in Je for some edge eE[Vi,Vi] incident to vi. As the budget allows for k=(2)further deletions, this implies that there is exactly one edge in E[Vi,Vi] with a deletion for {i,i}([1,]2). Since this edge needs to be incident to both vi and vi, it follows that {vii[1,]} is a multicolored clique.

We first show that if G is a yes-instance to Multicolored Clique, then (T,𝒥,c,k) is yes-instance of WTR. Given a clique {v1,,v} with viVi, we construct a solution to the WTR instance as follows: We remove the unique occurrence of x in job Jvi for each i[1,] and in job J{vi1,vi2} for i1<i2[1,]. We call the resulting strings w~v for vV and w~e for eE. As we removed k occcurences of x in total, it remains to show that the capacity constraints are satisfied. For letter zi, we have a load of 1 at time 1 (from job Jvi) and a load of |V1|1 at time to (from jobs Jv for VVi{vi}). Thus, the capacities of zi are satisfied. For vV{v1,,v}, the maximum load at any time is deg(v)=r, so its capacity is satisfied. We continue with letter vi for i[1,]. At time 4(q1)+3, the load of vi is either 0 (if if(q) or f(q)={i,i} for some i<i) or r (if f(q)={i,i} for some i>i; in this case, job Jvi and r1 edge jobs (all edges from E[Vi,Vi]δ(vi) except for {vi,vi}). In both cases, the capacity constraint is satisfied. Symetrically, at time 4(q1)+5, the load is either 0 (if if(q) or f(q)={i,i} for some i>i) or r (if f(q)={i,i} for some i<i; in this case, job Jvi and r1 edge jobs (all edges from E[Vi,Vi]δ(vi) except for {vi,vi}). In both cases, the capacity constraint is satisfied.

Lastly, we show that if (T,𝒥,c,k) is a yes-instance to WTR, then G is a yes-instance to Multicolored Clique. Since (T,𝒥,c,k) is a yes-instance, there is a subset 𝒥𝒥 of jobs with |𝒥| such that removing the unique occurrence of x in the corresponding strings results in a plan that satisfies the capacity constraints c. As the capacity of zi for i[] at time step 2 is |V1|1, it follows that for each i[1,], there exists some viVi with Jvi𝒥. As c(zi,1)=1, this vi is unique.

We next show that {v1,,v} forms a multicolored clique. Consider any q[1,(2)] and let f(q)={i1,i2} for some i1<i2. At time 4(q1)+3, job Jvi1 contributes one to the load of vi1. For each edge eE[Vi1,Vi2]δ(vi1), job Je contributes 1 to the load of vi1 at time 4(q1)+3 if and only if Je𝒥. As |E[Vi1,Vi2]δ(vi1)|=r=c(vi,4(q1)+3), this implies that there is at least one edge eE[Vi1,Vi2]δ(vi1) with Je𝒥. Considering the letter vi2 and time 4(q1)+5, symmetric arguments show that there is at least one edge eE[Vi1,Vi2]δ(vi2) with Je𝒥. As there are deletions in jobs Jv for vV and k=(2), this implies that for each i1<i2[], there is a unique eE[Vi,Vi]δ(vi1)δ(vi2) with Je𝒥, and this edge is incident to both vi1 and vi2, i.e., e={vi1,vi2}. This implies that {v1,,v} is a multicolored clique.

5 Conclusion

In this work, we initiated theoretical studies of repairing production schedules by the deletion of waiting times in a job shop setting by introducing a mathematical formulation of the problem. While our results mostly highlight the hardness of this problem, there are also few tractability results, for example when there are only few different machine types and a short time horizon is considered. We conclude with several proposals for further research questions:

  • Instead of removing waiting times, there are also other natural operations to modify a schedule, for example swapping jobs. What is the complexity of the resulting problem?

  • Our scheduling model is quite general, being based on job shop scheduling. Does basing the problem on simpler scheduling models, e.g. flow shops, result in a more tractable problem?

  • Given that eliminating all capacity violations turned out to be quite hard, one might try to not directly search for a “perfect” schedule, but instead only a “better” one (i.e., one with fewer capacity violations). Can such a schedule be found using the framework of parameterized local search?

  • While we study how to remove pre-planned waiting times, it might be interesting to study how these waiting times can be established in the first place. How can we create schedules that are robust in the sense that they can be repaired by removing waiting times?

References

  • [1] Alessandro Agnetis, Jean-Charles Billaut, Michael Pinedo, and Dvir Shabtay. Fifty years of research in scheduling – Theory and applications. European Journal of Operational Research, 2025. doi:10.1016/j.ejor.2025.01.034.
  • [2] Susanne Albers and Günter Schmidt. Scheduling with unexpected machine breakdowns. Discret. Appl. Math., 110(2-3):85–99, 2001. doi:10.1016/S0166-218X(00)00266-3.
  • [3] Kaja Balzereit, Niels Grüttemeier, Nils Morawietz, Dennis Reinhardt, Stefan Windmann, and Petra Wolf. Scalable neighborhood local search for single-machine scheduling with family setup times. CoRR, abs/2409.00771, 2024. doi:10.48550/arXiv.2409.00771.
  • [4] Nicolas Boria and Federico Della Croce. Reoptimization in machine scheduling. Theor. Comput. Sci., 540:13–26, 2014. doi:10.1016/J.TCS.2014.04.004.
  • [5] Marek Cygan, Fedor V. Fomin, Lukasz Kowalik, Daniel Lokshtanov, Dániel Marx, Marcin Pilipczuk, Michal Pilipczuk, and Saket Saurabh. Parameterized Algorithms. Springer, 2015. doi:10.1007/978-3-319-21275-3.
  • [6] Michael Dom, Daniel Lokshtanov, and Saket Saurabh. Kernelization lower bounds through colors and ids. ACM Trans. Algorithms, 11(2):13:1–13:20, 2014. doi:10.1145/2650261.
  • [7] M. R. Garey, David S. Johnson, and Larry J. Stockmeyer. Some simplified NP-complete graph problems. Theor. Comput. Sci., 1(3):237–267, 1976. doi:10.1016/0304-3975(76)90059-1.
  • [8] Kevin D Glazebrook. Scheduling tasks with exponential service times on parallel processors. J. Appl. Probab., 16(3):685–689, 1979.
  • [9] Hendrik W. Lenstra Jr. Integer programming with a fixed number of variables. Math. Oper. Res., 8(4):538–548, 1983. doi:10.1287/MOOR.8.4.538.
  • [10] Bala Kalyanasundaram and Kirk Pruhs. Fault-tolerant scheduling. SIAM J. Comput., 34(3):697–719, 2005. doi:10.1137/S0097539794261799.
  • [11] Richard M. Karp. Reducibility among combinatorial problems. Complexity of Computer Problems, pages 85–103, 1972. doi:10.1007/978-1-4684-2001-2_9.
  • [12] Bernhard Korte and Jens Vygen. Combinatorial Optimization. Springer, 6th edition, 2018.
  • [13] Patricio Lamas, Marcos Goycoolea, Bernardo K. Pagnoncelli, and Alexandra M. Newman. A target-time-windows technique for project scheduling under uncertainty. Eur. J. Oper. Res., 314(2):792–806, 2024. doi:10.1016/J.EJOR.2023.10.027.
  • [14] Arthur Müller and Lukas Vollenkemper. Reinforcement learning as an improvement heuristic for real-world production scheduling. In Proceedings of the 23rd IEEE International Conference on Machine Learning and Applications (ICMLA ’21). IEEE, 2024. doi:10.1109/ICMLA52953.2021.00085.
  • [15] Krzysztof Pietrzak. On the parameterized complexity of the fixed alphabet shortest common supersequence and longest common subsequence problems. J. Comput. Syst. Sci., 67(4):757–771, 2003. doi:10.1016/S0022-0000(03)00078-3.
  • [16] George A. Rovithakis, Stelios E. Perrakis, and Manolis A. Christodoulou. Application of a neural-network scheduler on a real manufacturing system. IEEE Trans. Control. Syst. Technol., 9(2):261–270, 2001. doi:10.1109/87.911378.
  • [17] Markus W. Schäffter. Scheduling with forbidden sets. Discret. Appl. Math., 72(1-2):155–166, 1997. doi:10.1016/S0166-218X(96)00042-X.
  • [18] Dvir Shabtay and Miri Gilenson. A state-of-the-art survey on multi-scenario scheduling. Eur. J. Oper. Res., 310(1):3–23, 2023. doi:10.1016/J.EJOR.2022.11.014.
  • [19] Martin Skutella, Maxim Sviridenko, and Marc Uetz. Unrelated machine scheduling with stochastic processing times. Math. Oper. Res., 41(3):851–864, 2016. doi:10.1287/MOOR.2015.0757.
  • [20] Stephen F Smith. Reactive scheduling systems. Intelligent scheduling systems, pages 155–192, 1995.
  • [21] V Subramaniam, AS Raheja, and K Rama Bhupal Reddy. Reactive repair tool for job shop schedules. Int. J. Prod. Res., 43(1):1–23, 2005.
  • [22] E. Szelke and G. Markus. A blackboard based perspective of reactive scheduling. In Artificial Intelligence in Reactive Scheduling, pages 60–77. Springer US, 1995. doi:10.1007/978-0-387-34928-2_6.
  • [23] Vadim G Vizing. On an estimate of the chromatic class of a p-graph. Diskret. analiz., 3:25–30, 1964.