Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Collections

Scheduling Theory

Machine, flowshop, jobshop and project scheduling: optimality of classic rules, complexity reductions, and approximation guarantees.

70 missions

Missions

21–40 of 70
OpenCompletedAll
🏆Completed
Algorithmic Game TheoryMechanism DesignOperations Research+1·Captain: mikedeng1

Algorithmic Mechanism Design I: MinWork Is a Strongly Truthful n-Approximation Mechanism for Task Scheduling on Unrelated MachinesResearch Paper

Motivation

Algorithmic mechanism design asks for algorithms whose inputs are held by self-interested parties. Each party reports its private data, the algorithm computes an outcome, and payments are arranged so that no party gains by misreporting. Nisan and Ronen introduced the field in Algorithmic Mechanism Design (Games Econ. Behav. 35, 2001). Their running example is task scheduling on unrelated machines: kkk tasks are distributed among nnn machines owned by different agents, each agent knows only its own processing times, and the designer wants to minimize the make-span.

Without incentives the problem is classical: minimizing make-span on unrelated machines is NP-hard and admits a polynomial 2-approximation (Lenstra, Shmoys, Tardos, 1990). With selfish agents the question changes: which approximation ratios can a truthful mechanism guarantee? This mission formalizes the paper's upper bound, the MinWork mechanism, which is the benchmark every later lower bound for truthful scheduling is compared with.

Timeline.

  • 1961: Vickrey introduces the second-price auction (J. Finance 16).
  • 1971–1973: Clarke and Groves generalize it to the VCG family of truthful mechanisms for utilitarian objectives (Groves, Econometrica 41, 1973).
  • 1999/2001: Nisan and Ronen show MinWork is a strongly truthful nnn-approximation, and that no truthful mechanism beats ratio 2.
  • 2007: Christodoulou, Koutsoupias and Vidali raise the deterministic lower bound to 1+21+\sqrt21+2​ for n≥3n \ge 3n≥3; Koutsoupias and Vidali later raise it to 1+φ≈2.6181+\varphi \approx 2.6181+φ≈2.618.
  • 2023: Christodoulou, Koutsoupias and Kovács prove the Nisan–Ronen conjecture: no deterministic truthful mechanism achieves a ratio below nnn (STOC 2023, arXiv:2301.11905), so MinWork is optimal among deterministic truthful mechanisms.

Setting

There are nnn agents and kkk tasks. Agent iii's type is the vector ti=(t1i,…,tki)t^i = (t^i_1,\dots,t^i_k)ti=(t1i​,…,tki​) of positive times, tji>0t^i_j > 0tji​>0 being the time agent iii needs to perform task jjj. A type vector is t=(t1,…,tn)t = (t^1,\dots,t^n)t=(t1,…,tn). An allocation xxx sends each task jjj to one agent; xix^ixi is the set of tasks agent iii receives. The make-span of xxx is

g(x,t)=max⁡i∑j∈xitji,g(x,t) = \max_{i} \sum_{j \in x^i} t^i_j ,g(x,t)=imax​j∈xi∑​tji​,

and agent iii's valuation is vi(x,ti)=−∑j∈xitjiv^i(x,t^i) = -\sum_{j \in x^i} t^i_jvi(x,ti)=−∑j∈xi​tji​.

A direct mechanism asks every agent to declare a type, computes an allocation x(d)x(d)x(d) from the declared vector ddd, and hands agent iii a payment pi(d)p^i(d)pi(d). Agent iii's utility is pi(d)+vi(x(d),ti)p^i(d) + v^i(x(d), t^i)pi(d)+vi(x(d),ti), with tit^iti its true type. The mechanism is truthful if declaring tit^iti maximizes agent iii's utility for every declaration of the others, and strongly truthful if truth-telling is the only such dominant strategy. An allocation rule is a ccc-approximation if g(x(t),t)≤c⋅g(y,t)g(x(t),t) \le c \cdot g(y,t)g(x(t),t)≤c⋅g(y,t) for every type vector ttt and every allocation yyy.

The MinWork mechanism allocates each task to an agent with minimal declared time for it, breaking ties arbitrarily. For each task it wins, an agent receives the second-best declared time min⁡i′≠idji′\min_{i' \ne i} d^{i'}_jmini′=i​dji′​:

pi(d)=∑j∈xi(d)min⁡i′≠idji′.p^i(d) = \sum_{j \in x^i(d)} \min_{i' \neq i} d^{i'}_j .pi(d)=j∈xi(d)∑​i′=imin​dji′​.

The Lean development uses the same names: load, makespan, IsTruthful, IsStronglyTruthful, IsApprox, IsMinWorkAlloc, secondBest, minTime, minWorkPay.

Formalization targets

Goal: Theorem 4.1

For n≥2n \ge 2n≥2 and every MinWork allocation rule xxx with payments ppp as above,

(x,p) is strongly truthfulandg(x(t),t)≤n⋅g(y,t)  for all positive t and all allocations y.(x,p)\ \text{is strongly truthful} \quad\text{and}\quad g(x(t),t) \le n \cdot g(y,t)\ \ \text{for all positive } t \text{ and all allocations } y .(x,p) is strongly truthfulandg(x(t),t)≤n⋅g(y,t)  for all positive t and all allocations y.

Milestones

  1. Theorem 3.1 (Groves): a VGC mechanism is truthful. This is an existing platform theorem, used as a reference.
  2. MinWork belongs to the VGC family. Its allocation maximizes ∑ivi(ti,x)\sum_i v^i(t^i,x)∑i​vi(ti,x), and its payment is ∑i′≠ivi′(ti′,x(t))+h−i\sum_{i'\ne i} v^{i'}(t^{i'},x(t)) + h^{-i}∑i′=i​vi′(ti′,x(t))+h−i with h−i=∑jmin⁡i′≠itji′h^{-i} = \sum_j \min_{i'\ne i} t^{i'}_jh−i=∑j​mini′=i​tji′​.
  3. Claim 4.2: MinWork is strongly truthful.
  4. g(x(t),t)≤∑jmin⁡itjig(x(t),t) \le \sum_{j} \min_i t^i_jg(x(t),t)≤∑j​mini​tji​.
  5. g(y,t)≥1n∑jmin⁡itjig(y,t) \ge \frac1n \sum_j \min_i t^i_jg(y,t)≥n1​∑j​mini​tji​ for every allocation yyy.
  6. Claim 4.3: MinWork is an nnn-approximation.

Significance

The theorem gives the first positive result for truthful scheduling: a mechanism that is truthful in the strongest sense and is within a factor nnn of optimal, whatever the tie-breaking rule. Every lower bound in the paper (Theorems 4.6, 4.10 and 4.12) and in the later literature measures itself against this ratio. Since the 2023 resolution of the Nisan–Ronen conjecture, the ratio nnn is known to be tight for deterministic truthful mechanisms.

The result is proved in the paper; it is not known to be formalized in any proof assistant. The platform already has Groves' theorem in an abstract form (AGT.vcg_incentive_compatible). This mission connects that abstract statement to a concrete combinatorial mechanism, and it adds the strict part of strong truthfulness for any number of tasks and agents, which the paper proves only for one task and two agents. The vocabulary (make-span over unrelated machines, direct scheduling mechanisms, strong truthfulness) is shared with the seven later missions of this series.

Difficulty

Truthfulness follows from Groves' theorem once MinWork is identified as a VGC mechanism. The identification requires the payment identity at every declared vector and under every tie-breaking rule, including ties at the winning time. The main difficulty is the strict part of strong truthfulness. A misreport that differs from the truth only on one task must still be shown to lose strictly for some declarations of the others. Those declarations must stay positive, and on every other task they must leave the outcome unchanged. The paper's proof covers only one task and two agents and leaves the general case as "similar". Its printed inequality also has the two utilities in the wrong order (see below), so it cannot be transcribed directly.

Formalization scope

  • Agents are Fin n and tasks are Fin k. An allocation is a function Fin k → Fin n, and an agent may receive no task. Types are positive reals, and every truthfulness and approximation quantifier ranges over positive true types, positive misreports and positive declarations of the others.
  • Payments are handed to the agent, so utility is the payment minus the true time spent. Payments are computed from the declared vector, never from true types.
  • The allocation rule is a parameter satisfying the MinWork specification (IsMinWorkAlloc). Every result holds for every tie-breaking rule, including rules that depend on the whole declared vector. No particular argmin is fixed.
  • n≥2n \ge 2n≥2 is a hypothesis of the goal and of the truthfulness items: with a single agent the paper's second-best minimum is undefined. The approximation items need only n≥1n \ge 1n≥1. There is no hypothesis on kkk.
  • The make-span and both minima are Finset.sup' / Finset.inf' over nonempty finite sets, so they are true maxima and minima with no default values.
  • Strong truthfulness is formalized as truthfulness plus: every misreport di≠tid^i \ne t^idi=ti is strictly worse than the truth for some positive declarations of the others. Given truthfulness this is equivalent to Definition 5. A formalization that states only that truth-telling is dominant, or proves strictness only for single-task instances, does not meet the goal. Neither does an existential ratio in place of nnn.
  • Printed slip: in the proof of Claim 4.2 (p. 177) the case di>tid^i > t^idi>ti reads "the utility for agent iii is ti−di<0t^i - d^i < 0ti−di<0, instead of 0 in the case of truth-telling". With the Definition 11 payments the misreporting agent loses the task (utility 0), and the truthful agent wins it with utility d3−i−ti>0d^{3-i} - t^i > 0d3−i−ti>0. The milestone text keeps the paper's words; the Lean statements assert what the argument establishes.
  • Out of scope: running time ("polynomial time"), and the paper's general revelation-principle framework (Proposition 2.1).
  • Welcome contributions: proofs of the milestones, and a reusable lemma connecting the local VGC milestone to AGT.vcg_incentive_compatible.

Selected references

  • N. Nisan, A. Ronen, Algorithmic Mechanism Design, Games and Economic Behavior 35 (2001) 166–196. https://doi.org/10.1006/game.1999.0790
  • T. Groves, Incentives in Teams, Econometrica 41 (1973) 617–631. https://doi.org/10.2307/1914085
  • W. Vickrey, Counterspeculation, Auctions, and Competitive Sealed Tenders, Journal of Finance 16 (1961) 8–37. https://doi.org/10.1111/j.1540-6261.1961.tb02789.x
  • J. K. Lenstra, D. B. Shmoys, É. Tardos, Approximation algorithms for scheduling unrelated parallel machines, Mathematical Programming 46 (1990) 259–271. https://doi.org/10.1007/BF01585745
  • G. Christodoulou, E. Koutsoupias, A. Kovács, A Proof of the Nisan-Ronen Conjecture, STOC 2023. https://arxiv.org/abs/2301.11905
9 thms4 active usersReviewed
🏆Completed
Algorithmic Game TheoryMechanism DesignOperations Research+1·Captain: mikedeng1

Algorithmic Mechanism Design II: A Lower Bound for Truthful Task SchedulingResearch Paper

Motivation

Algorithms deployed on the Internet often take their inputs from parties who own them and who may lie when lying pays. Nisan and Ronen's Algorithmic Mechanism Design (Games and Economic Behavior 35, 2001) proposed studying optimization problems in this setting: the algorithm designer may hand out payments, and must guarantee that the intended output is produced when every participant acts in its own interest. The paper's central test case is scheduling on unrelated machines, a standard problem of combinatorial optimization, in which the machines are the selfish participants and only they know how long each job takes them.

For this problem the paper shows that incentives cost a factor of two at least: with two or more machines, no mechanism can guarantee a make-span below twice the optimum. This was the first lower bound separating what incentive-compatible mechanisms can achieve from what ordinary approximation algorithms can achieve, and it started a line of work on the "Nisan–Ronen conjecture" (that the right factor for nnn machines is nnn), with improved lower bounds by Christodoulou, Koutsoupias and Vidali (Algorithmica, 2009) and by Koutsoupias and Vidali (Algorithmica, 2013), and a resolution announced by Christodoulou, Koutsoupias and Kovács (STOC 2023).

Setting

There are nnn agents (machines) i=1,…,ni = 1,\dots,ni=1,…,n and kkk tasks j=1,…,kj = 1,\dots,kj=1,…,k. Agent iii's private type is the vector ti=(t1i,…,tki)t^i = (t^i_1,\dots,t^i_k)ti=(t1i​,…,tki​) of positive real numbers, tjit^i_jtji​ being the time agent iii needs for task jjj; a type vector is t=(t1,…,tn)t = (t^1,\dots,t^n)t=(t1,…,tn). An allocation xxx assigns every task to one agent; xix^ixi is the set of tasks given to agent iii. For a set XXX of tasks write ti(X)=∑j∈Xtjit^i(X) = \sum_{j\in X} t^i_jti(X)=∑j∈X​tji​. The objective is the make-span

g(x,t)=max⁡iti(xi),g(x,t) = \max_{i} t^i(x^i),g(x,t)=imax​ti(xi),

and an allocation rule is a ccc-approximation if its make-span is at most ccc times that of every allocation, on every type vector.

A mechanism m=(o,p)m = (o,p)m=(o,p) gives each agent iii a set AiA^iAi of strategies. On a strategy profile a=(a1,…,an)a = (a^1,\dots,a^n)a=(a1,…,an) it outputs an allocation o(a)o(a)o(a) and hands agent iii a payment pi(a)p^i(a)pi(a). An agent of type tit^iti has utility pi(a)−ti(oi(a))p^i(a) - t^i(o^i(a))pi(a)−ti(oi(a)). A strategy is dominant if it maximizes the agent's utility whatever the others play. The mechanism implements a ccc-approximation if every agent of every type has a dominant strategy and every profile of dominant strategies yields a ccc-approximate allocation.

A direct mechanism (x,p)(x,p)(x,p) has AiA^iAi equal to the set of types, and is truthful if reporting the true type is dominant. For a truthful mechanism, the price pi(X,t−i)p^i(X,t^{-i})pi(X,t−i) is the payment agent iii receives when, against the others' reports t−it^{-i}t−i, some report of its own makes it receive exactly XXX (and 000 if none does); the price difference is Δi(A,B)=pi(A∪B,t−i)−pi(A,t−i)\Delta^i(A,B) = p^i(A\cup B,t^{-i}) - p^i(A,t^{-i})Δi(A,B)=pi(A∪B,t−i)−pi(A,t−i).

Formalization targets

Goal: Theorem 4.6

For every n≥2n\ge 2n≥2, k≥3k\ge3k≥3 and c<2c<2c<2, no mechanism with any strategy sets implements a ccc-approximation:

∀ (A,o,p):¬ Implements(o,p,c).\forall\, (A, o, p):\quad \neg\ \mathrm{Implements}(o,p,c).∀(A,o,p):¬ Implements(o,p,c).

Milestones

  1. Proposition 2.1 (revelation principle): a mechanism implementing a ccc-approximation yields a truthful direct mechanism whose allocation rule is a ccc-approximation.
  2. Theorem 4.6 for truthful mechanisms (§4.3): no truthful direct mechanism has a ccc-approximate allocation rule for c<2c<2c<2. With milestone 1 it gives the goal.
  3. Proposition 4.4 (independence): for a truthful mechanism, t1−i=t2−it_1^{-i}=t_2^{-i}t1−i​=t2−i​ and xi(t1)=xi(t2)x^i(t_1)=x^i(t_2)xi(t1​)=xi(t2​) imply pi(t1)=pi(t2)p^i(t_1)=p^i(t_2)pi(t1​)=pi(t2​).
  4. Proposition 4.5 (maximization): xi(t)x^i(t)xi(t) maximizes pi(X,t−i)−ti(X)p^i(X,t^{-i}) - t^i(X)pi(X,t−i)−ti(X) over attainable XXX.
  5. Lemma 4.7: the price-difference inequalities satisfied by xi(t)x^i(t)xi(t), and the uniqueness statement for sets satisfying them strictly.
  6. Claim 4.8: for two agents, all-ones types and 0<ε<10<\varepsilon<10<ε<1, moving agent 1's times to ε\varepsilonε on its own bundle and 1+ε1+\varepsilon1+ε elsewhere leaves the allocation unchanged.
  7. The even case of the ratio: at that perturbed instance the mechanism's make-span is ∣x2(t)∣|x^2(t)|∣x2(t)∣ while some allocation achieves 12∣x2(t)∣+kε\tfrac12|x^2(t)| + k\varepsilon21​∣x2(t)∣+kε.

Significance

The result. Theorem 4.6 shows that the requirement of dominant-strategy incentive compatibility, by itself, rules out approximation ratios below 222 for scheduling on unrelated machines, a problem for which polynomial-time 222-approximation algorithms that ignore incentives exist (Lenstra, Shmoys, Tardos 1990) and for which the exact optimum is computable in exponential time. Combined with the MinWork mechanism of the same paper (an nnn-approximation), it determines the optimal ratio for two machines. It is the base case of the Nisan–Ronen conjecture and the prototype of the "characterize truthful mechanisms by prices" technique used throughout later work on the conjecture.

Formalizing it. The theorem has been proved since 1999, but no machine-checked version is known to exist. The mission produces a formal account of general mechanisms with arbitrary strategy sets, dominant-strategy implementation, the revelation principle in that generality, and the price characterization of truthful mechanisms (independence and maximization). These are reusable for every other lower bound in this paper and for the later literature on the conjecture.

Difficulty

The statement quantifies over all mechanisms, with arbitrary strategy sets and arbitrary payment functions, so no finite search settles it. The revelation principle reduces to truthful direct mechanisms, but even these are an infinite-dimensional family: the allocation rule may break ties in any way, and prices may be any functions of the other agents' reports.

The printed argument also has two places that need care. Proposition 4.5 and Lemma 4.7, as printed, range over all sets of tasks, while Definition 12 gives unattainable sets price 000; the statements hold only over attainable sets, and are formalized that way. And the case where agent 2's bundle has odd size is dispatched in one sentence ("which still yields the same allocation"), which the preceding lemma does not justify when agent 2's best bundle at the perturbed prices is not unique. A complete formal proof of the goal must supply an argument for that case.

Formalization scope

  • Agents are Fin n, tasks Fin k; an allocation is a function Fin k → Fin n; bundles may be empty. The make-span is a finite maximum and assumes n≥1n\ge1n≥1 (NeZero n).
  • Types, declarations and misreports are strictly positive reals throughout (Definition 10). Utility is quasi-linear; payments are handed to the agent and may have either sign.
  • A general mechanism has strategy sets A : Fin n → Type u (any universe), output ooo and payments ppp on dependent strategy profiles. Implements requires both that every agent of every positive type has a dominant strategy and that every profile of dominant strategies yields a ccc-approximate allocation. Dominance is against every profile of the others, not only dominant ones. Without the existence clause, a mechanism with no dominant strategies would implement vacuously; the definition excludes that.
  • Thresholds made explicit: n≥2n\ge2n≥2 and k≥3k\ge3k≥3, both taken from the proof ("We prove the theorem for the case of two agents"; "Let k≥3k\ge3k≥3"). The goal holds for each fixed nnn and kkk and every c<2c<2c<2, for every mechanism, with no restriction on tie-breaking and no requirement of strong truthfulness. At n=1n=1n=1 the claim is false.
  • Proposition 2.1 is stated for task scheduling with the ccc-approximation specification; "truthful implementation" is read as truth-telling dominant and the truthful output ccc-approximate.
  • Printed slips: Proposition 4.5 and Lemma 4.7 are stated over attainable sets; the "Moreover" of Lemma 4.7 requires YYY attainable. The odd case of the ratio step is not a milestone.
  • The reduction from n>2n>2n>2 to two agents ("having the other agents be much slower") is not a separate milestone; the goal covers every n≥2n\ge2n≥2.
  • Running time ("polynomial-time computable") is out of scope and not modelled.

Contributions of any of the milestones are welcome, as are alternative proofs of the goal that avoid the terse odd case.

Selected references

  • N. Nisan, A. Ronen, Algorithmic Mechanism Design, Games and Economic Behavior 35 (2001) 166–196. https://doi.org/10.1006/game.1999.0790
  • A. Mas-Colell, M. D. Whinston, J. R. Green, Microeconomic Theory, Oxford University Press, 1995 (revelation principle, p. 871).
  • J. K. Lenstra, D. B. Shmoys, É. Tardos, Approximation algorithms for scheduling unrelated parallel machines, Mathematical Programming 46 (1990) 259–271. https://doi.org/10.1007/BF01585745
  • G. Christodoulou, E. Koutsoupias, A. Vidali, A lower bound for scheduling mechanisms, Algorithmica 55 (2009).
  • E. Koutsoupias, A. Vidali, A lower bound of 1+φ for truthful scheduling mechanisms, Algorithmica 66 (2013).
  • G. Christodoulou, E. Koutsoupias, A. Kovács, A proof of the Nisan–Ronen conjecture, STOC 2023.
11 thms3 active usersReviewed
🏆Completed
Algorithmic Game TheoryMechanism DesignOperations Research+1·Captain: mikedeng1

Algorithmic Mechanism Design III: No Additive Truthful Mechanism Achieves a c-Approximation for Task Scheduling for Any c < nResearch Paper

Motivation

Algorithmic mechanism design studies optimization problems whose inputs are held by self-interested agents: an algorithm must not only compute a good solution but also pay the agents so that reporting their data truthfully is in their own interest. Nisan and Ronen introduced the field in Algorithmic Mechanism Design (Games Econ. Behav. 35, 2001) with task scheduling on unrelated machines as the model problem. Each machine is owned by an agent who alone knows how long it takes for each task; the designer wants a schedule of small make-span.

The paper gives a truthful mechanism, MinWork, whose make-span is within a factor nnn of optimal, and a lower bound of 222 for every truthful mechanism. It conjectures that no truthful mechanism beats nnn (Conjecture 4.9) and proves the conjecture for two natural classes. This mission is about one of them, additive mechanisms (Theorem 4.10, p. 180).

Timeline of the gap between 222 and nnn:

  • 1999/2001: Nisan and Ronen prove the lower bound 222 for all truthful mechanisms and nnn for additive and for local mechanisms.
  • 2007: Christodoulou, Koutsoupias and Vidali raise the general lower bound to 1+21+\sqrt 21+2​ for n≥3n\ge 3n≥3 (SODA 2007; Algorithmica 2009).
  • 2008: Christodoulou, Koutsoupias and Vidali characterize the truthful mechanisms for two machines (ESA 2008; arXiv 0807.3427); in parallel, Dobzinski and Sundararajan (EC 2008) characterize them and show that for two machines no truthful mechanism beats 222.
  • 2023: Christodoulou, Koutsoupias and Kovács prove the Nisan–Ronen conjecture: no truthful mechanism beats nnn (STOC 2023; arXiv 2301.11905).

Setting

There are nnn agents i=1,…,ni=1,\dots,ni=1,…,n and kkk tasks j=1,…,kj=1,\dots,kj=1,…,k. The type of agent iii is a vector ti=(t1i,…,tki)t^i=(t^i_1,\dots,t^i_k)ti=(t1i​,…,tki​) of positive reals, tjit^i_jtji​ being the time agent iii needs for task jjj; a type vector t=(t1,…,tn)t=(t^1,\dots,t^n)t=(t1,…,tn) collects all types, and t−it^{-i}t−i denotes the types of the agents other than iii. An allocation x=(x1,…,xn)x=(x^1,\dots,x^n)x=(x1,…,xn) is a partition of the tasks, xix^ixi being the set given to agent iii. For a set XXX of tasks write ti(X)=∑j∈Xtjit^i(X)=\sum_{j\in X}t^i_jti(X)=∑j∈X​tji​. The make-span is

g(x,t)=max⁡iti(xi).g(x,t)=\max_i t^i(x^i).g(x,t)=imax​ti(xi).

A direct mechanism m=(x,p)m=(x,p)m=(x,p) maps every declared type vector ttt to an allocation x(t)x(t)x(t) and to payments pi(t)p^i(t)pi(t) handed to the agents. An agent with true type tit^iti gets utility pi(t)−ti(xi(t))p^i(t)-t^i(x^i(t))pi(t)−ti(xi(t)). The mechanism is truthful if, whatever the others declare, no agent gains by declaring a type other than its true one. It is a ccc-approximation if g(x(t),t)≤c g(y,t)g(x(t),t)\le c\,g(y,t)g(x(t),t)≤cg(y,t) for every positive type vector ttt and every allocation yyy.

The price offered for a set XXX to agent iii (Definition 12) is the payment pi(t′i,t−i)p^i(t'^i,t^{-i})pi(t′i,t−i) at any declaration t′it'^it′i for which the mechanism gives agent iii exactly XXX, and 000 if there is no such declaration. For truthful mechanisms this is well defined (Proposition 4.4, Independence). The mechanism is additive (Definition 13) if

pi(X,t−i)=∑j∈Xpi({j},t−i)p^i(X,t^{-i})=\sum_{j\in X}p^i(\{j\},t^{-i})pi(X,t−i)=j∈X∑​pi({j},t−i)

for every agent iii, type vector ttt and set XXX of tasks. MinWork, which gives each task to the fastest agent and pays it the second-fastest time, is additive.

Formalization targets

Goal: Theorem 4.10

For n≥1n\ge 1n≥1 agents and k≥n2k\ge n^2k≥n2 tasks, for every truthful additive mechanism (x,p)(x,p)(x,p) and every real c<nc<nc<n,

∃ t, ∃ y:g(x(t),t)>c⋅g(y,t).\exists\,t,\ \exists\,y:\qquad g\bigl(x(t),t\bigr)>c\cdot g(y,t).∃t, ∃y:g(x(t),t)>c⋅g(y,t).

The goal leaves the mechanism, its tie-breaking and ccc completely general. It says that the ratio nnn of MinWork is optimal among additive mechanisms.

Milestones

  1. Proposition 4.4 (Independence): the payment depends on agent iii's declaration only through its allocation.
  2. Proposition 4.5 (Maximization): xi(t)x^i(t)xi(t) maximizes pi(X,t−i)−ti(X)p^i(X,t^{-i})-t^i(X)pi(X,t−i)−ti(X) over the sets XXX agent iii can obtain against t−it^{-i}t−i.
  3. Pigeonhole: with k≥n2k\ge n^2k≥n2 tasks some agent receives at least nnn tasks.
  4. Claim 4.11: at the all-ones type vector ttt, lowering agent iii's times to 1−ϵ1-\epsilon1−ϵ on xi(t)x^i(t)xi(t) and ϵ\epsilonϵ elsewhere keeps all of xi(t)x^i(t)xi(t) with agent iii, provided the empty set is attainable for agent iii.
  5. The ratio step: at that perturbed type vector, an allocation giving agent iii a fixed set of nnn tasks has make-span at least (1−ϵ)n(1-\epsilon)n(1−ϵ)n, while some allocation has make-span at most 1+kϵ1+k\epsilon1+kϵ.

Significance

Theorem 4.10 shows that the gap between MinWork's ratio nnn and the general lower bound 222 cannot be closed by any mechanism that prices tasks separately, and so any better mechanism would have to couple the prices of different tasks. It was the first class-restricted confirmation of Conjecture 4.9, which was eventually proved for all truthful mechanisms (Christodoulou–Koutsoupias–Kovács 2023). The additive case is the cleanest entry point: its proof needs only the two basic properties of truthful mechanisms, Independence and Maximization, which every later lower bound also uses.

All results here are proved on paper, and none has a machine-checked proof on Prove2Me as of this mission's drafting. The mission produces a formal model of scheduling mechanisms and prices that other lower bounds can reuse, formal statements of Independence and Maximization, and a formal proof of Theorem 4.10. Along the way the formalization corrects two points of the printed argument (see Formalization scope).

Difficulty

The work is to extract prices from an arbitrary truthful mechanism. Prices are defined through the attainable sets of Definition 12. A natural first idea replaces the mechanism by per-task prices qji(t−i)q^i_j(t^{-i})qji​(t−i) that the agent maximizes against. That gives a different class, because Definition 13 constrains the price of every set of tasks, including sets the mechanism never allocates, whose price is 000.

Claim 4.11 is the critical step, and its printed argument does not go through for an arbitrary truthful additive mechanism. It needs the empty set to be attainable with price 000. For a mechanism with a bounded ratio this holds, because a very slow agent must receive nothing, but this has to be derived from the approximation hypothesis. From that, one has to show that every task of xi(t)x^i(t)xi(t) carries a single-task price of at least 111. The final step also needs care. The paper's "w.l.o.g. ∣x1∣=n|x^1|=n∣x1∣=n" is a further reduction, and the paper's bound g≥∣x1∣g\ge|x^1|g≥∣x1∣ has to be replaced by (1−ϵ)∣x1∣(1-\epsilon)|x^1|(1−ϵ)∣x1∣.

Formalization scope

  • Representation. Agents are Fin n, tasks Fin k, type vectors Fin n → Fin k → ℝ, and an allocation is a map Fin k → Fin n sending each task to its agent. taskSet x i is xix^ixi, and the make-span is a Finset.sup' over the nonempty set of agents ([NeZero n]). A mechanism is a pair alloc, pay of arbitrary functions; nothing about its tie-breaking is fixed.
  • Standing assumptions. Types are positive. Truthfulness, additivity and approximation quantify over positive type vectors only. Utility is quasi-linear, with payments handed to the agent.
  • Prices. price follows Definition 12 literally: the payment at a Classical.choose witness declaration when the set is attainable, and 000 otherwise. Additivity (IsAdditive) is required for every set of tasks, attainable or not, as Definition 13 states. It is a condition on prices, not on the payment function.
  • Explicit threshold. The goal assumes k≥n2k\ge n^2k≥n2, the value the paper's proof starts from; the printed theorem does not mention kkk. It is stated for every n≥1n\ge1n≥1; at n=1n=1n=1 it holds because all allocations coincide.
  • Printed slips, corrected. (1) Proposition 4.5 is stated over attainable sets: over all sets, with the price 000 of unattainable sets, it fails for truthful mechanisms that never give the agent nothing and charge it. (2) Claim 4.11 carries the added hypothesis that ∅\emptyset∅ is attainable for agent iii. Without it the claim is false: a mechanism that always gives agent iii its ∣x∣|x|∣x∣ cheapest tasks and pays nothing is truthful and additive. The claim is stated for an arbitrary agent iii instead of "agent 1 after relabelling". (3) The ratio step states g≥(1−ϵ)ng\ge(1-\epsilon)ng≥(1−ϵ)n where the paper prints g≥∣x1∣≥ng\ge|x^1|\ge ng≥∣x1∣≥n, and states the paper's w.l.o.g. ∣x1∣=n|x^1|=n∣x1∣=n as a hypothesis of the step, not of the goal.
  • Out of scope. Running time and the revelation principle are not modelled; the goal is stated for truthful direct mechanisms, as §4.3 fixes.
  • Ruled out. A trivializing encoding would define additivity through the payment function instead of the prices of Definition 12, fix nnn, prove the ratio for one c<nc<nc<n only, or drop truthfulness. The last makes the claim false: an optimal allocation rule with zero payments is additive and a 111-approximation. The goal here quantifies over every nnn, every c<nc<nc<n and every truthful additive mechanism.
  • Non-vacuity. Every hypothesis of the goal except the ratio is satisfiable (a constant allocation with zero payments is truthful and additive), and the bound is tight by MinWork.
  • Contributions welcome. Proofs of the milestones, a proof that a bounded-ratio mechanism makes ∅\emptyset∅ attainable for every agent, and reuse of the model for the local-mechanism bound (Theorem 4.12).

Selected references

  • N. Nisan, A. Ronen, Algorithmic Mechanism Design, Games and Economic Behavior 35 (2001) 166–196. https://doi.org/10.1006/game.1999.0790
  • G. Christodoulou, E. Koutsoupias, A. Vidali, A lower bound for scheduling mechanisms, SODA 2007; Algorithmica 55, 2009. https://doi.org/10.1007/s00453-008-9165-3
  • G. Christodoulou, E. Koutsoupias, A. Vidali, A characterization of 2-player mechanisms for scheduling, ESA 2008. https://arxiv.org/abs/0807.3427
  • S. Dobzinski, M. Sundararajan, On characterizations of truthful mechanisms for combinatorial auctions and scheduling, EC 2008, pp. 38–47.
  • G. Christodoulou, E. Koutsoupias, A. Kovács, A proof of the Nisan–Ronen conjecture, STOC 2023. https://doi.org/10.1145/3564246.3585176 (arXiv:2301.11905)
8 thms2 active usersReviewed
🏆Completed
Algorithmic Game TheoryMechanism DesignOperations Research+1·Captain: mikedeng1

Algorithmic Mechanism Design IV: No Local Truthful Mechanism Achieves a c-Approximation for Task Scheduling for Any c < nResearch Paper

Motivation

Nisan and Ronen's Algorithmic Mechanism Design (Games and Economic Behavior 35, 2001) asks how well a computational task can be carried out when its inputs are held by self-interested agents who may lie about them. Their test case is scheduling on unrelated machines: tasks must be assigned to agents (machines), each agent privately knows how long it needs for each task, and the planner wants to minimize the time at which the last agent finishes. The paper shows that the mechanism MinWork, which gives each task to the fastest agent and pays it the second-fastest time, is truthful and loses a factor of at most nnn against the optimum, and that no truthful mechanism can do better than a factor 222. It then conjectures (Conjecture 4.9) that the factor nnn cannot be improved by any truthful mechanism.

That conjecture became the Nisan–Ronen conjecture, one of the central questions of algorithmic mechanism design. A sequence of papers raised the general lower bound from 222 to 1+21 + \sqrt 21+2​ (Christodoulou, Koutsoupias and Vidali), to 1+φ≈2.6181 + \varphi \approx 2.6181+φ≈2.618 (Koutsoupias and Vidali) and to larger constants, and Christodoulou, Koutsoupias and Kovács (STOC 2023) finally proved the conjecture for all deterministic truthful mechanisms. In the original paper, Nisan and Ronen confirm the conjecture for two restricted classes of mechanisms, with short direct arguments. This mission concerns the second class, local mechanisms (Theorem 4.12).

Setting

There are kkk tasks j∈{1,…,k}j \in \{1, \dots, k\}j∈{1,…,k} and nnn agents i∈{1,…,n}i \in \{1, \dots, n\}i∈{1,…,n}. A type vector ttt records, for every agent iii and task jjj, the positive time tjit^i_jtji​ agent iii needs for task jjj. An allocation xxx assigns every task to one agent; xix^ixi is the set of tasks of agent iii. For a set XXX of tasks write ti(X)=∑j∈Xtjit^i(X) = \sum_{j \in X} t^i_jti(X)=∑j∈X​tji​. The make-span of xxx is g(x,t)=max⁡iti(xi)g(x, t) = \max_i t^i(x^i)g(x,t)=maxi​ti(xi).

A direct mechanism (x,p)(x, p)(x,p) asks every agent for its type, computes an allocation x(t)x(t)x(t) from the declarations, and hands agent iii the payment pi(t)p^i(t)pi(t). Agent iii's utility is pi(t)−ti(xi(t))p^i(t) - t^i(x^i(t))pi(t)−ti(xi(t)) measured with its true times. The mechanism is truthful if declaring the true type maximizes each agent's utility whatever the other agents declare. The allocation rule is a ccc-approximation if g(x(t),t)≤c⋅g(y,t)g(x(t), t) \le c \cdot g(y, t)g(x(t),t)≤c⋅g(y,t) for every type vector ttt and every allocation yyy.

For a truthful mechanism the payment to agent iii depends only on the set it receives and on the declarations t−it^{-i}t−i of the others (Proposition 4.4). This gives the price offered to agent iii for a set XXX (Definition 12):

pi(X,t−i)={pi(t′i,t−i)if some t′i gives xi(t′i,t−i)=X,0otherwise.p^i(X, t^{-i}) = \begin{cases} p^i(t'^i, t^{-i}) & \text{if some } t'^i \text{ gives } x^i(t'^i, t^{-i}) = X, \\ 0 & \text{otherwise.} \end{cases}pi(X,t−i)={pi(t′i,t−i)0​if some t′i gives xi(t′i,t−i)=X,otherwise.​

A mechanism is local (Definition 14) if pi(X,t−i)p^i(X, t^{-i})pi(X,t−i) depends only on the other agents' times {tjl:l≠i,j∈X}\{t^l_j : l \ne i, j \in X\}{tjl​:l=i,j∈X} on the tasks of XXX. MinWork is local: its price for XXX is ∑j∈Xmin⁡l≠itjl\sum_{j \in X} \min_{l \ne i} t^l_j∑j∈X​minl=i​tjl​.

Formalization targets

Goal: Theorem 4.12

For every n≥1n \ge 1n≥1, every k≥n2k \ge n^2k≥n2 and every real c<nc < nc<n, no truthful local mechanism is a ccc-approximation:

∀(x,p) truthful and local, ∀c<n:∃ t, yg(x(t),t)>c⋅g(y,t).\forall (x, p) \text{ truthful and local},\ \forall c < n:\quad \exists\, t,\ y \quad g(x(t), t) > c \cdot g(y, t).∀(x,p) truthful and local, ∀c<n:∃t, yg(x(t),t)>c⋅g(y,t).

The bound holds for every c<nc < nc<n, so together with MinWork it shows that nnn is the exact best ratio for local truthful mechanisms.

Milestones

  1. Proposition 4.4 (Independence). Payments depend only on the allocated set and on t−it^{-i}t−i.
  2. Proposition 4.5 (Maximization). xi(t)x^i(t)xi(t) maximizes pi(X,t−i)−ti(X)p^i(X, t^{-i}) - t^i(X)pi(X,t−i)−ti(X) over the sets XXX that agent iii can obtain.
  3. Lemma 4.13. Every type vector has type vectors arbitrarily close to it at which each agent's maximizing set is unique.
  4. Claim 4.14, first step. If xi(t)x^i(t)xi(t) is the unique maximizer, lowering agent iii's times on xi(t)x^i(t)xi(t) keeps xi(t)x^i(t)xi(t).
  5. Ratio step. An allocation that gives one agent nnn tasks of time about 111, while every other agent's own tasks are nearly free, has make-span about nnn, while splitting those nnn tasks gives make-span about 111.

Significance

The result. Theorem 4.12 settles the Nisan–Ronen conjecture for a natural class of mechanisms. Locality captures the mechanisms in which the price for a bundle of tasks is set only by the competition for those tasks. It includes MinWork and, more generally, every mechanism that prices tasks separately using the other agents' bids on them. The theorem says that for this class the trivial per-task auction is already optimal, so any improvement over the ratio nnn must use prices that depend on the other agents' times on tasks outside the bundle.

Formalizing it. The statement is not open: it follows from the 2023 proof of the Nisan–Ronen conjecture, and Nisan and Ronen's own argument is much shorter. That argument is a sketch, though. Lemma 4.13 rests on an informal measure-theoretic appeal, and the core claim relies on a maximization property stated over all sets of tasks. A machine-checked proof pins down exactly which properties of truthful mechanisms the short argument needs. None of these results is known to have been formalized. The definitions (type vectors, truthful mechanisms, prices, locality) are shared with the other missions of this series.

Difficulty

An argument that looks at one agent at a time does not go through. Changing one agent's declaration changes the prices offered to every other agent, so an allocation that is stable for one agent can shift for another. The argument needs a type vector at which every agent's choice is strict, and only then can it lower times agent by agent and follow the allocation. Producing such a type vector is Lemma 4.13. The printed argument for it applies a "for almost every type vector" statement to sets defined by the price functions of an arbitrary mechanism, which need not be measurable. A proof must therefore work without any regularity of the mechanism. A second difficulty is Definition 12's convention that a set the agent cannot obtain has price 000. Locality constrains these zero prices too, and the argument has to account for sets that are obtainable at one type vector and not at a nearby one.

Formalization scope

Agents are Fin n, tasks Fin k. An allocation is a function Fin k → Fin n, a type vector is Fin n → Fin k → ℝ, and a mechanism is a pair of functions alloc (declarations to allocation) and pay (declarations to the payment handed to each agent). Utilities are quasi-linear. All types, declarations and misreports are positive, and every truthfulness, locality and approximation quantifier ranges over positive type vectors. The make-span is a Finset.sup' over the nonempty set of agents ([NeZero n]).

Conventions and explicit thresholds:

  • k≥n2k \ge n^2k≥n2. The theorem is printed without a bound on the number of tasks, and its proof begins "Let k≥n2k \ge n^2k≥n2". The goal carries k≥n2k \ge n^2k≥n2 as a hypothesis.
  • Truthfulness is assumed. §4.3 assumes throughout that the mechanism is truthful (by the revelation principle this is no loss). The goal quantifies over all truthful local mechanisms.
  • Prices use Definition 12 literally, including the value 000 for sets the agent cannot obtain, and locality is Definition 14 applied to that price function over all sets XXX, not only single tasks. When several declarations give the same set, the price uses one chosen witness; by Proposition 4.4 the choice does not matter for truthful mechanisms.
  • Proposition 4.5 is stated over the sets the agent can obtain. As printed, over all subsets, it is false for a truthful mechanism that never leaves an agent idle and pays it negative amounts. Uniqueness of maximizers (Lemma 4.13, Claim 4.14) refers to the same family.
  • Lemma 4.13 uses Mathlib's norm on Fin n → Fin k → ℝ, the sup norm. No measurability of the mechanism is assumed.
  • Claim 4.14 is printed at tji=1t^i_j = 1tji​=1 with 0<ε<10 < \varepsilon < 10<ε<1. The first step is stated at any type vector, with 0<ε≤tji0 < \varepsilon \le t^i_j0<ε≤tji​ on the lowered tasks.
  • Running time and computability are out of scope.

Ruled-out trivializations: locality is not restricted to single tasks; the goal does not assume that maximizers are unique at every type vector (that is Lemma 4.13's conclusion at one point, not a hypothesis); and the bound holds for every c<nc < nc<n, not for some.

Needed infrastructure: finite sums over allocation fibres, sup norms on function spaces, and a genericity argument for finitely many affine functions (Lemma 4.13). The model file and the price and locality definitions are reusable in the other missions of the series. Proofs of individual milestones are welcome independently.

Selected references

  • N. Nisan, A. Ronen, Algorithmic Mechanism Design, Games and Economic Behavior 35 (2001) 166–196. https://doi.org/10.1006/game.1999.0790
  • A. Mas-Colell, M. D. Whinston, J. R. Green, Microeconomic Theory, Oxford University Press, 1995 (pp. 876–880, basic properties of truthful mechanisms).
  • G. Christodoulou, E. Koutsoupias, A. Vidali, A lower bound for scheduling mechanisms, Algorithmica 55 (2009).
  • E. Koutsoupias, A. Vidali, A lower bound of 1+φ for truthful scheduling mechanisms, Algorithmica 66 (2013).
  • G. Christodoulou, E. Koutsoupias, A. Kovács, A proof of the Nisan–Ronen conjecture, STOC 2023.
8 thms3 active usersReviewed
🏆Completed
Algorithmic Game TheoryMechanism DesignOperations Research+1·Captain: mikedeng1

Algorithmic Mechanism Design VIII: A Truthful Approximation Scheme for Bounded Scheduling with VerificationResearch Paper

Motivation

Algorithmic mechanism design asks for algorithms whose inputs are held by self-interested agents: the designer can pay the agents, and must choose payments so that each agent's own interest leads it to reveal what the algorithm needs. Nisan and Ronen introduced the framework with task scheduling on unrelated machines as the running example (Nisan–Ronen 2001). In the basic model, where payments depend only on what the agents declare, they showed that no truthful mechanism approximates the optimal make-span within a factor below 2, and that the natural mechanism only reaches a factor nnn.

Their Section 5 changes the information available: in a mechanism with verification the payments may also depend on the times in which the tasks were actually performed. With this extra information, an exact optimizer becomes a strongly truthful mechanism (Theorem 5.1, the Compensation-and-Bonus mechanism). Exact scheduling on unrelated machines is NP-hard, so the question is whether an approximation algorithm can take the optimizer's place. Theorem 5.6 of the paper shows that plugging a non-optimal algorithm into Compensation-and-Bonus destroys truthfulness in general. Theorem 5.9, the subject of this mission, shows that for the bounded problem a specific approximation scheme, the rounding algorithm of Horowitz and Sahni (1976), can be combined with a modified payment rule to give a truthful mechanism whose outcome is within a factor 1+ε1+\varepsilon1+ε of optimal.

Setting

There are nnn agents and kkk tasks. Agent iii needs time tjit^i_jtji​ for task jjj; the vector t=(tji)t = (t^i_j)t=(tji​) is the type vector, and agent iii alone knows its row tit^iti. In the bounded scheduling problem (Definition 33) there are fixed numbers 0<a<b0 < a < b0<a<b with a≤tji≤ba \le t^i_j \le ba≤tji​≤b for all i,ji, ji,j, and every declaration lies in the same range. An allocation xxx assigns each task to one agent; xix^ixi is the set of tasks of agent iii.

A strategy of agent iii has two parts: a declaration di∈[a,b]kd^i \in [a,b]^kdi∈[a,b]k, and an execution, which for each decision xxx of the mechanism specifies the actual time t~j≥tji\tilde t_j \ge t^i_jt~j​≥tji​ in which agent iii performs each task j∈xij \in x^ij∈xi. The mechanism chooses x=x(d)x = x(d)x=x(d) from the declarations alone and afterwards observes the actual times t~\tilde tt~. The objective is the make-span with actual times,

g(x,t~)=max⁡i∑j∈xit~j.g(x,\tilde t) = \max_i \sum_{j \in x^i} \tilde t_j .g(x,t~)=imax​j∈xi∑​t~j​.

Agent iii receives a payment pip^ipi and has utility pi−∑j∈xit~jp^i - \sum_{j \in x^i} \tilde t_jpi−∑j∈xi​t~j​.

The corrected time vector of agent iii keeps agent iii's actual times on its own tasks and the other agents' declarations elsewhere: corri(x,d,t~)j=t~j\mathrm{corr}^i(x,d,\tilde t)_j = \tilde t_jcorri(x,d,t~)j​=t~j​ for j∈xij \in x^ij∈xi and djld^l_jdjl​ for j∈xlj \in x^lj∈xl, l≠il \ne il=i. For a step δ>0\delta > 0δ>0, r^=δ⌈r/δ⌉\hat r = \delta\lceil r/\delta\rceilr^=δ⌈r/δ⌉ rounds rrr up to a multiple of δ\deltaδ, and g^(x,τ)=g(x,τ^)\hat g(x,\tau) = g(x,\hat\tau)g^​(x,τ)=g(x,τ^).

The rounding mechanism (Definition 34) allocates with an algorithm that exactly solves the problem with rounded declarations d^\hat dd^, and pays

pi=∑j∈xit~j  −  g^(x,corri(x,d,t~)).p^i = \sum_{j\in x^i}\tilde t_j \;-\; \hat g\big(x, \mathrm{corr}^i(x, d, \tilde t)\big).pi=j∈xi∑​t~j​−g^​(x,corri(x,d,t~)).

The first term, the compensation, uses exact actual times; the second, the bonus, uses rounded quantities.

A strategy is dominant if it is a best response to every declarations and executions of the others. The mechanism is truthful if every agent has a dominant strategy that declares its true type.

Formalization targets

Goal: Theorem 5.9 without running time

For every ε>0\varepsilon > 0ε>0, every 0<δ≤εa0 < \delta \le \varepsilon a0<δ≤εa and every allocation algorithm solving the rounded problem exactly, the rounding mechanism is truthful, and at every profile of dominant strategies from the class named in the proof (declarations with the true rounded values, executions whose rounded times equal the rounded true times),

g(x(d),t~)≤(1+ε) g(y,t)for every allocation y.g\big(x(d),\tilde t\big) \le (1+\varepsilon)\, g(y,t) \quad \text{for every allocation } y .g(x(d),t~)≤(1+ε)g(y,t)for every allocation y.

Milestones

  1. The solution of the rounded problem is a (1+ε)(1+\varepsilon)(1+ε)-approximation: g(x,t^)≤g(y,t^) ∀yg(x,\hat t) \le g(y,\hat t)\ \forall yg(x,t^)≤g(y,t^) ∀y implies g(x,t)≤(1+ε)g(y,t) ∀yg(x,t) \le (1+\varepsilon) g(y,t)\ \forall yg(x,t)≤(1+ε)g(y,t) ∀y.
  2. After rounding, g^\hat gg^​ is the make-span, g^(x,corr∗(x,d))=g(x,d^)\hat g(x,\mathrm{corr}^*(x,d)) = g(x,\hat d)g^​(x,corr∗(x,d))=g(x,d^), and each agent's utility equals its rounded bonus.
  3. Every strategy with the true rounded values is dominant.
  4. When all agents follow such strategies, the outcome is a (1+ε)(1+\varepsilon)(1+ε)-approximation.
  5. Truth-telling with minimal execution is dominant; hence the mechanism is truthful.

Significance

The result shows that verification does more than make exact optimization truthful: it lets a polynomial-time approximation scheme be implemented in dominant strategies, provided the bonus is computed on the same rounded instance the algorithm optimizes. This contrasts with Theorem 5.6, where an arbitrary approximation algorithm inside Compensation-and-Bonus is not truthful, and with the factor-2 lower bound of the basic model. The principle it illustrates is that the payments must reward exactly the objective the algorithm optimizes.

The paper gives only a proof sketch. Formalizing it makes the argument's hypotheses explicit: which rounding step suffices, what the allocation algorithm must satisfy, and over which strategy profiles the approximation guarantee holds. No machine-checked version of this theorem or of the Compensation-and-Bonus argument is known to exist.

Difficulty

The sketch reduces the theorem to "arguments similar to those in 5.1", but the rounded setting departs from Theorem 5.1 in two ways. Rounding is many-to-one, so an agent's declaration and execution are pinned down only up to their rounded values, and the algorithm's optimality holds only for the rounded instance. Consequently the claim that the strategies with the true rounded values are the only dominant ones does not survive arbitrary tie-breaking: an agent that is always favoured on ties can overstate its rounded time by one step without ever losing, and two such lies at one profile can push the make-span above the (1+ε)(1+\varepsilon)(1+ε) bound. The approximation guarantee therefore has to be stated for the strategy class the proof identifies, not derived from dominance alone. The remaining steps require exact bookkeeping of rounding across sums and of the corrected time vectors, which a proof sketch leaves implicit.

Formalization scope

  • Agents are Fin n with [NeZero n], tasks Fin k; allocations are functions Fin k → Fin n; the make-span is a Finset.sup' over agents. Types and declarations satisfy a ≤ t i j ≤ b with 0 < a < b; actual times are only bounded below by the true times.
  • roundUp δ r = δ * ⌈r / δ⌉. The statement holds for every δ∈(0,εa]\delta \in (0,\varepsilon a]δ∈(0,εa], which covers the intended choice δ=εa\delta = \varepsilon aδ=εa; the paper leaves δ\deltaδ as "a function of aaa and ε\varepsilonε".
  • The Horowitz–Sahni dynamic program is not formalized. The allocation algorithm is a parameter with the hypothesis that it solves the rounded problem exactly; ties are arbitrary, and the goal holds for every such algorithm. Running time ("polynomial time") is out of scope, and with it the role of the upper bound bbb, which is kept as part of the problem.
  • An execution is a function of the decision (Definition 18). Dominance quantifies over all declarations in [a,b][a,b][a,b] and all executions of the others.
  • The payment uses the allocation x(d)x(d)x(d) in the bonus. Definition 34 prints x(t^)x(\hat t)x(t^); since the rounding algorithm rounds the declarations itself, x(d)x(d)x(d) is the allocation actually computed. The hat on corr\mathrm{corr}corr is absorbed by g^\hat gg^​.
  • The goal's approximation part is restricted to dominant profiles of the class named in the proof, because the unrestricted form (Definition 3, every dominant profile) is false for some tie-breaking rules; an explicit two-agent, one-task instance is recorded in the goal's Formalization Note.
  • A formalization that measures the approximation with declared rather than actual times, lets the allocation read the true types, or states the approximation only at the truthful profile while claiming the general form, does not formalize this theorem.

Useful infrastructure: lemmas on Int.ceil rounding of finite sums and on Finset.sup' monotonicity, and a reusable model of mechanisms with verification. Proofs of the milestones in any order are welcome.

Selected references

  • N. Nisan, A. Ronen, Algorithmic Mechanism Design, Games and Economic Behavior 35 (2001) 166–196. https://doi.org/10.1006/game.1999.0790
  • E. Horowitz, S. Sahni, Exact and Approximate Algorithms for Scheduling Nonidentical Processors, Journal of the ACM 23 (1976) 317–327. https://doi.org/10.1145/321941.321951
8 thms2 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchOptimization·Captain: mikedeng1

Optimal Two- and Three-Stage Production Schedules with Setup Times Included 2: Johnson's Rule for Three MachinesResearch Paper

Motivation

Johnson's 1954 paper in Naval Research Logistics Quarterly is the starting point of machine scheduling theory. Its first section solves the two-machine flow shop: nnn items must pass through machine 1 and then machine 2, and an explicit ordering rule minimizes the total elapsed time. Its second section treats three machines. There the problem "loses some of the nice structure of the two-stage case" (p. 65), and the general three-machine problem was later shown to be strongly NP-hard (Garey, Johnson and Sethi, 1976). Johnson nevertheless identifies a restricted case, in which the middle machine is dominated by the first (or the last), where the two-machine rule still gives an optimal schedule. That case, and the structural facts behind it, are the content of this mission.

The three-machine results are still the reference point for polynomially solvable flow shops and for lower bounds in branch-and-bound methods for the general problem.

Timeline.

  • 1954: Johnson proves the two-machine rule (Theorem 1) and, for three machines, the reduction to a common ordering (Lemma 3), a closed form for the elapsed time, and optimality of the rule on Ai+BiA_i + B_iAi​+Bi​, Bi+CiB_i + C_iBi​+Ci​ when min⁡Ai≥max⁡Bj\min A_i \ge \max B_jminAi​≥maxBj​ (Theorem 2), with the mirror case min⁡Ci≥max⁡Bj\min C_i \ge \max B_jminCi​≥maxBj​ asserted.
  • 1976: Garey, Johnson and Sethi show that minimizing makespan in a three-machine flow shop is strongly NP-hard in general, so some restriction of Theorem 2's kind is unavoidable for an exact ordering rule.

Setting

There are nnn items and three machines. Item iii needs processing time Ai>0A_i > 0Ai​>0 on machine 1, Bi>0B_i > 0Bi​>0 on machine 2 and Ci>0C_i > 0Ci​>0 on machine 3, in that order. Each machine handles at most one item at a time, and processing is not interrupted.

A schedule assigns each item start times si1,si2,si3s^1_i, s^2_i, s^3_isi1​,si2​,si3​. It is feasible when all start times are at least 000 on machine 1, the processing intervals of distinct items on the same machine do not overlap, and si1+Ai≤si2s^1_i + A_i \le s^2_isi1​+Ai​≤si2​, si2+Bi≤si3s^2_i + B_i \le s^3_isi2​+Bi​≤si3​. The three machines may process the items in different orders. The total elapsed time (makespan) is max⁡i(si3+Ci)\max_i (s^3_i + C_i)maxi​(si3​+Ci​).

An ordering σ\sigmaσ lists the items, σ(k)\sigma(k)σ(k) being the item in position kkk. Its as-soon-as-possible schedule processes the items in the order σ\sigmaσ on every machine and starts each item on each machine as early as the rules allow. For an ordering, with positions 1,…,n1, \dots, n1,…,n, Johnson defines

Ku=∑i=1uAi−∑i=1u−1Bi,Hv=∑i=1vBi−∑i=1v−1Ci,K_u = \sum_{i=1}^{u} A_i - \sum_{i=1}^{u-1} B_i, \qquad H_v = \sum_{i=1}^{v} B_i - \sum_{i=1}^{v-1} C_i,Ku​=i=1∑u​Ai​−i=1∑u−1​Bi​,Hv​=i=1∑v​Bi​−i=1∑v−1​Ci​,

the sums running over the items in the first uuu (resp. vvv) positions.

Johnson's three-stage rule says that item iii definitely precedes item jjj when

min⁡(Ai+Bi, Cj+Bj)<min⁡(Aj+Bj, Ci+Bi)(IV)\min(A_i + B_i,\ C_j + B_j) < \min(A_j + B_j,\ C_i + B_i) \tag{IV}min(Ai​+Bi​, Cj​+Bj​)<min(Aj​+Bj​, Ci​+Bi​)(IV)

and calls them indifferent under equality. An ordering is consistent with (IV) when no item placed later is definitely preferred to an item placed earlier.

Formalization targets

Goal: Theorem 2 (p. 67)

If every AiA_iAi​ is at least every BjB_jBj​, then an ordering consistent with (IV) exists, and for every such ordering σ\sigmaσ the as-soon-as-possible schedule of σ\sigmaσ is feasible and satisfies

makespan⁡(as-soon-as-possible schedule of σ)≤makespan⁡(s)for every feasible schedule s.\operatorname{makespan}(\text{as-soon-as-possible schedule of } \sigma) \le \operatorname{makespan}(s) \quad \text{for every feasible schedule } s .makespan(as-soon-as-possible schedule of σ)≤makespan(s)for every feasible schedule s.

Milestones

  1. Lemma 3 (p. 65). Every feasible schedule is matched or beaten by the as-soon-as-possible schedule of some single ordering.
  2. Closed form (p. 66). For every ordering, the total idle time of machine 3 is ∑iYi=max⁡1≤u≤v≤n(Hv+Ku)\sum_i Y_i = \max_{1 \le u \le v \le n}(H_v + K_u)∑i​Yi​=max1≤u≤v≤n​(Hv​+Ku​), so that
makespan⁡=∑i=1nCi+max⁡1≤u≤v≤n(Ku+Hv),\operatorname{makespan} = \sum_{i=1}^{n} C_i + \max_{1 \le u \le v \le n} (K_u + H_v),makespan=i=1∑n​Ci​+1≤u≤v≤nmax​(Ku​+Hv​),

the "maximum walk" of p. 68. 3. Special case (p. 67). If min⁡Ai≥max⁡Bj\min A_i \ge \max B_jminAi​≥maxBj​ then max⁡u≤vKu=Kv\max_{u \le v} K_u = K_vmaxu≤v​Ku​=Kv​, so the makespan is ∑iCi+max⁡v(Hv+Kv)\sum_i C_i + \max_v (H_v + K_v)∑i​Ci​+maxv​(Hv​+Kv​). 4. (III) ⇔\Leftrightarrow⇔ (IV) (p. 67). Interchanging the items in positions j,j+1j, j+1j,j+1 changes HHH and KKK only at j,j+1j, j+1j,j+1, and the interchange is strictly worse for the diagonal terms exactly when (IV) holds. 5. Lemma 4 (p. 67). Relation (IV) is transitive, except when the middle item is indifferent to both others. 6. Mirror case (p. 68). The conclusion of Theorem 2 also holds when every CiC_iCi​ is at least every BjB_jBj​.

Significance

The result. Theorem 2 gives an O(nlog⁡n)O(n \log n)O(nlogn) exact method for a class of three-machine flow shops, in a problem that is strongly NP-hard in general. Lemma 3 says that, for three machines, permutation schedules are dominant; Johnson's example on p. 65 shows this fails for four machines. The closed form of milestone 2 expresses the makespan of any ordering as a longest path in a grid, the device behind most later flow-shop lower bounds.

Formalizing it. All results are proved on paper, some tersely: Lemma 3's proof is two lines and cites the wrong lemma, Lemma 4 is proved by reference to Lemma 2, and the mirror case is asserted without proof. A search of Mathlib and of the platform catalog found no machine-checked proof of any of them. The mission produces a checked account of the three-machine flow shop, including the comparison against all feasible schedules rather than only permutation schedules, and pins down the exact form of the hypotheses (see below).

Difficulty

The interchange argument of the two-machine case does not transfer directly. For a general ordering the makespan involves max⁡u≤v(Hv+Ku)\max_{u \le v}(H_v + K_u)maxu≤v​(Hv​+Ku​), and interchanging adjacent items changes terms that depend on everything placed earlier; the page notes that "the decision is not independent of what precedes the interchanged elements". The hypothesis min⁡A≥max⁡B\min A \ge \max BminA≥maxB is what makes KKK nondecreasing along the ordering, collapsing the double maximum to the diagonal. A second obstacle is that (IV) is not a total preorder: ties break transitivity, so passing from "no adjacent pair can be improved" to "optimal" needs the all-pairs consistency and the tie exception of Lemma 4. Finally, Lemma 3 is a statement about arbitrary start-time schedules, so the reduction to orderings must handle machines whose orders differ.

Formalization scope

Items are Fin n; processing times are real-valued functions A B C : Fin n → ℝ, assumed positive in each theorem that is about schedules (the paper's standing assumption, p. 61). A schedule is three start-time functions; feasibility is spelled out as above with non-overlap written as a disjunction of inequalities. The makespan is the maximum of the machine-3 completion times together with 000, so the empty instance has makespan 000. An ordering is an Equiv.Perm (Fin n) with σ k the item in position k; positions are 0-based, so the Lean K u, H v are the paper's Ku+1K_{u+1}Ku+1​, Hv+1H_{v+1}Hv+1​. Statements with maxima over positions assume n≥1n \ge 1n≥1.

Hypotheses made explicit or corrected:

  • min⁡Ai≥max⁡Bi\min A_i \ge \max B_iminAi​≥maxBi​ is read globally, Bj≤AiB_j \le A_iBj​≤Ai​ for all i,ji, ji,j, as in the section heading. The pointwise reading Bi≤AiB_i \le A_iBi​≤Ai​ makes Theorem 2 false (an instance with five items is recorded in the Formalization Note of the goal).
  • Consistency with (IV) is required for all pairs of positions, not only adjacent ones.
  • Lemma 4 carries Lemma 2's exception for an item indifferent to both others; without it the statement is false.
  • Lemma 3's proof cites "Lemma 2" where Lemma 1 is meant.
  • The interchange equivalence (milestone 4) is stated for arbitrary reals, which is stronger than the page needs.

Optimality in the goal is against every feasible schedule. A formalization that compares only orderings with each other, or that defines the objective as the closed form ∑C+max⁡(Ku+Hv)\sum C + \max(K_u + H_v)∑C+max(Ku​+Hv​), would drop Lemma 3's content and is ruled out: the makespan is the latest completion time of a start-time schedule. The existence clause keeps the optimality clause from being vacuous.

A complete development needs finite sums over initial segments of Fin n, Finset.sup', and permutation manipulations (adjacent transpositions, bubble-sort arguments). The feasibility model and the closed form are reusable for other flow-shop results; contributions of general lemmas on adjacent interchanges of permutations are welcome.

Selected references

  • S. M. Johnson, Optimal two- and three-stage production schedules with setup times included, Naval Research Logistics Quarterly 1(1):61–68, 1954. https://doi.org/10.1002/nav.3800010110
  • M. R. Garey, D. S. Johnson, R. Sethi, The complexity of flowshop and jobshop scheduling, Mathematics of Operations Research 1(2):117–129, 1976. https://doi.org/10.1287/moor.1.2.117
13 thms2 active usersReviewed
🏆Completed
Graph TheoryLinear OptimizationOperations Research+1·Captain: mikedeng1

Project Scheduling with Time Windows and Scarce Resources I: A Time-Feasible Schedule Exists iff the Project Network Has No Cycle of Positive LengthTextbook

Motivation

Project scheduling assigns start times to the activities of a project subject to constraints between them. The classical critical path method (CPM) of Kelley and Walker (1959) and the program evaluation and review technique (PERT) allow only minimum time lags: activity jjj may start no earlier than a given time after activity iii starts. Practice also needs maximum time lags: activity jjj must start no later than a given time after iii. These express deadlines, release dates, time windows and "no wait" couplings. Once maximum time lags are allowed, the project network has cycles and negative arc weights, and even the existence of a schedule is no longer automatic.

This mission is the first of a series on Neumann, Schwindt and Zimmermann, Project Scheduling with Time Windows and Scarce Resources (2nd ed., Springer 2003), a standard reference for resource-constrained project scheduling with general temporal constraints. Chapter 1 contains the temporal part of the theory: feasibility, earliest and latest schedules, floats, and the distance order. Every later chapter adds resource constraints on top of this layer.

Timeline. Roy (1964) introduced the Metra Potential Method, which is scheduling on activity-on-node networks with minimum time lags. Neumann (1975, Sect. 6.4) treated time windows through potentials on networks with arbitrary arc weights. Bartusch, Möhring and Radermacher (1988, Annals of Operations Research 16) developed the general theory of scheduling project networks with resource constraints and time windows, including the feasibility criterion stated below. The book collects these results in Chapter 1.

Setting

A project consists of n≥1n\ge 1n≥1 real activities 1,…,n1,\dots,n1,…,n and two fictitious activities, 000 (project beginning) and n+1n+1n+1 (project completion), so the node set is V={0,1,…,n+1}V=\{0,1,\dots,n+1\}V={0,1,…,n+1}. Each activity iii has an integer duration pip_ipi​, with p0=pn+1=0p_0=p_{n+1}=0p0​=pn+1​=0 and pi>0p_i>0pi​>0 for real activities.

A minimum time lag dijmin⁡d^{\min}_{ij}dijmin​ between two different activities becomes an arc ⟨i,j⟩\langle i,j\rangle⟨i,j⟩ of weight δij=dijmin⁡\delta_{ij}=d^{\min}_{ij}δij​=dijmin​. A maximum time lag dijmax⁡d^{\max}_{ij}dijmax​ becomes a backward arc ⟨j,i⟩\langle j,i\rangle⟨j,i⟩ of weight δji=−dijmax⁡\delta_{ji}=-d^{\max}_{ij}δji​=−dijmax​. There is at most one arc per ordered pair, keeping the tightest lag. The result is the activity-on-node (AoN) network N=(V,E,δ)N=(V,E,\delta)N=(V,E,δ), whose integer weights may be positive, negative or zero and which in general contains cycles. The book establishes that for every node iii there is a path from 000 to iii of nonnegative length and a path from iii to n+1n+1n+1 of length at least pip_ipi​ (p. 8, from Definition 1.1.1 and Remarks 1.1.2). This is the standing assumption of the chapter.

A schedule is a vector S=(S0,…,Sn+1)S=(S_0,\dots,S_{n+1})S=(S0​,…,Sn+1​) of real start times with S0=0S_0=0S0​=0 and Si≥0S_i\ge 0Si​≥0. It is time-feasible if it satisfies the temporal constraints

Sj−Si ≥ δij(⟨i,j⟩∈E),S_j-S_i\ \ge\ \delta_{ij}\qquad(\langle i,j\rangle\in E),Sj​−Si​ ≥ δij​(⟨i,j⟩∈E),

and ST\mathcal S_TST​ is the set of time-feasible schedules. A time-feasible schedule minimizing the project duration Sn+1S_{n+1}Sn+1​ is time-optimal.

The length of a path or cycle is the sum of its arc weights. For an integer L=LSn+1L=LS_{n+1}L=LSn+1​, which is either a prescribed maximum project duration dˉ\bar ddˉ or the shortest project duration, the temporal scheduling network N+N^+N+ adds the backward arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩ with weight −L-L−L. The distance dijd_{ij}dij​ is the length of a longest path from iii to jjj in N+N^+N+, with dii=0d_{ii}=0dii​=0. The earliest and latest start times are ESi=d0iES_i=d_{0i}ESi​=d0i​ and LSi=−di0LS_i=-d_{i0}LSi​=−di0​, the earliest completion time is ECi=ESi+piEC_i=ES_i+p_iECi​=ESi​+pi​, and the total float is TFi=LSi−ESiTF_i=LS_i-ES_iTFi​=LSi​−ESi​. The distance order ≺D\prec_D≺D​ is defined for i≠ji\ne ji=j by: i≺Dji\prec_D ji≺D​j if dij>0d_{ij}>0dij​>0, or dij=0d_{ij}=0dij​=0 and dji<0d_{ji}<0dji​<0.

Formalization targets

Goal: Theorem 1.3.3 (p. 10)

ST≠∅⟺N contains no cycle of positive length.\mathcal S_T\ne\emptyset\quad\Longleftrightarrow\quad N\ \text{contains no cycle of positive length}.ST​=∅⟺N contains no cycle of positive length.

The statement contains no constants. It is the consistency criterion for the temporal constraints and the entry condition for everything else in the book.

Milestones

  1. Distances, §1.3, p. 11, Eq. (1.3.3). If N+N^+N+ has no cycle of positive length, then ddd satisfies dij≥δijd_{ij}\ge\delta_{ij}dij​≥δij​ on E+E^+E+ and the triangle inequality dij≥dih+dhjd_{ij}\ge d_{ih}+d_{hj}dij​≥dih​+dhj​, and it is the smallest family that does.
  2. Earliest and latest schedules, §1.3, p. 12. Under the same hypothesis and the standing assumption, ES=(d0i)iES=(d_{0i})_iES=(d0i​)i​ is time-feasible and lies below every time-feasible schedule. LS=(−di0)iLS=(-d_{i0})_iLS=(−di0​)i​ is time-feasible, satisfies LSn+1≤LLS_{n+1}\le LLSn+1​≤L, and lies above every time-feasible schedule SSS with Sn+1≤LS_{n+1}\le LSn+1​≤L.
  3. Remark 1.3.2 (p. 10). If ST≠∅\mathcal S_T\ne\emptysetST​=∅, there is an integer-valued time-optimal schedule.
  4. Proposition 1.3.8 (p. 15). For a real activity iii, [LSi,ECi[≠∅[LS_i,EC_i[\ne\emptyset[LSi​,ECi​[=∅ if and only if iii is critical (TFi=0TF_i=0TFi​=0) or near-critical (0<TFi<pi0<TF_i<p_i0<TFi​<pi​). A further item of the mission, not a milestone, states the claim of §1.4, p. 17 (after Definition 1.4.3): if N+N^+N+ has no cycle of positive length, ≺D\prec_D≺D​ is a strict order on VVV.

Significance

The result itself. Theorem 1.3.3 tells when the temporal constraints of a project can be met at all. Milestones 1 and 2 identify the earliest and latest schedules with longest path lengths, which makes temporal scheduling a pair of longest-path computations (a forward pass from 000 and a backward pass to 000). Remark 1.3.2 justifies working in integer time. The distance order and the base time intervals [LSi,ECi[[LS_i,EC_i[[LSi​,ECi​[ are the inputs of the resource-constrained methods in Chapters 2 and 3: priority rules schedule along ≺D\prec_D≺D​, and base time intervals give lower bounds on resource usage.

Formalizing it. These results are classical and proved, but the book does not prove Theorem 1.3.3; it points to Neumann (1975) and Bartusch et al. (1988). To our knowledge they have no machine-checked form in this generality, with arbitrary integer weights, cycles, fictitious start and end nodes, and the backward arc of N+N^+N+. The CPM results for acyclic event networks with nonnegative durations already on Prove2Me are a special case. The definitions of this mission (project, AoN network, schedule, N+N^+N+, distances) are intended as the shared substrate for the later missions of the series, which add renewable and cumulative resources.

Difficulty

The necessity direction of the goal is a telescoping sum around a cycle. The sufficiency direction needs a schedule, and the natural candidate Si=d0iS_i=d_{0i}Si​=d0i​ requires three things: longest path lengths must be well defined, they must be finite, and they must satisfy the temporal constraints. With negative weights and cycles, the maximum over walks is unbounded when a positive cycle exists, and a walk-based definition gives nothing. A path-based definition gives a finite maximum but loses the concatenation property, so the triangle inequality becomes a statement about removing nonpositive cycles from walks. S0=0S_0=0S0​=0 and Si≥0S_i\ge 0Si​≥0 further depend on the standing assumption: without it, an arc ⟨i,0⟩\langle i,0\rangle⟨i,0⟩ with positive weight makes ST\mathcal S_TST​ empty although no cycle is positive. The same combinatorics of walks, paths and cycles is behind milestones 1 and 2 and the distance-order item.

Formalization scope

Nodes are Fin (n + 2): 000 is the project beginning and Fin.last (n + 1) the project completion. The field one_le_n records n≥1n\ge 1n≥1. Durations are natural numbers with p0=pn+1=0p_0=p_{n+1}=0p0​=pn+1​=0 and pi>0p_i>0pi​>0 for real activities. Arc weights are arbitrary integers on a loop-free Finset of ordered pairs, so parallel arcs cannot occur. Start times are real; integrality appears only as Remark 1.3.2.

Walks are functions Fin (m + 1) → Fin (n + 2). A path is an injective walk, and a cycle is a closed walk with at least one arc and distinct nodes apart from the repeated endpoint. Distances are maxima over the finitely many paths, valued in WithBot ℤ with ⊥=−∞\bot=-\infty⊥=−∞ for unreachable pairs. No supremum over an unbounded set is taken. ESiES_iESi​ and LSiLS_iLSi​ convert these to integers with junk value 000 for −∞-\infty−∞, and every milestone that uses them carries the hypotheses under which the distances are finite. The backward arc of N+N^+N+ has weight −L-L−L for an integer parameter LLL. If NNN already contains an arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩, the two arcs merge into one carrying the larger weight, as in the book's convention for parallel lags. The standing assumption of p. 8 is a named predicate and a hypothesis of the goal and of milestones 2 and 4.

A trivializing formalization is ruled out: weights are signed integers and cycles are allowed, so the no-positive-cycle condition is not vacuous, and the standing assumption is satisfiable by projects with maximum time lags.

A complete development needs cycle removal from closed walks, the Bellman-type characterization of longest paths without positive cycles, and total unimodularity or a direct integrality argument for Remark 1.3.2. The walk, path and distance layer is reusable for any difference-constraint system. Proofs of milestones and lemmas on walk decomposition are welcome.

Selected references

  • K. Neumann, C. Schwindt, J. Zimmermann, Project Scheduling with Time Windows and Scarce Resources, 2nd ed., Springer, 2003. https://doi.org/10.1007/978-3-540-24800-2
  • M. Bartusch, R. H. Möhring, F. J. Radermacher, "Scheduling project networks with resource constraints and time windows", Annals of Operations Research 16 (1988), 201–240. https://doi.org/10.1007/BF02283745
  • K. Neumann, Operations Research Verfahren, Band III, Hanser, 1975, Sect. 6.4.
  • B. Roy, Les problèmes d'ordonnancement: applications et méthodes, Dunod, 1964.
  • J. E. Kelley, M. R. Walker, "Critical-path planning and scheduling", Proceedings of the Eastern Joint Computer Conference, 1959, 160–173. https://doi.org/10.1145/1460299.1460318
  • R. K. Ahuja, T. L. Magnanti, J. B. Orlin, Network Flows, Prentice Hall, 1993, Sect. 5.4 and 5.6.
8 thms3 active usersReviewed
CombinatoricsGraph TheoryOperations Research+1·Captain: mikedeng1

Project Scheduling with Time Windows and Scarce Resources II: A Time-Feasible Strict Order Is Feasible iff It Breaks Up Every Minimal Forbidden SetTextbook

Motivation

Resource-constrained project scheduling asks for start times of the activities of a project so that prescribed time lags between activities are respected and, at every moment, the activities in progress do not require more of any renewable resource (staff, machines, reactors) than is available. When the time lags include maximum time lags (deadlines relative to other activities), even finding a feasible schedule is NP-hard, and the feasible region is in general neither convex nor connected. Branch-and-bound methods for this problem (the problem PS∣temp∣Cmax⁡PS|temp|C_{\max}PS∣temp∣Cmax​ in the notation of Neumann, Schwindt & Zimmermann) do not search over schedules directly. They search over strict orders of the activities, that is, over sets of precedence constraints "jjj starts after iii has finished".

This mission formalizes the theory behind that search, as developed by Bartusch, Möhring & Radermacher (1988) and presented in §2.3 of Neumann, Schwindt & Zimmermann, Project Scheduling with Time Windows and Scarce Resources (2nd ed., Springer 2003). Its goal, Theorem 2.3.10, says when a strict order resolves every resource conflict.

Setting

A project has activities V={0,1,…,n+1}V = \{0, 1, \dots, n+1\}V={0,1,…,n+1}, where 000 and n+1n+1n+1 are fictitious activities marking the project's start and completion and 1,…,n1, \dots, n1,…,n are the real activities (n≥1n \ge 1n≥1). Activity iii has duration pi∈Z≥0p_i \in \mathbb Z_{\ge 0}pi​∈Z≥0​, with p0=pn+1=0p_0 = p_{n+1} = 0p0​=pn+1​=0 and pi>0p_i > 0pi​>0 for real activities. Time lags are encoded in the project network NNN: an arc ⟨i,j⟩∈E\langle i, j\rangle \in E⟨i,j⟩∈E with integer weight δij\delta_{ij}δij​ imposes Sj−Si≥δijS_j - S_i \ge \delta_{ij}Sj​−Si​≥δij​. The book's standing assumptions give, for every node iii, a path from 000 to iii of nonnegative length and a path from iii to n+1n+1n+1 of length at least pip_ipi​.

A schedule is a vector S∈Rn+2S \in \mathbb R^{n+2}S∈Rn+2 with S0=0S_0 = 0S0​=0 and Si≥0S_i \ge 0Si​≥0. It is time-feasible if Sj−Si≥δijS_j - S_i \ge \delta_{ij}Sj​−Si​≥δij​ for all arcs. The set of time-feasible schedules is ST\mathcal S_TST​.

Each renewable resource k∈Rk \in \mathcal Rk∈R has a capacity RkR_kRk​, and activity iii uses rik≤Rkr_{ik} \le R_krik​≤Rk​ units of it while in progress, with r0k=rn+1,k=0r_{0k} = r_{n+1,k} = 0r0k​=rn+1,k​=0. The active set at time ttt is A(S,t)={i∣Si≤t<Si+pi}\mathcal A(S,t) = \{ i \mid S_i \le t < S_i + p_i\}A(S,t)={i∣Si​≤t<Si​+pi​}, and SSS is resource-feasible if ∑i∈A(S,t)rik≤Rk\sum_{i \in \mathcal A(S,t)} r_{ik} \le R_k∑i∈A(S,t)​rik​≤Rk​ for all kkk and all t≥0t \ge 0t≥0. The feasible region S\mathcal SS consists of the schedules that are both time-feasible and resource-feasible.

A strict order O⊆V×VO \subseteq V \times VO⊆V×V is an asymmetric, transitive relation. Its order polyhedron is

ST(O)={S∈ST∣Sj≥Si+pi for all (i,j)∈O}.\mathcal S_T(O) = \{ S \in \mathcal S_T \mid S_j \ge S_i + p_i \ \text{for all } (i,j) \in O\}.ST​(O)={S∈ST​∣Sj​≥Si​+pi​ for all (i,j)∈O}.

OOO is time-feasible if ST(O)≠∅\mathcal S_T(O) \ne \emptysetST​(O)=∅, and feasible if moreover ST(O)⊆S\mathcal S_T(O) \subseteq \mathcal SST​(O)⊆S. The order network N(O)N(O)N(O) adds to NNN, for each (i,j)∈O(i,j) \in O(i,j)∈O, an arc ⟨i,j⟩\langle i,j\rangle⟨i,j⟩ of weight pip_ipi​, or raises the weight of an existing arc to max⁡(δij,pi)\max(\delta_{ij}, p_i)max(δij​,pi​). A schedule SSS induces the strict order O(S)={(i,j)∣i≠j, Sj≥Si+pi}O(S) = \{(i,j) \mid i \ne j,\ S_j \ge S_i + p_i\}O(S)={(i,j)∣i=j, Sj​≥Si​+pi​}.

A set F⊆VF \subseteq VF⊆V is forbidden if ∑i∈Frik>Rk\sum_{i \in F} r_{ik} > R_k∑i∈F​rik​>Rk​ for some resource kkk. It is a minimal forbidden set if no proper subset of it is forbidden. F\mathcal FF denotes the set of minimal forbidden sets.

Formalization targets

Goal: Theorem 2.3.10 (Bartusch et al. 1988)

For every time-feasible strict order OOO,

O feasible  ⟺  ∀F∈F ∃ i,j∈F: N(O) has a path from i to j of length≥pi.O \text{ feasible} \iff \forall F \in \mathcal F\ \exists\, i, j \in F:\ N(O) \text{ has a path from } i \text{ to } j \text{ of length} \ge p_i .O feasible⟺∀F∈F ∃i,j∈F: N(O) has a path from i to j of length≥pi​.

Milestones

  1. Proposition 2.3.3. A strict order OOO is time-feasible if and only if N(O)N(O)N(O) has no cycle of positive length.
  2. Bartusch et al.'s criterion (quoted in the proof of Theorem 2.3.10). A schedule SSS is resource-feasible if and only if every F∈FF \in \mathcal FF∈F contains distinct i,ji, ji,j with Sj≥Si+piS_j \ge S_i + p_iSj​≥Si​+pi​.
  3. Proposition 2.3.6. For time-feasible SSS, the strict order O(S)O(S)O(S) is feasible if and only if S∈SS \in \mathcal SS∈S.
  4. Theorem 2.3.7. S=⋃O∈OST(O)\mathcal S = \bigcup_{O \in \mathcal O} \mathcal S_T(O)S=⋃O∈O​ST​(O), where O\mathcal OO is the finite set of inclusion-minimal feasible strict orders.
  5. Remark 2.3.11. A time-feasible schedule partitions FFF if and only if every A(S,t)∩F\mathcal A(S,t) \cap FA(S,t)∩F, t≥0t \ge 0t≥0, is feasible. A time-feasible order is feasible if and only if it breaks up all (equivalently, all minimal) forbidden sets. A time-feasible schedule is feasible if and only if it partitions all forbidden sets.

Significance

Theorem 2.3.10 turns the feasibility of a strict order, which is a statement about infinitely many schedules and all times ttt, into a finite check: one longest-path computation in N(O)N(O)N(O) for each minimal forbidden set. Together with Proposition 2.3.3 and the structural Theorem 2.3.7, it shows that S\mathcal SS is a finite union of polyhedra indexed by feasible strict orders. This justifies the enumeration schemes of Chapter 2 of the book (branching on the pairs that break up a minimal forbidden set) and the notions of active and stable schedules developed in later sections.

All results in this mission are proved in the literature. For the resource-feasibility criterion, the book cites Bartusch et al. (1988) instead of proving it. To our knowledge, none of these results has been machine-checked. A formal development would give a verified foundation for the order-based description of the feasible region, on which later missions of this series (active schedules, delaying modes, stable schedules) build.

Difficulty

The sufficiency half of the goal is short once the criterion is available: a path of length ≥pi\ge p_i≥pi​ in N(O)N(O)N(O) forces Sj≥Si+piS_j \ge S_i + p_iSj​≥Si​+pi​ on the whole order polyhedron. The necessity half carries the content. If for some minimal forbidden set FFF no path in N(O)N(O)N(O) between elements of FFF reaches the required length, one must construct a schedule in ST(O)\mathcal S_T(O)ST​(O) in which all activities of FFF are simultaneously in progress. This means adding the reverse constraints Sj−Si<piS_j - S_i < p_iSj​−Si​<pi​ for all i,j∈Fi, j \in Fi,j∈F to the temporal system without creating a cycle of positive length, while keeping S0=0S_0 = 0S0​=0 and S≥0S \ge 0S≥0. The obvious reading "no single arc gives a precedence, so they can overlap" fails because maximum time lags combine into long paths through activities outside FFF. The standing assumption that every node is reachable from 000 by a path of nonnegative length is needed here: without it the equivalence is false.

Formalization scope

  • The activity set is Fin (n + 2): 0 is the project start and Fin.last (n + 1) the project completion. Durations and resource data are natural numbers, arc weights are integers, and start times are real numbers.
  • Strict orders are finite sets of pairs, Finset (Fin (n+2) × Fin (n+2)), required to be asymmetric and transitive.
  • Resource constraints hold for every t≥0t \ge 0t≥0. The book's (2.1.4) writes 0≤t≤dˉ0 \le t \le \bar d0≤t≤dˉ. In Chapter 2 schedules are not bounded by dˉ\bar ddˉ, and the book's proofs and Remark 2.3.11 use all t≥0t \ge 0t≥0. This is a convention of the whole series, not a strengthening.
  • A path is a walk (nodes may repeat) and its length is the sum of its arc weights. A cycle of positive length is a closed walk with at least one arc and positive length. For a time-feasible order, N(O)N(O)N(O) has no cycle of positive length. In that case "some path of length ≥pi\ge p_i≥pi​" coincides with the book's "longest path length ≥pi\ge p_i≥pi​", so no supremum over paths appears.
  • The standing assumptions of the book form a single predicate Project.StandingAssumptions, which is a hypothesis of every theorem: n≥1n \ge 1n≥1; p0=pn+1=0p_0 = p_{n+1} = 0p0​=pn+1​=0 and pi>0p_i > 0pi​>0 otherwise; no loops; r0k=rn+1,k=0r_{0k} = r_{n+1,k} = 0r0k​=rn+1,k​=0 and rik≤Rkr_{ik} \le R_krik​≤Rk​; and the two path conditions of p. 8.
  • Minimal forbidden sets and inclusion-minimal feasible orders use Mathlib's Minimal, taken among forbidden sets and among feasible strict orders respectively.
  • The goal is an equivalence, and both directions are required. Weakening it to sufficiency, or dropping the time-feasibility of OOO or the minimality of FFF, would change the theorem. Keeping the book's cut-off t≤dˉt \le \bar dt≤dˉ would also change it, because a schedule could then have an unresolved conflict after dˉ\bar ddˉ and still be called feasible.
  • Theorem 1.3.3 of Chapter 1 (a time-feasible schedule exists if and only if the network has no cycle of positive length) is needed for Proposition 2.3.3 and is restated here for N(O)N(O)N(O). Chapter 1's mission is drafted separately.
  • Useful infrastructure beyond this mission: longest-path potentials on integer-weighted digraphs without positive cycles (feasibility of difference constraints), and the walk and cycle API on Network. Contributions of this general lemma layer are welcome.

Selected references

  • K. Neumann, C. Schwindt, J. Zimmermann, Project Scheduling with Time Windows and Scarce Resources, 2nd ed., Springer, 2003. https://doi.org/10.1007/978-3-540-24800-2
  • M. Bartusch, R. H. Möhring, F. J. Radermacher, Scheduling project networks with resource constraints and time windows, Annals of Operations Research 16 (1988), 201–240. https://doi.org/10.1007/BF02283745
9 thms2 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchOptimization·Captain: mikedeng1

Project Scheduling with Time Windows and Scarce Resources IV: Every Feasible Schedule Obeys a Minimal Delaying Mode of Each Forbidden SetTextbook

Motivation

Resource-constrained project scheduling with general temporal constraints, written PS∣temp∣Cmax⁡PS|temp|C_{\max}PS∣temp∣Cmax​, asks for start times of the activities of a project that respect minimum and maximum time lags between activities and the capacities of renewable resources (staff, machines, reactors), and that minimize the project duration. Deciding whether a feasible schedule exists at all is already NP-complete (Bartusch, Möhring and Radermacher, 1988), so exact methods are branch-and-bound procedures. The dominant family, going back to De Reyck and Herroelen (1998) and presented in Chapter 2 of Neumann, Schwindt and Zimmermann's monograph, branches on resource conflicts: whenever the currently computed schedule overloads a resource at some time ttt, the set of activities in progress at ttt is a forbidden set, and the node is split into children, each of which adds precedence constraints that resolve the conflict.

Such a scheme is only correct if the children together retain every feasible schedule. Theorem 2.5.7 of the book is exactly this completeness guarantee, and it is the reason the enumeration can be restricted to the small family of minimal delaying modes instead of arbitrary ways of breaking up a conflict. The same section also contains the preprocessing results (§2.5.2) that exploit two-element forbidden sets before any branching happens. This mission formalizes both.

Setting

A project has activities V={0,1,…,n+1}V=\{0,1,\dots,n+1\}V={0,1,…,n+1} with n≥1n\ge1n≥1; activity 000 is the project beginning and n+1n+1n+1 the project completion, both of duration 000, and every real activity i∈{1,…,n}i\in\{1,\dots,n\}i∈{1,…,n} has an integer duration pi>0p_i>0pi​>0. The project network NNN has arc set EEE and integer arc weights δij\delta_{ij}δij​; the arc ⟨i,j⟩\langle i,j\rangle⟨i,j⟩ imposes the temporal constraint Sj−Si≥δijS_j-S_i\ge\delta_{ij}Sj​−Si​≥δij​. A finite set R\mathcal RR of renewable resources is given; resource kkk has capacity Rk∈NR_k\in\mathbb NRk​∈N and activity iii uses rik∈Z≥0r_{ik}\in\mathbb Z_{\ge0}rik​∈Z≥0​ units of it, with rik≤Rkr_{ik}\le R_krik​≤Rk​ and r0k=rn+1,k=0r_{0k}=r_{n+1,k}=0r0k​=rn+1,k​=0.

A schedule is a vector S=(Si)i∈VS=(S_i)_{i\in V}S=(Si​)i∈V​ of real start times with S0=0S_0=0S0​=0 and Si≥0S_i\ge0Si​≥0. The active set at time ttt is A(S,t)={i∈V∣Si≤t<Si+pi}\mathcal A(S,t)=\{i\in V\mid S_i\le t<S_i+p_i\}A(S,t)={i∈V∣Si​≤t<Si​+pi​}. The schedule is time-feasible if it satisfies all temporal constraints, resource-feasible if

∑i∈A(S,t)rik≤Rk(k∈R, t≥0),\sum_{i\in\mathcal A(S,t)}r_{ik}\le R_k\qquad(k\in\mathcal R,\ t\ge0),i∈A(S,t)∑​rik​≤Rk​(k∈R, t≥0),

and feasible if it is both; S\mathcal SS denotes the set of feasible schedules.

A set F⊆VF\subseteq VF⊆V is forbidden if ∑i∈Frik>Rk\sum_{i\in F}r_{ik}>R_k∑i∈F​rik​>Rk​ for some kkk, a feasible set otherwise, and minimal forbidden if no proper subset is forbidden. For a forbidden FFF, a set B⊆FB\subseteq FB⊆F is a delaying alternative if F∖BF\setminus BF∖B is feasible, and a minimal delaying alternative if no proper subset of BBB is one. A minimal delaying mode for FFF is a pair (i,B)(i,B)(i,B) with BBB a minimal delaying alternative for FFF and i∈F∖Bi\in F\setminus Bi∈F∖B.

For §2.5.2, fix an integer upper bound UBUBUB on the project duration. The temporal scheduling network N+N^+N+ adds to NNN the arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩ with weight δn+1,0=−UB\delta_{n+1,0}=-UBδn+1,0​=−UB, and dijd_{ij}dij​ is the longest path length from iii to jjj in N+N^+N+ (−∞-\infty−∞ if there is no path, dii=0d_{ii}=0dii​=0).

Formalization targets

Goal: Theorem 2.5.7 (p. 49)

For every forbidden set FFF and every feasible schedule S∈SS\in\mathcal SS∈S there is a minimal delaying mode (i,B)(i,B)(i,B) for FFF with

Sj≥Si+pi(j∈B).S_j\ge S_i+p_i\qquad(j\in B).Sj​≥Si​+pi​(j∈B).

FFF is arbitrary (not necessarily minimal); BBB must be a minimal delaying alternative and iii must lie outside BBB.

Milestones

  1. Eqs. (2.5.2)–(2.5.3), p. 46. BBB is a minimal delaying alternative for a forbidden FFF iff F∖BF\setminus BF∖B is a maximal feasible subset of FFF, iff B⊆FB\subseteq FB⊆F,
∑i∈F∖Brik≤Rk (k∈R)and∀j∈B ∃k: ∑i∈F∖Brik+rjk>Rk.\sum_{i\in F\setminus B}r_{ik}\le R_k\ (k\in\mathcal R)\quad\text{and}\quad\forall j\in B\ \exists k:\ \sum_{i\in F\setminus B}r_{ik}+r_{jk}>R_k.i∈F∖B∑​rik​≤Rk​ (k∈R)and∀j∈B ∃k: i∈F∖B∑​rik​+rjk​>Rk​.
  1. Bartusch et al.'s criterion (proof of Theorem 2.3.10, p. 35). A schedule is resource-feasible iff every minimal forbidden set FFF contains distinct i,ji,ji,j with Sj≥Si+piS_j\ge S_i+p_iSj​≥Si​+pi​.
  2. Lemma 2.5.5, p. 49. A minimal delaying alternative for FFF is an inclusion-minimal set meeting every minimal forbidden F′⊆FF'\subseteq FF′⊆F.
  3. Theorem 2.5.11, p. 55. If {i,j}\{i,j\}{i,j} is a two-element forbidden set with dij<pid_{ij}<p_idij​<pi​ and dij>−pjd_{ij}>-p_jdij​>−pj​, then every feasible SSS with Sn+1≤UBS_{n+1}\le UBSn+1​≤UB satisfies Sj≥Si+piS_j\ge S_i+p_iSj​≥Si​+pi​.
  4. Eq. (2.5.7), p. 55. If for a two-element forbidden set {i,j}\{i,j\}{i,j} neither dij>−pjd_{ij}>-p_jdij​>−pj​ nor dji>−pid_{ji}>-p_idji​>−pi​ holds, then for all h,l∈Vh,l\in Vh,l∈V and every feasible SSS with Sn+1≤UBS_{n+1}\le UBSn+1​≤UB,
Sl≥Sh+min⁡(dhi+pi+djl, dhj+pj+dil).S_l\ge S_h+\min\bigl(d_{hi}+p_i+d_{jl},\ d_{hj}+p_j+d_{il}\bigr).Sl​≥Sh​+min(dhi​+pi​+djl​, dhj​+pj​+dil​).

Significance

The result. Theorem 2.5.7 is the completeness statement of the De Reyck–Herroelen enumeration scheme (Algorithm 2.5.8): if every child of a conflict node imposes the precedence constraints i→ji\to ji→j (j∈Bj\in Bj∈B) of one minimal delaying mode (i,B)(i,B)(i,B), the children's order polyhedra together contain all feasible schedules of the parent. Proposition 2.5.9(a), the correctness of the whole branch-and-bound procedure, rests on it. Because the objective does not enter, the book reuses the theorem for the regular and nonregular objectives of Chapter 3. Theorem 2.5.11 and inequality (2.5.7) are the preprocessing rules that shrink the time-feasible region before enumeration: each adds temporal constraints that every feasible schedule within the bound already satisfies, which raises the lower bound ESn+1ES_{n+1}ESn+1​ and prunes the enumeration.

Formalizing it. All statements are proved in the book (Bartusch et al.'s criterion is quoted from their 1988 paper with the necessity argument sketched). None of them has a machine-checked proof; the Prove2Me catalog contains precedence-only scheduling models (Brucker–Knust) and acyclic event networks (Kelley–Walker) but no model with time windows and forbidden sets. The mission produces a reusable library of forbidden sets, delaying alternatives and longest-path distances in networks with maximum time lags.

Difficulty

The obvious idea — pick any two overlapping activities and delay one — does not give a minimal delaying alternative with a single delaying activity iii common to all of BBB. The proof has to pass from the pairwise separations that resource-feasibility guarantees in each minimal forbidden subset to a set BBB that is simultaneously minimal as a delaying alternative and ordered behind one activity outside BBB. This needs the correspondence between delaying alternatives and hitting sets of the minimal forbidden subsets (Lemma 2.5.5) and the positivity of real durations to keep iii outside BBB. For the preprocessing results, the delicate part is relating longest paths in N+N^+N+, including the backward arc carrying −UB-UB−UB, to the start-time differences of every feasible schedule within the bound.

Formalization scope

Activities are Fin (n + 2), with n+1n+1n+1 as Fin.last (n + 1). Start times are real; durations, capacities, requirements and time lags are integers (natural numbers where the book says so). The standing assumptions of the book (at least one real activity, zero-duration dummies, positive durations of real activities, no loops, r0k=rn+1,k=0r_{0k}=r_{n+1,k}=0r0k​=rn+1,k​=0, rik≤Rkr_{ik}\le R_krik​≤Rk​, and paths in NNN from 000 to every node and from every node to n+1n+1n+1) are one hypothesis P.StandingAssumptions of every theorem.

Resource constraints are imposed for every t≥0t\ge0t≥0, not only for 0≤t≤dˉ0\le t\le\bar d0≤t≤dˉ as (2.1.4) literally writes. The book's proofs and Remark 2.3.11 use the t≥0t\ge0t≥0 reading; with the literal cut-off, schedules running past dˉ\bar ddˉ could violate capacities after dˉ\bar ddˉ, and Bartusch et al.'s criterion would fail.

Longest path lengths are maxima over simple paths, with values in WithBot ℝ (⊥ for −∞-\infty−∞). If N+N^+N+ has a cycle of positive length, no schedule satisfies the temporal constraints with Sn+1≤UBS_{n+1}\le UBSn+1​≤UB, and the statements using dijd_{ij}dij​ are vacuous, as in the book. The arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩ of N+N^+N+ has weight −UB-UB−UB, or max⁡(δn+1,0,−UB)\max(\delta_{n+1,0},-UB)max(δn+1,0​,−UB) if NNN already has such an arc. UBUBUB is an integer.

The goal is not trivial: it quantifies over minimal delaying modes only. A variant without the minimality of BBB, or allowing i∈Bi\in Bi∈B, would be nearly empty (take B=F∖{i}B=F\setminus\{i\}B=F∖{i}), and the statement here rules both out. Maximality in milestone 1 is taken among subsets of FFF.

Welcome contributions: the hitting-set correspondence between delaying alternatives and minimal forbidden subsets, the telescoping bound Sj−Si≥dijS_j-S_i\ge d_{ij}Sj​−Si​≥dij​ for feasible schedules, and proofs of any milestone. The definitions restate the setup of the series' earlier missions (II: order polyhedra) locally, because those are still drafts.

Selected references

  • K. Neumann, C. Schwindt, J. Zimmermann, Project Scheduling with Time Windows and Scarce Resources, 2nd ed., Springer, 2003, §2.5. https://doi.org/10.1007/978-3-540-24800-2
  • M. Bartusch, R. H. Möhring, F. J. Radermacher, Scheduling project networks with resource constraints and time windows, Annals of Operations Research 16 (1988) 201–240. https://doi.org/10.1007/BF02283745
  • B. De Reyck, W. Herroelen, A branch-and-bound procedure for the resource-constrained project scheduling problem with generalized precedence relations, European Journal of Operational Research 111 (1998) 152–174. https://doi.org/10.1016/S0377-2217(97)00305-6
8 thms3 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchOptimization·Captain: mikedeng1

Project Scheduling with Time Windows and Scarce Resources V: A Schedule Is Inventory-Feasible iff It Resolves Every Minimal Surplus and Shortage SetTextbook

Motivation

In make-to-order production, chemical process industries and other manufacturing settings modelled as projects, activities do not only occupy machines for a while: they also consume intermediate products at their start and deposit products into storage facilities at their completion. Storage is bounded above by a tank or warehouse capacity and below by a safety stock. Resources of this kind are called cumulative resources (or inventory resources, reservoirs in the constraint-programming literature). They were introduced into resource-constrained project scheduling by Neumann and Schwindt (2002), and Chapter 2 of Neumann, Schwindt and Zimmermann, Project Scheduling with Time Windows and Scarce Resources (2nd ed., Springer 2003), develops their theory in §2.12.

A scheduler handling cumulative resources needs a finite combinatorial description of which schedules respect the inventory bounds at every instant, because the time axis is continuous and cannot be checked point by point in a search procedure. Theorem 2.12.4 of the book gives such a description, and it is the basis of the branch-and-bound procedure of Neumann and Schwindt for the problem PSc∣temp∣Cmax⁡PSc|temp|C_{\max}PSc∣temp∣Cmax​.

Setting

A project consists of activities V={0,1,…,n+1}V=\{0,1,\dots,n+1\}V={0,1,…,n+1} with n≥1n\ge 1n≥1, where 000 is the project beginning and n+1n+1n+1 the project completion. Activity iii has an integer duration pi≥0p_i\ge 0pi​≥0, with p0=pn+1=0p_0=p_{n+1}=0p0​=pn+1​=0 and pi>0p_i>0pi​>0 for the real activities.

For each cumulative resource kkk in a set Rγ\mathcal R^\gammaRγ, every activity iii has an integer demand rikr_{ik}rik​. If rik<0r_{ik}<0rik​<0, activity iii withdraws −rik-r_{ik}−rik​ units of kkk at its start; if rik>0r_{ik}>0rik​>0, it deposits rikr_{ik}rik​ units at its completion; rik=0r_{ik}=0rik​=0 means kkk is not used. The demand r0kr_{0k}r0k​ of the project beginning is the initial stock. Write Vk−={i∣rik<0}V_k^-=\{i\mid r_{ik}<0\}Vk−​={i∣rik​<0} and Vk+={i∣rik>0}V_k^+=\{i\mid r_{ik}>0\}Vk+​={i∣rik​>0}. Each resource has a safety stock R‾k∈Z\underline R_k\in\mathbb ZR​k​∈Z and a storage capacity R‾k∈Z\overline R_k\in\mathbb ZRk​∈Z.

A schedule is a vector S=(Si)i∈VS=(S_i)_{i\in V}S=(Si​)i∈V​ of real start times with S0=0S_0=0S0​=0 and Si≥0S_i\ge 0Si​≥0. The active set and the inventory of kkk at time t≥0t\ge 0t≥0 are

Ak(S,t)={i∈Vk−∣Si≤t}∪{i∈Vk+∣Si+pi≤t},rk(S,t)=∑i∈Ak(S,t)rik.\mathcal A_k(S,t)=\{i\in V_k^-\mid S_i\le t\}\cup\{i\in V_k^+\mid S_i+p_i\le t\},\qquad r_k(S,t)=\sum_{i\in\mathcal A_k(S,t)} r_{ik}.Ak​(S,t)={i∈Vk−​∣Si​≤t}∪{i∈Vk+​∣Si​+pi​≤t},rk​(S,t)=i∈Ak​(S,t)∑​rik​.

The schedule is inventory-feasible if R‾k≤rk(S,t)≤R‾k\underline R_k\le r_k(S,t)\le\overline R_kR​k​≤rk​(S,t)≤Rk​ for all kkk and all t≥0t\ge 0t≥0.

Two standing assumptions of the section are used throughout: (2.12.1) R‾k≤∑i∈Vrik≤R‾k\underline R_k\le\sum_{i\in V}r_{ik}\le\overline R_kR​k​≤∑i∈V​rik​≤Rk​, so the final inventory is admissible; and Remark 2.12.2, R‾k≤0≤R‾k\underline R_k\le 0\le\overline R_kR​k​≤0≤Rk​.

A nonempty F⊆VF\subseteq VF⊆V is a kkk-surplus set if ∑i∈Frik>R‾k\sum_{i\in F}r_{ik}>\overline R_k∑i∈F​rik​>Rk​, and a kkk-shortage set if ∑i∈Frik<R‾k\sum_{i\in F}r_{ik}<\underline R_k∑i∈F​rik​<R​k​. A kkk-surplus set FFF is minimal if no kkk-surplus set arises from FFF by removing a nonempty set of replenishing activities, and none arises by adding a nonempty set of depleting activities. Minimal kkk-shortage sets are defined with the roles of replenishing and depleting activities exchanged. Fk+\mathcal F_k^+Fk+​ and Fk−\mathcal F_k^-Fk−​ denote the minimal kkk-surplus and kkk-shortage sets.

Formalization targets

Goal: Theorem 2.12.4

A schedule SSS is inventory-feasible if and only if

∀k, ∀F∈Fk+ ∃j∈F, i∉F: rjk>0, rik<0, Sj+pj≥Si,\forall k,\ \forall F\in\mathcal F_k^+\ \exists j\in F,\ i\notin F:\ r_{jk}>0,\ r_{ik}<0,\ S_j+p_j\ge S_i,∀k, ∀F∈Fk+​ ∃j∈F, i∈/F: rjk​>0, rik​<0, Sj​+pj​≥Si​, ∀k, ∀F∈Fk− ∃j∈F, i∉F: rjk<0, rik>0, Sj≥Si+pi.\forall k,\ \forall F\in\mathcal F_k^-\ \exists j\in F,\ i\notin F:\ r_{jk}<0,\ r_{ik}>0,\ S_j\ge S_i+p_i.∀k, ∀F∈Fk−​ ∃j∈F, i∈/F: rjk​<0, rik​>0, Sj​≥Si​+pi​.

Milestones

  1. The invariance claim after Remark 2.12.2 (p. 131): adding the same integer aka_kak​ to r0kr_{0k}r0k​, R‾k\underline R_kR​k​ and R‾k\overline R_kRk​ does not change the set of inventory-feasible schedules.
  2. Lemma 2.12.3 (a): for every kkk-surplus set FFF there is a minimal kkk-surplus set F′F'F′ with ∅≠F′∩Vk+⊆F∩Vk+\emptyset\ne F'\cap V_k^+\subseteq F\cap V_k^+∅=F′∩Vk+​⊆F∩Vk+​ and F′∩Vk−⊇F∩Vk−F'\cap V_k^-\supseteq F\cap V_k^-F′∩Vk−​⊇F∩Vk−​.
  3. Lemma 2.12.3 (b): the shortage counterpart.
  4. Theorem 2.12.4 (a) on its own: the upper constraints rk(S,t)≤R‾kr_k(S,t)\le\overline R_krk​(S,t)≤Rk​ hold for all t≥0t\ge 0t≥0 iff condition (a) holds.
  5. Theorem 2.12.4 (b) on its own: the lower constraints hold for all t≥0t\ge 0t≥0 iff condition (b) holds.

Significance

The theorem turns a constraint over a continuum of time points into finitely many disjunctions, each a choice among precedence relations. An inventory excess caused by a minimal surplus set is removed by a start-to-completion relation Sj+pj≥SiS_j+p_j\ge S_iSj​+pj​≥Si​ (a replenishment is postponed until after a withdrawal starts, equivalently a maximum time lag), and a shortage by a completion-to-start relation Sj≥Si+piS_j\ge S_i+p_iSj​≥Si​+pi​. Consequences stated in the book: the feasible region of PSc∣temp∣Cmax⁡PSc|temp|C_{\max}PSc∣temp∣Cmax​ is a finite union of polyhedra; branching on these relations, organized as pairs of strict orders and reflexive relations, is a complete search scheme; and minimal delaying alternatives for surplus and shortage sets can be enumerated. Because every problem with renewable resources can be rewritten as one with cumulative resources (p. 130), the book also concludes that this union of polyhedra is in general disconnected.

The result is proved in the book (and in Neumann and Schwindt, 2002). To our knowledge it has no machine-checked proof. This mission produces a Lean formalization of the model, of the one-sided minimality notion, and of the two-sided characterization with its supporting lemma.

Difficulty

The combinatorial core is simple to state but easy to state wrongly. The natural first idea, to use inclusion-minimal surplus sets as for renewable resources, gives a different family Fk+\mathcal F_k^+Fk+​ and a false theorem: the book's minimality allows removing only replenishing activities and adding only depleting ones. The existence lemma needs Remark 2.12.2 to keep at least one replenishing activity in the minimal set, and the sufficiency direction needs (2.12.1) to guarantee a depleting activity outside the minimal set. Both membership conditions of the active set are closed at ttt, so activities that deplete or replenish exactly at the critical instant must be counted on the correct side; a half-open reading changes which schedules are feasible. The initial stock r0kr_{0k}r0k​ is handled by the same active-set rule as any other demand, which matters for the invariance claim.

Formalization scope

  • Activities are Fin (n + 2), activity n+1n+1n+1 is Fin.last (n + 1); resources are an arbitrary type K. Demands, safety stocks and capacities are integers (ℤ); start times are reals (ℝ); durations are natural numbers cast to ℝ.
  • The inventory constraints are required for every t≥0t\ge 0t≥0. The book prints (2.12.2) for 0≤t≤dˉ0\le t\le\bar d0≤t≤dˉ, but its proof of Theorem 2.12.4 works with an arbitrary t≥0t\ge 0t≥0 (the necessity half uses the last completion time of a replenishing activity, which need not be at most dˉ\bar ddˉ). The two readings coincide for schedules with Sn+1≤dˉS_{n+1}\le\bar dSn+1​≤dˉ whose activities all finish by Sn+1S_{n+1}Sn+1​.
  • A schedule satisfies S0=0S_0=0S0​=0 and Si≥0S_i\ge 0Si​≥0 and is not required to be time-feasible; time lags play no role in this section's results and are not part of the model.
  • (2.12.1) and Remark 2.12.2 are explicit hypotheses (TotalDemandWithinBounds, BoundsStraddleZero) wherever the book's proofs use them. Surplus and shortage sets are nonempty by definition, and minimality uses proper inclusions.
  • A formalization in which Fk+\mathcal F_k^+Fk+​ is empty or trivial (for instance, minimality with non-strict inclusions, which no set satisfies) makes condition (a) vacuous; the definitions here follow p. 131 exactly, and a concrete instance with a nonempty Fk+\mathcal F_k^+Fk+​ has been checked locally.

Reusable parts: the cumulative-resource model and inventory profile, which later missions on continuous cumulative resources (§2.12.2) or on the NP-completeness of PSc∣temp∣Cmax⁡PSc|temp|C_{\max}PSc∣temp∣Cmax​ (Theorem 2.12.1) can build on. Contributions welcome: proofs of the lemmas, of either half of the theorem, and finite-sum lemmas about Finset.filter that the proofs need.

Selected references

  • K. Neumann, C. Schwindt, J. Zimmermann, Project Scheduling with Time Windows and Scarce Resources, 2nd ed., Springer, 2003, §2.12.1, pp. 128–135. https://doi.org/10.1007/978-3-540-24800-2
  • K. Neumann, C. Schwindt, Project scheduling with inventory constraints, Mathematical Methods of Operations Research 56 (2003) 513–533 (cited in the book as 2002). https://doi.org/10.1007/s001860200251
7 thms2 active usersReviewed
Graph TheoryLinear OptimizationOperations Research+1·Captain: mikedeng1

Project Scheduling with Time Windows and Scarce Resources VI: Stable, Semistable, Pseudostable and Quasistable Schedules Are Extreme Points of the Feasible RegionTextbook

Motivation

Resource-constrained project scheduling with minimum and maximum time lags is the model behind make-to-order production, process-industry batch planning and large engineering projects. When the objective is the project duration or another regular function (nondecreasing in every start time), an optimum can be found among schedules that cannot be shifted to the left. Many objectives in practice are nonregular: net present value, earliness–tardiness costs, resource levelling and resource investment. For these, delaying an activity can pay, and "shift as far left as possible" no longer identifies a finite set of candidate schedules.

Neumann, Nübel and Schwindt (Math. Methods Oper. Res. 52, 2000) answered this with classes of schedules defined by the absence of pairs of opposite shifts: stable, semistable, pseudostable and quasistable schedules, the mirror image of active, semiactive, pseudoactive and quasiactive schedules. Section 3.2 of Neumann, Schwindt and Zimmermann, Project Scheduling with Time Windows and Scarce Resources (Springer 2003), shows that these classes are exactly the extreme points of the feasible region and of its natural convex pieces. The classification of objective functions in §3.3, and every enumeration scheme of the later chapter, rests on that correspondence.

Setting

A project has activities V={0,1,…,n+1}V=\{0,1,\dots,n+1\}V={0,1,…,n+1} with n≥1n\ge1n≥1. Activity 000 is the project beginning and n+1n+1n+1 the project completion. Activity iii has an integer duration pip_ipi​, with p0=pn+1=0p_0=p_{n+1}=0p0​=pn+1​=0 and pi>0p_i>0pi​>0 otherwise. The project network NNN has node set VVV and arcs ⟨i,j⟩∈E\langle i,j\rangle\in E⟨i,j⟩∈E with integer weights δij\delta_{ij}δij​, each encoding a temporal constraint Sj−Si≥δijS_j-S_i\ge\delta_{ij}Sj​−Si​≥δij​. A prescribed deadline dˉ∈N\bar d\in\mathbb Ndˉ∈N is included as the backward arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩ of weight −dˉ-\bar d−dˉ. Renewable resources kkk have capacities RkR_kRk​, and activity iii uses rik≤Rkr_{ik}\le R_krik​≤Rk​ units while it runs.

A schedule is a vector S∈Rn+2S\in\mathbb R^{n+2}S∈Rn+2 of start times. The time-feasible region ST\mathcal S_TST​ collects the schedules with S0=0S_0=0S0​=0, S≥0S\ge0S≥0 and Sj−Si≥δijS_j-S_i\ge\delta_{ij}Sj​−Si​≥δij​ on every arc; it is a polyhedron, and a polytope when every activity precedes n+1n+1n+1 as in Remarks 1.1.2. A schedule is resource-feasible if at every time t≥0t\ge0t≥0 the running activities A(S,t)={i∣Si≤t<Si+pi}\mathcal A(S,t)=\{i\mid S_i\le t<S_i+p_i\}A(S,t)={i∣Si​≤t<Si​+pi​} use at most RkR_kRk​ units of every resource. The feasible region is S=ST∩SR\mathcal S=\mathcal S_T\cap\mathcal S_RS=ST​∩SR​. It is in general neither convex nor connected.

A schedule induces the strict order O(S)={(i,j)∣i≠j, Sj≥Si+pi}O(S)=\{(i,j)\mid i\ne j,\ S_j\ge S_i+p_i\}O(S)={(i,j)∣i=j, Sj​≥Si​+pi​}. For a strict order OOO, the order polytope is ST(O)={S∈ST∣Sj≥Si+pi ((i,j)∈O)}\mathcal S_T(O)=\{S\in\mathcal S_T\mid S_j\ge S_i+p_i\ ((i,j)\in O)\}ST​(O)={S∈ST​∣Sj​≥Si​+pi​ ((i,j)∈O)}. The order OOO is feasible if ∅≠ST(O)⊆S\emptyset\ne\mathcal S_T(O)\subseteq\mathcal S∅=ST​(O)⊆S. The schedule polytope of SSS is ST(O(S))\mathcal S_T(O(S))ST​(O(S)).

A shift moves a schedule SSS to S′≠SS'\neq SS′=S. It is global if both are feasible, local if in addition a continuous path inside S\mathcal SS joins them, order-preserving if O(S)⊆O(S′)O(S)\subseteq O(S')O(S)⊆O(S′), and order-monotone if O(S)O(S)O(S) and O(S′)O(S')O(S′) are comparable. Two shifts from SSS to S′S'S′ and S′′S''S′′ are opposite if S′′−S=λ(S′−S)S''-S=\lambda(S'-S)S′′−S=λ(S′−S) with λ<0\lambda<0λ<0. A feasible schedule is stable, semistable, pseudostable or quasistable if no pair of opposite global, local, order-monotone or order-preserving shifts, respectively, starts at it. It is antiactive if no global right-shift starts at it.

Formalization targets

Goal: Theorem 3.2.10

For every feasible schedule SSS:

(a) S antiactive  ⟺  S maximal in S,(b) S stable  ⟺  S∈ext⁡S,(c) S semistable  ⟺  S∈ext⁡CS, CS the component of S containing S,(d) S pseudostable  ⟺  S∈ext⁡ST(O) for all feasible O⊆O(S),(e) S quasistable  ⟺  S∈ext⁡ST(O(S)).\begin{aligned} &\text{(a) } S\text{ antiactive}\iff S\text{ maximal in }\mathcal S, \qquad \text{(b) } S\text{ stable}\iff S\in\operatorname{ext}\mathcal S,\\ &\text{(c) } S\text{ semistable}\iff S\in\operatorname{ext}C_S,\ C_S\text{ the component of }\mathcal S\text{ containing }S,\\ &\text{(d) } S\text{ pseudostable}\iff S\in\operatorname{ext}\mathcal S_T(O)\ \text{for all feasible }O\subseteq O(S),\\ &\text{(e) } S\text{ quasistable}\iff S\in\operatorname{ext}\mathcal S_T(O(S)). \end{aligned}​(a) S antiactive⟺S maximal in S,(b) S stable⟺S∈extS,(c) S semistable⟺S∈extCS​, CS​ the component of S containing S,(d) S pseudostable⟺S∈extST​(O) for all feasible O⊆O(S),(e) S quasistable⟺S∈extST​(O(S)).​

Milestones

  • Lemma 3.2.4: opposite order-preserving or order-monotone shifts can be taken uniform (all moved activities move by one common amount).
  • Lemma 3.2.8: pseudostable schedules are the local extreme points of S\mathcal SS, the points on no segment that lies entirely in S\mathcal SS.
  • Lemma 3.2.9: when SSS is not pseudostable, a segment through SSS can be found inside one order polytope ST(O)\mathcal S_T(O)ST​(O) with O⊆O(S)O\subseteq O(S)O⊆O(S) feasible.
  • Proposition 3.2.13: the quasistable schedules, and every class below them in Fig. 3.2.6, form finite sets.
  • Proposition 3.2.16: every vertex of ST\mathcal S_TST​ is the unique solution of S0=0S_0=0S0​=0, Sj−Si=δijS_j-S_i=\delta_{ij}Sj​−Si​=δij​ on the arcs of a spanning tree of NNN; for the minimal point, an outtree rooted at 000.
  • Theorem 3.2.18: SSS is quasistable iff it is the unique solution of such a tree system in the schedule network N(O(S))N(O(S))N(O(S)).
  • Remark 3.2.7: every activity of a quasistable schedule is tied to another one by a tight duration or time lag, so quasistable schedules are integer-valued.

Significance

The theorem makes four shift-defined classes computable objects: extreme points of explicit polytopes, or of a finite union of them. Together with Proposition 3.2.13, it gives each class of nonregular objective functions in §3.3 a finite candidate set of schedules among which an optimum can be sought (§3.2, p. 207). Theorem 3.2.18 gives the certificate for quasistable schedules: a spanning tree of the schedule network, which the later sections use to enumerate vertices.

The results are proved in the book, except Lemma 3.2.9, whose proof is cited to Neumann, Nübel and Schwindt (2000). As far as a search of the platform shows, none of them has been formalized. A formalization supplies the missing details, among them that connected and path components of S\mathcal SS coincide and the degenerate vertices behind the tree description. It also produces a reusable library of schedule classes on real-valued start times.

Difficulty

Part (b) is close to the definition, since a pair of opposite global shifts is a segment through SSS with feasible endpoints. The content is elsewhere. In (c) the definition speaks of continuous trajectories and the right-hand side of connected components, so the proof needs local path-connectedness of a finite union of polytopes. In (d) the feasible region is not convex: an order-monotone shift keeps SSS and S′S'S′ in a common order polytope, but S′S'S′ and S′′S''S′′ may lie in different ones. The segment through SSS has to be moved into a single order polytope ST(O)\mathcal S_T(O)ST​(O) with O⊆O(S)O\subseteq O(S)O⊆O(S), and that is Lemma 3.2.9. Proposition 3.2.16 and Theorem 3.2.18 need the passage from n+2n+2n+2 linearly independent tight constraints to a spanning tree. They must allow degenerate vertices, where several trees describe the same point, and must represent the nonnegativity constraints Si≥0S_i\ge0Si​≥0 by arcs of the network.

Formalization scope

Activities are Fin (n + 2); start times are real vectors Fin (n + 2) → ℝ with the pointwise order. Durations, capacities and requirements are natural numbers, and time lags integers. The deadline is the arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩ of weight −dˉ-\bar d−dˉ, which is always present, as §3.1 prescribes. Resource constraints are imposed for every t≥0t\ge0t≥0, not only for 0≤t≤dˉ0\le t\le\bar d0≤t≤dˉ as (3.1.2) writes; the proofs use the first reading. Extreme points are Mathlib's Set.extremePoints ℝ, maximal points are Maximal for the pointwise order, and components are connectedComponentIn. A local shift carries an explicit continuous map from unitInterval into S\mathcal SS. Strict orders are asymmetric, transitive relations on VVV. A spanning tree is an arc set of size n+1n+1n+1 whose underlying simple graph is connected. Its arcs must be arcs of NNN, resp. of N(O(S))N(O(S))N(O(S)), with their network weights, so an arbitrary equation system does not count.

The schedule classes are defined through shifts and nothing else. Defining "stable" as "extreme point", or "pseudostable" as "local extreme point", would make the goal and Lemma 3.2.8 tautologies, and such encodings are ruled out. Proposition 3.2.16 carries the book's standing convention (§1.2, p. 8) that every node is reached from 000 by a walk of nonnegative length. Without it the statement is false.

The definitions duplicate, under this mission's namespace, the model of the book's Chapter 2 missions (order polytopes, shifts, active classes). They are written to be merged with those once published. Contributions on the geometry of finite unions of polytopes, and on spanning-tree bases of difference constraint systems, are reusable beyond this mission.

Selected references

  • K. Neumann, C. Schwindt, J. Zimmermann, Project Scheduling with Time Windows and Scarce Resources, 2nd ed., Springer, 2003, §3.1–3.2. https://doi.org/10.1007/978-3-540-24800-2
  • K. Neumann, H. Nübel, C. Schwindt, Active and stable project scheduling, Mathematical Methods of Operations Research 52 (2000), 441–465. https://doi.org/10.1007/s001860000092
  • M. Bartusch, R. H. Möhring, F. J. Radermacher, Scheduling project networks with resource constraints and time windows, Annals of Operations Research 16 (1988), 199–240. https://doi.org/10.1007/BF02283745
12 thms1 active userReviewed
Complexity TheoryOperations ResearchOptimization+1·Captain: mikedeng1

Project Scheduling with Time Windows and Scarce Resources IX: Deciding Feasibility with Cumulative Resources Is NP-Complete Even for Acyclic Project NetworksTextbook

Motivation

Project scheduling with cumulative resources models production and logistics projects in which activities fill and empty storage: an activity withdraws material from an inventory when it starts and deposits its output when it completes, and every inventory must stay between a safety stock and a storage capacity. Neumann, Schwindt and Zimmermann treat this model in §2.12 of Project Scheduling with Time Windows and Scarce Resources (2nd ed., Springer 2003, doi:10.1007/978-3-540-24800-2) and use it in their process-industry applications.

Before any optimization, a scheduler must know whether a feasible schedule exists at all. Theorem 2.12.1 of the book answers the complexity of this question: it is NP-complete, and it stays NP-complete when the project network has no cycles. The contrast with renewable resources (machines, workers) is the point of the theorem: with renewable resources, the feasibility problem is NP-complete as well (Theorem 2.3.13, after Bartusch, Möhring and Radermacher, 1988), but an acyclic network always admits a feasible schedule when every requirement is within capacity.

The mission also collects the two other reductions the book proves in full: Proposition 2.5.4 (recognizing whether an activity lies in some minimal delaying alternative, the branching object of the book's branch-and-bound procedures, is NP-complete) and Proposition 3.4.2 (maximizing weighted start-time deviations, a resource-levelling objective, is NP-hard without any resource constraints).

Setting

A project has activities V={0,1,…,n+1}V=\{0,1,\dots,n+1\}V={0,1,…,n+1} with n≥1n\ge1n≥1; activity 000 is the project beginning and n+1n+1n+1 the project completion. Activity iii has a duration pi∈Np_i\in\mathbb Npi​∈N, with p0=pn+1=0p_0=p_{n+1}=0p0​=pn+1​=0. The project network NNN has arc set EEE; an arc ⟨i,j⟩\langle i,j\rangle⟨i,j⟩ with integer weight δij\delta_{ij}δij​ imposes the temporal constraint Sj−Si≥δijS_j-S_i\ge\delta_{ij}Sj​−Si​≥δij​ on the start times. A schedule is a real vector S=(Si)i∈VS=(S_i)_{i\in V}S=(Si​)i∈V​ with S0=0S_0=0S0​=0 and Si≥0S_i\ge0Si​≥0; it is time-feasible if it meets every temporal constraint.

Cumulative resources k∈Rγk\in\mathcal R^\gammak∈Rγ carry integer demands rikr_{ik}rik​: rik<0r_{ik}<0rik​<0 depletes −rik-r_{ik}−rik​ units at the start SiS_iSi​, rik>0r_{ik}>0rik​>0 replenishes rikr_{ik}rik​ units at the completion Si+piS_i+p_iSi​+pi​, and r0kr_{0k}r0k​ is the initial stock. The inventory at time ttt is

rk(S,t)=∑i: rik<0, Si≤trik+∑i: rik>0, Si+pi≤trik.r_k(S,t)=\sum_{i:\ r_{ik}<0,\ S_i\le t} r_{ik}+\sum_{i:\ r_{ik}>0,\ S_i+p_i\le t} r_{ik}.rk​(S,t)=i: rik​<0, Si​≤t∑​rik​+i: rik​>0, Si​+pi​≤t∑​rik​.

With safety stock R‾k\underline R_kR​k​ and storage capacity R‾k\overline R_kRk​ (integers, R‾k≤∑i∈Vrik≤R‾k\underline R_k\le\sum_{i\in V}r_{ik}\le\overline R_kR​k​≤∑i∈V​rik​≤Rk​ by (2.12.1)), SSS is feasible if it is time-feasible and R‾k≤rk(S,t)≤R‾k\underline R_k\le r_k(S,t)\le\overline R_kR​k​≤rk​(S,t)≤Rk​ for every kkk and every t≥0t\ge0t≥0. The decision problem of PSc∣temp∣Cmax⁡PSc|temp|C_{\max}PSc∣temp∣Cmax​ asks whether a feasible schedule exists.

For renewable resources k∈Rk\in\mathcal Rk∈R with capacities RkR_kRk​ and requirements rik∈Nr_{ik}\in\mathbb Nrik​∈N, a set F⊆VF\subseteq VF⊆V is forbidden if ∑i∈Frik>Rk\sum_{i\in F}r_{ik}>R_k∑i∈F​rik​>Rk​ for some kkk. A delaying alternative for FFF is a set B⊆FB\subseteq FB⊆F such that F∖BF\setminus BF∖B is not forbidden; it is minimal if no proper subset of BBB is one.

In PS∞∣temp,dˉ∣fPS\infty|temp,\bar d|fPS∞∣temp,dˉ∣f there are no resources, schedules must also satisfy Sn+1≤dˉS_{n+1}\le\bar dSn+1​≤dˉ, and the objective here is f(S)=−∑i∈V∑j>iwij∣Sj−Si∣f(S)=-\sum_{i\in V}\sum_{j>i}w_{ij}|S_j-S_i|f(S)=−∑i∈V​∑j>i​wij​∣Sj​−Si​∣ with weights wij≥0w_{ij}\ge0wij​≥0.

NP and NP-completeness are taken in the sense of Cook's Turing-machine formulation, with instances written in binary.

Formalization targets

Goal: Theorem 2.12.1

Assuming PARTITION is NP-complete,

L={codes of instances of PSc∣temp∣Cmax⁡ with a feasible schedule}  and  Lacyc=L∩{N acyclic}L=\{\text{codes of instances of }PSc|temp|C_{\max}\text{ with a feasible schedule}\}\ \text{ and }\ L_{\mathrm{acyc}}=L\cap\{N\text{ acyclic}\}L={codes of instances of PSc∣temp∣Cmax​ with a feasible schedule}  and  Lacyc​=L∩{N acyclic}

are both NP-complete.

Milestones

  1. Membership (proof of Theorem 2.12.1): L∈NPL\in\mathrm{NP}L∈NP and Lacyc∈NPL_{\mathrm{acyc}}\in\mathrm{NP}Lacyc​∈NP.
  2. Reduction correctness (proof of Theorem 2.12.1): for sizes s(1),…,s(ν)s(1),\dots,s(\nu)s(1),…,s(ν) with even sum, the project with r0=rn+1=−∑s(i)/2r_0=r_{n+1}=-\sum s(i)/2r0​=rn+1​=−∑s(i)/2, ri=s(i)r_i=s(i)ri​=s(i), R‾=R‾=0\underline R=\overline R=0R​=R=0, d0,n+1min⁡=1d^{\min}_{0,n+1}=1d0,n+1min​=1 has an acyclic network, and it has a feasible schedule iff the sizes split into two parts of equal sum.
  3. Polynomial transformation: PARTITION≤pLacyc\mathrm{PARTITION}\le_p L_{\mathrm{acyc}}PARTITION≤p​Lacyc​.
  4. Proof of Proposition 2.5.4, one resource: for j∗∈B⊆Fj^*\in B\subseteq Fj∗∈B⊆F, BBB is a minimal delaying alternative iff R−min⁡j∈Brj<∑i∈F∖Bri≤RR-\min_{j\in B}r_j<\sum_{i\in F\setminus B}r_i\le RR−minj∈B​rj​<∑i∈F∖B​ri​≤R.
  5. Proof of Proposition 2.5.4, with rj∗=1r_{j^*}=1rj∗​=1: a minimal delaying alternative contains j∗j^*j∗ iff some A⊆F∖{j∗}A\subseteq F\setminus\{j^*\}A⊆F∖{j∗} has ∑i∈Ari=R\sum_{i\in A}r_i=R∑i∈A​ri​=R.
  6. Proposition 2.5.4: assuming SUBSET SUM is NP-complete, deciding whether some minimal delaying alternative for a forbidden set FFF contains j∗∈Fj^*\in Fj∗∈F is NP-complete.
  7. Proof of Proposition 3.4.2: a graph has a cut of at least MMM edges iff the constructed instance has a schedule with Si∈{0,1}S_i\in\{0,1\}Si​∈{0,1} and ∑i<jwij∣Sj−Si∣≥M\sum_{i<j}w_{ij}|S_j-S_i|\ge M∑i<j​wij​∣Sj​−Si​∣≥M.
  8. Proposition 3.4.2: assuming SIMPLE MAX CUT is NP-complete, the decision version of PS∞∣temp,dˉ∣−∑∑wij∣Sj−Si∣PS\infty|temp,\bar d|-\sum\sum w_{ij}|S_j-S_i|PS∞∣temp,dˉ∣−∑∑wij​∣Sj​−Si​∣ is NP-hard.

The goal follows from milestones 1 and 3 together with the transfer of NP-completeness along ≤p\le_p≤p​ (on the platform as CookPvsNP.npComplete_of_polyReducible).

Significance

The result. Theorem 2.12.1 explains why the book's methods for cumulative resources enumerate precedence relations between depleting and replenishing activities (minimal surplus and shortage sets, Theorem 2.12.4) instead of relying on a constructive feasibility test: unless P = NP, no polynomial algorithm decides feasibility, even for acyclic networks, where the renewable-resource case is trivial. Proposition 2.5.4 does the same for the branching scheme of §2.5, and Proposition 3.4.2 places the resource-levelling objectives of Chapter 3 among the hard ones.

Formalizing it. The three results are proved in the book, as short reductions whose delicate steps are left implicit: the polynomial size of a certificate for real-valued schedules, the handling of instances outside the construction (odd sums, empty index sets, oversized items), and the passage from an optimization problem to its decision version. None of the three reductions is machine-checked anywhere known. The mission states them against a single Turing-machine model and a single binary encoding, reusing the published definitions CookPvsNP_defs, so that the reductions compose with the Cook–Levin development already on the platform.

Difficulty

The mathematical content of the reductions is short; the difficulty is in the complexity-theoretic layer. Two steps resist the obvious argument.

First, NP membership. The book's certificate is a schedule, and a schedule is a real vector: it is not a string. A verifier needs a finite certificate of polynomial length, and it is not immediate that a feasible instance has a feasible schedule with small rational (or integer) start times, since the inventory constraints involve strict orderings between event times.

Second, polynomial-time computability in a concrete Turing-machine model. The transformation must compute, on a one-tape machine, binary codes of sums and halves of the input sizes, an arc list of quadratic length, and must map malformed strings to fixed no-instances. Informal "clearly polynomial" arguments have to become explicit machine constructions or a reusable library of closure properties.

Formalization scope

The Lean development fixes the following conventions.

  • Activities are Fin (n + 2), with the completion Fin.last (n + 1); resources are Fin m. Start times are real.
  • The inventory constraints hold for every t≥0t\ge0t≥0, not only for 0≤t≤dˉ0\le t\le\bar d0≤t≤dˉ as (2.12.2) is printed; the book's proofs use this reading.
  • Real activities may have duration 000 in PSc∣temp∣Cmax⁡PSc|temp|C_{\max}PSc∣temp∣Cmax​: the reduction of Theorem 2.12.1 uses only such activities.
  • Instances are coded as lists of integers written in binary over the alphabet {0,1,−,#}\{0,1,-,\#\}{0,1,−,#}; arc weights are listed for every ordered pair of activities together with an arc indicator. Well-formedness (standing assumptions such as n≥1n\ge1n≥1, p0=pn+1=0p_0=p_{n+1}=0p0​=pn+1​=0, no loops, (2.12.1), rik≤Rkr_{ik}\le R_krik​≤Rk​) is part of each language.
  • "Acyclic" means the nodes admit a numbering increasing along every arc.
  • The NP-completeness of PARTITION, SUBSET SUM and SIMPLE MAX CUT (Karp, 1972) enters as a hypothesis of the corresponding theorem; these are not results of the book.
  • Proposition 3.4.2 is stated for the decision version of the optimization problem, with natural-number weights and threshold.

A trivializing formalization is ruled out: the languages contain only codes of well-formed instances, the encoding is injective, and the hypotheses on the source problems are true theorems, so the goal cannot hold vacuously or by a degenerate encoding.

A complete development needs closure properties of polynomial-time computable functions in Cook's model (composition, binary arithmetic, list manipulation), transitivity of ≤p\le_p≤p​, and a small-certificate lemma for systems of difference constraints with strict and non-strict inequalities. These are reusable for every NP-hardness proof stated in the same framework. Contributions to any of them, to the instance-level milestones 2, 4, 5 and 7, or to the NP-completeness of PARTITION, SUBSET SUM and SIMPLE MAX CUT in this model, are welcome.

Selected references

  • K. Neumann, C. Schwindt, J. Zimmermann, Project Scheduling with Time Windows and Scarce Resources, 2nd ed., Springer, 2003. doi:10.1007/978-3-540-24800-2
  • M. Bartusch, R. H. Möhring, F. J. Radermacher, Scheduling project networks with resource constraints and time windows, Annals of Operations Research 16, 1988.
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, W. H. Freeman, 1979.
  • R. M. Karp, Reducibility among combinatorial problems, in Complexity of Computer Computations, Plenum, 1972. doi:10.1007/978-1-4684-2001-2_9
  • S. Cook, The P versus NP Problem, Clay Mathematics Institute problem description. claymath.org
14 thms3 active usersReviewed
CombinatoricsOperations ResearchOptimization+1·Captain: mikedeng1

Online Scheduling of a Single Machine to Minimize Total Weighted Completion Time: Delayed SWPT Has Competitive Ratio 2Research Paper

Motivation

A single machine must process nnn jobs that arrive over time. Job jjj is released at time rjr_jrj​, needs pjp_jpj​ units of uninterrupted processing, and has weight wj>0w_j > 0wj​>0; the goal is to minimize the total weighted completion time ∑jwjCj\sum_j w_j C_j∑j​wj​Cj​. Offline, with all release dates equal to zero, Smith's rule (sequence by nondecreasing pj/wjp_j/w_jpj​/wj​) is optimal (Smith 1956); with arbitrary release dates the problem 1 ∣ rj ∣ ∑wjCj1\,|\,r_j\,|\,\sum w_j C_j1∣rj​∣∑wj​Cj​ is strongly NP-hard (Lenstra, Rinnooy Kan and Brucker 1977).

In the online version the scheduler learns of job jjj only at time rjr_jrj​, and at each moment must either start a released job or keep the machine idle. Its quality is measured by its competitive ratio: the worst case, over all instances, of the ratio between the online schedule's cost and the offline optimum. Release-date scheduling is one of the basic test cases of online optimization.

Timeline:

  • 1996. Hoogeveen and Vestjens show that no online algorithm has competitive ratio below 2, even with equal weights, and give the 2-competitive algorithm Delayed SPT for equal weights.
  • 1997. Hall, Schulz, Shmoys and Wein give a (3+ε)(3+\varepsilon)(3+ε)-competitive algorithm for arbitrary weights, based on geometric intervals and linear programming.
  • 1998. Phillips, Stein and Wein give another 2-competitive algorithm for equal weights, which does not extend to arbitrary weights.
  • 2002. Goemans, Queyranne, Schulz, Skutella and Wang obtain a (1+2)(1+\sqrt2)(1+2​)-competitive deterministic algorithm from an LP relaxation.
  • 2004. Anderson and Potts show that Delayed SWPT has competitive ratio exactly 2 for arbitrary positive weights, matching the lower bound.

Setting

An instance has jobs j∈J={1,…,n}j \in J = \{1,\dots,n\}j∈J={1,…,n} with integer release dates rj≥0r_j \ge 0rj​≥0, integer processing times pj≥1p_j \ge 1pj​≥1 and real weights wj>0w_j > 0wj​>0. A schedule assigns each job an integer start time SjS_jSj​. It is feasible if Sj≥rjS_j \ge r_jSj​≥rj​ for every jjj and no two intervals [Sj,Sj+pj)[S_j, S_j + p_j)[Sj​,Sj​+pj​) overlap; idle time is allowed. Its cost is C(S)=∑jwj(Sj+pj)C(S) = \sum_j w_j (S_j + p_j)C(S)=∑j​wj​(Sj​+pj​).

Delayed SWPT runs over unit time slots [t,t+1)[t, t+1)[t,t+1). When the machine is available at time ttt, it looks at the jobs released by ttt and not yet started, and selects one with the smallest ratio pj/wjp_j/w_jpj​/wj​. Ties go to the smaller pjp_jpj​, then to the smaller index. If pj≤tp_j \le tpj​≤t, it starts jjj at ttt and the machine is busy until t+pjt + p_jt+pj​. Otherwise the machine stays idle and the rule is applied again at t+1t+1t+1. The resulting schedule is written π\piπ, or dswpt I in Lean. In particular no job starts before time pjp_jpj​.

The proof uses three auxiliary problems:

  • the doubled problem (2P), with data (2rj,2pj,wj)(2r_j, 2p_j, w_j)(2rj​,2pj​,wj​);
  • the extended problem (E), with release dates rj′=max⁡{pj,f(rj)}r'_j = \max\{p_j, f(r_j)\}rj′​=max{pj​,f(rj​)}, where f(t)f(t)f(t) is the first time at or after ttt at which π\piπ leaves the machine free;
  • one unit-length gap job gtg_tgt​ for each slot [t,t+1)[t,t+1)[t,t+1) in which Delayed SWPT idles although a job jjj is available. The gap job has release date f(rj)f(r_j)f(rj​) and weight wj/pjw_j/p_jwj​/pj​.

The schedule πE\pi_EπE​ of (E) runs the original jobs as in π\piπ and each gtg_tgt​ in [t,t+1)[t,t+1)[t,t+1).

Formalization targets

Goal: Theorem 8

min⁡{ρ  :  ∑jwjCj(π)≤ρ∑jwjCj(S) for every instance and every feasible schedule S}=2.\min\Bigl\{\rho \;:\; \sum_j w_j C_j(\pi) \le \rho \sum_j w_j C_j(S)\ \text{for every instance and every feasible schedule } S\Bigr\} = 2.min{ρ:j∑​wj​Cj​(π)≤ρj∑​wj​Cj​(S) for every instance and every feasible schedule S}=2.

Lean: IsLeast {ρ | ∀ n I S, IsFeasible I.r I.p S → cost I.w I.p (dswpt I) ≤ ρ * cost I.w I.p S} 2. Both halves are required: the upper bound 222 and the fact that no smaller constant is valid for this algorithm.

Milestones

  1. πj≥pj\pi_j \ge p_jπj​≥pj​ for every job (§2) and rj′≤max⁡{2rj,pj}r'_j \le \max\{2r_j, p_j\}rj′​≤max{2rj​,pj​} (§3.2).
  2. πE\pi_EπE​ is feasible for (E) (§3.2).
  3. Lemma 1. If π∗\pi^*π∗ and μ∗\mu^*μ∗ are optimal for (P) and (2P), then C(μ∗)=2 C(π∗)C(\mu^*) = 2\,C(\pi^*)C(μ∗)=2C(π∗).
  4. Lemma 2. πE\pi_EπE​ is optimal for (E).
  5. Lemma 3. If μ∗\mu^*μ∗ is optimal for (2P) and a feasible σE\sigma_EσE​ for (E) satisfies
∑j∈JwjCj(σE)+∑g∈GwgCg(σE)≤∑j∈JwjCj(μ∗)+∑g∈GwgCg(πE),(1)\sum_{j\in J} w_j C_j(\sigma_E) + \sum_{g\in G} w_g C_g(\sigma_E) \le \sum_{j\in J} w_j C_j(\mu^*) + \sum_{g\in G} w_g C_g(\pi_E), \tag{1}j∈J∑​wj​Cj​(σE​)+g∈G∑​wg​Cg​(σE​)≤j∈J∑​wj​Cj​(μ∗)+g∈G∑​wg​Cg​(πE​),(1)

then C(π)≤2 C(S)C(\pi) \le 2\,C(S)C(π)≤2C(S) for every feasible SSS. 6. Inequality (1) holds for some feasible σE\sigma_EσE​, for every optimal μ∗\mu^*μ∗ of (2P) (§§3.4–3.6).

Significance

The theorem shows that a deterministic online algorithm can match the lower bound of Hoogeveen and Vestjens for arbitrary positive weights. This settles the best competitive ratio for deterministic online algorithms for 1 ∣ rj ∣ ∑wjCj1\,|\,r_j\,|\,\sum w_j C_j1∣rj​∣∑wj​Cj​. The algorithm needs no linear program. The analysis also does not compare the algorithm with a lower bound on the optimum. Instead it shows that the online schedule is optimal for a modified problem (E), and it converts an optimal schedule of (2P) into a schedule of (E).

The result was proved on paper in 2004. Neither Mathlib nor the Prove2Me catalog contains a machine-checked proof of it, or of any competitive ratio for online scheduling with release dates. This mission provides several reusable pieces:

  • an executable, verified-terminating definition of an online scheduling rule;
  • the doubling lemma for release-date problems;
  • the optimality criterion behind Lemma 2;
  • the block-by-block exchange argument of §§3.3–3.6.

Difficulty

The obvious argument fails at Lemma 2. Delayed SWPT is far from optimal for (P) itself, and its idle time is unbounded in relative terms. The proof therefore has to show that the inserted gap jobs make every idle slot "justified", so that a preemptive best-available argument becomes valid for (E). That argument rests on an optimality criterion of Belouadah, Posner and Potts (1992), which is not in Mathlib.

The second difficulty is inequality (1). Once μ∗\mu^*μ∗ is doubled and the gap jobs are inserted, nongap jobs must be shifted, and the gain of each gap-generating job must be charged against the delay of the gap jobs in its block. That accounting (Lemmas 4–7 of the paper) is an induction over blocks with signed differences of completion times.

The natural first idea, plain online SWPT (start the available job with the smallest pj/wjp_j/w_jpj​/wj​ whenever the machine is free), has no finite competitive ratio (Example 1 of the paper), so the delay πj≥pj\pi_j \ge p_jπj​≥pj​ is essential to the bound and must be tracked through the whole argument.

Formalization scope

Conventions committed to in Lean:

  • Data. Jobs are Fin n (0-based, so "smallest index" is the order of Fin n). Times are natural numbers, the paper's standing integer-data assumption (p. 688), and weights are real. Every instance carries pj≥1p_j \ge 1pj​≥1 and wj>0w_j > 0wj​>0.
  • Schedules and optimality. Schedules are integer start times. Feasibility, cost and optimality are defined for any finite job type, so (E), with job type Fin n ⊕ gapTimes I, uses the same notions. "Optimal" means optimal among all feasible nonpreemptive schedules with integer start times.
  • The algorithm. Delayed SWPT is a def: a unit-time simulation that compares ratios by cross-multiplication and re-applies the rule at every slot. It runs to the horizon ∑j(rj+2pj)+1\sum_j (r_j + 2p_j) + 1∑j​(rj​+2pj​)+1. A sorry-free check (not uploaded) shows that every job has started by then, and that the simulation reproduces Examples 3 and 4 of the paper, including the gap times 0,2,3,4,5,60,2,3,4,5,60,2,3,4,5,6 of Table 2.
  • Completion times. In (2P) the completion time is μj∗+2pj\mu^*_j + 2p_jμj∗​+2pj​, and gap jobs have unit length.

The goal quantifies over every feasible schedule of every instance. It cannot be met by restricting the competitor to schedules without idle time or to list schedules, by dropping release-date feasibility, or by leaving jobs unscheduled.

Out of scope:

  • The general lower bound "no online algorithm beats 2" (Example 2 of the paper, due to Hoogeveen and Vestjens) is not part of the mission. The lower half of the goal concerns Delayed SWPT only.
  • The Belouadah–Posner–Potts optimality criterion is an external ingredient of Lemma 2. Solvers may formalize it as a supporting theorem.

Infrastructure that a complete development needs:

  • simulation invariants for the algorithm;
  • exchange and left-shift arguments for single-machine schedules;
  • the job-splitting relaxation behind the best-available criterion.

The schedule vocabulary and the criterion are reusable for other release-date scheduling results. Contributions toward the block lemmas of §§3.3–3.6 (Lemmas 4–7, the bound (9)) are welcome as supporting theorems.

Selected references

  • E. J. Anderson and C. N. Potts, Online Scheduling of a Single Machine to Minimize Total Weighted Completion Time, Mathematics of Operations Research 29(3), 686–697, 2004. https://doi.org/10.1287/moor.1040.0092
  • J. A. Hoogeveen and A. P. A. Vestjens, Optimal On-Line Algorithms for Single-Machine Scheduling, IPCO 1996, LNCS 1084, 404–414. https://doi.org/10.1007/3-540-61310-2_30
  • L. A. Hall, A. S. Schulz, D. B. Shmoys and J. Wein, Scheduling to Minimize Average Completion Time: Off-line and On-line Approximation Algorithms, Mathematics of Operations Research 22(3), 513–544, 1997. https://doi.org/10.1287/moor.22.3.513
  • C. Phillips, C. Stein and J. Wein, Minimizing Average Completion Time in the Presence of Release Dates, Mathematical Programming 82, 199–223, 1998. https://doi.org/10.1007/BF01585872
  • M. X. Goemans, M. Queyranne, A. S. Schulz, M. Skutella and Y. Wang, Single Machine Scheduling with Release Dates, SIAM Journal on Discrete Mathematics 15(2), 165–192, 2002. https://doi.org/10.1137/S089548019936223X
  • H. Belouadah, M. E. Posner and C. N. Potts, Scheduling with Release Dates on a Single Machine to Minimize Total Weighted Completion Time, Discrete Applied Mathematics 36(3), 213–231, 1992. https://doi.org/10.1016/0166-218X(92)90255-9
  • J. K. Lenstra, A. H. G. Rinnooy Kan and P. Brucker, Complexity of Machine Scheduling Problems, Annals of Discrete Mathematics 1, 343–362, 1977. https://doi.org/10.1016/S0167-5060(08)70743-X
  • W. E. Smith, Various Optimizers for Single-Stage Production, Naval Research Logistics Quarterly 3, 59–66, 1956. https://doi.org/10.1002/nav.3800030106
11 thms3 active usersReviewed
Control TheoryOperations ResearchProbability+1·Captain: mikedeng1

Dynamic Scheduling of a System with Two Parallel Servers in Heavy Traffic with Resource Pooling: The Threshold Policy Is Asymptotically OptimalResearch Paper

Motivation

Many service systems route several classes of work to servers with overlapping skills: call centers with cross-trained agents, manufacturing cells with flexible machines, computing clusters with heterogeneous processors. Choosing which server works on which class at each moment is a dynamic scheduling problem. Exact optimal policies are out of reach except in toy cases, so heavy-traffic theory replaces the queueing system by a Brownian control problem, solves that limit problem, and then asks for a policy in the original system whose performance converges to the Brownian optimum. This programme was proposed by Harrison (Harrison 1988), and the parallel server system studied here is the example Harrison used (Harrison, Ann. Appl. Probab. 1998) to show that the greedy static priority rule can be very inefficient.

Bell and Williams (2001) gave the first proof of asymptotic optimality of a continuous-review policy for this system, with renewal arrivals and general service times. Harrison (1998) had treated Poisson arrivals and deterministic service times with a discrete-review policy and a pathwise criterion. Harrison and López (Queueing Systems, 1999) identified the complete resource pooling condition for general parallel server systems. The threshold policy and the proof method of Bell and Williams were later extended to multiserver systems (Bell and Williams, Electron. J. Probab., 2005).

Setting

There are two job classes and two servers. Server 1 serves class 1 (activity 1); server 2 serves class 1 (activity 2) and class 2 (activity 3). A sequence of such systems is indexed by r→∞r\to\inftyr→∞. On a probability space, i.i.d. sequences uˇk(i)\check u_k(i)uˇk​(i) (k=1,2k=1,2k=1,2) and vˇj(i)\check v_j(i)vˇj​(i) (j=1,2,3j=1,2,3j=1,2,3), i≥1i\ge1i≥1, are fixed: strictly positive, mutually independent, with mean one and finite variances αk2,βj2\alpha_k^2,\beta_j^2αk2​,βj2​. In system rrr the interarrival times are ukr(i)=uˇk(i)/λkru_k^r(i)=\check u_k(i)/\lambda_k^rukr​(i)=uˇk​(i)/λkr​ and the service times are vjr(i)=vˇj(i)/μjrv_j^r(i)=\check v_j(i)/\mu_j^rvjr​(i)=vˇj​(i)/μjr​. The renewal processes Akr(t)A_k^r(t)Akr​(t) and Sjr(t)S_j^r(t)Sjr​(t) count arrivals and potential service completions.

A scheduling control policy is an allocation T=(T1,T2,T3)T=(T_1,T_2,T_3)T=(T1​,T2​,T3​), where Tj(t)T_j(t)Tj​(t) is the time devoted to activity jjj in [0,t][0,t][0,t]. Each Tj(t)T_j(t)Tj​(t) is a random variable, each TjT_jTj​ is continuous and nondecreasing from 000, and so are the idle times I1=t−T1I_1=t-T_1I1​=t−T1​ and I2=t−T2−T3I_2=t-T_2-T_3I2​=t−T2​−T3​. The queue lengths

Q1(t)=A1(t)−S1(T1(t))−S2(T2(t)),Q2(t)=A2(t)−S3(T3(t))Q_1(t)=A_1(t)-S_1(T_1(t))-S_2(T_2(t)),\qquad Q_2(t)=A_2(t)-S_3(T_3(t))Q1​(t)=A1​(t)−S1​(T1​(t))−S2​(T2​(t)),Q2​(t)=A2​(t)−S3​(T3​(t))

must be nonnegative. Policies may anticipate the future. The rates satisfy Assumption 3.1: λ1>μ1\lambda_1>\mu_1λ1​>μ1​, 1−(λ1−μ1)/μ2=λ2/μ31-(\lambda_1-\mu_1)/\mu_2=\lambda_2/\mu_31−(λ1​−μ1​)/μ2​=λ2​/μ3​, and the rates converge at rate 1/r1/r1/r to limits with second-order parameters θ1,θ2\theta_1,\theta_2θ1​,θ2​. Assumption 3.2 is h1μ2≥h2μ3h_1\mu_2\ge h_2\mu_3h1​μ2​≥h2​μ3​, and Assumption 3.3 gives finite exponential moments near 000. With Q^r(t)=r−1Qr(r2t)\hat Q^r(t)=r^{-1}Q^r(r^2t)Q^​r(t)=r−1Qr(r2t) the cost is

J^r(Tr)=E(∫0∞e−γt h⋅Q^r(t) dt).\hat J^r(T^r)=\mathbf E\Big(\int_0^\infty e^{-\gamma t}\,h\cdot\hat Q^r(t)\,dt\Big).J^r(Tr)=E(∫0∞​e−γth⋅Q^​r(t)dt).

The threshold policy with Lr=[clog⁡r]L^r=[c\log r]Lr=[clogr] works as follows. Server 1 works whenever it has a class 1 job available. Server 2 serves class 1 with preemptive-resume priority when more than LrL^rLr class 1 jobs are present, and otherwise serves class 2. The Brownian benchmark is built from a two-dimensional Brownian motion X~\tilde XX~ with drift θ\thetaθ and diagonal covariance, from y=(1,μ2/μ3)y=(1,\mu_2/\mu_3)y=(1,μ2​/μ3​), and from the reflected process W~∗=y⋅X~+V~∗\tilde W^*=y\cdot\tilde X+\tilde V^*W~∗=y⋅X~+V~∗ with V~∗(t)=−inf⁡s≤ty⋅X~(s)\tilde V^*(t)=-\inf_{s\le t}y\cdot\tilde X(s)V~∗(t)=−infs≤t​y⋅X~(s). Its cost is J∗=E∫0∞e−γth2 W~∗(t)/y2 dtJ^*=\mathbf E\int_0^\infty e^{-\gamma t}h_2\,\tilde W^*(t)/y_2\,dtJ∗=E∫0∞​e−γth2​W~∗(t)/y2​dt.

Formalization targets

Goal: Theorem 5.3

For ccc larger than a constant c0c_0c0​ that depends only on the model data, and for every sequence {Tr}\{T^r\}{Tr} of scheduling control policies,

lim inf⁡r→∞J^r(Tr) ≥ J∗ = lim⁡r→∞J^r(Tr,∗),J∗<∞.\liminf_{r\to\infty}\hat J^r(T^r)\ \ge\ J^*\ =\ \lim_{r\to\infty}\hat J^r(T^{r,*}),\qquad J^*<\infty .r→∞liminf​J^r(Tr) ≥ J∗ = r→∞lim​J^r(Tr,∗),J∗<∞.

Milestones

  • Proposition B.1: the one-dimensional Skorokhod problem, its explicit solution and its minimality.
  • Appendix A, (181) and (184): Cramér-type deviation bounds for delayed renewal processes.
  • Theorem 7.2: after first reaching LrL^rLr, the class 1 queue stays within Lr−1L^r-1Lr−1 of the threshold, with probability tending to one.
  • Theorem 7.1: (Q^1r,I^1r)⇒(0,0)(\hat Q_1^r,\hat I_1^r)\Rightarrow(0,0)(Q^​1r​,I^1r​)⇒(0,0) under the threshold policy.
  • Lemma 8.1: the fluid-scaled threshold allocations converge to Tˉ∗(t)=(t,λ1−μ1μ2t,λ2μ3t)\bar T^*(t)=(t,\frac{\lambda_1-\mu_1}{\mu_2}t,\frac{\lambda_2}{\mu_3}t)Tˉ∗(t)=(t,μ2​λ1​−μ1​​t,μ3​λ2​​t).
  • Theorem 5.2 (state-space collapse): (Q^1r,Q^2r,I^1r,I^2r)⇒(0,Q~2∗,0,I~2∗)(\hat Q_1^r,\hat Q_2^r,\hat I_1^r,\hat I_2^r)\Rightarrow(0,\tilde Q_2^*,0,\tilde I_2^*)(Q^​1r​,Q^​2r​,I^1r​,I^2r​)⇒(0,Q~​2∗​,0,I~2∗​).
  • Lemma 9.3: along a subsequence achieving a finite lim inf⁡\liminfliminf cost, the fluid-scaled processes converge to (0,λt,μt,Tˉ∗,0)(0,\lambda t,\mu t,\bar T^*,0)(0,λt,μt,Tˉ∗,0).

A further draft theorem states that Definition 5.1 determines an admissible allocation, unique pathwise, whenever Lr≥1L^r\ge1Lr≥1.

Significance

The theorem proves that a simple state-dependent rule, which sends server 2 to class 1 only when the class 1 queue exceeds a logarithmic safety stock, is asymptotically optimal among all policies, including those that anticipate the future. The limiting cost is the explicit optimum of the Brownian control problem. The proof gives a template for heavy-traffic asymptotic optimality under complete resource pooling: a lower bound valid for every policy, and state-space collapse under the proposed policy. The residual process analysis of Section 7 shows how a threshold of order log⁡r\log rlogr makes starvation of server 1 negligible on the diffusion time scale.

The paper's results are proved but not machine-checked; no formal proof exists in any proof assistant. The mission asks for formal statements of the paper's main theorem and its supporting lemmas, followed by formal proofs. Parts of the development are independent of the paper: the one-dimensional Skorokhod map, renewal large deviation bounds, and convergence encodings on path space.

Difficulty

The lower bound must hold for arbitrary, possibly anticipating, policies, so no Markov structure is available. The argument has to pass through fluid limits of an arbitrary cost-minimizing subsequence and a pathwise minimality property, and Fatou's lemma for the limit needs uniform control. For the upper bound, the obvious approach, a static priority rule, is known to fail: it starves server 1 and produces a large class 1 queue. With a threshold policy, the hard step is to show that the class 1 queue, once at the threshold, rarely moves Lr−1L^r-1Lr−1 away from it over a time interval of length r2tr^2tr2t. That requires large deviation estimates for renewal processes started at random, multiparameter stopping times. Showing that J^r(Tr,∗)\hat J^r(T^{r,*})J^r(Tr,∗) converges to J∗J^*J∗, rather than only that the processes converge in distribution, also requires uniform integrability of the scaled queue lengths.

Formalization scope

Classes and activities are indexed by Fin 2 and Fin 3. The i.i.d. sequences keep the paper's index base i≥1i\ge1i≥1, and the systems are indexed by n∈Nn\in\mathbb Nn∈N with r=rn∈[1,∞)r=r_n\in[1,\infty)r=rn​∈[1,∞), rn→∞r_n\to\inftyrn​→∞. Time is real, and every condition is imposed for t≥0t\ge0t≥0. Admissibility is exactly (11)–(14). Measurability in (11) is with respect to the completion of P\mathbf PP, since the paper's space is complete. Finiteness of the renewal processes everywhere on Ω\OmegaΩ, which the paper obtains by discarding a null set, is a hypothesis. Queue lengths are real, costs are lower Lebesgue integrals in [0,∞][0,\infty][0,∞], counting processes take values in N∪{∞}\mathbb N\cup\{\infty\}N∪{∞}, and Λ\LambdaΛ, Λ∗\Lambda^*Λ∗ take values in the extended reals.

The constant c0c_0c0​ is existential and is chosen after the model data and before ccc, the policies and the Brownian motions. The threshold relations are required only for the systems with Lr≥1L^r\ge1Lr≥1, which are all but finitely many. Each convergence to a deterministic limit (Theorem 7.1, Lemmas 8.1 and 9.3) is stated as u.o.c. convergence in probability, the paper's own equivalence (p. 633). Theorem 5.2 is stated in coupling form: there are copies of the processes on one probability space, with Skorokhod paths and the same laws, that converge almost surely uniformly on compacts. This is equivalent to weak convergence in D4\mathbf D^4D4 to a limit with continuous paths. J∗J^*J∗ is defined by (44) from an arbitrary pair of independent standard Brownian motions (Mathlib's IsBrownianReal), not by a closed form.

Two formalizations would make the goal trivial, and both are excluded. Leaving out the requirement that Tr,∗T^{r,*}Tr,∗ actually follow the policy would make the goal false or empty. Narrowing the class of competing policies, for example to non-anticipating ones, would weaken the theorem. A draft theorem also states that the threshold allocation exists and is unique pathwise, so the hypothesis on Tr,∗T^{r,*}Tr,∗ can be satisfied.

The development needs renewal theory (functional central limit theorems, Cramér bounds), multiparameter stopping times, tightness in D\mathbf DD, the Skorokhod representation theorem, the reflection map, and properties of reflected Brownian motion. Contributions are welcome at every level: proofs of milestones, reusable lemmas on renewal processes and the Skorokhod map, and further lemmas of the paper (Lemmas 7.5, 7.6 and 9.2 are not yet stated).

Selected references

  • S. L. Bell and R. J. Williams, Dynamic scheduling of a system with two parallel servers in heavy traffic with resource pooling: asymptotic optimality of a threshold policy, Ann. Appl. Probab. 11 (2001) 608–649. https://doi.org/10.1214/aoap/1015345343
  • J. M. Harrison, Heavy traffic analysis of a system with parallel servers: asymptotic optimality of discrete-review policies, Ann. Appl. Probab. 8 (1998) 822–848.
  • J. M. Harrison and M. J. López, Heavy traffic resource pooling in parallel-server systems, Queueing Systems 33 (1999) 339–368.
  • J. M. Harrison, Brownian models of queueing networks with heterogeneous customer populations, in Stochastic Differential Systems, Stochastic Control Theory and Their Applications, Springer (1988) 147–186.
  • S. L. Bell and R. J. Williams, Dynamic scheduling of a parallel server system in heavy traffic with complete resource pooling: asymptotic optimality of a threshold policy, Electron. J. Probab. 10 (2005) 1044–1115.
  • J. M. Harrison, Brownian Motion and Stochastic Flow Systems, Wiley (1985).
14 thms1 active userReviewed
Linear OptimizationOperations ResearchOptimization·Captain: mikedeng1

Project Scheduling with Time Windows and Scarce Resources III: Active, Semiactive, Pseudoactive and Quasiactive Schedules Are Minimal Points of the Feasible RegionTextbook

Motivation

Exact and heuristic methods for resource-constrained project scheduling do not search the whole continuum of start-time vectors. They enumerate a finite candidate set that is guaranteed to contain an optimal schedule. For machine scheduling and precedence-only project scheduling the classical candidate sets (semiactive and active schedules) are defined by shifting single activities earlier. With general time lags — minimum and maximum delays between the starts of activities — several activities can be rigidly tied together, and single-activity shifts no longer describe the right candidate sets.

Neumann, Nübel and Schwindt (Neumann et al. 2000) introduced shifts of sets of activities and four resulting classes of schedules: active, semiactive, pseudoactive and quasiactive. Section 2.4 of Neumann, Schwindt and Zimmermann, Project Scheduling with Time Windows and Scarce Resources (Springer 2003), characterizes each class geometrically as the minimal points of a subset of the feasible region. The branch-and-bound procedures of §2.5 of the book enumerate exactly these objects: each enumeration node is a strict order OOO together with the minimal point of its order polyhedron.

Setting

A project has activities V={0,1,…,n+1}V=\{0,1,\dots,n+1\}V={0,1,…,n+1}; 000 and n+1n+1n+1 are fictitious activities marking the project beginning and completion, and 1,…,n1,\dots,n1,…,n are the real activities. Activity iii has an integer duration pip_ipi​, with p0=pn+1=0p_0=p_{n+1}=0p0​=pn+1​=0 and pi>0p_i>0pi​>0 otherwise. Time lags are the arcs ⟨i,j⟩∈E\langle i,j\rangle\in E⟨i,j⟩∈E of the project network NNN with integer weights δij\delta_{ij}δij​. A schedule is a vector S=(S0,…,Sn+1)S=(S_0,\dots,S_{n+1})S=(S0​,…,Sn+1​) of real start times with S0=0S_0=0S0​=0 and Si≥0S_i\ge0Si​≥0. It is time-feasible if Sj−Si≥δijS_j-S_i\ge\delta_{ij}Sj​−Si​≥δij​ for all ⟨i,j⟩∈E\langle i,j\rangle\in E⟨i,j⟩∈E; these schedules form the polyhedron ST\mathcal S_TST​.

There are renewable resources k∈Rk\in\mathcal Rk∈R with capacity RkR_kRk​; activity iii uses rik≤Rkr_{ik}\le R_krik​≤Rk​ units while it is in progress. With the active set A(S,t)={i∣Si≤t<Si+pi}\mathcal A(S,t)=\{i\mid S_i\le t<S_i+p_i\}A(S,t)={i∣Si​≤t<Si​+pi​}, a schedule is resource-feasible if ∑i∈A(S,t)rik≤Rk\sum_{i\in\mathcal A(S,t)}r_{ik}\le R_k∑i∈A(S,t)​rik​≤Rk​ for every kkk and every t≥0t\ge0t≥0. The feasible region is S=ST∩SR\mathcal S=\mathcal S_T\cap\mathcal S_RS=ST​∩SR​. It is in general neither convex nor connected.

A schedule induces the strict order O(S)={(i,j)∣i≠j, Sj≥Si+pi}O(S)=\{(i,j)\mid i\ne j,\ S_j\ge S_i+p_i\}O(S)={(i,j)∣i=j, Sj​≥Si​+pi​}. For a strict order OOO, the order polyhedron is ST(O)={S∈ST∣Sj≥Si+pi ∀(i,j)∈O}\mathcal S_T(O)=\{S\in\mathcal S_T\mid S_j\ge S_i+p_i\ \forall (i,j)\in O\}ST​(O)={S∈ST​∣Sj​≥Si​+pi​ ∀(i,j)∈O}. OOO is feasible if ∅≠ST(O)⊆S\emptyset\ne\mathcal S_T(O)\subseteq\mathcal S∅=ST​(O)⊆S. ST(O(S))\mathcal S_T(O(S))ST​(O(S)) is the schedule polyhedron of SSS.

A left-shift from SSS to S′S'S′ means S′≤SS'\le SS′≤S componentwise and S′≠SS'\ne SS′=S. For feasible S≠S′S\ne S'S=S′, the shift is global; it is local if a continuous trajectory x:[0,1]→Sx:[0,1]\to\mathcal Sx:[0,1]→S joins SSS to S′S'S′; it is order-preserving if O(S)⊆O(S′)O(S)\subseteq O(S')O(S)⊆O(S′) and order-monotone if O(S)⊆O(S′)O(S)\subseteq O(S')O(S)⊆O(S′) or O(S)⊇O(S′)O(S)\supseteq O(S')O(S)⊇O(S′). A feasible schedule is active, semiactive, pseudoactive or quasiactive if no global, local, order-monotone or order-preserving left-shift, respectively, starts at it. A minimal point of M⊆Rn+2\mathcal M\subseteq\mathbb R^{n+2}M⊆Rn+2 is a point S∈MS\in\mathcal MS∈M such that no S′∈MS'\in\mathcal MS′∈M satisfies S′≤SS'\le SS′≤S, S′≠SS'\ne SS′=S.

Formalization targets

Goal: Theorem 2.4.9

For a feasible schedule SSS:

(a) S active  ⟺  S is a minimal point of S,(b) S semiactive  ⟺  S is a minimal point of a component of S,(c) S pseudoactive  ⟺  S is the minimal point of ST(O) for every feasible strict order O⊆O(S),(d) S quasiactive  ⟺  S is the minimal point of ST(O(S)).\begin{aligned} &\text{(a) } S \text{ active} &&\iff S \text{ is a minimal point of } \mathcal S,\\ &\text{(b) } S \text{ semiactive} &&\iff S \text{ is a minimal point of a component of } \mathcal S,\\ &\text{(c) } S \text{ pseudoactive} &&\iff S \text{ is the minimal point of } \mathcal S_T(O) \text{ for every feasible strict order } O\subseteq O(S),\\ &\text{(d) } S \text{ quasiactive} &&\iff S \text{ is the minimal point of } \mathcal S_T(O(S)). \end{aligned}​(a) S active(b) S semiactive(c) S pseudoactive(d) S quasiactive​​⟺S is a minimal point of S,⟺S is a minimal point of a component of S,⟺S is the minimal point of ST​(O) for every feasible strict order O⊆O(S),⟺S is the minimal point of ST​(O(S)).​

Part (a) is close to a restatement of the definitions. The content lies in (b), which passes from trajectories to connected components; in (c), which replaces a condition on shifts by a condition on finitely many polyhedra; and in (d), which reduces quasiactivity to a single polyhedron.

Milestones

  1. Lemma 2.4.7: for a strict order OOO with ST(O)≠∅\mathcal S_T(O)\ne\emptysetST​(O)=∅, lb ST(O)lb\,\mathcal S_T(O)lbST​(O) is the unique minimal point of ST(O)\mathcal S_T(O)ST​(O).
  2. §2.4, p. 39: an order-monotone shift is local.
  3. §2.4, p. 42: AS⊆SAS⊆PAS⊆QAS\mathcal{AS}\subseteq\mathcal{SAS}\subseteq\mathcal{PAS}\subseteq\mathcal{QAS}AS⊆SAS⊆PAS⊆QAS.
  4. §2.4, p. 44: the pseudoactive schedules are exactly the local minimal points of S\mathcal SS in the Euclidean metric.
  5. Remark 2.4.10 (a): if S≠∅\mathcal S\ne\emptysetS=∅, some minimal point of S\mathcal SS is an optimal schedule.
  6. Remark 2.4.10 (b): quasiactive schedules are integer-valued, and S≠∅\mathcal S\neq\emptysetS=∅ iff an integer-valued optimal schedule exists.
  7. Proposition 2.10.2: Sn+1≤dˉ=∑i∈Vmax⁡(pi,max⁡⟨i,j⟩∈Eδij)S_{n+1}\le\bar d=\sum_{i\in V}\max(p_i,\max_{\langle i,j\rangle\in E}\delta_{ij})Sn+1​≤dˉ=∑i∈V​max(pi​,max⟨i,j⟩∈E​δij​) for every quasiactive SSS.

Significance

The characterization makes each schedule class checkable and enumerable. By (d), deciding quasiactivity is a longest-path computation in the schedule network. Deciding activeness is NP-hard (Neumann et al. 2000); the same holds for semiactive and pseudoactive schedules, which is why exact algorithms enumerate the quasiactive schedules. Remark 2.4.10 and Proposition 2.10.2 then give the two facts every such algorithm relies on: an optimal schedule lies among the (integer-valued) quasiactive schedules, and all of them fit into the horizon [0,dˉ][0,\bar d][0,dˉ]. Regular objective functions other than the project duration (§2.10) inherit the same candidate sets.

On the formal side, the mission produces a reusable model of PS∣temp∣Cmax⁡PS|temp|C_{\max}PS∣temp∣Cmax​ with real start times and general time lags: time-feasible and resource-feasible schedules, schedule-induced orders, order polyhedra and the four schedule classes. No part of this material is formalized on Prove2Me or, as far as is known, anywhere else. The results are all proved in the literature (Neumann et al. 2000; the book gives proofs or calls them obvious); the work here is to formalize them.

Difficulty

The obvious argument for (b) says "a trajectory stays in one component, so local shifts move within components". The converse needs that two schedules in the same connected component of S\mathcal SS are joined by a path inside S\mathcal SS. That is false for general sets and has to come from the structure of S\mathcal SS as a finite union of order polyhedra (the basic structural theorem of Bartusch, Möhring and Radermacher), which is not part of this mission's statements and must be proved on the way.

For (c), the difficulty is that an order-monotone shift may shrink the order O(S)O(S)O(S). The proof has to produce, from a feasible sub-order O⊆O(S)O\subseteq O(S)O⊆O(S) whose polyhedron has a smaller minimal point, a shift that is short enough to keep every overlap of SSS. This requires the resource feasibility of whole order polyhedra, i.e. that ST(O(S))⊆S\mathcal S_T(O(S))\subseteq\mathcal SST​(O(S))⊆S for feasible SSS. Resource feasibility is a condition on all times t≥0t\ge0t≥0, while the orders only record pairwise relations between activities.

Formalization scope

Activities are Fin (n + 2), with 0 and Fin.last (n + 1) the fictitious ones. Durations are natural numbers, arc weights integers, and start times real. Resource requirements and capacities are natural numbers. The resource constraints hold for every t≥0t\ge0t≥0; (2.1.4) writes 0≤t≤dˉ0\le t\le\bar d0≤t≤dˉ, but the book's proofs use the unrestricted form. Minimal points are Mathlib's Minimal for the componentwise order on Fin (n + 2) → ℝ. Components in (b) are connected components (connectedComponentIn), while local shifts are defined by continuous trajectories from the unit interval, as in Definition 2.4.3. The theorem is stated for feasible SSS, since the schedule classes consist of feasible schedules by definition. The lower bound lblblb is a vector of real infima and is used only for nonempty order polyhedra.

Defining "active" as "minimal point of S\mathcal SS", or any class through its right-hand side, would make the goal trivial. That is ruled out: every class is defined through the shifts of Definitions 2.4.1–2.4.6, including the trajectory condition and the orders O(S)O(S)O(S).

Remark 2.4.8 (the minimal point of ST(O)\mathcal S_T(O)ST​(O) is the vector of longest path lengths in N(O)N(O)N(O)) is not stated, since it needs path lengths and the reachability conventions of Remarks 1.1.2. Contributions of that network layer, and of the structural theorem S=⋃OST(O)\mathcal S=\bigcup_O\mathcal S_T(O)S=⋃O​ST​(O) (Theorem 2.3.7), are welcome as supporting lemmas.

Selected references

  • K. Neumann, C. Schwindt, J. Zimmermann, Project Scheduling with Time Windows and Scarce Resources, 2nd ed., Springer, 2003. https://doi.org/10.1007/978-3-540-24800-2
  • K. Neumann, H. Nübel, C. Schwindt, Active and stable project scheduling, Mathematical Methods of Operations Research 52 (2000), 441–465. https://doi.org/10.1007/s001860000092
  • M. Bartusch, R. H. Möhring, F. J. Radermacher, Scheduling project networks with resource constraints and time windows, Annals of Operations Research 16 (1988), 201–240. https://doi.org/10.1007/BF02283745
  • A. Sprecher, R. Kolisch, A. Drexl, Semi-active, active, and non-delay schedules for the resource-constrained project scheduling problem, European Journal of Operational Research 80 (1995), 94–102. https://doi.org/10.1016/0377-2217(93)E0294-8
10 thms2 active usersReviewed
Convex OptimizationOperations ResearchOptimization·Captain: mikedeng1

Project Scheduling with Time Windows and Scarce Resources VII: A Locally Quasiconcave Objective Always Has a Quasistable Optimal ScheduleTextbook

Motivation

Resource-constrained project scheduling asks for start times of the activities of a project that respect precedence-type time lags and the capacities of renewable resources (machines, crews, equipment). Classical project scheduling minimizes the project duration, a regular objective: delaying an activity never helps. Many objectives met in practice are not regular. The resource investment problem minimizes the cost of the resource capacities that must be procured; resource levelling problems minimize fluctuations of resource usage over time; the resource renting problem trades fixed procurement against time-dependent renting costs; net present value and earliness–tardiness objectives reward late as well as early starts. For such objectives the familiar fact that "some active schedule is optimal" fails, and algorithms need another finite set of candidate schedules that is guaranteed to contain an optimum.

Chapter 3 of Neumann, Schwindt and Zimmermann, Project Scheduling with Time Windows and Scarce Resources (2nd ed., Springer 2003, doi:10.1007/978-3-540-24800-2), organizes the objective functions of project scheduling into seven classes and pairs each class with a class of schedules that contains an optimal schedule. This mission formalizes §3.3 of that chapter. The classification goes back to Neumann, Nübel and Schwindt (2000) and Zimmermann (2001); the two locally defined classes, and the matching schedule classes of quasiactive and quasistable schedules, are the book's device for covering discontinuous resource-based objectives.

Setting

A project consists of activities V={0,1,…,n+1}V=\{0,1,\dots,n+1\}V={0,1,…,n+1}, n≥1n\ge 1n≥1, where 000 and n+1n+1n+1 are fictitious activities marking the project beginning and completion. Activity iii has an integer duration pip_ipi​ (p0=pn+1=0p_0=p_{n+1}=0p0​=pn+1​=0, pi>0p_i>0pi​>0 otherwise). The project network has an arc set EEE with integer weights δij\delta_{ij}δij​; a schedule is a vector S=(S0,…,Sn+1)S=(S_0,\dots,S_{n+1})S=(S0​,…,Sn+1​) of real start times with S0=0S_0=0S0​=0, S≥0S\ge 0S≥0, and it is time-feasible if Sj−Si≥δijS_j-S_i\ge\delta_{ij}Sj​−Si​≥δij​ for all ⟨i,j⟩∈E\langle i,j\rangle\in E⟨i,j⟩∈E. A maximum project duration dˉ∈N\bar d\in\mathbb Ndˉ∈N is prescribed through a backward arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩ of weight −dˉ-\bar d−dˉ, so Sn+1≤dˉS_{n+1}\le\bar dSn+1​≤dˉ. Each renewable resource kkk has capacity RkR_kRk​, activity iii uses rikr_{ik}rik​ units while in progress, and rk(S,t)r_k(S,t)rk​(S,t) is the total usage at time ttt. The feasible region S\mathcal SS consists of the time-feasible schedules with rk(S,t)≤Rkr_k(S,t)\le R_krk​(S,t)≤Rk​ for all kkk and ttt.

For an objective function f:R≥0n+2→Rf:\mathbb R^{n+2}_{\ge 0}\to\mathbb Rf:R≥0n+2​→R, problem PS∣temp,dˉ∣fPS|temp,\bar d|fPS∣temp,dˉ∣f asks for an optimal schedule: some S∈SS\in\mathcal SS∈S with f(S)≤f(S′)f(S)\le f(S')f(S)≤f(S′) for all S′∈SS'\in\mathcal SS′∈S.

A schedule induces the strict order O(S)={(i,j)∣i≠j, Sj≥Si+pi}O(S)=\{(i,j)\mid i\ne j,\ S_j\ge S_i+p_i\}O(S)={(i,j)∣i=j, Sj​≥Si​+pi​} of precedences it realizes. The equal-order set of SSS is

ST=(O(S))={S′ time-feasible∣Sj′≥Si′+pi ∀(i,j)∈O(S), O(S′)=O(S)},\mathcal S_T^{=}(O(S))=\{S'\text{ time-feasible}\mid S'_j\ge S'_i+p_i\ \forall (i,j)\in O(S),\ O(S')=O(S)\},ST=​(O(S))={S′ time-feasible∣Sj′​≥Si′​+pi​ ∀(i,j)∈O(S), O(S′)=O(S)},

a polytope with part of its boundary removed. The distinct equal-order sets partition S\mathcal SS into finitely many pieces.

Schedule classes are defined through shifts. A shift from a feasible SSS to a feasible S′≠SS'\ne SS′=S is order-preserving if O(S)⊆O(S′)O(S)\subseteq O(S')O(S)⊆O(S′); it is a left-shift if S′≤SS'\le SS′≤S. Two shifts from SSS to S′S'S′ and S′′S''S′′ are opposite if S′′−S=λ(S′−S)S''-S=\lambda(S'-S)S′′−S=λ(S′−S) with λ<0\lambda<0λ<0. A feasible schedule is active if no feasible left-shift exists, quasiactive if no order-preserving left-shift exists, stable if no pair of opposite shifts to feasible schedules exists, and quasistable if no pair of opposite order-preserving shifts exists.

Objective classes: fff is regular if S≤S′S\le S'S≤S′ implies f(S)≤f(S′)f(S)\le f(S')f(S)≤f(S′); quasiconcave on a set MMM if f(λS+(1−λ)S′)≥min⁡[f(S),f(S′)]f(\lambda S+(1-\lambda)S')\ge\min[f(S),f(S')]f(λS+(1−λ)S′)≥min[f(S),f(S′)] for S,S′∈MS,S'\in MS,S′∈M, λ∈[0,1]\lambda\in[0,1]λ∈[0,1]; lower semicontinuous if f(S)≤lim inf⁡S′→Sf(S′)f(S)\le\liminf_{S'\to S}f(S')f(S)≤liminfS′→S​f(S′) on R≥0n+2\mathbb R^{n+2}_{\ge 0}R≥0n+2​. Then fff is locally regular (class 6) if it is lower semicontinuous and regular on every equal-order set ST=(O(S))\mathcal S_T^{=}(O(S))ST=​(O(S)), S∈SS\in\mathcal SS∈S, and locally quasiconcave (class 7) if it is lower semicontinuous and quasiconcave on every such set.

Formalization targets

Goal: Theorem 3.3.13

For every locally quasiconcave fff,

S≠∅ ⟹ ∃ S quasistable with f(S)=min⁡S′∈Sf(S′).\mathcal S\ne\emptyset\ \Longrightarrow\ \exists\,S\ \text{quasistable with}\ f(S)=\min_{S'\in\mathcal S}f(S').S=∅ ⟹ ∃S quasistable with f(S)=S′∈Smin​f(S′).

Milestones

  • Class 1 (§3.3.2): every regular fff has an active optimal schedule when S≠∅\mathcal S\ne\emptysetS=∅.
  • Class 5 (§3.3.6): every quasiconcave fff has a stable optimal schedule when S≠∅\mathcal S\ne\emptysetS=∅.
  • Eq. (3.3.11): the equal-order sets form a finite partition of S\mathcal SS.
  • Propositions 3.3.5 and 3.3.6: the resource investment objective ∑kckmax⁡trk(S,t)\sum_k c_k\max_t r_k(S,t)∑k​ck​maxt​rk​(S,t) with ck≥0c_k\ge 0ck​≥0 is constant on each equal-order set and lower semicontinuous, hence locally regular.
  • Theorem 3.3.9: every locally regular fff has a quasiactive optimal schedule when S≠∅\mathcal S\ne\emptysetS=∅.

Significance

Quasiactive and quasistable schedules are finite in number: they are the minimal points and the vertices of the finitely many schedule polytopes. Theorem 3.3.13 therefore turns the minimization of any locally quasiconcave objective over a disconnected, non-convex feasible region into a finite search. Class 7 contains the resource levelling objectives ∑ck∑rkt2\sum c_k\sum r_{kt}^2∑ck​∑rkt2​ and ∑ck∑okt\sum c_k\sum o_{kt}∑ck​∑okt​, the total variation of the resource profiles, and the resource renting objective (Propositions 3.3.10 and 3.3.12, and Nübel 2001). The enumeration schemes and decision sets of §3.5–3.7 rest on this result, and Theorem 3.3.9 plays the same role for class 6 (resource investment, changeover times).

The results are proved in the book and the cited papers. As far as a search of the platform shows, none of them, and none of the schedule classes, has a machine-checked formalization; Mathlib supplies lower semicontinuity and quasiconcavity but nothing about schedules. The mission produces a checked version of the classification theorems in the book's exact generality: general time lags (cycles in the network allowed), real start times, and arbitrary objectives given only by their class.

Difficulty

The optimum need not exist a priori: objectives of classes 6 and 7 are discontinuous, and the feasible region is a finite union of polytopes that is in general disconnected. Existence of a minimizer needs compactness of S\mathcal SS (which depends on the deadline arc and the network's path structure) together with lower semicontinuity.

The main obstacle is that the objective is only controlled piecewise. Quasiconcavity holds on each equal-order set separately, and an equal-order set is not closed: a schedule polytope ST(O(S))\mathcal S_T(O(S))ST​(O(S)) also contains schedules inducing strictly larger orders, where the hypothesis on fff says nothing about its relation to the values on ST=(O(S))\mathcal S_T^{=}(O(S))ST=​(O(S)). The obvious argument, taking an optimal schedule and invoking quasiconcavity along the segment of a pair of opposite order-preserving shifts, only relates fff at points of one equal-order set, and it does not by itself produce a schedule that admits no such pair at all. The same issue arises for Theorem 3.3.9 with order-preserving left-shifts, which may cross from one equal-order set into another.

Formalization scope

Activities are Fin (n + 2), with 0 and Fin.last (n + 1) fictitious. Start times are real; objective functions are total functions (Fin (n + 2) → ℝ) → ℝ whose regularity, quasiconcavity and lower semicontinuity are required only on the nonnegative orthant (lower semicontinuity is Mathlib's LowerSemicontinuousOn on the orthant). The deadline Sn+1≤dˉS_{n+1}\le\bar dSn+1​≤dˉ is the network's backward arc, as in §3.1. The project structure records the book's standing property (p. 8) that from each node iii there is a path to n+1n+1n+1 of length at least pip_ipi​; this bounds every activity by dˉ\bar ddˉ. The resource constraints are imposed for all t≥0t\ge 0t≥0, which under that property is the book's 0≤t≤dˉ0\le t\le\bar d0≤t≤dˉ. The peak max⁡trk(S,t)\max_t r_k(S,t)maxt​rk​(S,t) in the resource investment objective is a supremum in N\mathbb NN over t≥0t\ge 0t≥0 of a nonempty finite set, hence attained.

"Optimal" always means minimizing fff over the whole feasible region S\mathcal SS, and the theorems quantify over every function in the class; a formalization with a fixed objective, or with optimality over a single polytope or a single equal-order set, would be a different and weaker statement. The schedule classes are defined through shifts, never as minimal or extreme points, so no statement is true by definition. The only hypothesis besides the class of fff is S≠∅\mathcal S\ne\emptysetS=∅.

The mission restates locally the project model, the induced orders and the shift classes also drafted by the companion missions on schedule classes of this series. Useful contributions beyond the milestones: compactness of S\mathcal SS and closedness of the schedule polytopes, the representation of S\mathcal SS as a finite union of feasible order polytopes, and the finiteness of the sets of quasiactive and quasistable schedules.

Selected references

  • K. Neumann, C. Schwindt, J. Zimmermann, Project Scheduling with Time Windows and Scarce Resources, 2nd ed., Springer, 2003, §3.3. doi:10.1007/978-3-540-24800-2
  • K. Neumann, H. Nübel, C. Schwindt, Active and stable project scheduling, Mathematical Methods of Operations Research 52 (2000), cited in the book as Neumann et al. (2000).
  • J. Zimmermann, Ablauforientiertes Projektmanagement: Modelle, Verfahren und Anwendungen, Gabler, 2001.
11 thms2 active usersReviewed
Graph TheoryLinear OptimizationOperations Research+1·Captain: mikedeng1

Project Scheduling with Time Windows and Scarce Resources VIII: A Vertex Schedule Maximizes the Net Present Value iff Its Spanning-Tree Subprojects Have the Right SignsTextbook

Motivation

Long-running projects such as construction, plant engineering or software development involve payments to and from the contractor at many points in time: disbursements when activities are carried out, progress payments when milestones are reached. When the planning horizon is long, money received later is worth less, and the natural financial objective is the net present value of all cash flows. Scheduling a project to maximize its net present value subject to minimum and maximum time lags was studied by Russell (1970) and Grinold (1972), and the problem is the prototype of a nonregular objective: delaying an activity can be profitable, because disbursements lose value when they are postponed.

This mission follows Chapter 3 of Neumann, Schwindt and Zimmermann, Project Scheduling with Time Windows and Scarce Resources (2nd ed., Springer 2003). The book shows that the net present value objective belongs to the class of binary-monotone objective functions (§3.3.5), and it uses this in §3.9.1 to give a combinatorial optimality criterion for the resource-free problem: a vertex schedule is optimal exactly when the subprojects cut off by the arcs of a spanning tree have net present values of the right sign (Proposition 3.9.2). That criterion drives the book's parametric analysis of the net present value as a function of the discount rate and the deadline.

Setting

A project consists of activities V={0,1,…,n+1}V=\{0,1,\dots,n+1\}V={0,1,…,n+1}, n≥1n\ge1n≥1, where 000 is the project beginning and n+1n+1n+1 the project completion. Activity iii has an integer duration pip_ipi​, with p0=pn+1=0p_0=p_{n+1}=0p0​=pn+1​=0 and pi>0p_i>0pi​>0 otherwise. Temporal constraints are the arcs of a project network N=⟨V,E;δ⟩N=\langle V,E;\delta\rangleN=⟨V,E;δ⟩: an arc ⟨i,j⟩\langle i,j\rangle⟨i,j⟩ with integer weight δij\delta_{ij}δij​ requires Sj−Si≥δijS_j-S_i\ge\delta_{ij}Sj​−Si​≥δij​ for the start times SiS_iSi​. A maximum project duration dˉ\bar ddˉ is the arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩ with weight −dˉ-\bar d−dˉ. The time-feasible region is

ST={S∈R≥0n+2∣S0=0, Sj−Si≥δij (⟨i,j⟩∈E)}.\mathcal S_T=\{S\in\mathbb R^{n+2}_{\ge0}\mid S_0=0,\ S_j-S_i\ge\delta_{ij}\ (\langle i,j\rangle\in E)\}.ST​={S∈R≥0n+2​∣S0​=0, Sj​−Si​≥δij​ (⟨i,j⟩∈E)}.

Let 0<β≤10<\beta\le10<β≤1 be the discount rate (β=1/(1+I)\beta=1/(1+I)β=1/(1+I) for an interest rate III) and ciF∈Rc_i^F\in\mathbb RciF​∈R the cash flow of activity iii, paid at its completion time Ci=Si+piC_i=S_i+p_iCi​=Si​+pi​. The problem (3.9.1) is

minimize f(S)=−∑i∈VciFβSi+pisubject to S∈ST,\text{minimize } f(S)=-\sum_{i\in V}c_i^F\beta^{S_i+p_i}\quad\text{subject to } S\in\mathcal S_T,minimize f(S)=−i∈V∑​ciF​βSi​+pi​subject to S∈ST​,

and a minimizer is a time-optimal schedule. A vertex of ST\mathcal S_TST​ is an extreme point. A spanning tree G=⟨V,EG⟩G=\langle V,E^G\rangleG=⟨V,EG⟩ is associated with SSS if EG⊆EE^G\subseteq EEG⊆E, EGE^GEG has n+1n+1n+1 arcs and a connected underlying undirected graph, and SSS is the unique solution of S0=0S_0=0S0​=0, Sj−Si=δijS_j-S_i=\delta_{ij}Sj​−Si​=δij​ for ⟨i,j⟩∈EG\langle i,j\rangle\in E^G⟨i,j⟩∈EG. Deleting a tree arc ⟨i,j⟩\langle i,j\rangle⟨i,j⟩ splits GGG into two subtrees; VijV_{ij}Vij​ is the node set of the one not containing 000. The arc is forward if the tree path from 000 passes it from iii to jjj and backward otherwise, and

npvij(S)=∑h∈VijchFβSh+phnpv^{ij}(S)=\sum_{h\in V_{ij}}c_h^F\beta^{S_h+p_h}npvij(S)=h∈Vij​∑​chF​βSh​+ph​

is the net present value of the subproject VijV_{ij}Vij​. Finally, fff is binary-monotone if it is monotone on every line {S+λz≥0∣λ∈R}\{S+\lambda z\ge0\mid\lambda\in\mathbb R\}{S+λz≥0∣λ∈R} with direction z∈{0,1}n+2z\in\{0,1\}^{n+2}z∈{0,1}n+2 (Definition 3.3.2).

Formalization targets

Goal: Proposition 3.9.2, pinned reading

Assume every node is reached from 000 by a path of nonnegative length (the standing convention of §1.2) and let SSS be a vertex of ST\mathcal S_TST​.

(sufficiency)G associated with S,  npvij(S)≥0 on forward arcs, npvij(S)≤0 on backward arcs ⟹ S time-optimal;\text{(sufficiency)}\quad G \text{ associated with } S,\ \ npv^{ij}(S)\ge0 \text{ on forward arcs},\ npv^{ij}(S)\le0 \text{ on backward arcs}\ \Longrightarrow\ S \text{ time-optimal};(sufficiency)G associated with S,  npvij(S)≥0 on forward arcs, npvij(S)≤0 on backward arcs ⟹ S time-optimal; (necessity, β<1)S time-optimal ⟹ ∃ G associated with S satisfying the sign conditions.\text{(necessity, } \beta<1)\quad S \text{ time-optimal}\ \Longrightarrow\ \exists\, G \text{ associated with } S \text{ satisfying the sign conditions}.(necessity, β<1)S time-optimal ⟹ ∃G associated with S satisfying the sign conditions.

The book states "if and only if … for each arc of the corresponding spanning tree", where the corresponding tree is chosen using optimality. The two directions above are the reading that makes the statement well defined: sufficiency for every associated tree, necessity for some associated tree.

Milestones

  1. §3.3.5: the net present value objective is binary-monotone and sum-separable.
  2. §3.9.1: if ST\mathcal S_TST​ is nonempty and bounded, some vertex of ST\mathcal S_TST​ is time-optimal.
  3. Proposition 3.2.16: every vertex of ST\mathcal S_TST​ has an associated spanning tree, an outtree rooted at 000 if the vertex is a minimal point.
  4. Proposition 3.5.4: a directed forest with at least one node has a source with at most one successor or a sink with exactly one predecessor.

Significance

Proposition 3.9.2 turns a nonconvex continuous optimization problem into a finite check on a spanning tree. Read as an economic statement, it says that at an optimal schedule no subproject with positive net present value can be started earlier and no subproject with negative net present value can be postponed. The book builds on it the parametric procedure of §3.9.1, which tracks the optimal tree as the discount rate or the deadline varies (Propositions 3.9.3 and 3.9.4), and the steepest descent method of §3.5.2 terminates exactly when the criterion holds.

The results are proved in the book, partly by reference to network optimization (Ahuja et al., 1993) and to Schwindt and Zimmermann (2001, 2002). None of them is formalized on the platform or, as far as is known, anywhere else. A formal proof would give the first machine-checked optimality certificate for a nonregular project scheduling objective, and the spanning-tree description of vertices (Proposition 3.2.16) is shared with Mission VI of this series.

Difficulty

The objective fff is neither convex nor concave when cash flows of both signs occur, so local optimality at a vertex does not imply global optimality by a convexity argument, and a first-order check along the edges of ST\mathcal S_TST​ is not obviously enough. The criterion is also not a statement about one tree: a degenerate vertex, where more than n+1n+1n+1 temporal constraints are binding, has several associated trees, and the sign conditions may hold on some and fail on others. Necessity therefore requires producing a suitable tree, not checking a given one. Finally, the combinatorial objects (the subtree VijV_{ij}Vij​, forward and backward orientation relative to the root) have to be connected to the geometry of ST\mathcal S_TST​ through Proposition 3.2.16, whose proof in the book is a citation.

Formalization scope

Activities are Fin (n + 2) with 0 the project beginning and Fin.last (n+1) the project completion; start times are real; durations are natural numbers and arc weights integers. The deadline is a structure field together with the backward arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩ of weight −dˉ-\bar d−dˉ. βx\beta^xβx is Real.rpow, and every statement assumes 0<β≤10<\beta\le10<β≤1 as the book does (p. 203). Vertices are Set.extremePoints ℝ. A spanning tree is a Finset of n+1n+1n+1 arcs whose SimpleGraph.fromRel is connected; VijV_{ij}Vij​ is the set of nodes not reachable from 000 once the arc is deleted.

Three readings are committed and disclosed in the item statements. Necessity is stated only for β<1\beta<1β<1: at β=1\beta=1β=1 the objective is constant, every schedule is optimal, and the sign conditions can fail on every tree. The standing convention of §1.2 (a path of nonnegative length from 000 to every node) is a hypothesis of Proposition 3.2.16 and of the goal; without it a vertex can be fixed by Si≥0S_i\ge0Si​≥0 rather than by arcs of NNN, and necessity fails. The existence of an optimal vertex assumes ST\mathcal S_TST​ nonempty and bounded, which the book asserts in §3.1. Chapter 3's resource constraints do not occur in this mission, which concerns PS∞∣temp,dˉ∣fPS\infty|temp,\bar d|fPS∞∣temp,dˉ∣f only.

The goal cannot be discharged by choosing the tree freely: associated trees must consist of arcs of NNN that are binding at SSS and determine SSS uniquely, and sufficiency must hold for every such tree. Contributions welcome beyond the milestones: a proof of Proposition 3.2.16 reusable by Mission VI, and a general lemma relating binding spanning trees of difference constraints to extreme points.

Selected references

  • K. Neumann, C. Schwindt, J. Zimmermann, Project Scheduling with Time Windows and Scarce Resources, 2nd ed., Springer, 2003, §3.1 (p. 203), §3.3.5 (pp. 224–225), §3.5.2 (p. 252), §3.9.1 (pp. 333–334). https://doi.org/10.1007/978-3-540-24800-2
  • A. H. Russell, "Cash flows in networks", Management Science 16 (1970), 357–373. https://doi.org/10.1287/mnsc.16.5.357
  • R. C. Grinold, "The payment scheduling problem", Naval Research Logistics Quarterly 19 (1972), 123–136.
  • C. Schwindt, J. Zimmermann, "A steepest ascent approach to maximizing the net present value of projects", Mathematical Methods of Operations Research 53 (2001), 435–450.
  • C. Schwindt, J. Zimmermann, "Parametrische Optimierung als Instrument zur Bewertung von Investitionsprojekten", Zeitschrift für Betriebswirtschaft 72 (2002), 593–617.
  • R. K. Ahuja, T. L. Magnanti, J. B. Orlin, Network Flows, Prentice Hall, 1993.
  • C. Berge, Graphs and Hypergraphs, North-Holland, Amsterdam, 1976.
9 thms2 active usersReviewed
Control TheoryOperations ResearchProbability+1·Captain: mikedeng1

Scheduling a Multi Class Queue with Many Exponential Servers: Asymptotic Optimality in Heavy Traffic: The HJB-Based Preemptive Policy Is Asymptotically Optimal Among Work-Conserving PoliciesResearch Paper

Motivation

Large call centers route several types of customers to a common pool of agents. When the pool is large and highly utilized, the relevant asymptotic regime is the quality-and-efficiency-driven (QED) or Halfin–Whitt regime (Halfin & Whitt 1981). The number of servers nnn grows while the offered load stays within O(n)O(\sqrt n)O(n​) of nnn. Waiting is then neither negligible nor overwhelming (Gans, Koole & Mandelbaum 2003).

Which class should a freed agent serve next? Exact optimization of a multi-class many-server queue with abandonment is intractable. The standard route is to solve a limiting diffusion control problem and translate its optimal control back into a policy for the queue. Atar, Mandelbaum and Reiman (Ann. Appl. Probab. 2004) carried this out for kkk customer classes, exponential service and abandonment, general renewal arrivals and general convex-type holding costs. They proved that the translated policy is asymptotically optimal. This mission formalizes that result for the preemptive policy.

Context:

  • Harrison & Zeevi (2004) studied the same multi-class many-server problem.
  • Bell & Williams (2001) proved asymptotic optimality of a threshold policy for a two-server system in conventional heavy traffic.
  • The present paper is the first to cover the QED regime with general costs and abandonment.

Setting

There are k≥1k\ge1k≥1 customer classes and nnn identical servers.

Primitives.

  • Arrivals. Class-iii customers arrive according to a renewal process AinA^n_iAin​ with interarrival times Uˇi(j)/λin\check U_i(j)/\lambda^n_iUˇi​(j)/λin​. Here the Uˇi(j)\check U_i(j)Uˇi​(j) are i.i.d., positive, of mean one and squared coefficient of variation CU,i2C^2_{U,i}CU,i2​.
  • Service. Service times are exponential with rate μin\mu^n_iμin​, represented by Poisson processes SinS^n_iSin​.
  • Abandonment. Waiting customers abandon at rate θin≥0\theta^n_i\ge0θin​≥0, represented by Poisson processes RinR^n_iRin​.

State. Xin(t)X^n_i(t)Xin​(t) is the number of class-iii customers in the system, Ψin(t)\Psi^n_i(t)Ψin​(t) the number in service and Φin=Xin−Ψin\Phi^n_i=X^n_i-\Psi^n_iΦin​=Xin​−Ψin​ the number waiting. The dynamics are

Xin(t)=Xi0,n+Ain(t)−Rin(∫0tΦin)−Sin(∫0tΨin),Ψn,Φn∈Z+k,∑iΨin≤n.X^n_i(t)=X^{0,n}_i+A^n_i(t)-R^n_i\Big(\int_0^t\Phi^n_i\Big)-S^n_i\Big(\int_0^t\Psi^n_i\Big),\qquad \Psi^n,\Phi^n\in\mathbb Z^k_+,\quad \textstyle\sum_i\Psi^n_i\le n .Xin​(t)=Xi0,n​+Ain​(t)−Rin​(∫0t​Φin​)−Sin​(∫0t​Ψin​),Ψn,Φn∈Z+k​,∑i​Ψin​≤n.

Policies.

  • A scheduling control policy (SCP) is the process Ψn\Psi^nΨn.
  • It is admissible if it does not anticipate the future beyond the time of the next arrival: past information is independent of future primitive increments.
  • It is work-conserving if no server idles while customers wait: (1⋅Xn−n)+=1⋅Φn(\mathbb 1\cdot X^n-n)^+=\mathbb 1\cdot\Phi^n(1⋅Xn−n)+=1⋅Φn.

Scaling and cost. In the QED scaling n−1λin→λin^{-1}\lambda^n_i\to\lambda_in−1λin​→λi​ with ∑iλi/μi=1\sum_i\lambda_i/\mu_i=1∑i​λi​/μi​=1. With ρi=λi/μi\rho_i=\lambda_i/\mu_iρi​=λi​/μi​ the centred processes are X^n=n−1/2(Xn−ρn)\hat X^n=n^{-1/2}(X^n-\rho n)X^n=n−1/2(Xn−ρn), Φ^n=n−1/2Φn\hat\Phi^n=n^{-1/2}\Phi^nΦ^n=n−1/2Φn and Ψ^n=n−1/2(Ψn−ρn)\hat\Psi^n=n^{-1/2}(\Psi^n-\rho n)Ψ^n=n−1/2(Ψn−ρn). The cost is

Cn=E∫0∞e−γtL~(Φ^n(t),Ψ^n(t)) dt.C^n=E\int_0^\infty e^{-\gamma t}\tilde L(\hat\Phi^n(t),\hat\Psi^n(t))\,dt .Cn=E∫0∞​e−γtL~(Φ^n(t),Ψ^n(t))dt.

The limiting control problem. It controls

X(t)=x+rW(t)+∫0tb(X(s),u(s)) ds,b(x,u)=ℓ+(μ−θ)(1⋅x)+u−μx,X(t)=x+rW(t)+\int_0^t b(X(s),u(s))\,ds,\qquad b(x,u)=\ell+(\mu-\theta)(\mathbb 1\cdot x)^+u-\mu x,X(t)=x+rW(t)+∫0t​b(X(s),u(s))ds,b(x,u)=ℓ+(μ−θ)(1⋅x)+u−μx,

where the control uuu takes values in the simplex Sk\mathbb S^kSk and WWW is a kkk-dimensional Brownian motion. The data are ri=(λiCU,i2+λi)1/2r_i=(\lambda_iC^2_{U,i}+\lambda_i)^{1/2}ri​=(λi​CU,i2​+λi​)1/2 and ℓi=λ^i−ρiμ^i\ell_i=\hat\lambda_i-\rho_i\hat\mu_iℓi​=λ^i​−ρi​μ^​i​. Its value V(x)V(x)V(x) is the infimum of E∫0∞e−γtL(X,u) dtE\int_0^\infty e^{-\gamma t}L(X,u)\,dtE∫0∞​e−γtL(X,u)dt, with L(x,u)=L~((1⋅x)+u,x−(1⋅x)+u)L(x,u)=\tilde L((\mathbb 1\cdot x)^+u,x-(\mathbb 1\cdot x)^+u)L(x,u)=L~((1⋅x)+u,x−(1⋅x)+u).

HJB equation and the proposed policy. The HJB equation is 12∑iri2∂iif+H(x,Df)−γf=0\tfrac12\sum_ir_i^2\partial_{ii}f+H(x,Df)-\gamma f=021​∑i​ri2​∂ii​f+H(x,Df)−γf=0 with H(x,p)=inf⁡u∈Sk[b(x,u)⋅p+L(x,u)]H(x,p)=\inf_{u\in\mathbb S^k}[b(x,u)\cdot p+L(x,u)]H(x,p)=infu∈Sk​[b(x,u)⋅p+L(x,u)]. Let hhh be a measurable selection of its minimizers. The proposed preemptive policy (P-SCP) sets the queue vector to Θ[(1⋅Xn−n)+h(X^n)]\Theta[(\mathbb 1\cdot X^n-n)^+h(\hat X^n)]Θ[(1⋅Xn−n)+h(X^n)], an integer rounding, and falls back to a static priority rule when that is infeasible.

Formalization targets

Goal: Theorem 2(i)

For a Cpol2C^2_{\mathrm{pol}}Cpol2​ solution fff of the HJB equation, a measurable minimizer selection hhh, and initial states with X^0,n→x\hat X^{0,n}\to xX^0,n→x:

lim⁡n→∞E∫0∞e−γtL~(Φ^tn,∗,Ψ^tn,∗) dt  ≤  lim inf⁡n→∞E∫0∞e−γtL~(Φ^tn,Ψ^tn) dt\lim_{n\to\infty}E\int_0^\infty e^{-\gamma t}\tilde L(\hat\Phi^{n,*}_t,\hat\Psi^{n,*}_t)\,dt\;\le\;\liminf_{n\to\infty}E\int_0^\infty e^{-\gamma t}\tilde L(\hat\Phi^n_t,\hat\Psi^n_t)\,dtn→∞lim​E∫0∞​e−γtL~(Φ^tn,∗​,Ψ^tn,∗​)dt≤n→∞liminf​E∫0∞​e−γtL~(Φ^tn​,Ψ^tn​)dt

This holds for every sequence of work-conserving admissible SCPs, and the left-hand limit exists and is finite. No constants are hard-coded.

Milestones

The milestones follow the proof:

  • on the diffusion side, Proposition 2 (well-posedness), Proposition 4 (stability and moment bounds), Proposition 5(i)–(ii) (growth and continuity of VVV) and Theorem 3 (VVV is the unique Cpol2C^2_{\mathrm{pol}}Cpol2​ HJB solution, and an optimal Markov policy exists);
  • on the queueing side, Proposition 1 (feedback rules give admissible SCPs), Lemmas 2–3 (moment bounds), Lemma 4(i)–(ii) (FCLT for the primitives and the fluid limit (Ψˉn,Φˉn)⇒(ρ,0)(\bar\Psi^n,\bar\Phi^n)\Rightarrow(\rho,0)(Ψˉn,Φˉn)⇒(ρ,0)), and Theorem 4(i)–(ii): lim inf⁡≥V(x)\liminf\ge V(x)liminf≥V(x) always, and lim sup⁡≤V(x)\limsup\le V(x)limsup≤V(x) under condition (49).

Significance

The result. Theorem 2(i) justifies using the diffusion control problem as a design tool for multi-class many-server systems. The policy is explicit given hhh, and it is optimal in the limit against all non-anticipating work-conserving policies, including those that use the full history and the time of the next arrival. The proof also identifies the limit cost with V(x)V(x)V(x).

Formalizing it. The result is proved on paper, with some steps (Proposition 1, the principle of optimality, the time-change and martingale limit theorems) given as sketches or citations. No part of it is machine-checked. A formalization requires:

  • a counting-process model of the queue;
  • a careful definition of non-anticipation;
  • a pathwise controlled SDE;
  • classical solvability of a semilinear elliptic HJB equation on Rk\mathbb R^kRk;
  • a weak-convergence argument in Skorokhod space.

Each of these is reusable well beyond this paper.

Difficulty

The obvious argument would show that X^n\hat X^nX^n converges to the controlled diffusion and pass the costs to the limit. This fails for two reasons:

  • the comparison class contains arbitrary non-Markov, history-dependent policies, so the queue does not converge to a single controlled diffusion;
  • the optimal selector hhh is in general discontinuous (for linear costs it is), so the proposed policy is not a continuous function of the state.

The proof instead compares every policy with the HJB solution through Itô's formula on the prelimit processes. This needs:

  • uniform moment bounds;
  • tightness of the integral processes;
  • the convergence of stochastic integrals of Kurtz and Protter;
  • and, for the proposed policy, the fact that the rounding Θ\ThetaΘ and the priority fallback perturb the minimizer by O(n−1/2)O(n^{-1/2})O(n−1/2).

Existence of a classical HJB solution on all of Rk\mathbb R^kRk, with only Hölder-continuous costs and polynomial growth, rests on a bounded-domain existence theorem for fully nonlinear elliptic equations.

Formalization scope

The Lean development commits to the following conventions.

  • Indexing and norms. Classes are Fin k with k≥1k\ge1k≥1; paper class iii is index i−1i-1i−1, so "class kkk" (highest priority, rounding remainder of Θ\ThetaΘ) is the last index. Vectors are Fin k → ℝ and ∥⋅∥\|\cdot\|∥⋅∥ is the paper's ℓ1\ell^1ℓ1 norm; the paper's ∣⋅∣|\cdot|∣⋅∣ on vectors is read the same way.
  • Probability space and paths. All systems share one complete probability space. Time is real and every condition is for t≥0t\ge0t≥0. The paper's "without loss" path regularity (finite arrival counts, Poisson paths Z+\mathbb Z_+Z+​-valued, nondecreasing and càdlàg) holds for every ω\omegaω.
  • Poisson processes are defined by independent Poisson increments; rate 000 gives the zero process.
  • Policies. A policy is a pair of real processes (Ψn,Xn)(\Psi^n,X^n)(Ψn,Xn) with integer values. Admissibility is Definition 2 verbatim, with the future σ\sigmaσ-field built from the next arrival time τin(t)\tau^n_i(t)τin​(t). Work conservation is (18).
  • Costs and value are lower Lebesgue integrals in [0,∞][0,\infty][0,∞], and lim⁡\limlim/lim inf⁡\liminfliminf are taken there. The integrands are nonnegative under work conservation.
  • Admissible systems range over sample spaces Ω : Type (universe 0). "Complete filtered probability space" means PPP complete with all null sets in F0\mathcal F_0F0​. Brownian motion is Mathlib's IsBrownianReal per coordinate, with independence and the (Ft)(\mathcal F_t)(Ft​)-Brownian property stated explicitly. VVV is the infimum over systems and their controlled processes.
  • Discount rate. γ>0\gamma>0γ>0 is a hypothesis; the paper leaves it implicit.
  • Initial states are integer vectors X0,n∈Z+kX^{0,n}\in\mathbb Z^k_+X0,n∈Z+k​ with n−1/2(X0,n−ρn)→xn^{-1/2}(X^{0,n}-\rho n)\to xn−1/2(X0,n−ρn)→x. The literal "X^0,n∈n−1/2Zk\hat X^{0,n}\in n^{-1/2}\mathbb Z^kX^0,n∈n−1/2Zk" would require ρin∈Z\rho_in\in\mathbb Zρi​n∈Z. Assumption 1(ii) is not imposed: each policy chooses its own initial split.
  • Lemma 3 is stated for all nnn beyond a threshold that depends on the sequence, with constants c,mˉc,\bar mc,mˉ chosen before xxx and the sequence. The printed all-nnn bound with ccc independent of xxx fails when the early terms X^0,n\hat X^{0,n}X^0,n are large.
  • Weak convergence to a continuous limit uses the coupling form CouplingConverges of the published BellWilliams2001.ThresholdPolicy.Paths; convergence to a deterministic limit is UocInProb.

The goal hypothesizes fff and hhh with the pointwise identity b(x,h(x))⋅Df(x)+L(x,h(x))=H(x,Df(x))b(x,h(x))\cdot Df(x)+L(x,h(x))=H(x,Df(x))b(x,h(x))⋅Df(x)+L(x,h(x))=H(x,Df(x)) for all xxx. An arbitrary "optimal Markov control policy" may differ from a minimizer selection on the Lebesgue-null lattice where X^n\hat X^nX^n lives, and that formalization would make the goal false. Restricting the comparators to feedback, Markov or nonpreemptive policies, fixing kkk, dropping abandonment, specializing to Poisson arrivals or linear costs, or imposing a common initial split would each trivialize or weaken the statement and is ruled out.

Not formalized:

  • Lemma 4(iii) (tightness);
  • Lemma 5 (Kurtz–Protter, which needs semimartingale theory absent from Mathlib);
  • Lemma 6 (convergence of Stieltjes integrals at limit points);
  • Proposition 5(iii);
  • the nonpreemptive results, Theorem 2(ii)–(iii).

Contributions are welcome on any milestone, and especially on infrastructure: Poisson and renewal processes, functional central limit theorems in Skorokhod space, classical solvability of elliptic HJB equations, and measurable selection of minimizers.

Selected references

  • R. Atar, A. Mandelbaum, M. I. Reiman, Scheduling a multi class queue with many exponential servers: asymptotic optimality in heavy traffic, Ann. Appl. Probab. 14(3), 2004. https://arxiv.org/abs/math/0407058
  • S. Halfin, W. Whitt, Heavy-traffic limits for queues with many exponential servers, Oper. Res. 29(3), 1981. https://doi.org/10.1287/opre.29.3.567
  • N. Gans, G. Koole, A. Mandelbaum, Telephone call centers: tutorial, review, and research prospects, Manuf. Serv. Oper. Manag. 5(2), 2003. https://doi.org/10.1287/msom.5.2.79.16071
  • J. M. Harrison, A. Zeevi, Dynamic scheduling of a multiclass queue in the Halfin–Whitt heavy traffic regime, Oper. Res. 52(2), 2004. https://doi.org/10.1287/opre.1040.0109
  • S. L. Bell, R. J. Williams, Dynamic scheduling of a system with two parallel servers in heavy traffic with resource pooling: asymptotic optimality of a threshold policy, Ann. Appl. Probab. 11(3), 2001. https://doi.org/10.1214/aoap/1015345343
  • T. G. Kurtz, P. Protter, Weak limit theorems for stochastic integrals and stochastic differential equations, Ann. Probab. 19(3), 1991. https://doi.org/10.1214/aop/1176990334
16 thms1 active userReviewed
Graph TheoryOperations ResearchTheoretical Computer Science·Captain: mikedeng1

On the Approximability of Single-Machine Scheduling with Precedence Constraints 5: An r-Approximate Vertex Cover of the Variable-Cost Graph Yields an (r + ε)-Approximate Vertex Cover of GResearch Paper

Why the variable cost matters

The problem 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_jC_j1∣prec∣∑wj​Cj​ asks for an order in which to process jobs on one machine, respecting precedence constraints, so as to minimize the weighted sum of completion times. It is strongly NP-hard, and for decades the best approximation ratio known has been 222, achieved by several unrelated algorithms (LP relaxations, Sidney decompositions, primal–dual methods).

Correa and Schulz (2005) and Ambühl and Mastrolilli (2009) showed that the problem is a special case of weighted vertex cover: its objective splits into a fixed cost, the same for every feasible solution, and a variable cost, which equals the weight of a vertex cover in an auxiliary graph GPSG^S_{\mathbf P}GPS​. Approximating vertex cover in GPSG^S_{\mathbf P}GPS​ within a factor α\alphaα therefore approximates the scheduling problem within α\alphaα. Uhan observed that the classical 2-approximations owe their guarantee to the fixed cost and can be arbitrarily bad on the variable cost alone.

Section 8 of Ambühl, Mastrolilli, Mutsanas and Svensson, Math. Oper. Res. 36(4) (2011) (DOI), proves the converse: approximating the variable cost is as hard as approximating vertex cover itself. A better-than-2 algorithm for 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_jC_j1∣prec∣∑wj​Cj​ must therefore either exploit the fixed cost or improve on the best known approximation for vertex cover, a long-standing open question.

Setting

Scheduling instance. A finite set NNN of jobs, a partial order P=(N,P)\mathbf P = (N,P)P=(N,P) (reflexive; (i,j)∈P(i,j) \in P(i,j)∈P, i≠ji \ne ji=j, means iii precedes jjj), processing times pj≥0p_j \ge 0pj​≥0 and weights wj≥0w_j \ge 0wj​≥0.

Incomparable pairs. Jobs x,yx,yx,y are incomparable, x∥yx \parallel yx∥y, if neither (x,y)(x,y)(x,y) nor (y,x)(y,x)(y,x) lies in PPP. The set inc⁡(P)\operatorname{inc}(\mathbf P)inc(P) consists of the ordered pairs (x,y)(x,y)(x,y) with x∥yx \parallel yx∥y.

The vertex cover graph GPSG^S_{\mathbf P}GPS​. One node per incomparable pair (i,j)(i,j)(i,j), of weight w(i,j)=piwjw_{(i,j)} = p_i w_jw(i,j)​=pi​wj​. Distinct nodes (i,j)(i,j)(i,j) and (k,ℓ)(k,\ell)(k,ℓ) are adjacent when, in one of the two orders, j=kj=kj=k and i=ℓi=\elli=ℓ, or j=kj=kj=k and (i,ℓ)∈P(i,\ell)\in P(i,ℓ)∈P, or (i,ℓ),(k,j)∈P(i,\ell),(k,j) \in P(i,ℓ),(k,j)∈P. For a set CCC of nodes, w(C)=∑u∈Cwuw(C) = \sum_{u\in C} w_uw(C)=∑u∈C​wu​; for a vertex cover CCC this is the variable cost, and τw(GPS)\tau_w(G^S_{\mathbf P})τw​(GPS​) is its minimum over all vertex covers.

The instance S(G,k)S(G,k)S(G,k). Given a graph G=(V,E)G=(V,E)G=(V,E) with V={v1,…,vn}V = \{v_1,\dots,v_n\}V={v1​,…,vn​} and k>0k > 0k>0, the instance has jobs vi′v'_ivi′​ (processing time k−ik^{-i}k−i, weight 000) and vi′′v''_ivi′′​ (processing time 000, weight kik^{i}ki), and precedence constraints vi′<vj′′v'_i < v''_jvi′​<vj′′​ and vj′<vi′′v'_j < v''_ivj′​<vi′′​ for each edge {vi,vj}∈E\{v_i,v_j\} \in E{vi​,vj​}∈E, plus vi′<vj′′v'_i < v''_jvi′​<vj′′​ for all i<ji<ji<j. The nodes (vi′,vi′′)(v'_i, v''_i)(vi′​,vi′′​) of GPSG^S_{\mathbf P}GPS​ have weight 111 and are called heavy; all others are light. For a set CCC of nodes, CG={vi:(vi′,vi′′)∈C}C_G = \{v_i : (v'_i,v''_i)\in C\}CG​={vi​:(vi′​,vi′′​)∈C}. The vertex cover number of GGG is τ(G)\tau(G)τ(G).

Formalization targets

Goal: Theorem 8.1

For every graph GGG on nnn vertices, every r≥1r \ge 1r≥1, ε>0\varepsilon>0ε>0 and every k≥1k \ge 1k≥1 with k>n2r/εk > n^2r/\varepsilonk>n2r/ε: if CCC is a vertex cover of GPSG^S_{\mathbf P}GPS​ for S=S(G,k)S = S(G,k)S=S(G,k) with w(C)≤r τw(GPS)w(C) \le r\,\tau_w(G^S_{\mathbf P})w(C)≤rτw​(GPS​), then CGC_GCG​ is a vertex cover of GGG,

∣CG∣≤r(τ(G)+n2k),|C_G| \le r\Bigl(\tau(G) + \frac{n^2}{k}\Bigr),∣CG​∣≤r(τ(G)+kn2​),

and, when E≠∅E \ne \emptysetE=∅,

∣CG∣≤r(1+n2k)τ(G)<(r+ε) τ(G).|C_G| \le r\Bigl(1+\frac{n^2}{k}\Bigr)\tau(G) < (r+\varepsilon)\,\tau(G).∣CG​∣≤r(1+kn2​)τ(G)<(r+ε)τ(G).

The paper words the theorem as "approximating the variable cost of 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_jC_j1∣prec∣∑wj​Cj​ is as hard as approximating vertex cover"; the statement above is the mathematical content its proof establishes.

Milestones (§8, p. 664)

  1. In GPSG^S_{\mathbf P}GPS​, every heavy node has weight 111, every light node has weight at most 1/k1/k1/k, and the light nodes have total weight at most n2/kn^2/kn2/k (for k≥1k \ge 1k≥1).
  2. Heavy nodes (vi′,vi′′)(v'_i,v''_i)(vi′​,vi′′​) and (vj′,vj′′)(v'_j,v''_j)(vj′​,vj′′​) are adjacent if and only if {vi,vj}∈E\{v_i,v_j\}\in E{vi​,vj​}∈E; for k>1k>1k>1 the subgraph induced by the weight-1 nodes is isomorphic to GGG via (vi′,vi′′)↦vi(v'_i,v''_i) \mapsto v_i(vi′​,vi′′​)↦vi​.

Significance

The result. Theorem 8.1 is one half of an equivalence: by Theorem 2.1 (Correa–Schulz, Ambühl–Mastrolilli), minimizing the variable cost is a special case of weighted vertex cover; by Theorem 8.1, it is also as hard to approximate. Any hardness of approximation for vertex cover (NP-hardness of factor 1.361.361.36 by Dinur and Safra; factor 2−δ2-\delta2−δ under the unique games conjecture by Khot and Regev) transfers to the variable cost. It also explains why the known 2-approximations must rely on the fixed cost, and it frames the later result of Bansal and Khot that the full objective is hard to approximate within 2−δ2-\delta2−δ under a variant of the unique games conjecture.

Formalizing it. The theorem is proved in the paper; no machine-checked version exists. The mission produces a checked account of the reduction: the vertex cover graph of an arbitrary precedence-constrained instance, the adjacency-poset instance built from a graph, and the quantitative transfer of approximation ratios. The definition of GPSG^S_{\mathbf P}GPS​ is shared with the other missions of this series.

Difficulty

The construction is short; the care is in the bookkeeping. One must check that the precedence relation is a partial order, determine exactly which ordered pairs are incomparable, verify that two heavy nodes are adjacent only through the third clause of the adjacency rule and only when the corresponding vertices are adjacent in GGG, and bound the weights of all remaining nodes, including the many nodes of weight 000. The transfer then compares an approximate cover of GPSG^S_{\mathbf P}GPS​ with an optimal one whose heavy part comes from an optimal cover of GGG; the additive error n2/kn^2/kn2/k must be converted into a multiplicative one, which requires τ(G)≥1\tau(G)\ge 1τ(G)≥1.

A first reading of the page suggests that GPSG^S_{\mathbf P}GPS​ has at most n2n^2n2 nodes; it does not. The pairs (vi′,vj′)(v'_i,v'_j)(vi′​,vj′​) and (vi′′,vj′′)(v''_i,v''_j)(vi′′​,vj′′​) with i≠ji\ne ji=j are incomparable nodes of weight 000, so there can be up to 4n2−2n4n^2-2n4n2−2n nodes. Only nodes of positive weight are few.

Formalization scope

  • Model. Jobs form a finite type; precedence constraints are an explicit reflexive partial-order relation P : N → N → Prop. Processing times and weights are nonnegative reals. The vertex cover graph is a SimpleGraph on the subtype of incomparable ordered pairs, using the symmetric closure of the printed adjacency rule without loops. Vertex covers are Mathlib's SimpleGraph.IsVertexCover; τ(G)\tau(G)τ(G) is Mathlib's vertexCoverNum, finite for a finite graph and converted with toNat; τw(GPS)\tau_w(G^S_{\mathbf P})τw​(GPS​) is a minimum over finite vertex covers.
  • The instance. The graph is a SimpleGraph (Fin n); i : Fin n stands for vi+1v_{i+1}vi+1​, so exponents are i+1i+1i+1 and the order i<ji<ji<j is that of Fin n. Jobs are Fin n ⊕ Fin n (v′v'v′ left, v′′v''v′′ right). The parameter kkk is in R≥0\mathbb R_{\ge 0}R≥0​.
  • Added hypotheses. k≥1k \ge 1k≥1, implicit in the page ("k>n2r/εk > n^2r/\varepsilonk>n2r/ε" does not imply it when ε\varepsilonε is large, and for k<1k<1k<1 the light nodes outweigh the heavy ones). The isomorphism with GGG is stated for k>1k > 1k>1, since at k=1k=1k=1 some light nodes also have weight 111. The multiplicative bound requires E≠∅E \ne \emptysetE=∅; the additive bound holds for every graph.
  • Not formalized. The phrases "approximation algorithm", "polynomial time" and "as hard as"; the passage from vertex covers of GPSG^S_{\mathbf P}GPS​ to schedules (Theorem 2.1, cited from Correa–Schulz and Ambühl–Mastrolilli); the fixed cost. What is stated instead is the explicit map C↦CGC \mapsto C_GC↦CG​ and the ratio it achieves, r(1+n2/k)<r+εr(1+n^2/k) < r+\varepsilonr(1+n2/k)<r+ε. The false count "at most n2n^2n2 vertices" is not stated.
  • Ruled out. The goal is not a statement about an arbitrary graph or an assumed cover of GGG: it concerns the specific instance S(G,k)S(G,k)S(G,k) and every CCC that is an rrr-approximate vertex cover of its graph, and the fact that CGC_GCG​ covers GGG is a conclusion, not a hypothesis.
  • Welcome contributions. Proofs of the two milestones and of the goal; general lemmas on GPSG^S_{\mathbf P}GPS​ (weights of vertex covers, behaviour under induced subgraphs) are reusable across the series.

Selected references

  • C. Ambühl, M. Mastrolilli, N. Mutsanas, O. Svensson, On the Approximability of Single-Machine Scheduling with Precedence Constraints, Mathematics of Operations Research 36(4):653–669, 2011. https://doi.org/10.1287/moor.1110.0512
  • J. R. Correa, A. S. Schulz, Single-Machine Scheduling with Precedence Constraints, Mathematics of Operations Research 30(4):1005–1021, 2005. https://doi.org/10.1287/moor.1050.0158
  • C. Ambühl, M. Mastrolilli, Single Machine Precedence Constrained Scheduling Is a Vertex Cover Problem, Algorithmica 53(4):488–503, 2009. https://doi.org/10.1007/s00453-008-9251-6
  • I. Dinur, S. Safra, On the Hardness of Approximating Minimum Vertex Cover, Annals of Mathematics 162(1):439–485, 2005. https://doi.org/10.4007/annals.2005.162.439
  • S. Khot, O. Regev, Vertex Cover Might Be Hard to Approximate to within 2 − ε, Journal of Computer and System Sciences 74(3):335–349, 2008. https://doi.org/10.1016/j.jcss.2007.06.019
  • N. Bansal, S. Khot, Optimal Long Code Test with One Free Bit, FOCS 2009, 453–462. https://doi.org/10.1109/FOCS.2009.23
5 thms1 active userReviewed
CombinatoricsGraph TheoryOperations Research+1·Captain: mikedeng1

On the Approximability of Single-Machine Scheduling with Precedence Constraints 4: Vertex Cover in Connected Graphs of Degree ≤ 3 Reduces to Weighted Vertex Cover for Interval-Order InstancesResearch Paper

Motivation

In the single-machine scheduling problem 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_jC_j1∣prec∣∑wj​Cj​, a set NNN of nnn jobs, each with a processing time pj≥0p_j\ge 0pj​≥0 and a weight wj≥0w_j\ge 0wj​≥0, is processed on one machine without interruption, subject to precedence constraints given by a partial order PPP on NNN. The aim is to minimize the weighted sum of completion times ∑jwjCj\sum_j w_jC_j∑j​wj​Cj​. The problem is strongly NP-hard for general precedence constraints (Lawler 1978; Lenstra and Rinnooy Kan 1978), and its approximability was a recurring open question in scheduling theory (Schuurman and Woeginger 1999).

A line of work by Chudak and Hochbaum, Correa and Schulz, and Ambühl and Mastrolilli showed that the problem is a special case of minimum weighted vertex cover in a graph GPSG^S_PGPS​ built from the instance. Many problems on partial orders become polynomial when the order is an interval order, so it is natural to ask whether this one does too. Section 7 of Ambühl, Mastrolilli, Mutsanas and Svensson (Math. Oper. Res. 2011) answers no: the problem stays NP-hard on interval orders. The proof is a reduction from vertex cover in connected graphs of maximum degree 3. This mission formalizes the correctness of that reduction.

Setting

A poset P=(N,P)P=(N,P)P=(N,P) is read as a reflexive relation: (x,y)∈P(x,y)\in P(x,y)∈P means x≤yx\le yx≤y. Jobs x,yx,yx,y are incomparable if neither (x,y)(x,y)(x,y) nor (y,x)(y,x)(y,x) is in PPP, and inc⁡(P)\operatorname{inc}(P)inc(P) is the set of ordered incomparable pairs. The vertex cover graph GPSG^S_PGPS​ has one node (i,j)(i,j)(i,j) for each incomparable pair, weighted piwjp_iw_jpi​wj​. Two nodes (i,j)(i,j)(i,j) and (k,ℓ)(k,\ell)(k,ℓ) are adjacent if j=kj=kj=k and i=ℓi=\elli=ℓ, or j=kj=kj=k and (i,ℓ)∈P(i,\ell)\in P(i,ℓ)∈P, or (i,ℓ)∈P(i,\ell)\in P(i,ℓ)∈P and (k,j)∈P(k,j)\in P(k,j)∈P. Write w(CI)w(C_I)w(CI​) for the minimum weight of a vertex cover of GISG^S_IGIS​.

A poset is an interval order if each element xxx can be assigned a closed real interval [ax,bx][a_x,b_x][ax​,bx​] such that x<yx<yx<y if and only if bx<ayb_x<a_ybx​<ay​.

The reduction starts from a graph G=(V,E)G=(V,E)G=(V,E) with vertices v1,…,vNv_1,\dots,v_Nv1​,…,vN​ and a spanning tree T=(V,ET)T=(V,E_T)T=(V,ET​) rooted at v1v_1v1​, numbered so that each parent comes before its children. The paper uses a breadth-first search tree.

  • Stage 1. The graph G′G'G′ is built from TTT. Each viv_ivi​ gets a pendant path vi−u1i−u2iv_i - u^i_1 - u^i_2vi​−u1i​−u2i​. Each non-tree edge {vi,vj}∈E∖ET\{v_i,v_j\}\in E\setminus E_T{vi​,vj​}∈E∖ET​ with i<ji<ji<j gets the path vi−e1ij−e2ij−u2jv_i - e^{ij}_1 - e^{ij}_2 - u^j_2vi​−e1ij​−e2ij​−u2j​. The non-tree edges themselves are not edges of G′G'G′.
  • Stage 2. The scheduling instance SSS has jobs s0s_0s0​, s1,…,sNs_1,\dots,s_Ns1​,…,sN​, m1,…,mNm_1,\dots,m_Nm1​,…,mN​, e1,…,eNe_1,\dots,e_Ne1​,…,eN​, and bijb_{ij}bij​ for each non-tree edge. Their intervals, processing times and weights are given in a table on p. 662. For example, sjs_jsj​ has interval [i,j][i,j][i,j], processing time 1/kj1/k^j1/kj and weight kik^iki, where viv_ivi​ is the parent of vjv_jvj​. The precedence constraints III are the interval order of these intervals. With nnn the number of jobs, the parameter is k=n2+1k=n^2+1k=n2+1.
  • The set DDD. It is {(s0,s1)}∪{(si,sj):vi parent of vj}∪{(si,mi),(mi,ei)}∪{(si,bij),(bij,mj)}\{(s_0,s_1)\}\cup\{(s_i,s_j): v_i \text{ parent of } v_j\}\cup\{(s_i,m_i),(m_i,e_i)\}\cup\{(s_i,b_{ij}),(b_{ij},m_j)\}{(s0​,s1​)}∪{(si​,sj​):vi​ parent of vj​}∪{(si​,mi​),(mi​,ei​)}∪{(si​,bij​),(bij​,mj​)}. The graph GI′G'_IGI′​ is the subgraph of GISG^S_IGIS​ induced by DDD.

Formalization targets

Goal: Theorem 7.1 (p. 661)

For every connected graph GGG of maximum degree at most 333, every parent-first spanning tree TTT and every m∈Nm\in\mathbb Nm∈N, the precedence constraints III of SSS form an interval order, and

G has a vertex cover of size≤m  ⟺  ⌊w(CI)⌋≤m+∣V∣+∣E∖ET∣.G \text{ has a vertex cover of size} \le m \iff \lfloor w(C_I)\rfloor \le m + |V| + |E\setminus E_T|.G has a vertex cover of size≤m⟺⌊w(CI​)⌋≤m+∣V∣+∣E∖ET​∣.

Milestones, in the order the proof uses them

  • Claim 1 (p. 662): τ(G′)=τ(G)+∣V∣+∣E∖ET∣\tau(G') = \tau(G)+|V|+|E\setminus E_T|τ(G′)=τ(G)+∣V∣+∣E∖ET​∣, where τ\tauτ is the vertex cover number.
  • Remark 7.1 (p. 662): for jobs with intervals [a,b][a,b][a,b] and [c,d][c,d][c,d] and a≤da\le da≤d, pi≤1/k⌈b⌉p_i\le 1/k^{\lceil b\rceil}pi​≤1/k⌈b⌉ and wj≤k⌈c⌉w_j\le k^{\lceil c\rceil}wj​≤k⌈c⌉. On incomparable pairs piwj∈{1}∪[0,1/k]p_iw_j\in\{1\}\cup[0,1/k]pi​wj​∈{1}∪[0,1/k]. Moreover, piwj≥kp_iw_j\ge kpi​wj​≥k forces b<cb<cb<c, and piwj=1p_iw_j=1pi​wj​=1 forces ⌈b⌉=⌈c⌉\lceil b\rceil=\lceil c\rceil⌈b⌉=⌈c⌉.
  • Claim 2 (p. 663): an incomparable pair (i,j)(i,j)(i,j) has piwj=1p_iw_j=1pi​wj​=1 if it is in DDD, and piwj≤1/kp_iw_j\le 1/kpi​wj​≤1/k otherwise.
  • Claim 3 (p. 663): GI′≅G′G'_I\cong G'GI′​≅G′.
  • §7, p. 664: with k=n2+1k=n^2+1k=n2+1, ∑(i,j)∈inc⁡(I)∖Dpiwj<1\sum_{(i,j)\in\operatorname{inc}(I)\setminus D}p_iw_j<1∑(i,j)∈inc(I)∖D​pi​wj​<1, and hence w(CI′)=⌊w(CI)⌋w(C'_I)=\lfloor w(C_I)\rfloorw(CI′​)=⌊w(CI​)⌋.

Significance

The result. Interval orders are a standard tractable class: several scheduling and order-theoretic problems that are hard in general become polynomial on them (Papadimitriou and Yannakakis 1979). Theorem 7.1 puts 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_jC_j1∣prec∣∑wj​Cj​ outside this pattern. Section 6 of the same paper shows that the problem nonetheless has a 3/23/23/2-approximation on interval orders, so hardness and approximability are separated on this class. The paper also remarks that the proof makes weighted vertex cover NP-hard to approximate within some factor r>1r>1r>1 on the graphs GISG^S_IGIS​ arising from interval orders.

Formalizing it. The theorem is proved in the paper. Nothing in this mission is open, and none of it has been machine-checked before. The work splits into the following parts:

  • a gadget argument on unweighted vertex covers (Claim 1, after Alimonti and Kann);
  • an exact case analysis of incomparable pairs in a concrete interval order (Remark 7.1, Claim 2);
  • a graph isomorphism (Claim 3);
  • a rounding argument that links weighted and unweighted optima.

The definitions of GPSG^S_PGPS​ and of minimum-weight vertex covers are shared with the other missions of this series.

Difficulty

The construction is explicit, and each step is elementary. The work is in the bookkeeping. Claim 2 requires classifying every incomparable pair of jobs, including pairs of different kinds such as (bij,sℓ)(b_{ij}, s_\ell)(bij​,sℓ​), by comparing ceilings of interval endpoints. Half-integer endpoints (mim_imi​, bijb_{ij}bij​) are exactly what separates weight-one pairs from comparable ones. Claim 3 requires checking adjacency in GISG^S_IGIS​ for all pairs of nodes of DDD in both directions. The paper writes out two cases in each direction and calls the rest similar.

Claim 1 has a direction that is not simply local. A vertex cover of G′G'G′ that misses both endpoints of a non-tree edge has to be repaired by swapping gadget vertices, and the repair must be repeated without increasing the size.

A natural first idea is to treat the light nodes (weight at most 1/k1/k1/k) as negligible one at a time. This does not suffice: the argument needs their total weight to stay below 111, which is what forces kkk to grow with n2n^2n2.

Formalization scope

  • Graph and tree. GGG is a SimpleGraph (Fin N); vertex vi+1v_{i+1}vi+1​ is i, and the root is index 0. The tree is a TreeLayout: a parent function returning none exactly at the root, with each parent of smaller index and adjacent in GGG. The statements hold for every such layout. This is stronger than the paper's breadth-first tree, and the proof uses only "parent before child".
  • Hypotheses of the goal. Connectivity and the degree bound ((G.neighborSet v).ncard ≤ 3) are kept as in the paper. They matter only for the NP-completeness of the source problem.
  • Jobs. The jobs form an inductive type with one constructor per row of the table. Their order is a PartialOrder instance: x≤yx\le yx≤y iff x=yx=yx=y or bx<ayb_x<a_ybx​<ay​. Processing times and weights are real numbers. Section 1 of the paper asks for nonnegative integers, but the instance uses 1/kj1/k^j1/kj and the formalization follows the instance as printed.
  • Constants. The constants are explicit: k=n2+1k=n^2+1k=n2+1 with nnn the cardinality of the job type, and c=∣V∣+∣E∖ET∣c=|V|+|E\setminus E_T|c=∣V∣+∣E∖ET​∣. Remark 7.1 and Claim 2 are stated for every real k>1k>1k>1.
  • Optimum values. w(CI)w(C_I)w(CI​) is a minimum over the finite family of vertex covers. Unweighted cover numbers are Mathlib's SimpleGraph.vertexCoverNum. The floor is Nat.floor, which agrees with the integer floor because w(CI)≥0w(C_I)\ge 0w(CI​)≥0.
  • Not formalized. The goal's wording ("NP-hard") is not formalized. Neither are the NP-completeness of degree-3 vertex cover (Garey, Johnson and Stockmeyer), the polynomial size of the construction, or Theorem 2.1 (cited), which turns a vertex cover of GISG^S_IGIS​ into a schedule. What is stated is the correctness of the reduction: the instance has interval-order constraints, and its optimum decides the vertex cover question.
  • No trivialization. The instance SSS is built from GGG and TTT exactly as in the table. The goal quantifies over all graphs and layouts, never over an instance SSS assumed to have the properties.
  • Contributions. Contributions are welcome on any milestone. Claims 1 and 3 are independent of the weights, and Claim 2 is independent of the graph theory.

Selected references

  • C. Ambühl, M. Mastrolilli, N. Mutsanas, O. Svensson, On the Approximability of Single-Machine Scheduling with Precedence Constraints, Math. Oper. Res. 36(4):653–669, 2011. https://doi.org/10.1287/moor.1110.0512
  • C. Ambühl, M. Mastrolilli, Single machine precedence constrained scheduling is a vertex cover problem, Algorithmica 53(4):488–503, 2009. https://doi.org/10.1007/s00453-008-9251-1
  • J. R. Correa, A. S. Schulz, Single machine scheduling with precedence constraints, Math. Oper. Res. 30(4):1005–1021, 2005. https://doi.org/10.1287/moor.1050.0158
  • P. Alimonti, V. Kann, Some APX-completeness results for cubic graphs, Theoret. Comput. Sci. 237(1–2):123–134, 2000. https://doi.org/10.1016/S0304-3975(98)00158-3
  • M. R. Garey, D. S. Johnson, L. Stockmeyer, Some simplified NP-complete graph problems, Theoret. Comput. Sci. 1(3):237–267, 1976. https://doi.org/10.1016/0304-3975(76)90059-1
  • C. H. Papadimitriou, M. Yannakakis, Scheduling interval-ordered tasks, SIAM J. Comput. 8(3):405–409, 1979. https://doi.org/10.1137/0208031
10 thms1 active userReviewed
PreviousNext

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me