Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Combinatorics

630 missions · 352 completed

The mathematics of finite and discrete structures — counting the arrangements of a set, deciding when a configuration meeting prescribed constraints can exist, and characterizing the patterns such structures are forced to contain. It encompasses enumerative and extremal combinatorics, graph theory, design theory, and additive combinatorics, with deep ties to algebra, probability, and computer science.

Missions

Open278Completed352All630
🏆Completed
Linear OptimizationOperations Research·Captain: mikedeng1

Some Polyhedra Related to Combinatorial Problems I: A Minimizing Vertex of the Corner Polyhedron Gives an Optimal Integer Solution When the Right-Hand Side Lies Deep in the Basis ConeResearch Paper

Motivation

An integer program in standard form, max⁡c⋅x\max c\cdot xmaxc⋅x subject to Ax=bAx=bAx=b, x≥0x\ge 0x≥0, xxx integer, is easy to state and hard to solve. Its linear programming relaxation is solved by the simplex method, which ends at an optimal basis BBB. The difference between the two problems is the integrality of xxx, and R. E. Gomory's 1969 paper Some polyhedra related to combinatorial problems (Linear Algebra Appl. 2, 451–558, doi:10.1016/0024-3795(69)90017-2) measures that difference with a finite Abelian group. Dropping only the nonnegativity of the basic variables leaves a problem whose feasible set has a simple structure: its convex hull is the corner polyhedron, and its integer points are the solutions of an equation in the group Zm/BZm\mathbb Z^m/B\mathbb Z^mZm/BZm.

The first section of the paper uses this to explain when the integer program reduces to the group problem. Its asymptotic theorem says that if bbb lies far enough inside the cone where BBB is LP-optimal, then an optimal solution of the group problem gives an optimal integer solution. Gomory first proved this in 1965 (On the relation between integer and noninteger solutions to linear programs, PNAS 53, 260–265, doi:10.1073/pnas.53.2.260) by a counting argument on group solutions. The 1969 paper proves it again through the geometry of the corner polyhedron. The vertices of that polyhedron are irreducible, and an irreducible point is small. This is the version formalized here.

Setting

Let A=(B,N)A=(B,N)A=(B,N) be an integer m×(m+n)m\times(m+n)m×(m+n) matrix, with BBB a nonsingular m×mm\times mm×m matrix (the first mmm columns) and NNN the m×nm\times nm×n matrix of the remaining columns N1,…,NnN_1,\dots,N_nN1​,…,Nn​. As on p. 452 of the paper, AAA contains an m×mm\times mm×m unit matrix. Let b∈Zmb\in\mathbb Z^mb∈Zm and c=(cB,cN)∈Rm+nc=(c_B,c_N)\in\mathbb R^{m+n}c=(cB​,cN​)∈Rm+n. The integer program (2) is

max⁡ c⋅xsubject toAx=b, x≥0, x integer.\max\ c\cdot x\quad\text{subject to}\quad Ax=b,\ x\ge 0,\ x\ \text{integer}.max c⋅xsubject toAx=b, x≥0, x integer.

The relative prices are cN∗=−cN+cBB−1Nc_N^*=-c_N+c_BB^{-1}NcN∗​=−cN​+cB​B−1N, with components cm+i∗c^*_{m+i}cm+i∗​. BBB is an optimal LP basis when B−1b≥0B^{-1}b\ge0B−1b≥0 and every cm+i∗≥0c^*_{m+i}\ge0cm+i∗​≥0.

The group is G=M(I)/M(B)\mathcal G=M(I)/M(B)G=M(I)/M(B), where M(I)=ZmM(I)=\mathbb Z^mM(I)=Zm and M(B)M(B)M(B) is the lattice of integer combinations of the columns of BBB. Let fff be the quotient map. G\mathcal GG is finite of order D=∣det⁡B∣D=|\det B|D=∣detB∣. Let N\mathcal NN be the set of nonzero elements among fN1,…,fNnfN_1,\dots,fN_nfN1​,…,fNn​, and put g0=fbg_0=fbg0​=fb. The group equation is

∑g∈Nt(g)⋅g=g0,t(g)∈Z≥0,\sum_{g\in\mathcal N}t(g)\cdot g=g_0,\qquad t(g)\in\mathbb Z_{\ge0},g∈N∑​t(g)⋅g=g0​,t(g)∈Z≥0​,

and P(G,N,g0)P(\mathcal G,\mathcal N,g_0)P(G,N,g0​) is the convex hull of its solutions in RN\mathbb R^{\mathcal N}RN. A nonnegative integer ttt is irreducible if any two integer vectors s,rs,rs,r with 0≤s,r≤t0\le s,r\le t0≤s,r≤t and ∑s(g)⋅g=∑r(g)⋅g\sum s(g)\cdot g=\sum r(g)\cdot g∑s(g)⋅g=∑r(g)⋅g are equal.

The group costs are c∗(g)=min⁡i: fNi=gcm+i∗c^*(g)=\min_{i:\,fN_i=g}c^*_{m+i}c∗(g)=mini:fNi​=g​cm+i∗​. Problem (8) minimizes ∑gc∗(g)t(g)\sum_g c^*(g)t(g)∑g​c∗(g)t(g) over the solutions of the group equation. A corresponding vertex of a point t∗t^*t∗ (Remark 1) is a nonnegative integer xN∗x_N^*xN∗​ with three properties: t∗(g)=∑i: fNi=gxm+i∗t^*(g)=\sum_{i:\,fN_i=g}x^*_{m+i}t∗(g)=∑i:fNi​=g​xm+i∗​; at most one positive xm+i∗x^*_{m+i}xm+i∗​ per group element; and xm+i∗=0x^*_{m+i}=0xm+i∗​=0 when fNi=0ˉfN_i=\bar0fNi​=0ˉ. It uses only least cost columns when xm+i∗>0x^*_{m+i}>0xm+i∗​>0 only if cm+i∗=c∗(fNi)c^*_{m+i}=c^*(fN_i)cm+i∗​=c∗(fNi​).

The cone KB={y∈Rm:B−1y≥0}K_B=\{y\in\mathbb R^m:B^{-1}y\ge0\}KB​={y∈Rm:B−1y≥0} is where BBB is LP-optimal. KB(d)K_B(d)KB​(d) is the set of its points at Euclidean distance at least ddd from its frontier, and lmax⁡l_{\max}lmax​ is the largest Euclidean length of a column NiN_iNi​.

Formalization targets

Goal: THEOREM 4 (p. 462)

If b∈KB(lmax⁡(D−1))b\in K_B\bigl(l_{\max}(D-1)\bigr)b∈KB​(lmax​(D−1)), t∗t^*t∗ is a vertex of P(G,N,fb)P(\mathcal G,\mathcal N,fb)P(G,N,fb) minimizing (8), and xN∗x_N^*xN∗​ is a corresponding vertex using only least cost columns, then

x∗=(B−1(b−NxN∗), xN∗) is an optimal integer solution of (2).x^*=\bigl(B^{-1}(b-Nx_N^*),\ x_N^*\bigr)\ \text{is an optimal integer solution of (2)}.x∗=(B−1(b−NxN∗​), xN∗​) is an optimal integer solution of (2).

Milestones

  1. THEOREM 1 (p. 459): an irreducible ttt satisfies ∏g∈N(1+t(g))≤∣G∣\prod_{g\in\mathcal N}(1+t(g))\le|\mathcal G|∏g∈N​(1+t(g))≤∣G∣.
  2. p. 459: ∣M(I)/M(B)∣=∣det⁡B∣|M(I)/M(B)|=|\det B|∣M(I)/M(B)∣=∣detB∣.
  3. THEOREM 2 (p. 460): every vertex of P(G,N,g0)P(\mathcal G,\mathcal N,g_0)P(G,N,g0​) is an irreducible integer point.
  4. Proof of THEOREM 4 (p. 462): an irreducible ttt has ∑gt(g)≤∣G∣−1\sum_g t(g)\le|\mathcal G|-1∑g​t(g)≤∣G∣−1.
  5. THEOREM 3 (p. 461): a corresponding vertex of a minimizing vertex is optimal for (2) whenever B−1(b−NxN∗)≥0B^{-1}(b-Nx_N^*)\ge0B−1(b−NxN∗​)≥0.
  6. Proof of THEOREM 4: ∥NxN∗∥≤lmax⁡(D−1)\|Nx_N^*\|\le l_{\max}(D-1)∥NxN∗​∥≤lmax​(D−1).
  7. Proof of THEOREM 4: y∈KB(d)y\in K_B(d)y∈KB​(d) and ∥v∥≤d\|v\|\le d∥v∥≤d imply y−v∈KBy-v\in K_By−v∈KB​.
  8. Display (10), p. 464: the optimal x∗x^*x∗ of THEOREM 4 satisfies ∏i=1n(1+xm+i∗)≤D\prod_{i=1}^n(1+x^*_{m+i})\le D∏i=1n​(1+xm+i∗​)≤D.

Significance

THEOREM 4 explains why integer programs with large right-hand sides are often easy. For bbb outside a band of width lmax⁡(D−1)l_{\max}(D-1)lmax​(D−1) along the boundary of the cone KBK_BKB​, the integer optimum is the LP solution B−1bB^{-1}bB−1b plus a correction that depends only on the class of bbb modulo the lattice of BBB. That correction is periodic in bbb and can be tabulated once for the DDD group elements. The vertex form adds a structural statement: the correction can be read off a vertex of the corner polyhedron, so the faces and vertices of corner polyhedra (the subject of the rest of the paper) bear directly on integer programming. Display (10) shows such optimal solutions are small: their nonbasic parts satisfy ∏(1+xm+i)≤D\prod(1+x_{m+i})\le D∏(1+xm+i​)≤D.

All results here are proved in the paper. None of them is formalized yet. The 1965 version of the asymptotic theorem, posed as THEOREM 1 of the mission on Gomory (1965), is in the same family of statements. Its hypothesis is that the group solution is optimal; here the hypothesis is that it is a vertex of P(G,N,g0)P(\mathcal G,\mathcal N,g_0)P(G,N,g0​). Its counting lemma (a short optimal group solution) and its bound ∥Ny∥≤(D−1)l\|Ny\|\le(D-1)l∥Ny∥≤(D−1)l have counterparts in milestones 4 and 6.

Difficulty

An optimal solution of the group problem need not be short. Group solutions can be made longer without changing their cost (add a zero-cost cycle), and a long nonbasic part xNx_NxN​ can push b−NxNb-Nx_Nb−NxN​ out of the cone KBK_BKB​, making xBx_BxB​ negative. The goal must therefore use the vertex hypothesis. Vertices are irreducible (THEOREM 2), and irreducibility bounds their size through a pigeonhole count in the finite group (THEOREM 1). The formal work is spread across three settings: the convex geometry of the hull of an infinite integer set, the arithmetic of the quotient Zm/BZm\mathbb Z^m/B\mathbb Z^mZm/BZm and its order, and the linear algebra connecting B−1B^{-1}B−1, the cone KBK_BKB​ and its frontier.

Formalization scope

Vectors on N\mathcal NN are functions on the subtype of a Finset of the group. Integer solutions are N\mathbb NN-valued, and P(G,N,g0)P(\mathcal G,\mathcal N,g_0)P(G,N,g0​) is convexHull ℝ of their real images. THEOREMS 1 and 2 are stated for an arbitrary finite Abelian group and any finite N\mathcal NN not containing 000, as in the paper. The integer program uses the basis as the first mmm columns, Fin m ⊕ Fin n. B−1B^{-1}B−1 is the real matrix inverse, and every statement assumes det⁡B≠0\det B\ne0detB=0. Optimality of xxx means feasibility plus c⋅y≤c⋅xc\cdot y\le c\cdot xc⋅y≤c⋅x for every feasible integer yyy; no supremum is taken. The loose phrases of the paper are read as follows.

  • "Vertex" means an extreme point.
  • "Every vertex is irreducible" means every vertex is the image of an integer solution, and that solution is irreducible.
  • "Optimal linear programming basis" means B−1b≥0B^{-1}b\ge0B−1b≥0 and c∗≥0c^*\ge0c∗≥0.
  • "Corresponding vertex x∗x^*x∗ of Px(B,N,b)P_x(B,N,b)Px​(B,N,b)" means conditions (i)–(iii) of Remark 1, which the paper calls easily verified, taken as the definition.
  • "Minimizing (8)" means minimizing over the integer solutions of the group equation.
  • "Euclidean distance from the frontier" is explicit: ∑kvk2\sqrt{\sum_k v_k^2}∑k​vk2​​, frontier in the usual topology, so KB(0)=KBK_B(0)=K_BKB​(0)=KB​.

The printed slip on p. 462 (∑(1+t∗(g))\sum(1+t^*(g))∑(1+t∗(g)) for THEOREM 1's product) is corrected in the Lean. Dropping the vertex hypothesis from the goal is not a valid formalization (the statement becomes false in general), and neither is defining irreducibility with real r,sr,sr,s (THEOREM 1 then fails). The existence of a minimizing vertex, which the paper takes for granted, is not posed.

A complete development needs the finiteness and order of Zm/BZm\mathbb Z^m/B\mathbb Z^mZm/BZm (Mathlib's AddSubgroup.index_eq_natAbs_det), extreme points of convex hulls (extremePoints_convexHull_subset), and an argument that a closed cone with yyy deep inside it contains the ball around yyy. The group-polyhedron definitions and THEOREMS 1–2 apply to any finite Abelian group and are reusable beyond this mission. Proofs of any milestone are welcome independently.

Selected references

  • R. E. Gomory, Some polyhedra related to combinatorial problems, Linear Algebra and Its Applications 2 (1969), 451–558. https://doi.org/10.1016/0024-3795(69)90017-2
  • R. E. Gomory, On the relation between integer and noninteger solutions to linear programs, Proc. Nat. Acad. Sci. USA 53 (1965), 260–265. https://doi.org/10.1073/pnas.53.2.260
  • R. E. Gomory, Faces of an integer polyhedron, Proc. Nat. Acad. Sci. USA 57 (1967), 16–18. https://doi.org/10.1073/pnas.57.1.16
  • B. L. van der Waerden, Modern Algebra, English translation of the 2nd German edition, Ungar, New York, 1949–1950 (computation of fff, cited p. 456).
11 thms2 active usersReviewed
🏆Completed
Operations ResearchTheoretical Computer Science·Captain: mikedeng1

Bounds on Multiprocessing Timing Anomalies 1: Changing the List, Relaxing the Order, Shortening the Tasks and Using n′ Processors Multiplies the Finishing Time by at Most 1 + (n − 1)/n′Research Paper

Motivation

Scheduling work on several identical processors creates a plausible expectation: completing tasks sooner, removing constraints, or adding a processor should not delay the overall finish. A fixed priority-list rule defeats that expectation. When a task becomes ready earlier, it can occupy a processor that would otherwise have run another task; that decision changes later availability. Graham gives explicit schedules in which changing only the list raises the finishing time from 12 to 14, removing two precedence constraints raises it to 16, shortening every task raises it to 13, or adding one processor raises it to 15. These examples establish that the effect is a property of the scheduling rule, rather than of a change in the total set of tasks. See Graham, §2, pp. 417–419.

This mission concerns the upper limit of that effect. An operations researcher comparing a list schedule before and after a change in available processors or task data needs a bound that survives all four changes at once. The result is also a compact benchmark for formalizing algorithms whose behavior depends on changing availability: the model must describe when tasks become ready, when processors are idle, and which ready task a processor takes next. Graham's Theorem 1, proved in the 1969 paper, supplies the bound.

Setting

There are rrr tasks T1,…,TrT_1,\ldots,T_rT1​,…,Tr​ and nnn identical processors. Each task TjT_jTj​ has a positive processing time μ(Tj)\mu(T_j)μ(Tj​). A strict precedence order ≺\prec≺ specifies tasks that must be completed first: Ti≺TjT_i\prec T_jTi​≺Tj​ means TjT_jTj​ cannot start before TiT_iTi​ ends. A priority list LLL orders all tasks once each. It may place a successor before a predecessor, because the processor skips any task that is not ready. These conventions are those of Graham's system in §2, p. 416.

A schedule gives each task a processor and a start time Sj≥0S_j\ge0Sj​≥0. Once started, the task occupies that processor continuously on [Sj,Sj+μ(Tj))[S_j,S_j+\mu(T_j))[Sj​,Sj​+μ(Tj​)). Two tasks on one processor cannot overlap, while tasks may meet at an endpoint. Precedence means Si+μ(Ti)≤SjS_i+\mu(T_i)\le S_jSi​+μ(Ti​)≤Sj​ whenever Ti≺TjT_i\prec T_jTi​≺Tj​. The finishing time ω\omegaω is the largest completion time Sj+μ(Tj)S_j+\mu(T_j)Sj​+μ(Tj​). The paper's list rule starts the first currently ready task in LLL whenever a processor can take a task. A ready task that has not started cannot coexist with an idle processor. If one task starts while another ready task waits, the started task is earlier in LLL.

The same task set is run a second time, with processing times μ′\mu'μ′, order ≺′\prec'≺′, list L′L'L′, processor count n′n'n′, and finishing time ω′\omega'ω′. The data satisfy μ′(Tj)≤μ(Tj)\mu'(T_j)\le\mu(T_j)μ′(Tj​)≤μ(Tj​) for every jjj and ≺′⊆≺\prec'\subseteq\prec≺′⊆≺. Thus the second run may have shorter tasks and fewer precedence relations; both lists may differ arbitrarily. Both runs obey the same list-scheduling rule Graham, §3, p. 419.

Formalization targets

General anomaly bound

For r,n,n′≥1r,n,n'\ge1r,n,n′≥1, positive processing times in both runs, and the two list schedules just described, the target is Graham's Theorem 1:

ω′ω≤1+n−1n′.\frac{\omega'}{\omega}\le 1+\frac{n-1}{n'}.ωω′​≤1+n′n−1​.

The numerator and denominator refer to the actual runs, with potentially different lists and processor counts. The conclusion does not require the new schedule to finish earlier. Its factor permits the anomalies shown in §2 while preventing an arbitrarily large increase for fixed n,n′n,n'n,n′. The supporting targets follow the proof's numbered displays: a precedence chain covering idle times in the second run (1), an idle-time estimate for that chain (2), a comparison of its length with the first run (3)–(4), and the work-volume bound used in (5) Graham, pp. 419–420.

Significance

The theorem places a quantitative limit on the damage caused by a list rule after simultaneous changes in task durations, precedence, ordering, and capacity. When n′=nn'=nn′=n, the factor is 2−1/n2-1/n2−1/n; when the first run has one processor, it is 111. The statement applies to every finite task set and every pair of lists, so it can be reused when a scheduling method produces a list without additional structural guarantees. Graham states that the bound is best possible, citing earlier examples; this mission formalizes the inequality in this paper and does not claim to formalize those external examples Graham, p. 420.

The result is proved in the source paper. The work here is to give its objects precise Lean definitions and to supply machine-checked proofs of the bound and its supporting claims. A formal development needs to account for continuous real start times, processor assignments, half-open execution intervals, and the behavior of the list rule when several tasks or processors become available together. The resulting model of a nonpreemptive list schedule and the basic volume and chain estimates can serve later formalizations of parallel scheduling. The present proposal contains statements and definitions; its theorem bodies are open proof obligations.

Difficulty

Total work alone does not control the anomaly: processors can be idle while precedence keeps available work from starting, and changing precedence can alter which tasks occupy the processors later. A direct comparison of corresponding task completion times between the two runs also fails in Graham's §2 examples. The central difficulty is expressing the second run's idle periods through work on a single precedence chain and then relating that chain to the first run's finishing time. In Lean, the chain and the time intervals must be represented without silently changing endpoint conventions or treating a ready but unstarted task as running.

Formalization scope

Tasks are Fin r and processors are Fin n, with zero-based Lean indices corresponding to Graham's one-based subscripts. Processing times and start times are real. A priority list is a bijection between positions and tasks, so it contains every task exactly once and need not respect precedence. The order is strict. Feasible schedules are nonpreemptive and use half-open running intervals. The finishing time is a finite maximum of completion times; for no tasks it is defined as zero, while the ratio theorem assumes r>0r>0r>0 so its denominator is positive. Both processor counts are positive. The strict positivity of both time functions is inherited from the paper's μ:T→(0,∞)\mu:T\to(0,\infty)μ:T→(0,∞) convention.

The model captures Graham's rule by no unforced idleness and list priority. It does not encode the smaller-index tie rule for identical processors, since the tie affects only processor labels and not start times or ω\omegaω. The covering-chain milestone represents Graham's set BBB as times before ω′\omega'ω′ when not all processors are busy; all-idle times before the finish cannot occur under the list rule. The idle-time milestone writes the total duration of empty tasks as n′ω′−∑jμ′(Tj)n'\omega'-\sum_j\mu'(T_j)n′ω′−∑j​μ′(Tj​) using display (5). Natural-number subtraction is not used in the constants: n−1n-1n−1 and n′−1n'-1n′−1 are real differences. The companion existence statement and a concrete two-task sanity check rule out a vacuous list-schedule predicate. Contributions can address that existence theorem, the chain construction, interval accounting, the volume bound, and the final ratio.

Selected references

  • R. L. Graham, Bounds on Multiprocessing Timing Anomalies, SIAM Journal on Applied Mathematics 17(2), 416–429, 1969. DOI: 10.1137/0117039.
7 thms2 active usersReviewed
🏆Completed
Operations ResearchOptimization·Captain: mikedeng1

Optimization and Approximation in Deterministic Sequencing and Scheduling: A Survey 1: The Optimal Two-Machine Open Shop Makespan Is max{T₁, T₂, maxⱼ(aⱼ + bⱼ)}Research Paper

Motivation

Shop scheduling asks how to sequence jobs that each need processing on several machines. In an open shop, a job's operations may be executed in any order, as in testing stations, repair bays, or classroom and examination timetables, where the order in which a candidate visits the stations is irrelevant. The objective studied here is the makespan Cmax⁡C_{\max}Cmax​, the time at which the last operation finishes.

The survey of Graham, Lawler, Lenstra and Rinnooy Kan (Ann. Discrete Math. 5, 1979) introduced the three-field notation α∣β∣γ\alpha|\beta|\gammaα∣β∣γ that the scheduling literature still uses, and classified the complexity of the problems it names. For the open shop, its §5.2.1 presents a simplified exposition of the result of Gonzalez and Sahni (J. ACM 23, 1976): with two machines and no preemption, the obvious lower bound on the makespan is always achieved. The same page records that the three-machine case O3∥Cmax⁡O3\|C_{\max}O3∥Cmax​ is binary NP-hard, so two machines is exactly where the problem is easy.

Timeline. Gonzalez and Sahni (1976) gave the linear-time algorithm for O2∥Cmax⁡O2\|C_{\max}O2∥Cmax​, proved NP-hardness for O3∥Cmax⁡O3\|C_{\max}O3∥Cmax​, and gave a polynomial algorithm for the preemptive problem O∣pmtn∣Cmax⁡O|pmtn|C_{\max}O∣pmtn∣Cmax​. Graham et al. (1979, §5.2.1) gave the shorter construction formalized here, and observed (§5.2.2) that it implies preemption brings no advantage for m=2m = 2m=2. Lenstra (cited as forthcoming in the survey) showed O2∣rj∣Cmax⁡O2|r_j|C_{\max}O2∣rj​∣Cmax​, O2∣tree∣Cmax⁡O2|tree|C_{\max}O2∣tree∣Cmax​ and O∥Cmax⁡O\|C_{\max}O∥Cmax​ unary NP-hard.

Setting

There are nnn jobs J1,…,JnJ_1, \dots, J_nJ1​,…,Jn​ and two machines M1M_1M1​, M2M_2M2​. Job JjJ_jJj​ has an operation on M1M_1M1​ of length aj≥0a_j \ge 0aj​≥0 and an operation on M2M_2M2​ of length bj≥0b_j \ge 0bj​≥0. There is no preemption, and every job is available at time 000. A schedule assigns start times s1(j)s_1(j)s1​(j) and s2(j)s_2(j)s2​(j) to the two operations of JjJ_jJj​, which then occupy [s1(j),s1(j)+aj)[s_1(j), s_1(j)+a_j)[s1​(j),s1​(j)+aj​) on M1M_1M1​ and [s2(j),s2(j)+bj)[s_2(j), s_2(j)+b_j)[s2​(j),s2​(j)+bj​) on M2M_2M2​.

A schedule is feasible if

  1. all start times are nonnegative;
  2. each machine processes at most one job at a time: the intervals of distinct jobs on the same machine do not overlap;
  3. each job is processed on at most one machine at a time: the two intervals of the same job do not overlap, in either order.

Write T1=∑jajT_1 = \sum_j a_jT1​=∑j​aj​ and T2=∑jbjT_2 = \sum_j b_jT2​=∑j​bj​ for the two machine loads. The survey's construction uses the sets

A={Jj∣aj≥bj},B={Jj∣aj<bj},A = \{J_j \mid a_j \ge b_j\}, \qquad B = \{J_j \mid a_j < b_j\},A={Jj​∣aj​≥bj​},B={Jj​∣aj​<bj​},

two distinct jobs JrJ_rJr​, JlJ_lJl​ with ar≥max⁡Jj∈Abja_r \ge \max_{J_j \in A} b_jar​≥maxJj​∈A​bj​ and bl≥max⁡Jj∈Bajb_l \ge \max_{J_j \in B} a_jbl​≥maxJj​∈B​aj​, and A′=A−{Jr,Jl}A' = A - \{J_r, J_l\}A′=A−{Jr​,Jl​}, B′=B−{Jr,Jl}B' = B - \{J_r, J_l\}B′=B−{Jr​,Jl​}.

Formalization targets

Goal: the optimal makespan

Cmax⁡∗=max⁡{T1, T2, max⁡j (aj+bj)},C^*_{\max} = \max\Big\{T_1,\ T_2,\ \max_j\,(a_j + b_j)\Big\},Cmax∗​=max{T1​, T2​, jmax​(aj​+bj​)},

and the optimum is attained. Formally, for every T≥0T \ge 0T≥0: a feasible schedule completing every operation by TTT exists if and only if T1≤TT_1 \le TT1​≤T, T2≤TT_2 \le TT2​≤T and aj+bj≤Ta_j + b_j \le Taj​+bj​≤T for all jjj. The goal mentions neither AAA, BBB, JrJ_rJr​, JlJ_lJl​ nor the case analysis; those are the milestones.

Milestones (in the order of the argument)

  1. Two distinct jobs JrJ_rJr​, JlJ_lJl​ with the required bounds exist when n≥2n \ge 2n≥2.
  2. Fig. 5.1: the blocks B′∪{Jl}B' \cup \{J_l\}B′∪{Jl​} and A′∪{Jr}A' \cup \{J_r\}A′∪{Jr​}, with A′A'A′ and B′B'B′ in arbitrary order, have feasible staircase schedules without idle time.
  3. Fig. 5.2: if T1−al≥T2−brT_1 - a_l \ge T_2 - b_rT1​−al​≥T2​−br​, the blocks combine into a feasible schedule of all jobs ending by T1+brT_1 + b_rT1​+br​.
  4. Case (1): if moreover ar≤T2−bra_r \le T_2 - b_rar​≤T2​−br​, some feasible schedule has length at most max⁡{T1,T2}\max\{T_1, T_2\}max{T1​,T2​}.
  5. Case (2): if moreover ar>T2−bra_r > T_2 - b_rar​>T2​−br​, some feasible schedule has length at most max⁡{T1,ar+br}\max\{T_1, a_r + b_r\}max{T1​,ar​+br​}.
  6. The symmetric case T1−al<T2−brT_1 - a_l < T_2 - b_rT1​−al​<T2​−br​: some feasible schedule has length at most max⁡{T1,T2,al+bl}\max\{T_1, T_2, a_l + b_l\}max{T1​,T2​,al​+bl​}.
  7. The lower bound: every feasible schedule has Cmax⁡≥max⁡{T1,T2,max⁡j(aj+bj)}C_{\max} \ge \max\{T_1, T_2, \max_j(a_j + b_j)\}Cmax​≥max{T1​,T2​,maxj​(aj​+bj​)}.

Significance

The result. The theorem gives a closed form for the optimal makespan of a two-machine open shop, together with a linear-time construction of an optimal schedule. The survey uses it immediately: since the bound is also a lower bound for preemptive schedules, O2∣pmtn∣Cmax⁡O2|pmtn|C_{\max}O2∣pmtn∣Cmax​ is solved by the same schedules, so preemption gives no advantage on two machines (§5.2.2). It is the standard example of a shop problem whose trivial lower bound is tight, and the contrast with the binary NP-hard O3∥Cmax⁡O3\|C_{\max}O3∥Cmax​ marks the complexity boundary for nonpreemptive open shops.

Formalizing it. The result is classical and proved on paper; no machine-checked proof is known to exist on this platform. The mission produces a reusable model of nonpreemptive two-machine open-shop schedules with real processing times, a verified lower bound, and a verified constructive argument. Unlike many existence-of-schedule results, the survey's construction is explicit (orders and start times), so it can be formalized directly rather than through an abstract existence argument.

Difficulty

The lower bound is the easy half. The difficulty is the construction: a schedule must meet the bound simultaneously on both machines and for every job. The natural first idea, running every job on M1M_1M1​ then M2M_2M2​ in some order (a flow-shop schedule), fails: Johnson's rule then gives a makespan that can exceed max⁡{T1,T2,max⁡j(aj+bj)}\max\{T_1, T_2, \max_j(a_j+b_j)\}max{T1​,T2​,maxj​(aj​+bj​)}, because the open shop needs some job to visit M2M_2M2​ first. The survey's construction moves exactly one job, JrJ_rJr​, to the front of M2M_2M2​, and its correctness depends on the choice of JrJ_rJr​ and JlJ_lJl​ and on the case split on T1−alT_1 - a_lT1​−al​ versus T2−brT_2 - b_rT2​−br​. Figures 5.3 and 5.4 are drawn for T1≥T2T_1 \ge T_2T1​≥T2​; in the other subcases the start times shown in the figures need adjustment, and the formal statements of the cases assert lengths rather than the figures' exact start times.

Formalization scope

All declarations live in the namespace SchedSurvey.O2. Jobs are Fin n, 0-based (JjJ_jJj​ is index j−1j-1j−1). Processing times and start times are real numbers with aj,bj≥0a_j, b_j \ge 0aj​,bj​≥0; the survey's integer data are a special case, and every claim of §5.2.1 holds over the reals. A schedule is a pair of start-time functions s₁ s₂ : Fin n → ℝ. Interval non-overlap is s + p ≤ s' ∨ s' + p' ≤ s, so touching intervals are allowed and a zero-length operation occupies nothing. Feasibility (IsFeasible) contains both disjointness constraints of §2.1, the machine constraint and the job constraint, with the two operations of a job in either order. The job constraint is essential: without it the optimum would be max⁡{T1,T2}\max\{T_1, T_2\}max{T1​,T2​} and the goal false.

"Length at most LLL" is CompletesBy a b S L: every operation ends by LLL. Optimality is stated in threshold form, which avoids taking a supremum or infimum over a possibly empty set and is equivalent to "Cmax⁡∗C^*_{\max}Cmax∗​ equals the maximum and is attained". A maximum over AAA or BBB appears as a bound on every member, which is also correct for empty AAA or BBB. The orders of A′A'A′ and B′B'B′ are duplicate-free lists whose members are exactly those sets; back-to-back start times are given by contigStart.

The milestone on the choice of JrJ_rJr​, JlJ_lJl​ assumes n≥2n \ge 2n≥2, which the page presupposes; the goal does not, and covers n≤1n \le 1n≤1 as well. The lower bound assumes T≥0T \ge 0T≥0, which matters only for n=0n = 0n=0.

A trivializing formalization is ruled out: feasibility includes both disjointness constraints and nonnegative start times, the threshold is quantified over all T≥0T \ge 0T≥0, and no constant is fixed.

Contributions welcome: proofs of the lower bound (a sum of disjoint intervals inside [0,T][0, T][0,T]), of the list-based block lemmas, of the case lemmas, and of the goal from them. The interval and back-to-back-schedule lemmas are reusable for other shop problems.

Selected references

  • R.L. Graham, E.L. Lawler, J.K. Lenstra, A.H.G. Rinnooy Kan, Optimization and approximation in deterministic sequencing and scheduling: a survey, Annals of Discrete Mathematics 5 (1979) 287–326. https://doi.org/10.1016/S0167-5060(08)70356-X
  • T. Gonzalez, S. Sahni, Open shop scheduling to minimize finish time, Journal of the ACM 23(4) (1976) 665–679. https://doi.org/10.1145/321978.321985
  • S.M. Johnson, Optimal two- and three-stage production schedules with setup times included, Naval Research Logistics Quarterly 1 (1954) 61–68. https://doi.org/10.1002/nav.3800010110
9 thms2 active usersReviewed
🏆Completed
Operations ResearchTheoretical Computer Science·Captain: mikedeng1

Bounds on Multiprocessing Timing Anomalies 2: Scheduling the k Longest Independent Tasks Optimally First Gives ω(k)/ω₀ ≤ 1 + (1 − 1/n)/(1 + ⌊k/n⌋)Research Paper

Motivation

Parallel processing poses a basic allocation question: given tasks of known lengths and several identical processors, how much time is lost when tasks are assigned by a simple list rule instead of an optimal allocation? The finishing time is the moment the last processor completes its work. For independent tasks, a list rule starts each next task on a processor that becomes free first. Such a rule is easy to execute, but the resulting allocation can be worse than the best partition of tasks among processors. Graham's 1969 paper studies precise worst-case ratios for this model. Its Theorem 3 asks what guarantee follows when the first kkk tasks of the list are chosen among the longest and scheduled optimally before the remaining tasks are appended.

The paper also notes two endpoints of this family. With k=0k=0k=0, the rule is ordinary list scheduling and has ratio at most 2−1/n2-1/n2−1/n. With k=nk=nk=n, the nnn longest tasks can start one per processor, giving ratio at most 3/2−1/(2n)3/2-1/(2n)3/2−1/(2n). Theorem 3 expresses these as instances of one bound, including the integer jump at each multiple of nnn Graham, pp. 427–428.

Setting

There are r>0r>0r>0 independent tasks, indexed by jin{0,…,r−1}jin\{0,\ldots,r-1\}jin{0,…,r−1}, and n>0n>0n>0 identical processors, indexed by p∈{0,…,n−1}p\in\{0,\ldots,n-1\}p∈{0,…,n−1}. Task jjj takes a strictly positive time μj\mu_jμj​. A priority list LLL contains every task exactly once. Its first kkk tasks are kkk of the longest tasks; equal lengths may be ordered either way. As a processor becomes available, it takes the next task in the list. Once started, a task runs without interruption. With no precedence constraints, every unstarted task is ready, so a processor runs its assigned tasks consecutively from time zero.

An assignment σ\sigmaσ records the processor that executes each task. The load ℓp(σ)\ell_p(\sigma)ℓp​(σ) of processor ppp is the sum of its assigned task lengths. The finishing time is the largest load, ω(k)=max⁡pℓp(σ)\omega(k)=\max_p\ell_p(\sigma)ω(k)=maxp​ℓp​(σ). For a list assignment, each task is assigned to a processor whose load from earlier list tasks is smallest at that step. If two processors are tied, either can be selected without changing the finishing time. For any assignment τ\tauτ, let ωk(τ)\omega_k(\tau)ωk​(τ) be the largest load contributed by only the first kkk list tasks. The algorithm requires its own first kkk assignments to minimize ωk(τ)\omega_k(\tau)ωk​(τ) over all assignments τ\tauτ; the remaining tasks may occur in any list order.

Write ω0\omega_0ω0​ for the smallest finishing time achievable for all rrr tasks. The paper introduces it as the minimum over possible lists and later describes the same problem as minimizing the largest part sum over partitions of the task lengths Graham, pp. 421, 428. The formalization uses the partition or assignment form of ω0\omega_0ω0​.

Formalization targets

Theorem 3

For 0≤k≤r0\le k\le r0≤k≤r, when the first kkk longest tasks have an optimal prefix schedule, Graham's bound is

ω(k)ω0≤1+1−1/n1+⌊k/n⌋.\frac{\omega(k)}{\omega_0} \le 1+\frac{1-1/n}{1+\lfloor k/n\rfloor}.ω0​ω(k)​≤1+1+⌊k/n⌋1−1/n​.

Here ⌊k/n⌋\lfloor k/n\rfloor⌊k/n⌋ is the greatest integer at most k/nk/nk/n. The theorem also says the bound is best possible when kkk is divisible by nnn Graham, Theorem 3, p. 427. The equality examples are a separate companion statement in this mission.

Source milestones

The milestones state the finite-volume lower bound on ω0\omega_0ω0​, the case in which completing all tasks takes no longer than completing the first kkk, and the paper's numbered inequalities (14), (15), and (16). They use α∗\alpha^*α∗ for the greatest task length after the first kkk positions. These are individual mathematical claims from the proof on p. 427, with their original wording and formulas recorded alongside the formal statements.

Significance

The theorem gives a quantitative tradeoff between work spent optimizing a prefix and the worst-case finishing time of the full list. Its denominator is 1+⌊k/n⌋1+\lfloor k/n\rfloor1+⌊k/n⌋, so the guarantee improves when the optimized prefix contains another full processor's worth of long tasks. The example at multiples of nnn shows that the stated coefficient cannot be uniformly reduced for the algorithm as specified Graham, p. 427.

The result is proved in the paper. This mission's remaining work is a machine-checked proof of its exact statement and the source's intermediate inequalities. The published identical-machine makespan definition is reused, while the list-assignment rule, optimal prefix, and longest-remaining-task quantity are made explicit here. Those definitions can also support later work on list scheduling without precedence constraints. No machine-checked proof of Theorem 3 is claimed by this proposal.

Difficulty

An optimal schedule for the first kkk long tasks does not make the entire list optimal: short tasks appended later can affect which processor finishes last. A bound based only on average total work ignores this final imbalance; a bound based only on the largest remaining task ignores how many long tasks every assignment must place together. The proof must relate both constraints to one finishing time while preserving the floor ⌊k/n⌋\lfloor k/n\rfloor⌊k/n⌋. Replacing that floor with the real quotient changes the claim, and allowing an arbitrary prefix arrangement loses the algorithm's required optimality.

Formalization scope

Tasks and processors are finite indexed sets Fin r and Fin n; their zero-based indices correspond to the paper's one-based TjT_jTj​ and PiP_iPi​. Task lengths are positive real numbers, and both rrr and nnn are positive. The list is a permutation, not merely a sequence that might omit or repeat a task. The condition k≤rk\le rk≤r is explicit because choosing kkk tasks from rrr requires it; the paper treats r>kr>kr>k inside the proof and the case k=rk=rk=r is covered by the theorem. There is no precedence relation in this mission.

The list rule is a predicate on assignments: each next task goes to a processor with least current load. It omits the paper's smaller-index tie convention, since tied identical processors can be interchanged without changing the finishing time. The first kkk positions must contain kkk longest tasks, and their assignment must minimize the prefix finishing time among all assignments of those tasks. Those conditions are part of the goal; the inequalities (14)–(16) are conclusions to establish, not assumptions of the goal. A separate existence statement and a checked concrete instance ensure the list predicate is satisfiable.

All maxima and minima range over finite nonempty processor or assignment sets. The optimum is a minimum over assignments, equivalent to the paper's minimum over lists in the independent-task model. The quantity α∗\alpha^*α∗ is the largest duration outside the prefix when k<rk<rk<r, and is defined as zero when none remains; statements using it require k<rk<rk<r. Division is by n>0n>0n>0 and by ω0>0\omega_0>0ω0​>0. Lean's natural-number quotient represents ⌊k/n⌋\lfloor k/n\rfloor⌊k/n⌋; real subtraction is used for coefficients such as n−1n-1n−1. Contributions toward the finite load identities, the optimum-over-lists equivalence, and proofs of the listed bounds are within scope.

Selected references

  • R. L. Graham, Bounds on Multiprocessing Timing Anomalies, SIAM Journal on Applied Mathematics 17(2), 416–429, 1969. DOI: 10.1137/0117039.
8 thms2 active usersReviewed
🏆Completed
Linear OptimizationOperations Research·Captain: mikedeng1

A Linear Programming Approach to the Cutting Stock Problem—Part II 1: The Lexicographic Knapsack Method Terminates, and M_1 Equals max(c_1, the Knapsack Optimum for the Longest Stock Length L_1)Research Paper

Motivation

The cutting stock problem asks how to cut standard stock rolls of length LLL into pieces of demanded lengths l1,…,lml_1,\dots,l_ml1​,…,lm​ so that the demands are met at least cost. Gilmore and Gomory solved its linear programming relaxation by column generation (Part I, Opns. Res. 9 (1961)): every cutting pattern is a column, and the column to bring into the basis is found by solving an integer knapsack problem whose objective is given by the current dual prices. In Part II (Opns. Res. 11 (1963)) they replaced the dynamic-programming knapsack solver of Part I by a lexicographic enumeration with a bound test, the Knapsack Method, and reported it to be about five times as fast as dynamic programming on their problem set (p. 866). The method is an early enumeration-with-bounds (branch-and-bound type) procedure for the integer knapsack, and the same pricing problem remains the inner loop of column generation for cutting stock, bin packing and their many variants.

Setting

There are mmm demanded lengths l1,…,lm>0l_1,\dots,l_m>0l1​,…,lm​>0 with linear programming prices b1,…,bm≥0b_1,\dots,b_m\ge0b1​,…,bm​≥0, and kkk standard stock lengths L1,…,LkL_1,\dots,L_kL1​,…,Lk​ with costs c1,…,ckc_1,\dots,c_kc1​,…,ck​. A column for LjL_jLj​ is a vector (α)m=(a1,…,am)(\alpha)_m=(a_1,\dots,a_m)(α)m​=(a1​,…,am​) of nonnegative integers with ∑iliai≤Lj\sum_i l_ia_i\le L_j∑i​li​ai​≤Lj​. The knapsack problem (1) for LjL_jLj​ is

Mˉj=max⁡{∑i=1mbiai:ai∈Z≥0, ∑i=1mliai≤Lj},\bar M_j=\max\Big\{\sum_{i=1}^m b_ia_i : a_i\in\mathbb Z_{\ge0},\ \sum_{i=1}^m l_ia_i\le L_j\Big\},Mˉj​=max{i=1∑m​bi​ai​:ai​∈Z≥0​, i=1∑m​li​ai​≤Lj​},

and the quantity to compute is Mj=max⁡(cj,Mˉj)M_j=\max(c_j,\bar M_j)Mj​=max(cj​,Mˉj​): a column improves the linear program exactly when Mˉj>cj\bar M_j>c_jMˉj​>cj​.

For a prefix (α)s=(a1,…,as)(\alpha)_s=(a_1,\dots,a_s)(α)s​=(a1​,…,as​) write λ⋅(α)s=∑i≤sliai\lambda\cdot(\alpha)_s=\sum_{i\le s}l_ia_iλ⋅(α)s​=∑i≤s​li​ai​ and β⋅(α)s=∑i≤sbiai\beta\cdot(\alpha)_s=\sum_{i\le s}b_ia_iβ⋅(α)s​=∑i≤s​bi​ai​. Step (1) sorts the data, b1/l1≥⋯≥bm/lmb_1/l_1\ge\dots\ge b_m/l_mb1​/l1​≥⋯≥bm​/lm​ and L1>⋯>LkL_1>\dots>L_kL1​>⋯>Lk​, and appends a slack item am+1a_{m+1}am+1​ with bm+1=0b_{m+1}=0bm+1​=0, lm+1=1l_{m+1}=1lm+1​=1. The method then generates vectors in decreasing lexicographic order:

  • Steps (2) and (7) fill the coefficients after a prefix greedily, ai=[(L−used)/li]a_i=[(L-\text{used})/l_i]ai​=[(L−used)/li​], for L=L1L=L_1L=L1​ at the start and L=LtL=L_tL=Lt​ later.
  • Step (3) records β⋅(α)m\beta\cdot(\alpha)_mβ⋅(α)m​ as the new MjM_jMj​ for each j≥tj\ge tj≥t with Lj≥λ⋅(α)mL_j\ge\lambda\cdot(\alpha)_mLj​≥λ⋅(α)m​ and β⋅(α)m>Mj\beta\cdot(\alpha)_m>M_jβ⋅(α)m​>Mj​.
  • Step (4) finds the last nonzero coefficient asa_sas​.
  • Step (5) lowers asa_sas​ by one and looks for the smallest jjj with Lj≥λ⋅(α)sL_j\ge\lambda\cdot(\alpha)_sLj​≥λ⋅(α)s​ and (Lj−λ⋅(α)s)bs+1>(Mj−β⋅(α)s)ls+1(L_j-\lambda\cdot(\alpha)_s)b_{s+1}>(M_j-\beta\cdot(\alpha)_s)l_{s+1}(Lj​−λ⋅(α)s​)bs+1​>(Mj​−β⋅(α)s​)ls+1​. This is the bound obtained by relaxing the integrality of as+1a_{s+1}as+1​.
  • Step (6) backs up to the previous nonzero coefficient when no such jjj exists, and stops when there is none.

In Lean the method is a transition system Step on states (a, s, t, M, phase). Reachable means reachable from the start state of Step (2).

Formalization targets

Goal: the method terminates and computes M1M_1M1​

Under Step (1)'s ordering, with li>0l_i>0li​>0 and bi≥0b_i\ge0bi​≥0:

every run is finite,a non-final state has a successor,M1final=max⁡(c1,Mˉ1),\text{every run is finite},\quad \text{a non-final state has a successor},\quad M_1^{\text{final}}=\max\big(c_1,\bar M_1\big),every run is finite,a non-final state has a successor,M1final​=max(c1​,Mˉ1​),

and at the end every MjM_jMj​ is cjc_jcj​ or the price of a column for LjL_jLj​, with cj≤Mjc_j\le M_jcj​≤Mj​. With one stock length (k=1k=1k=1) this is the paper's correctness claim in full: the method solves the knapsack problem (1).

Milestones

  1. Steps (2)/(7): the greedy completion is the lexicographically largest extension of a prefix that fits the capacity.
  2. Step (4): the lexicographic predecessor of a fitting vector extends (α1)s(\alpha^1)_s(α1)s​.
  3. The relaxation bound β⋅(α′)m≤β⋅(α)s+bs+1(L−λ⋅(α)s)/ls+1\beta\cdot(\alpha')_m\le\beta\cdot(\alpha)_s+b_{s+1}(L-\lambda\cdot(\alpha)_s)/l_{s+1}β⋅(α′)m​≤β⋅(α)s​+bs+1​(L−λ⋅(α)s​)/ls+1​ for every fitting extension.
  4. After the inequality of Step (5)'s test fails, no vector between (α)s(\alpha)_s(α)s​ and the successor of Step (6) improves MMM.
  5. The vectors tested in Step (3) strictly decrease lexicographically.
  6. The invariant on entry to Step (3).
  7. Identical prices: among lengths with equal prices only the shortest is needed (p. 869).

Significance

The goal certifies the pricing oracle of the Gilmore–Gomory column-generation method: the knapsack solver it calls returns the true optimum. The pricing step is exact only if that holds, and so is the conclusion "no improving column exists, the LP is solved". Milestones 1–4 are reusable on their own. They are the textbook ingredients of lexicographic and depth-first branch-and-bound for the integer knapsack: the greedy completion, the density-ordered LP bound, and the pruning rule.

The paper's argument is a short paragraph, "clear from the order in which the mmm-vectors are generated". Writing it out shows that the claim is false for several stock lengths as the steps are printed. With m=1m=1m=1, l1=b1=1l_1=b_1=1l1​=b1​=1, L=(5,2)L=(5,2)L=(5,2) and c=(0,0)c=(0,0)c=(0,0), the method tests (5)(5)(5). The test of Step (5) at (4)(4)(4) then fails for j=2j=2j=2 only because 4>L24>L_24>L2​, Step (6) stops, and the run ends with M2=0M_2=0M2​=0 while the column (2)(2)(2) has price 2. This mission therefore states what the printed algorithm does achieve. The result is the full claim for k=1k=1k=1 and for the longest stock length in general. A sorry-free Lean proof of the counterexample run accompanies the mission. No part of the method or its correctness has been formalized before, as far as a search of the platform shows.

Difficulty

The method skips large parts of the lexicographic order in three ways:

  • Step (7) jumps to the greedy completion for LtL_tLt​ rather than L1L_1L1​.
  • Step (6) jumps past every smaller value of asa_sas​ once one test fails.
  • Step (3) updates only j≥tj\ge tj≥t.

Each skip has to be justified by a bound that holds for every vector skipped, and the justifications interact. A skipped vector may fit a longer stock length but not LtL_tLt​, and the bound for L1L_1L1​ must cover exactly those vectors. The invariant that survives all three moves is not the naive "every vector above the current one has been examined": the vectors are not examined, only bounded. Termination needs the lexicographic decrease together with finiteness of the fitting vectors, and progress needs the reachability invariants 1≤s≤m1\le s\le m1≤s≤m and as≠0a_s\ne0as​=0 that make Step (5) meaningful.

Formalization scope

  • Items are Fin m from 0. The paper's level sss is a natural number, and the prefix (α)s(\alpha)_s(α)s​ is the coordinates with index < s. Stock lengths are Fin k, and L1L_1L1​ is index 0. The slack item enters only through bs+1b_{s+1}bs+1​ and ls+1l_{s+1}ls+1​ at s=ms=ms=m. The lexicographic order is Mathlib's Pi.Lex, and [x][x][x] is Nat.floor.

  • Mj=max⁡(cj,Mˉj)M_j=\max(c_j,\bar M_j)Mj​=max(cj​,Mˉj​) is the predicate IsKnapsackMax, never a real supremum, which would be 0 on an empty set. The algorithm's output is never defined as the maximum. The goal includes termination and progress, so its correctness clause is not vacuous. Every invariant is stated for reachable states only.

  • The goal adds two disclosed hypotheses:

    • li>0l_i>0li​>0: demanded lengths.
    • bi≥0b_i\ge0bi​≥0: with a negative price the method is wrong, and the paper applies it to cutting-stock prices, which are nonnegative.

    Step (1) is a hypothesis on the data, not a step.

  • The goal also makes these corrections:

    • The claim is narrowed. For every jjj it states cj≤Mjc_j\le M_jcj​≤Mj​ and attainment, and the upper bound for L1L_1L1​ only (see Significance).
    • The lexicographic order is repaired. The page's definition misses vectors that differ in their first coefficient.
    • Step (7) reads λ·(α)_m. The page's "satisfies Lt≥λ⋅(α)sL_t\ge\lambda\cdot(\alpha)_sLt​≥λ⋅(α)s​" is read as λ⋅(α)m\lambda\cdot(\alpha)_mλ⋅(α)m​.
    • Step (4) gets a stopping case. On the zero vector, where the page has no case, the method stops.
    • Milestone 4 is stated for the failed inequality, not for a failed test, because the failed test is what makes the page's successor claim false.
  • Not formalized:

    • Step (3′), the variant computing max⁡j(Mj−cj)\max_j(M_j-c_j)maxj​(Mj​−cj​).
    • The record of patterns.
    • The dynamic-programming recursion and the timing estimates.

Contributions welcome: proofs of the milestones, and a corrected multi-length variant of the method with a proof that it computes every MjM_jMj​.

Selected references

  • P. C. Gilmore, R. E. Gomory, A linear programming approach to the cutting stock problem—Part II, Operations Research 11(6) (1963), 863–888. https://doi.org/10.1287/opre.11.6.863
  • P. C. Gilmore, R. E. Gomory, A linear programming approach to the cutting-stock problem, Operations Research 9(6) (1961), 849–859. https://doi.org/10.1287/opre.9.6.849
  • G. B. Dantzig, Discrete-variable extremum problems, Operations Research 5(2) (1957), 266–288. https://doi.org/10.1287/opre.5.2.266
10 thms2 active usersReviewed
🏆Completed
Operations ResearchOptimization·Captain: mikedeng1

One-Processor Scheduling with Symmetric Earliness and Tardiness Penalties 3: For Maximum Discrepancy, Some Optimum Schedule Is in Standard Form and Has the Release Time PropertyResearch Paper

Motivation

In just-in-time scheduling each job has a preferred time, and finishing it early is as undesirable as finishing it late. Garey, Tarjan and Wilfong (Math. Oper. Res. 13 (1988), 330–348) study the one-processor version in which task TiT_iTi​ has a length li≥0l_i \ge 0li​≥0 and a preferred starting time ai≥0a_i \ge 0ai​≥0, and the discrepancy of a task started at time sis_isi​ is ∣si−ai∣|s_i - a_i|∣si​−ai​∣. The penalty is symmetric: earliness and tardiness cost the same. The paper shows that minimizing the total discrepancy is NP-complete, gives an O(Nlog⁡N)O(N \log N)O(NlogN) algorithm for a fixed task order, and, in §3, gives an efficient algorithm for minimizing the maximum discrepancy max⁡i∣si−ai∣\max_i |s_i - a_i|maxi​∣si​−ai​∣: a linear-time test for a given bound γ\gammaγ after an O(Nlog⁡N)O(N \log N)O(NlogN) sort, combined with a search over γ\gammaγ.

This mission formalizes the structural core of §3: the normal-form theorems (Theorem 4, Corollaries 1 and 2) on which the paper's algorithm for maximum discrepancy is built. The two companion missions of the series cover the NP-completeness result (THEOREM 1) and the fixed-order algorithm (THEOREMS 2 and 3).

Setting

Fix a bound γ≥0\gamma \ge 0γ≥0. A schedule has maximum discrepancy at most γ\gammaγ exactly when every task starts no earlier than its release time ri=max⁡{0,ai−γ}r_i = \max\{0, a_i - \gamma\}ri​=max{0,ai​−γ} and finishes no later than its deadline di=ai+li+γd_i = a_i + l_i + \gammadi​=ai​+li​+γ (p. 342). These release times and deadlines are in special form: with α=2γ\alpha = 2\gammaα=2γ, for every task

di−ri=li+αor(ri=0 and di<li+α).d_i - r_i = l_i + \alpha \quad\text{or}\quad \big(r_i = 0 \text{ and } d_i < l_i + \alpha\big).di​−ri​=li​+αor(ri​=0 and di​<li​+α).

From here on the data are arbitrary reals li≥0l_i \ge 0li​≥0, ri≥0r_i \ge 0ri​≥0, did_idi​ in special form for some α≥0\alpha \ge 0α≥0.

A schedule is an execution order σ\sigmaσ (σ(p)\sigma(p)σ(p) is the task in position ppp) with starting times sis_isi​, such that a task finishes no later than any later-positioned task starts: sσ(p)+lσ(p)≤sσ(q)s_{\sigma(p)} + l_{\sigma(p)} \le s_{\sigma(q)}sσ(p)​+lσ(p)​≤sσ(q)​ for p<qp < qp<q. It is feasible if ri≤sir_i \le s_iri​≤si​ and si+li≤dis_i + l_i \le d_isi​+li​≤di​ for every iii. Its makespan is max⁡i(si+li)\max_i (s_i + l_i)maxi​(si​+li​), and it is optimum if it is feasible and no feasible schedule has a smaller makespan.

Let TNT_NTN​ be a task of largest deadline. A schedule has the release time property if every task executed after TNT_NTN​ has release time strictly later than the starting time of TNT_NTN​. When the tasks are indexed with d1≤⋯≤dNd_1 \le \dots \le d_Nd1​≤⋯≤dN​, a schedule is in standard form with split index jjj, 0≤j≤N−10 \le j \le N-10≤j≤N−1, if it begins with an optimum schedule of T1,…,TjT_1, \dots, T_jT1​,…,Tj​, followed by TNT_NTN​ started at the maximum of rNr_NrN​ and the completion time of that first part, followed by Tj+1,…,TN−1T_{j+1}, \dots, T_{N-1}Tj+1​,…,TN−1​ in this order without idle time.

Formalization targets

Goal: Corollary 2 (p. 344)

∃ a feasible schedule  ⟹  ∃ an optimum schedule that is in standard form and has the release time property.\exists\ \text{a feasible schedule} \;\Longrightarrow\; \exists\ \text{an optimum schedule that is in standard form and has the release time property.}∃ a feasible schedule⟹∃ an optimum schedule that is in standard form and has the release time property.

Milestones

  1. Theorem 4 (pp. 342–343): if a feasible schedule exists, some optimum schedule has the release time property with respect to any task of largest deadline.
  2. The deadline split (proof of Corollary 1, p. 344): in every feasible schedule with the release time property, a task Ti≠TNT_i \neq T_NTi​=TN​ precedes TNT_NTN​ if and only if di≤sN+αd_i \le s_N + \alphadi​≤sN​+α.
  3. Corollary 1 (p. 344): some optimum schedule executes before TNT_NTN​ exactly the tasks of deadline at most β\betaβ, for some β≥0\beta \ge 0β≥0.
  4. Release before finish (p. 344): every task executed after TNT_NTN​ has ri≤sN+lNr_i \le s_N + l_Nri​≤sN​+lN​.
  5. Normalization after TNT_NTN​ (p. 344): an optimum schedule with the release time property can be changed, without moving TNT_NTN​ later and without touching the tasks before it, so that there is no idle time from the start of TNT_NTN​ on and the later tasks run in deadline order.

Two companion items state the reformulation of the maximum-discrepancy bound as release times and deadlines, and the special form with α=2γ\alpha = 2\gammaα=2γ.

Significance

Corollary 2 reduces the search for an optimum schedule of T1,…,TnT_1, \dots, T_nT1​,…,Tn​ to nnn candidates, one per split index, each assembled from an optimum schedule of a shorter prefix. This is the dynamic program of §3.2, which the paper implements in O(N)O(N)O(N) time after an O(Nlog⁡N)O(N \log N)O(NlogN) sort; combined with a search over γ\gammaγ it minimizes the maximum discrepancy. Without the special form, deciding whether one processor can meet arbitrary release times and deadlines is NP-complete (reference [4] of the paper, Garey and Johnson 1979), so the normal form is what separates the tractable case from the general one.

The results are proved in the paper; no machine-checked proof of them is known on Prove2Me. A formal proof of Corollary 2 would certify the correctness of the split-index recursion and, together with the companions, of the reduction from maximum discrepancy to this release-time/deadline problem.

Difficulty

The obvious argument fails in Theorem 4. Exchanging a straggler (a task after TNT_NTN​ released no later than TNT_NTN​ starts) with TNT_NTN​ shifts the tasks between them, and for general release times and deadlines those tasks can become infeasible. The paper's argument uses the special form at every step: a task released after TNT_NTN​ starts has ri>0r_i > 0ri​>0, hence di−ri=li+αd_i - r_i = l_i + \alphadi​−ri​=li​+α exactly, and this equality is what bounds how far tasks may move. It also needs an extremal choice of the optimum schedule (fewest tasks after TNT_NTN​, then fewest tasks between TNT_NTN​ and the first straggler), which requires showing that optimum schedules exist over the reals. The passage to standard form then combines this with an earliest-deadline exchange for the tasks after TNT_NTN​ and the replacement of the first part by an optimum sub-schedule without losing feasibility of the later tasks.

Formalization scope

Tasks are indexed by Fin N (0-based: the paper's T1,…,TNT_1, \dots, T_NT1​,…,TN​ are 0,…,N−10, \dots, N-10,…,N−1; the paper's TNT_NTN​ in Corollary 2 is the index N−1N-1N−1). All data are real. A schedule is a permutation σ : Fin N ≃ Fin N with starting times s : Fin N → ℝ; "executed before/after" refers to σ, not to a comparison of starting times, because zero-length tasks may share a starting time. Execution intervals meet at most at endpoints. The makespan is a supremum over Fin N; "optimum" quantifies over all feasible schedules, not only standard-form ones.

The standing hypotheses on every structural item are α≥0\alpha \ge 0α≥0, ri≥0r_i \ge 0ri​≥0, li≥0l_i \ge 0li​≥0 (from γ≥0\gamma \ge 0γ≥0, the max⁡{0,⋅}\max\{0, \cdot\}max{0,⋅} in rir_iri​, and nonnegative lengths); since si≥ris_i \ge r_isi​≥ri​, start times are nonnegative, as the paper assumes throughout. The special form is a hypothesis of every structural item and cannot be dropped. The paper's dummy task T0T_0T0​ (r0=d0=l0=0r_0 = d_0 = l_0 = 0r0​=d0​=l0​=0) is not a task; an empty first part completes at time 000. Corollary 1's "the set of all tasks with deadline β\betaβ or less" is read as excluding TNT_NTN​, the only reading under which it is true. The page's "(which can be assumed optimum)" is part of the standard form. The deadline-split and release-before-finish milestones are stated for every feasible schedule, as their arguments allow.

A trivializing formalization is excluded: the standard form pins TNT_NTN​'s start to max⁡(rN,C)\max(r_N, C)max(rN​,C) and the later tasks to consecutive positions in index order, so the goal is not Corollary 1 restated with an unconstrained split.

Running times (O(Nlog⁡N)O(N \log N)O(NlogN), O(N)O(N)O(N)) and the algorithm of §3.2 are out of scope. The development needs only finite permutations, finite suprema and the exchange arguments of §3.1; the existence of an optimum schedule (minimum makespan over finitely many orders, earliest-start schedules) is reusable for other single-machine problems with release times and deadlines. Proofs of any milestone, and a sorry-free proof of the existence of optimum schedules, are welcome.

Selected references

  • M. R. Garey, R. E. Tarjan, G. T. Wilfong, One-Processor Scheduling with Symmetric Earliness and Tardiness Penalties, Mathematics of Operations Research 13(2):330–348, 1988. https://doi.org/10.1287/moor.13.2.330
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, W. H. Freeman, 1979 (reference [4] of the paper, cited on p. 342 for the NP-completeness of one-processor scheduling with release times and deadlines). ISBN 0-7167-1045-5
7 thms2 active usersReviewed
🏆Completed
Operations ResearchOptimization·Captain: mikedeng1

One-Processor Scheduling with Symmetric Earliness and Tardiness Penalties 2: For a Fixed Task Order, the Block-Shifting Algorithm Minimizes the Total DiscrepancyResearch Paper

Motivation

A task may have a preferred execution time because a material delivery, external event, or downstream operation is timed to it. Starting early can be as costly as starting late. Garey, Tarjan, and Wilfong study this situation on one processor, where tasks cannot overlap but idle time between tasks is permitted. Their total-discrepancy problem is hard when the execution order is free; fixing the order leaves a substantial timing problem because each start can still move and can force changes to earlier starts. The paper gives an explicit scheduling procedure for that case and proves that it minimizes the sum of deviations from preferred starts (Garey, Tarjan, and Wilfong, 1988, §§1–2).

This mission targets that procedure and its correctness theorem. It also records the paper's equal-length result: when task lengths are identical, a minimum-cost schedule exists in preferred-start order, so the fixed-order procedure applies after sorting. These are two precise claims about total discrepancy, distinct from the paper's separate maximum-discrepancy problem.

Setting

There are nnn tasks, indexed i=0,…,n−1i=0,\ldots,n-1i=0,…,n−1. Task iii has a nonnegative length lil_ili​, a nonnegative preferred starting time aia_iai​, and an actual starting time sis_isi​. Once started, it occupies the one processor until si+lis_i+l_isi​+li​. Each start must be nonnegative. For the fixed-order problem, task iii must finish before task i+1i+1i+1 starts, so si+li≤si+1s_i+l_i\le s_{i+1}si​+li​≤si+1​. This permits idle time when the inequality is strict. The task's discrepancy is ∣si−ai∣|s_i-a_i|∣si​−ai​∣; because its preferred completion is ai+lia_i+l_iai​+li​, this is also its absolute completion-time discrepancy. The total discrepancy is

cost⁡n(s)=∑i=0n−1∣si−ai∣.\operatorname{cost}_n(s)=\sum_{i=0}^{n-1}|s_i-a_i|.costn​(s)=i=0∑n−1​∣si​−ai​∣.

A block is a maximal consecutive set of tasks with no idle time between neighboring tasks. Within a block, Decrease counts tasks whose actual starts are later than preferred, and Increase counts tasks whose actual starts are no later than preferred. These names describe the effect of moving the entire block earlier: the discrepancy of a Decrease task falls initially, whereas that of an Increase task rises. The algorithm's state is a schedule SnS_nSn​ for the first nnn tasks. It inserts the next task at its preferred start when the previous task has finished, and at that finish time otherwise. In the latter case it may move the final block earlier until one of the paper's stopping events occurs: the block reaches time zero, a late task reaches its preferred start, or the block meets its predecessor (§2.2, p. 337).

For the equal-length result, schedules may execute tasks in any order. Two distinct tasks are feasible together when one completes before the other starts; meeting at endpoints is allowed. The objective remains the same sum of absolute discrepancies.

Formalization targets

Fixed-order optimality

The main target is the paper's Theorem 2. For all nonnegative task data and every nnn, the actual schedule SnS_nSn​ produced by the block-shifting algorithm is feasible and satisfies

cost⁡n(Sn)≤cost⁡n(s)for every feasible fixed-order schedule s.\operatorname{cost}_n(S_n)\le\operatorname{cost}_n(s) \qquad\text{for every feasible fixed-order schedule }s.costn​(Sn​)≤costn​(s)for every feasible fixed-order schedule s.

This includes the empty schedule and schedules with idle time. The milestone list follows the paper's own assertions: a balanced final block can move earlier without changing cost; every block of SnS_nSn​ has more Increase tasks than Decrease tasks or begins at zero (Lemma 7); and the two claims in Case 2 of the proof record the cost of adding a late task and the comparison with schedules that place it earlier (Theorem 2 and Lemma 7, p. 338).

Equal-length schedules

The companion target is Theorem 3. If every task has a common length L≥0L\ge0L≥0 and the preferred starts are indexed so that ai≤ai+1a_i\le a_{i+1}ai​≤ai+1​, then among all feasible schedules, including those with another task order, at least one minimum-cost schedule starts the tasks in index order:

∃s  [cost⁡n(s)≤cost⁡n(t) for every feasible t]with si≤si+1 whenever i+1<n.\exists s\;\bigl[ \operatorname{cost}_n(s)\le\operatorname{cost}_n(t) \text{ for every feasible }t \bigr] \quad\text{with }s_i\le s_{i+1}\text{ whenever }i+1<n.∃s[costn​(s)≤costn​(t) for every feasible t]with si​≤si+1​ whenever i+1<n.

The paper uses this statement to connect free-order equal-length scheduling to its fixed-order procedure (Theorem 3, p. 340).

Significance

Theorem 2 certifies an explicit schedule, not just the existence of an optimum. It fixes the objective value for a prescribed execution order and gives a baseline against which any other legal timing of the same tasks can be compared. Theorem 3 supplies the ordering fact needed to use that result when all lengths agree. Together they explain why a problem that is difficult for unrestricted task lengths still has these structured solvable cases (Garey, Tarjan, and Wilfong, 1988, abstract and §2.5).

The paper proves both theorems on paper. The work here is to formalize its schedule construction, cost, block boundaries, and comparison classes in Lean, then obtain machine-checked proofs of the stated targets. The draft theorem declarations compile with proof placeholders; this proposal does not claim that they are already machine-checked results. A completed development would make the model and the algorithm available for later formal work on scheduling with idle time and symmetric earliness and tardiness penalties.

Difficulty

Moving a task closer to its preferred start can move neighboring tasks farther from theirs. A simple task-by-task choice therefore does not establish global optimality. Nor does a count of late and early tasks in one block, by itself, justify arbitrary earlier movements of individual tasks. The difficult comparison in the paper is between the algorithm's schedule and a competing feasible schedule whose new task starts earlier: the whole final block and its constrained relative movements matter. The proof also has to account for the boundary at time zero and for blocks that merge when shifted. These are the reasons the explicit procedure and its block invariant need careful statements (§§2.2–2.3, pp. 337–338).

Formalization scope

The Lean development represents task data and starts as functions from natural-number indices to real numbers; only indices below nnn count. Task iii in Lean is the paper's Ti+1T_{i+1}Ti+1​. Preferred times, lengths, and legal starting times are nonnegative. The start-time condition is a standing convention used by the paper's algorithm, especially its zero-boundary stopping rule. Fixed-order feasibility requires successive tasks to be separated by at least the earlier task's length; unrestricted feasibility uses pairwise nonoverlap. An empty schedule is legal and has cost zero. The representation allows zero-length tasks, as does the paper's li≥0l_i\ge0li​≥0 model.

The algorithm is defined by its stated insertion and single-block-shift operations. Defining SnS_nSn​ as a chosen minimizer would erase the content of Theorem 2; the goal also explicitly asserts feasibility so that an illegal low-cost function cannot qualify. Theorem 3 compares with every feasible unrestricted-order schedule, not just schedules already in preferred-start order. The core definitions, the block invariant, and the cost comparisons are useful beyond this particular proof. Formalizing the paper's running-time bounds and heap implementation is outside this mission.

Selected references

  • Michael R. Garey, Robert E. Tarjan, and Gordon T. Wilfong, One-Processor Scheduling with Symmetric Earliness and Tardiness Penalties, Mathematics of Operations Research 13(2), 330–348, 1988. DOI: 10.1287/moor.13.2.330.
6 thms2 active usersReviewed
🏆Completed
Theoretical Computer Science·Captain: marwahaha

Truly Subquadratic 3SUM on a Word RAMResearch Paper

3SUM asks whether three distinct positions in a list of integers contain values summing to zero. It is a central problem in fine-grained complexity.

This mission formalizes a deterministic O(n1.9992)O(n^{1.9992})O(n1.9992) algorithm for 3SUM on nnn polynomially bounded integers, on a word RAM with O(log⁡n)O(\log n)O(logn)-bit words. Alman and Vassilevska Williams's 2026 breakthrough gives a truly subquadratic bound, refuting the integer 3SUM hypothesis in this model.

The mission uses the accompanying Lean specification and shares the Exact Triangle milestone with Truly Subcubic APSP on a Word RAM.

References

  1. Josh Alman and Virginia Vassilevska Williams, Truly Subquadratic 3SUM and Truly Subcubic APSP via Triangles in Sparse Lopsided Graphs, 2026, Theorem 22.
  2. Anthropic, formal-math / 3sum-apsp: accompanying Lean formalization, Theorem_22_3SUM.
51 thms2 active usersReviewed
🏆Completed
Operations ResearchTheoretical Computer Science·Captain: mikedeng1

Fibonacci Heaps and Their Uses in Improved Network Optimization Algorithms: O(log n) Amortized Time for Delete Min and Delete, O(1) for Every Other F-Heap OperationResearch Paper

Why priority queues with cheap decrease key matter

Many network optimization algorithms spend most of their time in a heap (priority queue): a set of items with real keys that supports inserting an item, finding and deleting an item of minimum key, melding two heaps, decreasing the key of an item, and deleting an arbitrary item. Dijkstra's shortest-path algorithm performs at most one delete min per vertex and may perform a decrease key for each edge, so on a graph with nnn vertices and mmm edges its running time is governed by how cheap decrease key is. Binomial queues support the heap operations in O(log⁡n)O(\log n)O(logn) worst-case time.

Fredman and Tarjan (J. ACM 34 (1987)) introduced Fibonacci heaps (F-heaps), in which delete min and delete take O(log⁡n)O(\log n)O(logn) amortized time and every other operation, decrease key included, takes O(1)O(1)O(1) amortized time. This brings Dijkstra's algorithm to O(nlog⁡n+m)O(n \log n + m)O(nlogn+m) and improves the best bounds then known for all-pairs shortest paths, weighted bipartite matching (the assignment problem) and minimum spanning trees. F-heaps descend from the binomial queues of Vuillemin (1978) and use the potential technique of amortized analysis of Sleator and Tarjan. This mission formalizes the paper's §2 analysis, culminating in THEOREM 1 on p. 604.

Setting

An F-heap is a list of heap-ordered rooted trees: every child's key is at least its parent's key. Each node holds an item, the item's real key and a mark bit; the children of a node are kept in the order in which they were linked to it, earliest first. The rank r(x)r(x)r(x) of a node is its number of children. A collection is a list of heaps, named by position; a run starts with no heaps.

Linking two trees with roots xxx and yyy makes yyy the last child of xxx if the key of xxx is smaller, and otherwise makes xxx the last child of yyy; the root that becomes a child is unmarked. The operations act as follows:

  • make heap adds an empty heap; find min acts on a nonempty heap and returns a root of minimum key; insert adds a one-node tree for a new item; meld concatenates two root lists.
  • delete min removes a root xxx of minimum key, adds the children of xxx to the roots, then repeats the linking step, "find any two trees whose roots have the same rank, and link them", in any order, until all root ranks are distinct.
  • decrease key (Δ,i,h)(\Delta, i, h)(Δ,i,h), Δ≥0\Delta \ge 0Δ≥0, lowers the key of iii and, if the node xxx of iii is not a root, cuts xxx from its parent, making its subtree a new tree. delete (i,h)(i, h)(i,h) cuts xxx the same way, destroys it and adds its children to the roots; a minimum root is deleted as in delete min.
  • After every cut of a child of ppp, if ppp is not a root it is marked if unmarked and cut from its own parent if marked: a cascading cut, which may repeat upwards.

The actual time of an operation is one unit, plus one per linking step, plus one per cut, plus for delete min the scan of the rank array (the largest rank involved). The potential Φ\PhiΦ of a collection is the number of trees plus twice the number of marked nonroot nodes, and the amortized time of an operation is its actual time plus the increase of Φ\PhiΦ.

Formalization targets

Goal: THEOREM 1 (p. 604)

There is an absolute constant C>0C > 0C>0 such that for every sequence of TTT operations started from no heaps, with ntn_tnt​ the number of items in the heap that operation ttt acts on, just before it acts,

∑t<Tcostt  ≤  ∑t<Tbt,bt={C(log⁡2nt+1)delete min or delete,Cotherwise.\sum_{t<T} \mathrm{cost}_t \;\le\; \sum_{t<T} b_t, \qquad b_t = \begin{cases} C(\log_2 n_t + 1) & \text{delete min or delete},\\ C & \text{otherwise.}\end{cases}t<T∑​costt​≤t<T∑​bt​,bt​={C(log2​nt​+1)C​delete min or delete,otherwise.​

The paper's wording is "the total time is at most the total amortized time, where the amortized time is O(log⁡n)O(\log n)O(logn) for each delete min or delete operation and O(1)O(1)O(1) for each of the other operations"; the goal pins the O(⋅)O(\cdot)O(⋅) to one constant and does not name the potential.

Milestones

  1. LEMMA 1: the iiith child of a node, in linking order, has rank at least i−2i - 2i−2.
  2. Fk+2≥φkF_{k+2} \ge \varphi^kFk+2​≥φk, with FkF_kFk​ the Fibonacci numbers and φ=(1+5)/2\varphi = (1+\sqrt5)/2φ=(1+5​)/2.
  3. COROLLARY 1: a node of rank kkk has at least Fk+2≥φkF_{k+2} \ge \varphi^kFk+2​≥φk descendants, itself included.
  4. A node of rank kkk in an nnn-item heap has φk≤n\varphi^k \le nφk≤n, hence k≤log⁡n/log⁡φk \le \log n / \log\varphik≤logn/logφ.
  5. make heap, find min and meld leave Φ\PhiΦ unchanged; insert raises it by one.
  6. delete min raises Φ\PhiΦ by at most log⁡n/log⁡φ\log n / \log \varphilogn/logφ minus the number of linking steps.
  7. decrease key raises Φ\PhiΦ by at most three minus the number of cascading cuts.
  8. From no heaps, the total amortized time bounds the total actual time.

Significance

THEOREM 1 is the data-structure fact behind the paper's algorithmic results: Dijkstra's algorithm in O(nlog⁡n+m)O(n\log n + m)O(nlogn+m), all-pairs shortest paths and the assignment problem in O(n2log⁡n+nm)O(n^2 \log n + nm)O(n2logn+nm), and minimum spanning trees in O(m β(m,n))O(m\,\beta(m,n))O(mβ(m,n)).

The result has been proved since 1987. What the mission adds is a machine-checked version of the paper's own analysis on an operational model with arbitrary linking order, unmarking on link, marking or cascading on every cut, and the paper's unit costs. The proposal states the theorem and its supporting claims; their proofs remain to be formalized.

Difficulty

The obvious argument, that every tree is a binomial tree and so every rank is at most log⁡2n\log_2 nlog2​n, fails as soon as cuts are allowed: decrease key and delete remove subtrees arbitrarily, and the trees of an F-heap need not be binomial. The rank bound has to come from an invariant of the whole history (LEMMA 1), which only holds because a node loses at most one child before it is cut itself; without cascading cuts a node of rank kkk can be left with only k+1k+1k+1 descendants, and no logarithmic rank bound holds. A second difficulty is that the actual cost of one decrease key is unbounded, since a single operation may trigger arbitrarily many cascading cuts; a per-operation worst-case argument cannot give THEOREM 1, and the bound holds only for whole sequences started from no heaps.

Formalization scope

All declarations sit in the namespace FibHeap.Amort. Trees are a nested inductive FTree with item ℕ, key ℝ, a mark bit and a list of children in linking order; a heap is List FTree (its roots) and a collection is List Heap. Operations are a step relation Step s op s' d whose data d records linking steps, cuts, cascading cuts and the scan, so that cost d = 1 + links + cuts + scan is determined by what the operation did. The linking loop allows every order of linking; a link tie makes the first root a child of the second, while delete min may choose any minimum-key root. Runs start from no heaps; reachable collections consist of heap-ordered trees with each item in at most one heap. LEMMA 1, COROLLARY 1 and the delete-min bound are stated for reachable collections only, as the paper's "a node in an F-heap" intends. The paper's standing assumption that an item is in only one heap at a time is enforced by the precondition of insert and is not a hypothesis. The minimum pointer and returned items are not stored in the step relation; find min is a state-preserving step charged one unit. A heap destroyed by meld remains as an empty heap at its index.

Logarithms in the goal are base 2 (footnote 1, p. 600), with a +1+1+1 so that one-item heaps have positive budget. The paper's constant "≤1.4404log⁡n\le 1.4404 \log n≤1.4404logn" on p. 604 is a rounding slip (1/log⁡2φ=1.44042…1/\log_2 \varphi = 1.44042\ldots1/log2​φ=1.44042…); milestones use the exact log⁡φn\log_\varphi nlogφ​n.

The goal is not trivialized: the cost counts every linking step and every cut as fixed by the step relation, CCC is quantified before the run, and a companion theorem (step_exists) shows that every valid operation can be performed, so runs containing delete min, decrease key and cascading cuts exist. A further supporting statement gives an explicit form of the delete potential bound. A complete development needs facts about the nested inductive (sizes, subtree lists, marked counts under list edits), an invariant for LEMMA 1 recording when each child was linked, and the Fibonacci inequality; contributions of any of these as separate lemmas are welcome.

Selected references

  • M. L. Fredman and R. E. Tarjan, Fibonacci heaps and their uses in improved network optimization algorithms, J. ACM 34(3):596–615, 1987. https://doi.org/10.1145/28869.28874
  • J. Vuillemin, A data structure for manipulating priority queues, Comm. ACM 21(4):309–315, 1978. https://doi.org/10.1145/359460.359478
  • R. E. Tarjan, Amortized computational complexity, SIAM J. Algebraic Discrete Methods 6(2):306–318, 1985. https://doi.org/10.1137/0606031
  • E. W. Dijkstra, A note on two problems in connexion with graphs, Numer. Math. 1:269–271, 1959. https://doi.org/10.1007/BF01386390
10 thms2 active usersReviewed
🏆Completed
Graph TheoryOperations Research·Captain: mikedeng1

Matching, Euler Tours and the Chinese Postman II: The Parity Problem Is a Minimum 1-Matching of the Odd Nodes under Shortest-Path LengthsResearch Paper

Motivation

A postal route that must traverse every street and return to its starting point is a Chinese postman tour. In a street network, crossings are nodes, streets are edges, and the cost of a route is the sum of the lengths of the streets it traverses, counted with repetition. The practical question is which streets must be repeated, and how often. Edmonds and Johnson's 1973 paper separates that question from constructing the resulting tour: choose the repeated edges first, then find an Euler tour of the resulting multigraph. This mission concerns their reduction of the first question to a minimum-weight matching of the odd-degree nodes.

The reduction matters because it turns a route problem with arbitrary repeated edge traversals into a finite pairing problem. The graph of odd nodes is complete, with the weight of a pair given by its shortest-path distance in the original network. A matching chooses which odd nodes to pair, and paths realizing those pairs identify edges to repeat. The paper also explains why overlapping paths can be replaced so that no edge needs to be repeated by this construction more than once (§3, pp. 92–93).

Setting

A graph G=(N,E)G=(N,E)G=(N,E) has finite node and edge sets. Each edge has two distinct endpoints. Different edges may join the same pair of nodes, so parallel streets remain distinct. The graph is connected when every pair of nodes is joined by an edge-labelled walk. A walk records both its node sequence and its edge sequence; a tour is a closed walk. An Euler tour uses each edge exactly once, while a postman tour uses each edge at least once. Edge lengths cec_ece​ are real and nonnegative, as in §2, p. 89.

The degree of vvv is the number of named edges meeting it. Let TTT be the set of nodes of odd degree. For every edge eee, let xe∈Nx_e\in\mathbb Nxe​∈N be the number of additional traversals. The parity equations require a nonnegative integer wvw_vwv​ at each node such that

∑e∋vxe=2wv+bv,bv={1v∈T,0v∉T.\sum_{e\ni v}x_e=2w_v+b_v,\qquad b_v=\begin{cases}1&v\in T,\\0&v\notin T.\end{cases}e∋v∑​xe​=2wv​+bv​,bv​={10​v∈T,v∈/T.​

Their objective is z(x)=∑e∈Ecexez(x)=\sum_{e\in E}c_ex_ez(x)=∑e∈E​ce​xe​. These are the paper's (3.1)–(3.4), with natural-number values encoding integrality and nonnegativity (§3, pp. 90–91). The equations express the even-degree condition after adding xex_exe​ copies of each original edge.

For nodes u,vu,vu,v, write d(u,v)d(u,v)d(u,v) for the length of a shortest path in GGG. This distance is backed by an attaining edge-simple path and is a lower bound on the length of every walk from uuu to vvv. Connectivity and nonnegative edge lengths make that specification meaningful. Let GpG_pGp​ be the complete graph on TTT, with edge {u,v}\{u,v\}{u,v} of weight d(u,v)d(u,v)d(u,v). A 1-matching MMM pairs every node of TTT with exactly one distinct node. Its length is L(M)=∑{u,v}∈Md(u,v)L(M)=\sum_{\{u,v\}\in M}d(u,v)L(M)=∑{u,v}∈M​d(u,v), equivalently half the sum of d(v,f(v))d(v,f(v))d(v,f(v)) over v∈Tv\in Tv∈T for the pairing involution fff (§3, p. 92).

Formalization targets

Odd-node matching reduction

The goal asserts the existence of a 1-matching of the odd nodes and both directions of the cost comparison:

∀M  ∃x satisfying the parity equations:xe∈{0,1}, z(x)≤L(M),\forall M\;\exists x\text{ satisfying the parity equations}:\quad x_e\in\{0,1\},\ z(x)\le L(M),∀M∃x satisfying the parity equations:xe​∈{0,1}, z(x)≤L(M), ∀x satisfying the parity equations  ∃M:L(M)≤z(x).\forall x\text{ satisfying the parity equations}\;\exists M:\quad L(M)\le z(x).∀x satisfying the parity equations∃M:L(M)≤z(x).

The existence conclusion and these comparisons identify the two attained minimum values. The goal does not prescribe a particular matching algorithm. Its supporting targets are the paper's uncrossing of overlapping shortest paths, parity of their union, reduction of multiplicities to zero or one, and decomposition of a parity solution into paths between odd nodes (§3, pp. 92–93).

Postman-tour reduction

The companion target relates a tour's edge counts to the parity equations. A postman tour traverses edge eee exactly 1+xe1+x_e1+xe​ times, and its length is

∑e∈Ece(1+xe)=∑e∈Ece+z(x).\sum_{e\in E}c_e(1+x_e) =\sum_{e\in E}c_e+z(x).e∈E∑​ce​(1+xe​)=e∈E∑​ce​+z(x).

Conversely, positive integer edge counts with even incidence at every node can be ordered into a postman tour of a connected graph. The supporting Euler theorem handles a connected even multigraph (§2, p. 89).

Significance

The matching result characterizes the optimal additional cost of an undirected postman route. It says that the parity correction can be chosen from paths pairing odd nodes, and that no feasible choice of repeated traversals can beat the best such pairing. The companion identifies how that correction becomes an actual tour. The distinction is useful when one wants to change the matching method while retaining the mathematical guarantee on route length.

A formal development would provide a reusable account of finite multigraph walks whose edge names distinguish parallel edges, of shortest-path distances with explicit attainment, and of parity subgraphs decomposed into odd-to-odd paths and even-degree residual tours. These pieces support later formal work on route inspection, TTT-joins, and matching-based graph algorithms. The paper proves the mathematical result. The exact multigraph and cost statements drafted here are open formalization targets; the local prior-art search found related results for simple graphs and metric traveling-salesman models, but no published statement with this paper's objects and conclusion.

Difficulty

A direct shortest-path construction can assign the same original edge to two matching pairs. Simply setting xe=1x_e=1xe​=1 for the union then loses the additivity of path lengths, and setting xe=2x_e=2xe​=2 loses the promised zero-or-one form. The uncrossing claim must control both the matching cost and the number of times each original edge appears. In the converse direction, an arbitrary feasible xxx may contain repeated edges and cycles; the argument must account for these without assuming that every supported edge belongs to a path between odd nodes (§3, pp. 92–93).

Formalization scope

Lean uses the published EdmondsMatching65.Polyhedron.Graph for GGG: finite node and edge types and an endpoint map into unordered node pairs. A proof field excludes loops; distinct edge names may have identical endpoints. Walks use lists of nodes and named edges with zero-based indices, and the one-node, zero-edge walk is allowed. Edge multiplicities and parity witnesses are natural numbers. Costs and shortest-path lengths are real. The distance relation requires both an attaining edge-simple path and a bound against all walks; no real infimum over a possibly empty family is used. The matching is a fixed-point-free involution on the odd nodes, extended by the identity elsewhere, and each pair contributes half of its two directed distance terms.

The goal assumes a connected graph and nonnegative edge lengths, exactly as the section does. The tour-only results require a nonempty node type because a list-based tour must have a starting node; without it an edgeless graph on no nodes would satisfy connectivity and even degree vacuously. The path-decomposition statement permits unused even-degree cycles. No independent set of “odd nodes” is accepted as input: it is computed from the degree of GGG.

The development needs finite multigraph incidence, edge-labelled walks, shortest paths, finite matchings, parity arithmetic, path uncrossing, and Euler tours. Proofs or useful intermediate results for these objects are welcome, including results applicable to multigraphs beyond the postman setting.

Selected references

  • Jack Edmonds and Ellis L. Johnson, Matching, Euler tours and the Chinese postman, Mathematical Programming 5 (1973), 88–124. DOI: 10.1007/BF01580113.
8 thms2 active usersReviewed
🏆Completed
Graph Theory·Captain: mikedeng1

Matching, Euler Tours and the Chinese Postman III: Next-Node Lists Describe an Euler Tour iff Their Last-Exit Edges Form an ArborescenceResearch Paper

Motivation

The Chinese postman problem asks for a shortest closed walk in a graph that traverses every edge at least once. Edmonds and Johnson (Math. Programming 5 (1973) 88–124) reduce it to two subproblems: choosing, by matching, which edges to duplicate so that every node has even degree, and then actually finding an Euler tour of the resulting even graph, a closed walk using every edge exactly once. Euler's theorem guarantees such a tour exists in a connected even graph, but an algorithm has to produce one and represent it compactly.

Section 5 of the paper gives two such representations and two algorithms. This mission formalizes the next-node representation: for each node nnn, a list of the nodes the tour goes to when leaving nnn. Theorem 5.1, which the paper attributes to its reference [1] (van Aardenne-Ehrenfest and de Bruijn, 1951) and to Berge's textbook, characterises exactly which lists describe an Euler tour; Theorem 5.2 shows the paper's next-node algorithm produces such lists. The same correspondence between Euler tours and arborescences is the combinatorial core of the BEST theorem on counting Euler tours of directed graphs (van Aardenne-Ehrenfest and de Bruijn, 1951).

Setting

A graph GGG is a finite set of nodes and a finite set EEE of edges, each edge meeting two distinct nodes; several edges may meet the same pair of nodes. The degree deg⁡(n)\deg(n)deg(n) counts the edges meeting nnn, and #(n,m)\#(n,m)#(n,m) counts the edges meeting both nnn and mmm. A tour is an alternating sequence (n1,e1,n2,…,el,nl+1=n1)(n_1, e_1, n_2, \dots, e_l, n_{l+1} = n_1)(n1​,e1​,n2​,…,el​,nl+1​=n1​) in which eie_iei​ meets nin_ini​ and ni+1n_{i+1}ni+1​; an Euler tour contains every edge exactly once. GGG is connected when every two nodes are joined by such a sequence, and even when every degree is even.

Next-node lists assign to each node nnn a list Ln(1),…,Ln(kn)L_n(1), \dots, L_n(k_n)Ln​(1),…,Ln​(kn​) of nodes. Given a starting node rrr, the lists specify a walk: start at rrr; the first time a node nnn is reached (for rrr, at the start), leave it to Ln(kn)L_n(k_n)Ln​(kn​); the second time, to Ln(kn−1)L_n(k_n - 1)Ln​(kn​−1); and so on, so the lists are read from the last entry to the first. The walk stops at the first node whose list is used up. The lists describe an Euler tour when this walk uses every entry of every list and its node sequence is the node sequence of an Euler tour of GGG. With parallel edges, the lists name nodes only, and the edge used at each step is chosen so that every edge is used once.

An arborescence with root rrr is a tree whose edges are directed so that every node n≠rn \ne rn=r has exactly one edge directed away from it, and rrr has none.

Formalization targets

Goal: Theorem 5.1 (p. 112)

For any node rrr of an even, connected graph GGG, next-node lists describe an Euler tour from rrr if and only if

(i) kn=12deg⁡(n)  ∀n,(ii) #(n,m)=∣{i:Lm(i)=n}∣+∣{i:Ln(i)=m}∣  ∀n,m,(iii) { n→Ln(1):n≠r } is an arborescence with root r.\text{(i)}\ k_n = \tfrac12\deg(n)\ \ \forall n,\qquad \text{(ii)}\ \#(n,m) = |\{i : L_m(i) = n\}| + |\{i : L_n(i) = m\}|\ \ \forall n, m,\qquad \text{(iii)}\ \{\,n \to L_n(1) : n \ne r\,\}\ \text{is an arborescence with root } r .(i) kn​=21​deg(n)  ∀n,(ii) #(n,m)=∣{i:Lm​(i)=n}∣+∣{i:Ln​(i)=m}∣  ∀n,m,(iii) {n→Ln​(1):n=r} is an arborescence with root r.

The printed condition (iii) names the edges (n,Ln(kn))(n, L_n(k_n))(n,Ln​(kn​)); the paper's own proofs of Theorems 5.1 and 5.2 use (n,Ln(1))(n, L_n(1))(n,Ln​(1)), the last exit from nnn, and the printed version is false on a five-edge example. The goal states Ln(1)L_n(1)Ln​(1).

Milestones

The proof on pp. 112–113 splits into five claims, each a milestone:

  1. an Euler tour's lists satisfy (i) and (ii);
  2. every node set SSS with r∉Sr \notin Sr∈/S has a last-exit edge n→Ln(1)n \to L_n(1)n→Ln​(1) leaving it;
  3. a map with one outgoing edge per node n≠rn \ne rn=r and an edge leaving every set avoiding rrr is an arborescence;
  4. under (i) and (ii) the walk from rrr returns to rrr having traversed every edge meeting rrr;
  5. at that point every node nnn incident to an untraversed edge differs from rrr and still has its last exit Ln(1)L_n(1)Ln​(1) unused.

Companion: Theorem 5.2 (p. 112)

The lists produced by any complete run of the next-node algorithm of §5.2 satisfy (i), (ii) and (iii); together with Theorem 5.1, they describe an Euler tour.

Significance

Theorem 5.1 replaces a global property, that a walk uses each edge exactly once, by three local checks: list lengths, pair counts, and a tree condition on one distinguished entry per node. That is what makes the next-node representation usable. An algorithm only has to build lists with (i)–(iii), and Theorem 5.2 shows the paper's algorithm does so for every nondeterministic choice it makes. The same arborescence condition counts Euler tours in the BEST theorem: each arborescence and each ordering of the remaining exits gives a distinct tour.

The results have been proved since 1951. What this mission proposes is a Lean formalization on multigraphs with parallel edges, which is the setting the Chinese postman reduction needs, since duplicated edges are parallel. A search of Mathlib and of the platform finds no formalization of Theorem 5.1, of next-node lists, or of arborescences in this form. Mathlib has Eulerian walks only for simple graphs, and the existing platform statements of Euler's theorem (for example SupplyChainTheory.euler_tour_iff) are about simple graphs and existence only.

Difficulty

The necessity of (i) and (ii) is counting. The substance is in two places.

First, necessity of (iii): the last exits must be shown to leave every set avoiding rrr. The tempting choice of the first exits Ln(kn)L_n(k_n)Ln​(kn​), which is what the theorem's printed text says, is wrong: on an Euler tour the first exits can form a cycle, so any argument that does not distinguish first from last exits cannot work.

Second, the converse: a walk that follows lists satisfying only (i) and (ii) can return to rrr early and leave edges untraversed, so (iii) is indispensable, and the statement depends on exactly which list entries have been used when the walk stops. Making that bookkeeping precise for a walk defined by a counter per node, and turning the used entries into an actual edge sequence in a multigraph, is the main formal work. For Theorem 5.2 the difficulty is the invariant of a nondeterministic algorithm: every trail it starts closes at its starting node because the graph is even.

Formalization scope

All declarations live in the namespace ChinesePostman.NextNode. A graph is a local ChinesePostman.NextNode.Graph V E structure on finite types: ends : E → Sym2 V and a proof that no edge is a loop, so parallel edges are distinct elements of E with the same ends. Walks are pairs of a node list and an edge list. A list LnL_nLn​ is a Lean List V, so Ln(i)L_n(i)Ln​(i) is the entry at index i−1i-1i−1 and knk_nkn​ is the length. The walk specified by the lists is computed by follow L r, a fuel-bounded iteration whose state is the current node and the number of departures made from each node. It returns the node sequence and the final counters. DescribesEulerTour G r L requires every counter to reach the list length and some edge list to make the node sequence an Euler tour. Condition (ii) is quantified over all pairs including n=mn = mn=m, where it forbids a node in its own list. Condition (iii) requires every list of a node n≠rn \ne rn=r to be nonempty, its chosen exit to be an edge of GGG, and requires the iterates of the parent map n↦Ln(1)n \mapsto L_n(1)n↦Ln​(1) (with r↦rr \mapsto rr↦r) to reach rrr from every node. The next-node algorithm is an inductive step relation on states, and "produced" means reachable termination from the Step 0 state.

Added hypotheses: Theorem 5.2 and its corollary assume an edge e0e_0e0​ meeting rrr, which Step 0 needs; a one-node graph without edges satisfies Theorem 5.1 trivially and is not excluded. Milestones drop the hypotheses their claims do not use (connectivity, evenness), which makes them stronger.

A trivializing formalization is ruled out explicitly: the lists are read backwards, from Ln(kn)L_n(k_n)Ln​(kn​) to Ln(1)L_n(1)Ln​(1), and must all be used up. Reading them forwards, dropping the "uses every entry" clause, or stating (iii) with an empty-list default for the parent would each make a different, and in some cases false, statement.

The development needs walks on multigraphs, counting arguments for closed walks (in-degree equals out-degree along the walk), and a finite pigeonhole/orbit argument for the arborescence. The arborescence lemma (milestone 3) is graph-free and reusable. Proofs of any milestone, and of the multigraph walk lemmas they need, are welcome.

Selected references

  • J. Edmonds and E. L. Johnson, Matching, Euler tours and the Chinese postman, Mathematical Programming 5 (1973) 88–124. https://doi.org/10.1007/BF01580113
  • T. van Aardenne-Ehrenfest and N. G. de Bruijn, Circuits and trees in oriented linear graphs, Simon Stevin 28 (1951) 203–217 (reference [1] of the paper, cited there as "Circuits and trees in oriented graphs").
  • C. Berge, The Theory of Graphs and its Applications, Wiley, New York, 1962 (reference [3] of the paper).
8 thms2 active usersReviewed
🏆Completed
Optimization·Captain: mikedeng1

Minimum Partition of a Matroid Into Independent Subsets: A Matroid Splits Into k Independent Sets Iff Every Subset A Has |A| ≤ k·r(A)Research Paper

Motivation

A matroid abstracts the dependence structure of the columns of a matrix: which sets of columns are linearly independent, forgetting the numbers themselves. Graphs are the special case of incidence matrices over the integers mod 2, where a set of edges is independent exactly when it contains no cycle. A natural combinatorial question in this setting is the minimum partition (or covering) problem: how few independent sets are needed to cover all the elements? For a graph this is the least number of forests that cover its edges, its arboricity.

J. Edmonds answered this question for every finite matroid in Minimum partition of a matroid into independent subsets (J. Res. NBS 69B, 1965, doi:10.6028/jres.069b.004). The answer is a good characterization in Edmonds's sense: there is a short certificate both that a partition into kkk independent sets exists (the partition itself) and that it does not (a single subset that is too large for its rank). The paper is also the source of the matroid partition algorithm, which together with matroid intersection became a basic tool of combinatorial optimization.

Timeline:

  • 1935 — H. Whitney introduces matroids through rank and independence postulates (doi:10.2307/2371182).
  • 1961–1964 — C. St. J. A. Nash-Williams and W. T. Tutte characterize graphs whose edges split into kkk forests, and graphs containing kkk edge-disjoint spanning trees (doi:10.1112/jlms/s1-39.1.12).
  • 1965 — Edmonds proves the partition theorem for all matroids, with Nash-Williams's forest theorem as a corollary; the companion paper "Lehman's switching game and a theorem of Tutte and Nash-Williams" treats the packing side.

Setting

Let MMM be a finite set of elements and call some subsets of MMM independent. Edmonds's two axioms are:

  • Axiom 1. Every subset of an independent set is independent. (A system satisfying this is an independence system.)
  • Axiom 2. For any subset AAA of the elements, all maximal independent sets contained in AAA contain the same number of elements.

A matroid is a system satisfying both. The rank r(A)r(A)r(A) of A⊆MA\subseteq MA⊆M is the largest cardinality of an independent subset of AAA; in a matroid every maximal independent subset of AAA has exactly r(A)r(A)r(A) elements.

A circuit is a minimal dependent set. A span (closed set) is a set SSS such that no circuit has exactly one element outside SSS, and the span S(A)S(A)S(A) of AAA is the minimal span containing AAA. A partition into kkk independent sets is a family I1,…,IkI_1,\dots,I_kI1​,…,Ik​ of mutually disjoint independent sets whose union is MMM; any number of the IiI_iIi​ may be empty.

Formalization targets

Goal: THEOREM 1 (p. 69)

For every finite matroid MMM and every integer k≥0k\ge 0k≥0, the elements of MMM can be partitioned into kkk independent sets if and only if there is no A⊆MA\subseteq MA⊆M with

∣A∣>k⋅r(A).|A|>k\cdot r(A).∣A∣>k⋅r(A).

The condition ranges over all subsets AAA of MMM, as printed.

Milestones

The milestones follow the paper's own lemmas and the claims of its main proof, in order.

  1. PROPOSITION 1. Axioms 1 and 2 are equivalent to Axioms 1 and 2', where Axiom 2' says that the union of an independent set and one element contains at most one circuit.
  2. PROPOSITION 2. Axioms 1 and 2' are equivalent to the circuit axioms 1c (no circuit contains another) and 2c (if distinct circuits C1,C2C_1,C_2C1​,C2​ share eee, then C1∪C2−eC_1\cup C_2-eC1​∪C2​−e contains a circuit).
  3. PROPOSITION 3 (Lehman). If e∈C1∩C2e\in C_1\cap C_2e∈C1​∩C2​ and a∈C1−C2a\in C_1-C_2a∈C1​−C2​ for circuits C1,C2C_1,C_2C1​,C2​, some circuit CCC satisfies a∈C⊆C1∪C2−ea\in C\subseteq C_1\cup C_2-ea∈C⊆C1​∪C2​−e.
  4. PROPOSITION 4. e∈S(A)e\in S(A)e∈S(A) if and only if e∈Ae\in Ae∈A or some circuit CCC has C−A={e}C-A=\{e\}C−A={e}.
  5. PROPOSITION 5. S(A)S(A)S(A) is the unique maximal set containing AAA with rank r(A)r(A)r(A); in particular r(S(I))=∣I∣r(S(I))=|I|r(S(I))=∣I∣ for independent III.
  6. "Only if" (§1.3). For any independence system, a partition into kkk independent sets gives ∣A∣≤k⋅r(A)|A|\le k\cdot r(A)∣A∣≤k⋅r(A) for all AAA.
  7. Counting step (§1.6). If xxx lies in none of kkk disjoint sets IiI_iIi​, and x∈Sx\in Sx∈S with ∣S∣≤k r(S)|S|\le k\,r(S)∣S∣≤kr(S), then some IiI_iIi​ has ∣Ii∩S∣<r(S)|I_i\cap S|<r(S)∣Ii​∩S∣<r(S).
  8. Rank drop (§1.6). For a span SSS and an independent III with ∣I∩S∣<r(S)|I\cap S|<r(S)∣I∩S∣<r(S): S(I∩S)⊆SS(I\cap S)\subseteq SS(I∩S)⊆S and r(S(I∩S))<r(S)r(S(I\cap S))<r(S)r(S(I∩S))<r(S).
  9. One-element augmentation (§1.6). If ∣S∣≤k r(S)|S|\le k\,r(S)∣S∣≤kr(S) for every span SSS, any kkk disjoint independent sets missing xxx can be rearranged into kkk disjoint independent sets covering exactly the old elements and xxx.
  10. "If" from spans (§1.6). If ∣S∣≤k⋅r(S)|S|\le k\cdot r(S)∣S∣≤k⋅r(S) holds for every span SSS, then MMM partitions into kkk independent sets.

A further item states the paper's COROLLARY (Nash-Williams): the edges of a loopless finite graph GGG can be colored with kkk colors so that no circuit is one color if and only if no nonempty node set UUU has ∣EU∣>k(∣U∣−1)|E_U|>k(|U|-1)∣EU​∣>k(∣U∣−1), where EUE_UEU​ is the set of edges with both ends in UUU.

Significance

THEOREM 1 is the base case of matroid union: the sets coverable by kkk independent sets form a matroid whose rank is min⁡A(∣M∖A∣+k r(A))\min_{A}\bigl(|M\setminus A|+k\,r(A)\bigr)minA​(∣M∖A∣+kr(A)), and the paper's proof is the origin of the polynomial-time matroid partition algorithm and of Edmonds's later matroid intersection theory. Its graphic case gives the arboricity formula of Nash-Williams, used for sparse-graph decompositions and orientation problems. The criterion is a certificate: a single set AAA with ∣A∣>k r(A)|A|>k\,r(A)∣A∣>kr(A) proves that no partition exists.

The theorem has been proved since 1965 and appears in every matroid textbook. As far as the platform and Mathlib record, it has no machine-checked proof: Mathlib has matroids, circuits, closure and strong circuit elimination for its own Matroid structure, but no matroid union or partition theorem. This mission produces a formal proof from Edmonds's own axioms (Axiom 2 as an equal-cardinality condition, not an exchange axiom), together with the equivalences between his axiom systems, which are themselves standard results.

Difficulty

The "only if" direction is a three-line count. The difficulty lies in the converse: a greedy procedure that fills the kkk sets one after another fails even for matroids, so an element xxx that fits in no current set must be accommodated by moving other elements between sets. Such a chain of exchanges must keep all kkk sets independent and disjoint at every step and must terminate; stating it with sets whose contents change while their labels persist is the part that needs careful bookkeeping in a formal proof. A second, smaller difficulty is the translation between the paper's equal-cardinality Axiom 2 and the circuit and span language its proof uses (PROPOSITIONS 1–5).

Formalization scope

  • Elements. A finite type α\alphaα with decidable equality; MMM is Finset.univ, subsets are Finset α, and the independent sets are a predicate Indep : Finset α → Prop.
  • Axiom 1 and the rank are the published Whitney definitions IndepI1 and rankOfIndep (WhitneyMatroid.RankIndep.Postulates). The rank is integer-valued, and every inequality ∣A∣|A|∣A∣ versus k r(A)k\,r(A)kr(A) is compared in Z\mathbb ZZ with k∈Nk\in\mathbb Nk∈N.
  • Axiom 2 is stated exactly as printed: maximal independent subsets of each AAA have equal cardinality. It is not replaced by an augmentation axiom or by Mathlib's Matroid.
  • Added hypothesis. Every matroid statement assumes that ∅\emptyset∅ is independent. The paper's rank presupposes it; for the empty family Axioms 1 and 2 hold vacuously and THEOREM 1 fails (an empty ground set and k=1k=1k=1).
  • Spans are formalized from the paper's words, ∣C∖S∣≠1|C\setminus S|\ne1∣C∖S∣=1 for every circuit CCC; the formula printed beside them, ∣S∩C∣≠1|S\cap C|\ne1∣S∩C∣=1, is a misprint. S(A)S(A)S(A) is the intersection of all spans containing AAA.
  • Partitions are indexed families Fin k → Finset α of pairwise disjoint independent sets covering MMM, with empty parts allowed and k=0k=0k=0 included.
  • The COROLLARY uses a graph given by ends : E → Sym2 V (parallel edges allowed, loops excluded), with independence as linear independence of incidence columns over Z/2Z\mathbb Z/2\mathbb ZZ/2Z. It quantifies over nonempty UUU: as printed, U=∅U=\emptysetU=∅ makes the statement false.
  • Not a trivialization. The goal mentions neither spans nor circuits and quantifies over all subsets AAA; the spans-only condition is a separate milestone, and the augmentation step requires the new family to cover exactly the old elements plus xxx.

Contributions welcome: proofs of the axiom equivalences (reusable for any work with Edmonds-style or circuit-axiom matroids), the span lemmas, the augmentation step, and bridges to Mathlib's Matroid that would let its circuit and closure library be reused.

Selected references

  • J. Edmonds, Minimum partition of a matroid into independent subsets, J. Res. Nat. Bur. Standards 69B (1965), 67–72. doi:10.6028/jres.069b.004
  • J. Edmonds, Lehman's switching game and a theorem of Tutte and Nash-Williams, J. Res. Nat. Bur. Standards 69B (1965), 73–77. doi:10.6028/jres.069b.005
  • C. St. J. A. Nash-Williams, Decomposition of finite graphs into forests, J. London Math. Soc. 39 (1964), 12. doi:10.1112/jlms/s1-39.1.12
  • H. Whitney, On the abstract properties of linear dependence, Amer. J. Math. 57 (1935), 509–533. doi:10.2307/2371182
17 thms2 active usersReviewed
🏆Completed
Theoretical Computer Science·Captain: wurtle

Truly Subcubic APSP on a Word RAMResearch Paper

All-pairs shortest paths (APSP) computes the shortest distance between every pair of vertices in a weighted graph. It is a fundamental problem in graph algorithms and a fundamental problem in fine-grained complexity.

This mission formalizes a deterministic O(n2.99942)O(n^{2.99942})O(n2.99942) algorithm for APSP on directed nnn-vertex graphs with polynomially bounded integer edge weights and no negative cycles, on a word RAM with O(log⁡n)O(\log n)O(logn)-bit words. This is a truly subcubic bound and refutes the integer-weight APSP hypothesis in this model.

References:

  1. Josh Alman and Virginia Vassilevska Williams, Truly Subquadratic 3SUM and Truly Subcubic APSP via Triangles in Sparse Lopsided Graphs, 2026, Theorem 22.
  2. Anthropic, formal-math / 3sum-apsp: accompanying Lean formalization, Theorem_22_APSP.
51 thms2 active usersReviewed
🏆Completed
Operations ResearchTheoretical Computer Science·Captain: mikedeng1

Approximation Algorithms for Scheduling Unrelated Parallel Machines 3: On the 3-Dimensional Matching Instances, a Schedule within ρ < 3/2 of Optimal Has Makespan ≤ 2 Iff a Matching ExistsResearch Paper

Motivation

Minimum makespan scheduling on unrelated parallel machines, written R ∣∣ Cmax⁡R\,||\,C_{\max}R∣∣Cmax​, is one of the basic models of scheduling theory: nnn jobs must each be run on one of mmm machines, and the time a job takes depends arbitrarily on the machine. Lenstra, Shmoys and Tardos (CWI Report OS-R8714, 1987; journal version Math. Programming 46, 1990) gave a polynomial 2-approximation algorithm for this problem and, in the same paper, a matching limit from below: no polynomial algorithm can guarantee a factor smaller than 3/23/23/2 unless P=NPP = NPP=NP (their Corollary 2).

The two bounds have stood for more than three decades. Closing the gap between 3/23/23/2 and 222 for R ∣∣ Cmax⁡R\,||\,C_{\max}R∣∣Cmax​ is listed among the central open problems of approximation algorithms (Williamson and Shmoys, The Design of Approximation Algorithms, 2011, open problem on unrelated machines; Schuurman and Woeginger, Polynomial time approximation algorithms for machine scheduling: ten open problems, J. Scheduling 1999). The lower bound of 3/23/23/2 is the subject of this mission.

Timeline:

  • 1979: Graham, Lawler, Lenstra and Rinnooy Kan introduce the classification R ∣∣ Cmax⁡R\,||\,C_{\max}R∣∣Cmax​.
  • 1987/1990: Lenstra, Shmoys and Tardos prove NP-completeness of deciding makespan at most 3 (Theorem 4, giving the bound 4/34/34/3) and at most 2 (Theorem 5, giving 3/23/23/2), both by reduction from 3-dimensional matching.
  • Since then: the 3/23/23/2 hardness bound has not been improved for general R ∣∣ Cmax⁡R\,||\,C_{\max}R∣∣Cmax​; progress has concentrated on special cases such as the restricted-assignment problem.

Setting

3-dimensional matching. There are three disjoint sets A={a1,…,an}A = \{a_1,\dots,a_n\}A={a1​,…,an​}, B={b1,…,bn}B = \{b_1,\dots,b_n\}B={b1​,…,bn​}, C={c1,…,cn}C = \{c_1,\dots,c_n\}C={c1​,…,cn​} and a family F={T1,…,Tm}F = \{T_1,\dots,T_m\}F={T1​,…,Tm​} of triples, each with one element of AAA, one of BBB and one of CCC. A matching is a subfamily F′F'F′ with ∣F′∣=n|F'| = n∣F′∣=n whose union is A∪B∪CA \cup B \cup CA∪B∪C. In Lean an instance is a map T:Fin m→Fin n×Fin n×Fin nT : \mathrm{Fin}\,m \to \mathrm{Fin}\,n \times \mathrm{Fin}\,n \times \mathrm{Fin}\,nT:Finm→Finn×Finn×Finn, and HasMatching T says that a matching exists.

The scheduling model. Machines are Fin m\mathrm{Fin}\,mFinm, jobs are Fin N\mathrm{Fin}\,NFinN, and pijp_{ij}pij​ is the positive integer processing time of job jjj on machine iii. A schedule σ\sigmaσ assigns each job to one machine; the load of machine iii is ∑j:σ(j)=ipij\sum_{j:\sigma(j)=i} p_{ij}∑j:σ(j)=i​pij​ and the makespan Cmax⁡(σ)C_{\max}(\sigma)Cmax​(σ) is the largest load. These are the published definitions MatousekLP.Scheduling.load and makespan.

The instance of Theorem 5. The triples containing aja_jaj​ are the triples of type jjj; let tjt_jtj​ be their number. From TTT the paper builds a scheduling instance with mmm machines, machine iii corresponding to TiT_iTi​, and with

  1. 2n2n2n element jobs, one for each bkb_kbk​ and each clc_lcl​;
  2. tj−1t_j - 1tj​−1 dummy jobs of type jjj, for each jjj.

If Ti=(aj,bk,cl)T_i = (a_j, b_k, c_l)Ti​=(aj​,bk​,cl​), machine iii processes the element jobs of bkb_kbk​ and clc_lcl​ in time 111, each dummy job of type jjj in time 222, and every other job in time 333. In Lean the processing-time matrix is P T.

A ρ\rhoρ-approximate schedule of an instance is a schedule σ\sigmaσ with Cmax⁡(σ)≤ρ Cmax⁡(τ)C_{\max}(\sigma) \le \rho\, C_{\max}(\tau)Cmax​(σ)≤ρCmax​(τ) for every schedule τ\tauτ of the same instance.

Formalization targets

Goal: Corollary 2 (p. 8)

For every instance TTT, every ρ<3/2\rho < 3/2ρ<3/2, and every ρ\rhoρ-approximate schedule σ\sigmaσ of the instance of Theorem 5 built from TTT,

Cmax⁡(σ)≤2  ⟺  T has a matching.C_{\max}(\sigma) \le 2 \iff T \text{ has a matching}.Cmax​(σ)≤2⟺T has a matching.

This is the mathematical content of "for every ρ<3/2\rho < 3/2ρ<3/2 there is no polynomial ρ\rhoρ-approximation algorithm unless P=NPP = NPP=NP": a ρ\rhoρ-approximate answer on the reduced instance decides 3-dimensional matching.

Theorem 5 (p. 7)

For every instance TTT,

∃ σ: Cmax⁡(σ)≤2  ⟺  T has a matching.\exists\,\sigma:\ C_{\max}(\sigma) \le 2 \iff T \text{ has a matching}.∃σ: Cmax​(σ)≤2⟺T has a matching.

Milestones from the proof of Theorem 5 (pp. 7–8)

  1. If every tj≥1t_j \ge 1tj​≥1, then ∑j(tj−1)+n=m\sum_j (t_j - 1) + n = m∑j​(tj​−1)+n=m: there are m−nm - nm−n dummy jobs.
  2. A matching yields a schedule with makespan at most 222.
  3. A schedule with makespan at most 222 yields a matching.

Significance

The result. Theorem 5 shows that R ∣∣ Cmax⁡R\,||\,C_{\max}R∣∣Cmax​ is NP-hard already when every processing time lies in {1,2,3}\{1,2,3\}{1,2,3} and the target makespan is 222. Corollary 2 turns this into the best known inapproximability bound for the problem, 3/23/23/2, which sits opposite the paper's own factor-2 algorithm. The same reduction pattern (machines as triples, dummy jobs that pin machines down) is reused in many later hardness proofs for scheduling and assignment problems.

Formalizing it. The result is proved on paper; no machine-checked version is known to exist. A formal proof requires the full combinatorial argument behind the reduction, including the counting of machines per type and the case where some element of AAA lies in no triple, which the paper does not discuss. Together with mission 1 of this series (the factor-2 algorithm), it fixes both ends of the [3/2,2][3/2, 2][3/2,2] gap in one formal library.

Difficulty

The construction is short, but the converse direction carries the weight: from an arbitrary schedule with makespan at most 222 one must recover a matching, although nothing in the schedule singles out the triples that form it. The obvious first idea, to read the matching off the machines that receive unit-time jobs, fails as it stands: a machine may receive one unit-time job or none, and the argument must exclude this for every schedule, including the degenerate instances in which some aja_jaj​ lies in no triple or triples repeat. For Corollary 2, the passage from the real-valued approximation guarantee to the threshold 222 depends on the strict inequality ρ<3/2\rho < 3/2ρ<3/2; at ρ=3/2\rho = 3/2ρ=3/2 the statement is no longer implied by Theorem 5.

Formalization scope

  • Machines are Fin m; 3DM elements are Fin n; a 3DM instance is T : Fin m → Fin n × Fin n × Fin n, so a triple may repeat. A matching is a set of nnn triple indices covering every aja_jaj​, bkb_kbk​, clc_lcl​ (the covering form of the paper's definition, which forces disjointness).
  • The jobs of the reduced instance are the sum type Fin n⊕Fin n⊕Σj Fin(tj−1)\mathrm{Fin}\,n \oplus \mathrm{Fin}\,n \oplus \Sigma_j\,\mathrm{Fin}(t_j - 1)Finn⊕Finn⊕Σj​Fin(tj​−1). They are enumerated by a fixed bijection with Fin N\mathrm{Fin}\,NFinN so that the published MatousekLP.Scheduling.makespan applies. All statements quantify over every schedule, so the choice of bijection is irrelevant.
  • tj−1t_j - 1tj​−1 is natural-number subtraction. When some tj=0t_j = 0tj​=0 (a case the paper leaves aside), type jjj has no dummy job; the reduced instance then has neither a matching nor a schedule of makespan at most 222, so Theorem 5 and Corollary 2 hold as stated, with no hypothesis tj≥1t_j \ge 1tj​≥1. Only the dummy-count milestone assumes tj≥1t_j \ge 1tj​≥1.
  • Processing times are natural numbers cast to R\mathbb RR; makespans are real.
  • Not formalized: "polynomial" (running time of an algorithm and of the reduction), membership in NP, "NP-complete" and "unless P=NPP = NPP=NP". Theorem 5 is stated as the reduction's equivalence; Corollary 2 is stated as the fact that any ρ\rhoρ-approximate schedule of the reduced instance, ρ<3/2\rho < 3/2ρ<3/2, decides 3-dimensional matching. The reduction is evidently polynomial, but this is not stated.
  • The approximation hypothesis compares σ\sigmaσ with every schedule τ\tauτ of the same reduced instance; a goal in which σ\sigmaσ is an arbitrary schedule, or is bounded against an unrelated quantity, would be a different statement. The strict inequality ρ<3/2\rho < 3/2ρ<3/2 is essential and kept.
  • Contributions welcome: proofs of the three milestones, of Theorem 5 from them, and of Corollary 2; reusable lemmas about integer-valued makespans and about sum-type job sets.

Selected references

  • J. K. Lenstra, D. B. Shmoys, É. Tardos, Approximation algorithms for scheduling unrelated parallel machines, CWI Report OS-R8714, Centre for Mathematics and Computer Science, Amsterdam, 1987 (FOCS 1987); the version cited for every theorem number in this mission.
  • 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
  • R. L. Graham, E. L. Lawler, J. K. Lenstra, A. H. G. Rinnooy Kan, Optimization and approximation in deterministic sequencing and scheduling: a survey, Annals of Discrete Mathematics 5 (1979) 287–326. https://doi.org/10.1016/S0167-5060(08)70356-X
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, Freeman, 1979 (3-dimensional matching, problem [SP1]).
  • P. Schuurman, G. J. Woeginger, Polynomial time approximation algorithms for machine scheduling: ten open problems, Journal of Scheduling 2 (1999) 203–213. https://doi.org/10.1002/(SICI)1099-1425(199909/10)2:5<203::AID-JOS26>3.0.CO;2-5
  • D. P. Williamson, D. B. Shmoys, The Design of Approximation Algorithms, Cambridge University Press, 2011. https://doi.org/10.1017/CBO9780511921735
8 thms2 active usersReviewed
🏆Completed
Operations ResearchOptimization·Captain: mikedeng1

Assortment Optimization under Variants of the Nested Logit Model 6: For δ > 1, the Powers-of-δ LP Optimum Scaled by (δ^(2γ̄+1), δ^(γ̄+1)) Is Feasible for the Full LPResearch Paper

Motivation

Assortment optimization asks which products a retailer should offer when customers choose among the offered products according to a probabilistic choice model, so as to maximize expected revenue. Under the nested logit model the products are grouped into nests: a customer first picks a nest, then a product inside it. The model is the standard relaxation of the independence of irrelevant alternatives property of the multinomial logit, and it is used throughout revenue management and transportation demand modelling.

Davis, Gallego and Topaloglu (Oper. Res. 62(2), 2014) classify the complexity of the assortment problem under four variants of the nested logit model, according to whether the dissimilarity parameters are at most one and whether a customer who picks a nest always buys there. The problem is NP-hard as soon as some dissimilarity parameter exceeds one (their Theorem 5), and also when the nests have positive no-purchase weights (Theorem 8), so for the general variant approximation is the realistic aim. §6.2 of the paper gives an approximation scheme for the most general variant: for any δ>1\delta > 1δ>1 it restricts each nest to a short list of candidate assortments indexed by the powers of δ\deltaδ, solves a small linear program, and loses at most a factor δ2γˉ+1\delta^{2\bar\gamma+1}δ2γˉ​+1 of the optimal revenue. This mission formalizes the guarantee behind that scheme, Theorem 12.

Setting

There are nests MMM and products N={1,…,n}N = \{1, \dots, n\}N={1,…,n} in each nest. Product jjj of nest iii has revenue rij≥0r_{ij} \ge 0rij​≥0 and preference weight vij>0v_{ij} > 0vij​>0; the products are ordered so that ri1≥⋯≥rinr_{i1} \ge \dots \ge r_{in}ri1​≥⋯≥rin​. Nest iii has a no-purchase weight vi0≥0v_{i0} \ge 0vi0​≥0 and a dissimilarity parameter γi>0\gamma_i > 0γi​>0, and v0≥0v_0 \ge 0v0​≥0 is the weight of leaving without choosing a nest. Offering Si⊆NS_i \subseteq NSi​⊆N in nest iii gives

Vi(Si)=vi0+∑j∈Sivij,Ri(Si)=∑j∈SirijvijVi(Si),Π(S1,…,Sm)=∑iVi(Si)γiRi(Si)v0+∑iVi(Si)γi.V_i(S_i) = v_{i0} + \sum_{j \in S_i} v_{ij}, \qquad R_i(S_i) = \frac{\sum_{j \in S_i} r_{ij} v_{ij}}{V_i(S_i)}, \qquad \Pi(S_1, \dots, S_m) = \frac{\sum_i V_i(S_i)^{\gamma_i} R_i(S_i)}{v_0 + \sum_i V_i(S_i)^{\gamma_i}}.Vi​(Si​)=vi0​+j∈Si​∑​vij​,Ri​(Si​)=Vi​(Si​)∑j∈Si​​rij​vij​​,Π(S1​,…,Sm​)=v0​+∑i​Vi​(Si​)γi​∑i​Vi​(Si​)γi​Ri​(Si​)​.

The optimal expected revenue Z∗=max⁡ΠZ^* = \max \PiZ∗=maxΠ is the optimal value of the linear program (3): minimize xxx subject to v0x≥∑iyiv_0 x \ge \sum_i y_iv0​x≥∑i​yi​ and yi≥Vi(Si)γi(Ri(Si)−x)y_i \ge V_i(S_i)^{\gamma_i}(R_i(S_i) - x)yi​≥Vi​(Si​)γi​(Ri​(Si​)−x) for every nest iii and every Si⊆NS_i \subseteq NSi​⊆N. Problem (4) keeps the second family of constraints only for a candidate collection of assortments in each nest.

Let γˉ=max⁡iγi\bar\gamma = \max_i \gamma_iγˉ​=maxi​γi​, assumed >1> 1>1 throughout §6, and fix δ>1\delta > 1δ>1. Put viL=vi0+min⁡jvijv^L_i = v_{i0} + \min_j v_{ij}viL​=vi0​+minj​vij​, viU=vi0+∑jvijv^U_i = v_{i0} + \sum_j v_{ij}viU​=vi0​+∑j​vij​, and let liLl^L_iliL​, liUl^U_iliU​ be the least integers with δl≥viL\delta^{l} \ge v^L_iδl≥viL​, δl≥viU\delta^l \ge v^U_iδl≥viU​. For each level l=liL,…,liUl = l^L_i, \dots, l^U_il=liL​,…,liU​, problem (15) maximizes ∑j∈Srijvij\sum_{j \in S} r_{ij} v_{ij}∑j∈S​rij​vij​ over the assortments with δl−1≤Vi(S)≤δl\delta^{l-1} \le V_i(S) \le \delta^lδl−1≤Vi​(S)≤δl; its value is G^il\hat G_{il}G^il​. An assortment S^il\hat S_{il}S^il​ is feasible for (15) and satisfies δ∑j∈S^ilrijvij≥G^il\delta \sum_{j \in \hat S_{il}} r_{ij} v_{ij} \ge \hat G_{il}δ∑j∈S^il​​rij​vij​≥G^il​. The candidate collection of nest iii is {S^il:l=liL,…,liU}∪{∅}\{\hat S_{il} : l = l^L_i, \dots, l^U_i\} \cup \{\emptyset\}{S^il​:l=liL​,…,liU​}∪{∅}.

Formalization targets

Goal: Theorem 12 (p. 28)

If (x^,y^)(\hat x, \hat y)(x^,y^​) is an optimal solution of problem (4) over the candidate collections {S^il}∪{∅}\{\hat S_{il}\} \cup \{\emptyset\}{S^il​}∪{∅}, then

(δ2γˉ+1x^, δγˉ+1y^) is feasible for problem (3).\big(\delta^{2\bar\gamma+1}\hat x,\ \delta^{\bar\gamma+1}\hat y\big) \text{ is feasible for problem (3).}(δ2γˉ​+1x^, δγˉ​+1y^​) is feasible for problem (3).

Milestones (Appendix A.6, pp. 53–54)

  1. x^≥0\hat x \ge 0x^≥0.
  2. Every nonempty assortment of nest iii lies in some level liL≤l≤liUl^L_i \le l \le l^U_iliL​≤l≤liU​.
  3. If δl−1≤a≤δl\delta^{l-1} \le a \le \delta^lδl−1≤a≤δl, then aγ−1≥(δl)γ−1δ−[γ−1]+a^{\gamma-1} \ge (\delta^l)^{\gamma-1}\delta^{-[\gamma-1]^+}aγ−1≥(δl)γ−1δ−[γ−1]+.
  4. Under the same hypothesis, (δl)γ−1≥δ−[1−γ]+aγ−1(\delta^l)^{\gamma-1} \ge \delta^{-[1-\gamma]^+} a^{\gamma-1}(δl)γ−1≥δ−[1−γ]+aγ−1.
  5. δγˉδ−[γi−1]+δ−[1−γi]+≥1\delta^{\bar\gamma}\delta^{-[\gamma_i-1]^+}\delta^{-[1-\gamma_i]^+} \ge 1δγˉ​δ−[γi​−1]+δ−[1−γi​]+≥1 and δγˉ+γi+1≤δ2γˉ+1\delta^{\bar\gamma+\gamma_i+1} \le \delta^{2\bar\gamma+1}δγˉ​+γi​+1≤δ2γˉ​+1.
  6. The scaled pair satisfies δγˉ+1y^i≥Vi(Si)γi(Ri(Si)−δ2γˉ+1x^)\delta^{\bar\gamma+1}\hat y_i \ge V_i(S_i)^{\gamma_i}(R_i(S_i) - \delta^{2\bar\gamma+1}\hat x)δγˉ​+1y^​i​≥Vi​(Si​)γi​(Ri​(Si​)−δ2γˉ​+1x^) for every nonempty SiS_iSi​.
  7. The same inequality for Si=∅S_i = \emptysetSi​=∅.

Companions

  • The guarantee stated after Theorem 12: with v0>0v_0 > 0v0​>0, the assortment assembled from the candidates solving problem (5) earns at least Z∗/δ2γˉ+1Z^*/\delta^{2\bar\gamma+1}Z∗/δ2γˉ​+1.
  • Proposition 15 (p. 51): when (15) is feasible, one of the explicit assortments S^(JL,JS)\hat S(J_L, J_S)S^(JL​,JS​), built from at most q=⌈δ/(δ−1)⌉q = \lceil \delta/(\delta-1)\rceilq=⌈δ/(δ−1)⌉ large and at most qqq small products by a greedy continuous knapsack, is feasible for (15) and within a factor δ\deltaδ of G^il\hat G_{il}G^il​.
  • The count liU−liL≤1+log⁡δ(viU/viL)l^U_i - l^L_i \le 1 + \log_\delta(v^U_i/v^L_i)liU​−liL​≤1+logδ​(viU​/viL​) (p. 28).

Significance

Theorem 12 is the analytical half of the approximation scheme. Together with the paper's Theorem 1 it shows that a linear program with 1+m1 + m1+m variables and at most 1+m(2+log⁡δ(viU/viL))1 + m(2 + \log_\delta(v^U_i/v^L_i))1+m(2+logδ​(viU​/viL​)) constraints yields an assortment within a factor δ2γˉ+1\delta^{2\bar\gamma+1}δ2γˉ​+1 of the optimum, for the most general variant, which is NP-hard, and Proposition 15 makes the candidate assortments computable. Letting δ↓1\delta \downarrow 1δ↓1 trades accuracy for running time, so the result is the paper's answer to how well the general problem can be approximated by this LP approach.

The result is proved in the paper; the work here is to formalize the known proof. To the best of a search of the Prove2Me library (local index and platform mirror, October 2026), no nested logit approximation result, and no knapsack lemma matching Proposition 15, has a machine-checked statement or proof. The definitions of the shared model (the instance, ViV_iVi​, RiR_iRi​, Π\PiΠ, the linear programs (3) and (4)) are common to the six missions of this series.

Difficulty

The obvious argument compares an arbitrary assortment SiS_iSi​ with the candidate of its level and uses the candidate's constraint in (4). This fails as a direct comparison because the nest weight Vi(⋅)γiV_i(\cdot)^{\gamma_i}Vi​(⋅)γi​ enters twice with different exponents, γi−1\gamma_i - 1γi​−1 in front of the revenue and γi\gamma_iγi​ in front of x^\hat xx^, and within one level ViV_iVi​ may vary by a factor δ\deltaδ. When γi>1\gamma_i > 1γi​>1 and when γi≤1\gamma_i \le 1γi​≤1 the monotonicity of t↦tγi−1t \mapsto t^{\gamma_i - 1}t↦tγi​−1 goes in opposite directions, so the losses must be tracked separately by [γi−1]+[\gamma_i - 1]^+[γi​−1]+ and [1−γi]+[1 - \gamma_i]^+[1−γi​]+, and they are absorbed by the common factor δγˉ\delta^{\bar\gamma}δγˉ​ only because γˉ>1\bar\gamma > 1γˉ​>1. The other half, Proposition 15, concerns a knapsack with both a lower and an upper bound on the total weight, where a feasible solution must be produced as well as a good objective value.

Formalization scope

Products are Fin n (0,…,n−10, \dots, n-10,…,n−1 for 1,…,n1, \dots, n1,…,n), nests a finite type, powers VγV^{\gamma}Vγ real powers, and x/0=0x/0 = 0x/0=0, which gives Ri(∅)=0R_i(\emptyset) = 0Ri​(∅)=0. The shared model carries the standing assumptions of §1 with the disclosed pins vij>0v_{ij} > 0vij​>0, rij≥0r_{ij} \ge 0rij​≥0 and γi>0\gamma_i > 0γi​>0 (the page allows γi=0\gamma_i = 0γi​=0 and zero-weight padding products, under which its convention 0γi=00^{\gamma_i} = 00γi​=0 and its proofs fail). Every statement of the mission also assumes n≥1n \ge 1n≥1 (for viLv^L_iviL​), γˉ\bar\gammaγˉ​ is the greatest of the γi\gamma_iγi​ (so at least one nest exists) with γˉ>1\bar\gamma > 1γˉ​>1, and δ>1\delta > 1δ>1. Integer powers δl\delta^lδl are zpow, and liLl^L_iliL​, liUl^U_iliU​ are written ⌈log⁡δviL⌉\lceil \log_\delta v^L_i\rceil⌈logδ​viL​⌉, ⌈log⁡δviU⌉\lceil \log_\delta v^U_i\rceil⌈logδ​viU​⌉, which equal the minima of the paper. An optimal solution of (4) is a feasible pair whose xxx is minimal among feasible pairs. The companion guarantee adds v0>0v_0 > 0v0​>0, the pin of Theorem 1. In Proposition 15 the capacity row of (33) includes v0v_0v0​, correcting a printed slip, and the running-time claim is not stated.

The assortments S^il\hat S_{il}S^il​ enter Theorem 12 as an arbitrary family with the two properties of p. 28; the theorem is stated for every such family. A level at which (15) has no feasible assortment imposes nothing on S^il\hat S_{il}S^il​, as on the page. Requiring S^il\hat S_{il}S^il​ to lie in its level unconditionally would make the hypothesis unsatisfiable for most instances and the goal vacuous; that encoding, and any goal that assumes displays (34)–(36) or the exponent bounds, is ruled out.

Contributions welcome: proofs of the real-power milestones (3)–(5), which are self-contained, of the two cases (6)–(7), and of Proposition 15, whose fractional knapsack lemma (the greedy solution of a continuous knapsack sorted by ratio is optimal) is reusable beyond this mission.

Selected references

  • J. M. Davis, G. Gallego, H. Topaloglu, Assortment optimization under variants of the nested logit model, Operations Research 62(2), 2014; cited from the revised manuscript of June 18, 2013. https://doi.org/10.1287/opre.2014.1256
  • A. M. Frieze, M. R. B. Clarke, Approximation algorithms for the m-dimensional 0–1 knapsack problem: worst-case and probabilistic analyses, European Journal of Operational Research 15(1), 1984. (cited on p. 50 of the paper for the continuous knapsack (33); link not verified here)
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, W. H. Freeman, 1979.
11 thms2 active usersReviewed
🏆Completed
Operations ResearchTheoretical Computer Science·Captain: mikedeng1

Approximation Algorithms for Scheduling Unrelated Parallel Machines 4: With Times in {p, q} and gcd(p, q) = 1, a q-Dimensional Matching Exists Iff a Schedule Has Makespan ≤ pqResearch Paper

Motivation

Minimum makespan scheduling on unrelated parallel machines asks for an assignment of nnn jobs to mmm machines, where job jjj takes pijp_{ij}pij​ time units on machine iii, so that the largest machine load is as small as possible. In the three-field notation of Graham, Lawler, Lenstra and Rinnooy Kan (1979) it is R ∥ Cmax⁡R\,\|\,C_{\max}R∥Cmax​. It models load balancing across heterogeneous processors, workers or production lines, and it is a standard test case for linear-programming rounding.

Lenstra, Shmoys and Tardos (FOCS 1987; CWI Report OS-R8714; Mathematical Programming 46, 1990) gave a polynomial 2-approximation algorithm for R ∥ Cmax⁡R\,\|\,C_{\max}R∥Cmax​ and showed that no polynomial algorithm achieves a ratio below 3/23/23/2 unless P = NP. They also asked which restrictions on the processing times keep the problem hard. Their Section 5 answers this when only two distinct processing times occur:

  • all pij=1p_{ij}=1pij​=1: trivial;
  • all pij∈{1,∞}p_{ij}\in\{1,\infty\}pij​∈{1,∞}: bipartite cardinality matching;
  • all pij∈{1,2}p_{ij}\in\{1,2\}pij​∈{1,2}: polynomial by matching techniques (Theorem 6);
  • all pij∈{p,q}p_{ij}\in\{p,q\}pij​∈{p,q} with p<qp<qp<q, 2p≠q2p\ne q2p=q: NP-hard (Theorem 7).

Theorem 7 is the paper's last result and closes this classification. Its proof generalizes the reduction of Theorem 4, which handles the case {1,3}\{1,3\}{1,3}, from 3-dimensional matching to qqq-dimensional matching.

Setting

Machines are indexed by i∈{1,…,m}i\in\{1,\dots,m\}i∈{1,…,m} and jobs by jjj. A schedule σ\sigmaσ assigns every job to exactly one machine. The load of machine iii is ∑j:σ(j)=ipij\sum_{j:\sigma(j)=i}p_{ij}∑j:σ(j)=i​pij​, and the makespan Cmax⁡(σ)C_{\max}(\sigma)Cmax​(σ) is the largest load.

qqq-dimensional matching. An instance is a ground set UUU of qnqnqn elements and a family S1,…,SmS_1,\dots,S_mS1​,…,Sm​ of qqq-element subsets of UUU. A matching is a subfamily F′⊆{1,…,m}F'\subseteq\{1,\dots,m\}F′⊆{1,…,m} with ∣F′∣=n|F'|=n∣F′∣=n and ⋃i∈F′Si=U\bigcup_{i\in F'}S_i=U⋃i∈F′​Si​=U; its members are then pairwise disjoint.

The instance of Theorem 7. Fix natural numbers 0<p<q0<p<q0<p<q that are relatively prime. Build a scheduling instance with mmm machines, machine iii corresponding to SiS_iSi​, and two kinds of jobs:

  • qnqnqn element jobs, one per u∈Uu\in Uu∈U, with piu=pp_{iu}=ppiu​=p if u∈Siu\in S_iu∈Si​ and piu=qp_{iu}=qpiu​=q otherwise;
  • p(m−n)p(m-n)p(m−n) dummy jobs, each taking qqq time units on every machine.

Every processing time lies in {p,q}\{p,q\}{p,q}.

The instance of Theorem 4. For disjoint sets A={a1,…,an}A=\{a_1,\dots,a_n\}A={a1​,…,an​}, BBB, CCC of the same size and triples Ti=(aj,bk,cl)T_i=(a_j,b_k,c_l)Ti​=(aj​,bk​,cl​), i=1,…,mi=1,\dots,mi=1,…,m, there are 3n3n3n element jobs, one per element of A∪B∪CA\cup B\cup CA∪B∪C, and m−nm-nm−n dummy jobs. Machine iii processes the element jobs of aja_jaj​, bkb_kbk​, clc_lcl​ in one time unit and every other job in three time units. A 3-dimensional matching is a subfamily of nnn triples covering A∪B∪CA\cup B\cup CA∪B∪C.

Formalization targets

Goal: Theorem 7 (p. 8; proof pp. 8–9)

For relatively prime 0<p<q0<p<q0<p<q and a family of qqq-subsets S1,…,SmS_1,\dots,S_mS1​,…,Sm​ of a qnqnqn-element set, the instance above has all processing times in {p,q}\{p,q\}{p,q}, and

∃ σ: Cmax⁡(σ)≤pq⟺∃ F′⊆{1,…,m}: ∣F′∣=n, ⋃i∈F′Si=U.\exists\,\sigma:\ C_{\max}(\sigma)\le pq \quad\Longleftrightarrow\quad \exists\,F'\subseteq\{1,\dots,m\}:\ |F'|=n,\ \bigcup_{i\in F'}S_i=U .∃σ: Cmax​(σ)≤pq⟺∃F′⊆{1,…,m}: ∣F′∣=n, i∈F′⋃​Si​=U.

Milestones

  1. Theorem 4 (p. 7): on the 3-dimensional matching instance with times in {1,3}\{1,3\}{1,3}, a schedule with makespan at most 333 exists iff a matching exists.
  2. Matching gives a schedule (pp. 8–9): if SSS has a matching, some schedule has Cmax⁡≤pqC_{\max}\le pqCmax​≤pq.
  3. No idle time, two kinds of machine (p. 9): in every schedule with Cmax⁡≤pqC_{\max}\le pqCmax​≤pq, three things hold. Every load equals pqpqpq. Every element job runs at length ppp, on a machine whose tuple contains it. Each machine processes either exactly qqq element jobs and no dummy job, or exactly ppp dummy jobs and no element job.

Significance

The result. Theorems 6 and 7 together say exactly which two-valued restrictions of R ∥ Cmax⁡R\,\|\,C_{\max}R∥Cmax​ are tractable. Up to scaling, the polynomial case is {1,2}\{1,2\}{1,2}, and every other pair {p,q}\{p,q\}{p,q} with p<qp<qp<q is NP-hard. Theorem 4 is the base case. It also yields Corollary 1 of the paper: no polynomial ρ\rhoρ-approximation with ρ<4/3\rho<4/3ρ<4/3 exists unless P = NP. Reductions of this "no idle time" kind are a standard template for hardness of restricted-assignment and two-value scheduling, a line still active in the study of the restricted assignment problem and of the "graph balancing" special case.

Formalizing it. The result is proved, in two short paragraphs. The proof leaves to the reader the construction of the schedule from a matching and the "easy number theoretic argument" behind the converse, and both become explicit here. The development also produces reusable pieces: a qqq-dimensional matching definition over an arbitrary family of qqq-sets, and two reduction instances built on the published load and makespan of MatousekLP.Scheduling.Schedule. To our knowledge no machine-checked proof of these reductions exists.

Difficulty

The forward direction is a direct construction. The converse is where care is needed. A schedule with makespan at most pqpqpq may a priori place an element job on a machine whose tuple does not contain it, at length qqq, and may mix element and dummy jobs on one machine. Looking at one machine at a time does not exclude either: a single machine can carry such a mixed load below pqpqpq. What excludes them is a property of the whole schedule together with the arithmetic of ppp and qqq. The hypotheses are sharp for this step: for p>qp>qp>q element jobs are cheaper on foreign machines, and for gcd⁡(p,q)>1\gcd(p,q)>1gcd(p,q)>1 a load ap+bq=pqap+bq=pqap+bq=pq with a,b>0a,b>0a,b>0 becomes possible.

Formalization scope

  • Model. Machines are Fin m. Schedules are maps from jobs to machines. Load and makespan are those of the published MatousekLP.Scheduling.Schedule, with natural-number processing times cast to R\mathbb RR. The jobs are enumerated as Fin (q * n + p * (m - n)) (Theorem 7) and Fin (3 * n + (m - n)) (Theorem 4), element jobs first, through finSumFinEquiv.
  • Ground set and family. The ground set is Fin (q * n) and the family is indexed by machines, so repeated tuples are allowed. The qqq-partite structure of qqq-dimensional matching is not imposed: the reduction does not use it, and statements over all families of qqq-sets contain the qqq-partite case. Theorem 4 keeps the tripartite structure, with triples of indices in Fin n × Fin n × Fin n.
  • Threshold. The thresholds are exactly pqpqpq and 333, with "makespan at most".
  • Complexity wording not formalized. "NP-hard" (Theorem 7) and "NP-complete" (Theorem 4) are not formalized: membership in NP, polynomial size of the reductions, and the hardness of the matching problems are not stated. What is stated is the equivalence each reduction establishes.
  • Hypotheses. Theorem 7's goal assumes 0<p<q0<p<q0<p<q, gcd⁡(p,q)=1\gcd(p,q)=1gcd(p,q)=1 and ∣Si∣=q|S_i|=q∣Si​∣=q. The paper's general case gcd⁡(p,q)=g>1\gcd(p,q)=g>1gcd(p,q)=g>1 is its stated "without loss of generality": divide all times by ggg, which divides every makespan by ggg. It is not part of the formal statement. The paper's hypothesis 2p≠q2p\ne q2p=q is dropped because the reduction does not use it. With coprime p<qp<qp<q it excludes only (p,q)=(1,2)(p,q)=(1,2)(p,q)=(1,2), where the equivalence still holds but the source problem is bipartite matching. The formal goal is therefore stronger than the paper's.
  • Edge cases. For m<nm<nm<n, natural-number subtraction gives no dummy jobs; both sides of each equivalence are false when n≥1n\ge1n≥1. This stands in for the paper's "trivial 'no' instance". For n=0n=0n=0 the empty family is a matching and the dummy jobs fill the machines exactly.
  • No trivialization. The goal is an equivalence about a constructed instance whose processing times are given by explicit definitions. Neither side is assumed, the schedule is quantified over all maps from jobs to machines, and the instance is not a free parameter pinned by hypotheses.
  • Welcome contributions. Besides the milestones: counting lemmas relating a machine's load to the numbers of jobs of each length on it, and lemmas on the makespan of MatousekLP.Scheduling.Schedule (it bounds every load; it is attained when m≥1m\ge1m≥1), which are reusable across scheduling reductions.

Selected references

  • J. K. Lenstra, D. B. Shmoys, É. Tardos, Approximation algorithms for scheduling unrelated parallel machines, CWI Report OS-R8714, Amsterdam, 1987 (the version formalized here); journal version in Mathematical Programming 46 (1990) 259–271. https://doi.org/10.1007/BF01585745
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, Freeman, 1979 (3-dimensional matching, problem SP1).
  • R. L. Graham, E. L. Lawler, J. K. Lenstra, A. H. G. Rinnooy Kan, Optimization and approximation in deterministic sequencing and scheduling: a survey, Annals of Discrete Mathematics 5 (1979) 287–326. https://doi.org/10.1016/S0167-5060(08)70356-X
  • J. Matoušek, B. Gärtner, Understanding and Using Linear Programming, Springer, 2007, §8.3 (the schedule, load and makespan definitions reused here). https://doi.org/10.1007/978-3-540-30717-4
7 thms2 active usersReviewed
🏆Completed
Operations ResearchOptimizationTheoretical Computer Science·Captain: mikedeng1

Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments III: The Three Revenue-Ordered Approximation Bounds Are TightResearch Paper

Motivation

A retailer or an airline chooses which products to offer, and customers choose among what is offered, or buy nothing. Choosing the offer set that maximizes expected revenue is the assortment problem, central to revenue management (Talluri and van Ryzin, 2004). It is NP-hard even for simple mixtures of logit models, so practice relies on heuristics. The most common one is the revenue-ordered assortments strategy: offer the products whose price is above a threshold, and pick the best threshold.

Berbeglia and Joret (arXiv:1606.01371) analyse this heuristic under every regular choice model, the broad class in which enlarging the offer set never raises the probability of choosing any given alternative, including not buying. They prove three approximation guarantees, then show that none of the three can be improved. This mission formalizes that last result, Theorem 3.4.

Timeline:

  • Talluri and van Ryzin (2004) showed that revenue-ordered assortments are optimal under the multinomial logit model.
  • Rusmevichientong, Shmoys, Tong and Topaloglu (2014) showed the problem NP-hard under a mixture of two logit models, and proved that revenue-ordered assortments earn at least OPT/(e(1+ln⁡(rk/r1)))\mathrm{OPT}/(e(1+\ln(r_k/r_1)))OPT/(e(1+ln(rk​/r1​))) under mixed logit.
  • Aouad, Farias, Levi and Segev (2018) proved an Ω(1/ln⁡(rk/r1))\Omega(1/\ln(r_k/r_1))Ω(1/ln(rk​/r1​)) guarantee under random utility models, and that the assortment problem there is NP-hard to approximate within Ω(1/n1−ϵ)\Omega(1/n^{1-\epsilon})Ω(1/n1−ϵ) and Ω(1/log⁡1−ϵ(rk/r1))\Omega(1/\log^{1-\epsilon}(r_k/r_1))Ω(1/log1−ϵ(rk​/r1​)).
  • Berbeglia and Joret (2016–2020, Algorithmica) proved the guarantees 1/k1/k1/k, 1/∑i(ri−ri−1)/ri≥1/(1+ln⁡(rk/r1))1/\sum_{i}(r_i-r_{i-1})/r_i\ge1/(1+\ln(r_k/r_1))1/∑i​(ri​−ri−1​)/ri​≥1/(1+ln(rk​/r1​)) and a purchase-probability bound under any regular model, and an instance on which all three, in their sum forms, are attained in the limit.

Setting

The products form a finite set C\mathcal CC. A system of choice probabilities gives, for every offer set S⊆CS\subseteq\mathcal CS⊆C and product xxx, the probability P(x,S)\mathcal P(x,S)P(x,S) that a customer buys xxx. The no-purchase probability is P(0,S)=1−∑x∈SP(x,S)\mathcal P(0,S)=1-\sum_{x\in S}\mathcal P(x,S)P(0,S)=1−∑x∈S​P(x,S). The model is regular if

  1. P(x,S)≥0\mathcal P(x,S)\ge0P(x,S)≥0 for x∈C∪{0}x\in\mathcal C\cup\{0\}x∈C∪{0};
  2. P(x,S)=0\mathcal P(x,S)=0P(x,S)=0 for x∉Sx\notin Sx∈/S;
  3. ∑x∈SP(x,S)≤1\sum_{x\in S}\mathcal P(x,S)\le1∑x∈S​P(x,S)≤1;
  4. P(x,S)≥P(x,S′)\mathcal P(x,S)\ge\mathcal P(x,S')P(x,S)≥P(x,S′) for S⊆S′S\subseteq S'S⊆S′ and every x∈S∪{0}x\in S\cup\{0\}x∈S∪{0}.

Each product has a revenue r(x)>0r(x)>0r(x)>0. Offering SSS earns rev(S)=∑x∈SP(x,S) r(x)\mathrm{rev}(S)=\sum_{x\in S}\mathcal P(x,S)\,r(x)rev(S)=∑x∈S​P(x,S)r(x), and OPT=max⁡S⊆Crev(S)\mathrm{OPT}=\max_{S\subseteq\mathcal C}\mathrm{rev}(S)OPT=maxS⊆C​rev(S).

Let r1<⋯<rkr_1<\dots<r_kr1​<⋯<rk​ be the distinct revenues, with r0:=0r_0:=0r0​:=0, and Si={x:r(x)≥ri}S_i=\{x: r(x)\ge r_i\}Si​={x:r(x)≥ri​}. The heuristic earns

RO=max⁡i∈[k]rev(Si).\mathrm{RO}=\max_{i\in[k]}\mathrm{rev}(S_i).RO=i∈[k]max​rev(Si​).

For an optimal S∗S^*S∗ let Ni=∑x∈S∗, r(x)≥riP(x,S∗)N_i=\sum_{x\in S^*,\,r(x)\ge r_i}\mathcal P(x,S^*)Ni​=∑x∈S∗,r(x)≥ri​​P(x,S∗), Nk+1=0N_{k+1}=0Nk+1​=0, and ℓ\ellℓ the largest index with Nℓ>0N_\ell>0Nℓ​>0. Section 3 of the paper proves

  • (A) OPT≤k⋅RO\mathrm{OPT}\le k\cdot\mathrm{RO}OPT≤k⋅RO (Theorem 3.1);
  • (B) OPT≤Dr⋅RO\mathrm{OPT}\le D_r\cdot\mathrm{RO}OPT≤Dr​⋅RO, with Dr=∑i=1kri−ri−1ri≤1+ln⁡(rk/r1)D_r=\sum_{i=1}^{k}\frac{r_i-r_{i-1}}{r_i}\le 1+\ln(r_k/r_1)Dr​=∑i=1k​ri​ri​−ri−1​​≤1+ln(rk​/r1​) (Theorem 3.2);
  • (C) OPT≤DN(S∗)⋅RO\mathrm{OPT}\le D_N(S^*)\cdot\mathrm{RO}OPT≤DN​(S∗)⋅RO, with DN(S∗)=∑i=1ℓNi−Ni+1Ni≤1+ln⁡(N1/Nℓ)D_N(S^*)=\sum_{i=1}^{\ell}\frac{N_i-N_{i+1}}{N_i}\le 1+\ln(N_1/N_\ell)DN​(S∗)=∑i=1ℓ​Ni​Ni​−Ni+1​​≤1+ln(N1​/Nℓ​) (Theorem 3.3).

Formalization targets

Goal: Theorem 3.4

For every k≥1k\ge1k≥1 and every δ>0\delta>0δ>0 there are a finite nonempty product set, a regular P\mathcal PP, revenues r>0r>0r>0 with exactly kkk distinct values, and an optimal S∗S^*S∗ with N1>0N_1>0N1​>0, such that

k⋅RO<(1+δ) OPT,Dr⋅RO<(1+δ) OPT,DN(S∗)⋅RO<(1+δ) OPT.k\cdot\mathrm{RO}<(1+\delta)\,\mathrm{OPT},\qquad D_r\cdot\mathrm{RO}<(1+\delta)\,\mathrm{OPT},\qquad D_N(S^*)\cdot\mathrm{RO}<(1+\delta)\,\mathrm{OPT}.k⋅RO<(1+δ)OPT,Dr​⋅RO<(1+δ)OPT,DN​(S∗)⋅RO<(1+δ)OPT.

So none of the bounds (A), (B), (C) stays true when multiplied by 1+δ1+\delta1+δ, for any number kkk of distinct revenues.

The tight instance (milestones)

The paper's witness has products (i,j)(i,j)(i,j) with i∈[k]i\in[k]i∈[k] and j∈[i]j\in[i]j∈[i]. Product (i,j)(i,j)(i,j) has revenue ε−j\varepsilon^{-j}ε−j, and P((i,j),S)=εi\mathcal P((i,j),S)=\varepsilon^iP((i,j),S)=εi when (i,j)∈S(i,j)\in S(i,j)∈S and (i,1),…,(i,j−1)∉S(i,1),\dots,(i,j-1)\notin S(i,1),…,(i,j−1)∈/S, and 000 otherwise, for 0<ε≤120<\varepsilon\le\tfrac120<ε≤21​. The milestones follow the proof on pp. 10–11:

  • the axiom (iii) bound;
  • (7);
  • (8);
  • (9);
  • regularity;
  • the distinct revenues ri=ε−ir_i=\varepsilon^{-i}ri​=ε−i;
  • RO=rev(C)<1/(1−ε)\mathrm{RO}=\mathrm{rev}(\mathcal C)<1/(1-\varepsilon)RO=rev(C)<1/(1−ε);
  • OPT=k\mathrm{OPT}=kOPT=k, attained by {(i,i)}\{(i,i)\}{(i,i)};
  • the limits OPT/RO→k\mathrm{OPT}/\mathrm{RO}\to kOPT/RO→k and Dr→kD_r\to kDr​→k;
  • Ni=εi+⋯+εkN_i=\varepsilon^i+\dots+\varepsilon^kNi​=εi+⋯+εk and DN→kD_N\to kDN​→k as ε→0+\varepsilon\to0^+ε→0+.

Significance

With Theorems 3.1–3.3, Theorem 3.4 closes the analysis: in terms of the parameters kkk, DrD_rDr​ and DND_NDN​, revenue-ordered assortments are understood exactly under regular choice. It complements the hardness results of Aouad et al., which already show that no efficient strategy can do much better than (A) and (B) in general; Theorem 3.4 shows that the analysis of this particular heuristic is exact. The instance is also a concrete regular choice model in which the optimal offer set is far from every nested one.

On the formal side, the mission produces a reusable formal definition of regular discrete choice models, revenue-ordered assortments and the three bound quantities. These are shared, under other sub-namespaces, with the companion missions on Theorems 3.2 and 3.3. The result is proved in the paper. As far as is known it has not been machine-checked anywhere; the remaining work is formalizing the paper's proof, including the steps it leaves to the reader.

Difficulty

The proof is a construction, and the paper verifies most of it in a sentence each. The work lies in those sentences:

  • the regularity of the instance at the no-purchase option, (9), which needs the row decomposition (8);
  • the claim that {(i,i)}\{(i,i)\}{(i,i)} is optimal among all 2k(k+1)/22^{k(k+1)/2}2k(k+1)/2 offer sets, asserted without proof;
  • the claim that the full set is the best threshold set;
  • the evaluation of DN(S∗)D_N(S^*)DN​(S∗), which needs ℓ=k\ell=kℓ=k.

The naive witness, one instance per bound, does not help: the goal asks for one instance with exactly kkk revenues on which all three bounds are nearly attained, for every kkk.

Formalization scope

  • Products are an arbitrary finite type C. The no-purchase option is not a product: P(0,S)\mathcal P(0,S)P(0,S) is the derived quantity noPurchase P S. Offer sets are Finset C.
  • IsRegular carries axioms (i)–(iv), with (i) and (iv) stated both for products and for the no-purchase option. Regularity at x=0x=0x=0 is essential: without it the guarantees fail.
  • OPT\mathrm{OPT}OPT is Finset.sup' over all subsets, including ∅\emptyset∅. RO\mathrm{RO}RO is Finset.sup' over the kkk threshold sets only, and needs Nonempty C.
  • The distinct revenues are the sorted image of r, indexed by Fin k from 000: the Lean index iii is the paper's i+1i+1i+1, and r0=0r_0=0r0​=0 is a separate case. ℓ\ellℓ lives in WithBot (Fin k).
  • Bounds are multiplicative (OPT≤D⋅RO\mathrm{OPT}\le D\cdot\mathrm{RO}OPT≤D⋅RO); there is no ratio RO/OPT\mathrm{RO}/\mathrm{OPT}RO/OPT in the goal. Limits are along 𝓝[>] 0.
  • The tight instance keeps the paper's 1-based pairs (i,j)(i,j)(i,j) as a subtype of Fin (k+1) × Fin (k+1). Its regularity is a theorem to prove, never a field assumed.

Ruled out:

  • a fixed kkk, since "for every kkk" is the content;
  • tightness of the logarithmic forms 1/(1+ln⁡(rk/r1))1/(1+\ln(r_k/r_1))1/(1+ln(rk​/r1​)) and 1/(1+ln⁡ν)1/(1+\ln\nu)1/(1+lnν), which this instance does not show (there 1+ln⁡(rk/r1)→∞1+\ln(r_k/r_1)\to\infty1+ln(rk​/r1​)→∞ while OPT/RO→k\mathrm{OPT}/\mathrm{RO}\to kOPT/RO→k), so only the sum forms are claimed tight;
  • an instance whose regularity is assumed;
  • a heuristic maximizing over all subsets.

Contributions welcome: proofs of the milestones, especially (8), (9) and the optimality of {(i,i)}\{(i,i)\}{(i,i)}, and general lemmas on sorted distinct values of a finite function.

Selected references

  • G. Berbeglia and G. Joret, Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments, arXiv:1606.01371v3, 2019; published in Algorithmica, 2020. https://arxiv.org/abs/1606.01371
  • K. Talluri and G. van Ryzin, Revenue Management Under a General Discrete Choice Model of Consumer Behavior, Management Science 50(1), 2004. https://doi.org/10.1287/mnsc.1030.0147
  • P. Rusmevichientong, D. Shmoys, C. Tong and H. Topaloglu, Assortment Optimization under the Multinomial Logit Model with Random Choice Parameters, Production and Operations Management 23(11), 2014. https://doi.org/10.1111/poms.12191
  • A. Aouad, V. Farias, R. Levi and D. Segev, The Approximability of Assortment Optimization Under Ranking Preferences, Operations Research 66(6), 2018. https://doi.org/10.1287/opre.2018.1724
15 thms2 active usersReviewed
🏆Completed
Discrete GeometryGraph Theory·Captain: mikedeng1

Diameter of Polyhedra: Limits of Abstraction 3: Larman's Bound 2^(d−1)·n Holds in the Base AbstractionResearch Paper

Motivation

The diameter of a polyhedron is the largest number of edges needed to travel between two vertices of its graph, its 1-skeleton. Whether this number is bounded by a polynomial in the dimension ddd and the number of facets nnn is a long-standing open question in polyhedral combinatorics; edge walks are the moves of vertex-following methods such as the simplex method, so the question bounds what any such method can achieve in the worst case. For fixed dimension, Larman proved in 1970 that the diameter grows at most linearly in nnn, with a constant exponential in ddd (Larman 1970).

Eisenbrand, Hähnle, Razborov and Rothvoß ask how much of the geometry of polyhedra such upper bounds actually use. They isolate a single combinatorial path property shared by the graphs of non-degenerate polyhedra and announce that the known upper bounds, Larman's among them, "continue to hold in the base abstraction" defined by that property (Dagstuhl abstract, p. 3). The next sentence of the abstract reads: "While the first bound is merely an adaption of the proof in [13], our proof of the second bound is much simpler than the one which was proved for polyhedra in [16]."

The held source is the five-page Dagstuhl extended abstract dated June 15, 2010, which states the results and proves none of them. The full proofs appear in Mathematics of Operations Research 35 (2010), 786–794 (journal article). The goal is an asserted theorem of the paper, not an open conjecture; the formal proof remains to be supplied.

Setting

Let [n]={1,…,n}[n]=\{1,\ldots,n\}[n]={1,…,n} and write ([n]d)\binom{[n]}{d}(d[n]​) for the family of its ddd-element subsets. A graph G=(V,E)G=(V,E)G=(V,E) in the base abstraction has a nonempty vertex set V⊆([n]d)V\subseteq\binom{[n]}{d}V⊆(d[n]​): every vertex is a ddd-element set of labels, and distinct vertices carry distinct sets. The edges are arbitrary subject to condition i): for each u,v∈Vu,v\in Vu,v∈V there exists a path connecting uuu and vvv whose intermediate vertices all contain u∩vu\cap vu∩v. Condition i) makes GGG connected. The class of all such graphs is Bd,n\mathcal B_{d,n}Bd,n​; ddd is called the dimension and nnn the number of facets of the abstraction, and D(d,n)D(d,n)D(d,n) denotes the largest diameter of a graph in Bd,n\mathcal B_{d,n}Bd,n​ (p. 2).

The distance dist⁡G(u,v)\operatorname{dist}_G(u,v)distG​(u,v) is the least number of edges of a path from uuu to vvv. The motivating example is the graph of a non-degenerate polyhedron (each vertex on exactly ddd facets), each vertex labelled by its ddd facets: two vertices are joined by a path inside the minimal face containing both, and every vertex of that path lies on the facets the two endpoints share (pp. 2–3). The class Bd,n\mathcal B_{d,n}Bd,n​ is not restricted to such graphs.

Formalization targets

Larman's bound in the base abstraction

The abstract prints the bound as Δu(d,n)≤2d−1⋅n\Delta_u(d,n)\le 2^{d-1}\cdot nΔu​(d,n)≤2d−1⋅n and states that it continues to hold in the base abstraction. The target is the pointwise form of D(d,n)≤2d−1 nD(d,n)\le 2^{d-1}\, nD(d,n)≤2d−1n:

∀d≥1, ∀n, ∀G∈Bd,n, ∀u,v∈V(G):dist⁡G(u,v)≤2d−1 n.\forall d\ge 1,\ \forall n,\ \forall G\in\mathcal B_{d,n},\ \forall u,v\in V(G):\qquad \operatorname{dist}_G(u,v)\le 2^{d-1}\, n .∀d≥1, ∀n, ∀G∈Bd,n​, ∀u,v∈V(G):distG​(u,v)≤2d−1n.

In Lean this is the existence, for every pair u,vu,vu,v, of a walk from uuu to vvv with at most 2d−1n2^{d-1}n2d−1n edges. The statement quantifies over every graph satisfying the paper's definition, not over polyhedron graphs or a chosen family.

Significance

For fixed ddd the bound is linear in nnn, and it holds for a class defined only by finite sets and one path condition. Together with the almost-quadratic lower bound D(n/4,n)=Ω(n2/log⁡n)D(n/4,n)=\Omega(n^2/\log n)D(n/4,n)=Ω(n2/logn) announced in the same paper, it shows what this one property of polyhedra can and cannot deliver: linear growth in nnn for fixed dimension, but not a linear bound uniform in ddd and nnn. Since the graph of every non-degenerate polyhedron lies in Bd,n\mathcal B_{d,n}Bd,n​ (a separate mission of this series), the result also recovers a Larman-type bound for polyhedra, with constant 2d−12^{d-1}2d−1.

Prove2Me already has proved polytope versions: Hirsch.larman_bound (the sharper constant n 2d−3n\,2^{d-3}n2d−3 for nonempty bounded HHH-polytopes), Hirsch.larman_high_dimension, Hirsch.larman_layer_recursion, and Hirsch.tight_row_interval, the polytope form of the interval property that condition i) abstracts. They concern a different object, the vertex-edge graph of an HHH-polytope, and are credited as prior work; none states the bound for every graph of Bd,n\mathcal B_{d,n}Bd,n​. The platform definition Hirsch_clf (connected layer families) is a related abstraction in layered form and is not the graph class used here.

Difficulty

The class contains graphs with no geometric realization, so a proof cannot use convexity, faces as polyhedra, or linear-algebraic facts about vertices; only labels, intersections and walks are available. Condition i) guarantees a path preserving the common labels of two endpoints, but says nothing directly about the length of that path. Mere connectivity on ddd-subsets of [n][n][n] gives only the trivial bound ∣V∣−1≤(nd)−1|V|-1\le\binom{n}{d}-1∣V∣−1≤(dn​)−1, which is far from linear in nnn. A linear bound in nnn has to be extracted from condition i) alone, uniformly over all graphs of the class. The abstract describes its proof only as much simpler than Larman's for polyhedra and states no intermediate lemma, so the mission carries no milestones.

Formalization scope

  • [n][n][n] is represented by Fin n, labels 0,…,n−10,\ldots,n-10,…,n−1; only the number of labels matters.
  • A vertex set is V : Finset (Finset (Fin n)) and the graph is G : SimpleGraph V; vertices are therefore distinct label sets, as in the paper. InB d n V G requires V.Nonempty, every vertex of cardinality ddd, and CondI V G.
  • CondI asks for a walk all of whose vertices contain u∩vu\cap vu∩v. This is the source's path condition: a walk contains a path on a subset of its vertices, and the endpoints contain u∩vu\cap vu∩v automatically.
  • "Continue to hold in the base abstraction" is read as the bound for every graph of Bd,n\mathcal B_{d,n}Bd,n​ and every pair of its vertices. No numeric D(d,n)D(d,n)D(d,n) is defined, so no supremum with a default value enters; a bounded walk witnesses the distance bound.
  • d≥1d\ge 1d≥1 is an explicit hypothesis: ddd is the dimension, and it makes d−1d-1d−1 ordinary subtraction in N\mathbb NN. At d=0d=0d=0 the class consists only of the one-vertex graph on ∅\emptyset∅, so nothing is lost. For n<dn<dn<d the class is empty.
  • The constant is 2d−12^{d-1}2d−1, the paper's constant for the abstraction; the polytope constant 2d−32^{d-3}2d−3 of Hirsch.larman_bound is not claimed here.

Trivializing formalizations are ruled out: labels are the vertices themselves (not an arbitrary labelling that could repeat), and condition i) is required for all pairs, not for one chosen path or only adjacent vertices. A proof only for polyhedron graphs, for graphs with extra structure, or with a larger constant leaves the goal open. Useful contributions include general lemmas about CondI and walk-length bookkeeping for SimpleGraph on a finite vertex type; these are reusable by the other missions of this series.

Selected references

  • F. Eisenbrand, N. Hähnle, A. Razborov, T. Rothvoß, Diameter of Polyhedra: Limits of Abstraction, Dagstuhl Seminar Proceedings 10211 (Flexible Network Design), 2010. https://drops.dagstuhl.de/opus/volltexte/2010/2724
  • F. Eisenbrand, N. Hähnle, A. Razborov, T. Rothvoß, Diameter of Polyhedra: Limits of Abstraction, Mathematics of Operations Research 35(4), 786–794, 2010. https://doi.org/10.1287/moor.1100.0470
  • D. G. Larman, Paths on polytopes, Proceedings of the London Mathematical Society (3) 20, 161–178, 1970. https://doi.org/10.1112/plms/s3-20.1.161
2 thms2 active usersReviewed
🏆Completed
Operations ResearchProbabilityTheoretical Computer Science·Captain: mikedeng1

Online Stochastic Matching: Beating 1-1/e 2: The Suggested Matching Algorithm Achieves 1 − 1/e with High Probability, and This Is Tight Even in ExpectationResearch Paper

Motivation

Online bipartite matching models decisions that must be made as requests arrive: an advertiser can be assigned to a compatible request once, and an assignment cannot be revised when later requests reveal a better alternative. Display advertising is the motivating example in Feldman, Mehta, Mirrokni, and Muthukrishnan's 2009 preprint: an ad server knows which advertisers accept each kind of impression and has estimates of how frequently those kinds will arrive. The question is how much this advance distributional information helps when decisions still have to be immediate.

This mission studies the paper's first offline-guided algorithm. It computes a maximum matching for the expected traffic and follows that matching during the actual random run. The paper shows that this natural policy reaches the familiar 1−1/e1-1/e1−1/e performance level and that its own analysis cannot be improved for this policy, even if performance is averaged over runs. The result sets the baseline for the paper's later two-suggested-matchings algorithm, which improves on that level under its stated assumptions. Feldman et al., §§1 and 4.1

Setting

An instance consists of a finite advertiser set AAA, a finite impression-type set III, and allowed edges E⊆A×IE\subseteq A\times IE⊆A×I. For each type iii, the nonnegative integer eie_iei​ is its expected number of arrivals. There are n=∑i∈Iei>0n=\sum_{i\in I}e_i>0n=∑i∈I​ei​>0 arrivals. Each arrival independently has type iii with probability ei/ne_i/nei​/n. A type with ei=0e_i=0ei​=0 remains in the graph but has zero arrival probability. On an arrival of type iii, an online algorithm may assign it to an adjacent advertiser that has not been assigned before, or may leave it unassigned. The hindsight optimum, OPT(ω)\mathrm{OPT}(\omega)OPT(ω), is the largest matching in the realization graph with one vertex for every arrival position, including separate vertices for repeated types. Feldman et al., §2, pp. 3–4

The suggested matching algorithm first selects any maximum integral flow in the expected-instance network. Each advertiser has capacity one, and each type iii has capacity eie_iei​. Equivalently, the selected edges form a maximum degree-capped bipartite matching M⊆EM\subseteq EM⊆E. Let A∗A^*A∗ be the advertisers covered by MMM. When type iii arrives, the algorithm chooses each advertiser joined to iii by a selected edge with probability 1/ei1/e_i1/ei​; any remaining probability chooses no advertiser. It assigns the chosen advertiser if available and otherwise makes no assignment. Its number of assignments is ALG(ω)\mathrm{ALG}(\omega)ALG(ω). This rule includes the algorithm's random choice in addition to the random arrival types. Feldman et al., §4.1, p. 5

The canonical residual cut of MMM places each advertiser and impression type on the source side if it is reachable from the source by residual edges. Write ATA_TAT​ for advertisers on the sink side and ISI_SIS​ for types on the source side. These sets describe the comparison between the expected-instance matching and the optimum of a realized run. The mission also uses occupancy: after nnn independent uniform throws into nnn bins, count how many bins in a fixed subset receive at least one ball. Feldman et al., §§2.1 and 4.1

Formalization targets

Theorem 4's general-instance guarantee is expressed in the additive form established by its analysis. For every ε>0\varepsilon>0ε>0, there are δ>0\delta>0δ>0 and NNN, uniform across all finite instances, all maximum integral flows, and all valid ways of realizing the algorithm's random choice, such that n≥Nn\ge Nn≥N implies

Pr⁡ ⁣[ALG(ω)≥(1−e−1)OPT(ω)−εn]≥1−e−δn.\Pr\!\left[\mathrm{ALG}(\omega)\ge(1-e^{-1})\mathrm{OPT}(\omega)-\varepsilon n\right] \ge 1-e^{-\delta n}.Pr[ALG(ω)≥(1−e−1)OPT(ω)−εn]≥1−e−δn.

The complete bipartite family gives the tightness target. When A=I={0,…,n−1}A=I=\{0,\ldots,n-1\}A=I={0,…,n−1} and every ei=1e_i=1ei​=1, every maximum expected-instance matching is perfect. For every run, OPT=n\mathrm{OPT}=nOPT=n, and

E[ALG]=n(1−(1−1n)n),lim⁡n→∞E[ALG]n=1−e−1.\mathbb E[\mathrm{ALG}] =n\left(1-\left(1-\frac1n\right)^n\right), \qquad \lim_{n\to\infty}\frac{\mathbb E[\mathrm{ALG}]}{n}=1-e^{-1}.E[ALG]=n(1−(1−n1​)n),n→∞lim​nE[ALG]​=1−e−1.

The milestone list follows the paper's Fact 1 and the named passages “Bounding ALG,” “Bounding OPT,” and “Tightness of the Analysis” in §4.1. The cut identity ∣A∗∣=∣AT∣+∑i∈ISei|A^*|=|A_T|+\sum_{i\in I_S}e_i∣A∗∣=∣AT​∣+∑i∈IS​​ei​ is a separate deterministic milestone. Feldman et al., Theorem 4 and §4.1, pp. 5–6

Significance

The result establishes exactly what this one-matching policy achieves under integer-frequency independent arrivals. It gives a guarantee for the actual number of assignments relative to the best assignment made with hindsight, and a family on which the limiting expected ratio equals the guarantee. That tight family explains why the paper introduces a second suggested matching rather than seeking a stronger bound for the same policy. Feldman et al., §4

The mathematical theorem is proved in the paper; the mission asks for its machine-checked formalization. The development would also provide reusable finite models of repeated-type realizations, integral degree-capped bipartite matchings, uniform occupancy, and residual reachability cuts. No machine-checked proof of these mission items is being claimed by this draft.

Difficulty

The expected-instance matching is selected before the arrivals, but the hindsight optimum can exploit the actual multiplicities of every type. Counting only the ads selected online does not compare the algorithm with that hindsight optimum. Also, a type may have several selected advertisers when ei>1e_i>1ei​>1, so replacing the algorithm's random choice by a deterministic designated ad would change its law. The cut and concentration statements have to apply uniformly to every maximum integral flow, including flows chosen by different tie-breaking rules. Feldman et al., §4.1, pp. 5–6

Formalization scope

Advertisers and types are finite Lean types. The edge relation is a finite set, and eie_iei​ and nnn are natural numbers with n=∑iei>0n=\sum_i e_i>0n=∑i​ei​>0. An integral maximum flow is represented by a maximum cardinality edge set with advertiser degree at most one and type-iii degree at most eie_iei​. The theorem quantifies over every such set. This is the unit-advertiser, integer-type-capacity flow used in §4.1, without a separate real-valued flow object. The canonical cut is defined through residual reachability. The hindsight optimum maximizes over matchings of the realized graph, with each arrival position distinct.

The run uses nnn independent uniform draws from ∑i{0,…,ei−1}\sum_i\{0,\ldots,e_i-1\}∑i​{0,…,ei​−1}. A valid labelling assigns each selected advertiser at type iii to a distinct copy. The drawn copy determines the type and, if labelled, the ad selected by the algorithm. This gives type probability ei/ne_i/nei​/n, conditional ad probability 1/ei1/e_i1/ei​ on selected edges, and the remaining “no ad” probability. Counts, probabilities, and expectations use finite sums, so there is no integrability convention. Since n>0n>0n>0, the run sample space is nonempty; no value of ALG/OPT\mathrm{ALG}/\mathrm{OPT}ALG/OPT at OPT=0\mathrm{OPT}=0OPT=0 is needed. The complete-graph family has n≥1n\ge1n≥1.

The paper writes 1−e−Ω(n)1-e^{-\Omega(n)}1−e−Ω(n) in both bounding passages. The Lean statements spell this out as ∀ε>0, ∃δ>0, ∃N, ∀\forall\varepsilon>0,\ \exists\delta>0,\ \exists N,\ \forall∀ε>0, ∃δ>0, ∃N, ∀ instances with n≥Nn\ge Nn≥N, with δ,N\delta,Nδ,N preceding the instance. Theorem 4's ratio language is represented by the additive estimate its proof yields; a vanishing ratio error requires a separate lower bound on OPT/n\mathrm{OPT}/nOPT/n. The exact finite-nnn expectation and its limit make “tight, even in expectation” precise. The printed Fact 1 exponent is εn/2\varepsilon n/2εn/2, while its Appendix A proof yields ε2n/2\varepsilon^2n/2ε2n/2; the formalized concentration statement uses the proved exponent and the milestone retains the printed wording. Neither a ratio with a zero denominator nor a labelling that changes the algorithm's choice law is accepted as a shortcut.

Contributions toward the occupancy bound, the residual cut identity, the realized matching bound, and the complete-graph expectation are welcome. The finite occupancy and matching interfaces are intended for reuse beyond this specific algorithm.

Selected references

  • Jon Feldman, Aranyak Mehta, Vahab Mirrokni, and S. Muthukrishnan, Online Stochastic Matching: Beating 1-1/e, arXiv:0905.4100v1, 2009; FOCS 2009. Preprint
10 thms2 active usersReviewed
🏆Completed
Discrete GeometryGraph Theory·Captain: mikedeng1

Diameter of Polyhedra: Limits of Abstraction 2: The Kalai–Kleitman Bound n^(1+log d) Holds in the Base AbstractionResearch Paper

Motivation

The diameter of a polyhedron is the largest number of edges needed to travel between two vertices of its graph. Bounding that number from the dimension and the number of facets is a central question in polyhedral combinatorics because edge walks are the discrete moves available to many vertex-based optimization methods. The Kalai–Kleitman bound gives a quasi-polynomial upper bound in these parameters. Eisenbrand, Hähnle, Razborov and Rothvoß ask how much of the geometry is needed to retain such a bound. They isolate one path property of polyhedron graphs and show that the upper bound persists for every graph satisfying it, including graphs that need not arise from a polyhedron (Dagstuhl abstract, pp. 1–3).

This mission concerns that extension of the Kalai–Kleitman bound. The held source is the five-page Dagstuhl extended abstract dated June 15, 2010. It states the result without its proof. The full proofs appeared in Mathematics of Operations Research 35 (2010), 786–794 (journal article). The goal is an asserted result of the paper, not an open conjecture; the formal proof remains to be supplied.

Setting

Let [n]={1,…,n}[n]=\{1,\ldots,n\}[n]={1,…,n}, and let ([n]d)\binom{[n]}{d}(d[n]​) be the family of its ddd-element subsets. A graph G=(V,E)G=(V,E)G=(V,E) in the paper's base abstraction has a nonempty vertex set V⊆([n]d)V\subseteq\binom{[n]}{d}V⊆(d[n]​). Thus each graph vertex is labelled by a distinct ddd-element subset. The edges can be chosen freely subject to condition i): for every pair u,v∈Vu,v\in Vu,v∈V, there is a path from uuu to vwhoseintermediateverticesallcontainv whose intermediate vertices all contain vwhoseintermediateverticesallcontainu\cap v.Theconditionimpliesthatthegraphisconnected.Thepaperdenotestheclassofthesegraphsby. The condition implies that the graph is connected. The paper denotes the class of these graphs by .Theconditionimpliesthatthegraphisconnected.Thepaperdenotestheclassofthesegraphsby\mathcal B_{d,n}andcallsand callsandcallsdthedimensionandthe dimension andthedimensionandn$ the number of facets of the abstraction (definition and footnote, p. 2).

The distance between two vertices is the minimum number of edges in a path between them. A graph's diameter is the largest of these distances, and D(d,n)D(d,n)D(d,n) is the largest diameter among graphs in Bd,n\mathcal B_{d,n}Bd,n​. Condition i) is motivated by non-degenerate polyhedra: if a vertex lies on exactly ddd facets, its tight facets provide a ddd-element label, and two vertices can be connected inside the smallest face containing them. Every vertex on that path retains the facets common to the endpoints (pp. 2–3). This motivation does not restrict GGG to a realizable polyhedron graph.

Formalization targets

Kalai–Kleitman bound in the base abstraction

For positive dimension, the target is the pointwise form of the paper's assertion that the upper bound continues to hold in the base abstraction:

∀d≥1, ∀n, ∀G∈Bd,n, ∀u,v∈V(G),dist⁡G(u,v)≤n1+log⁡2d.\forall d\geq 1,\ \forall n,\ \forall G\in\mathcal B_{d,n},\ \forall u,v\in V(G),\quad \operatorname{dist}_G(u,v)\leq n^{1+\log_2 d}.∀d≥1, ∀n, ∀G∈Bd,n​, ∀u,v∈V(G),distG​(u,v)≤n1+log2​d.

Equivalently, every pair has a connecting path with at most n1+log⁡2dn^{1+\log_2 d}n1+log2​d edges. Consequently D(d,n)≤n1+log⁡2dD(d,n)\leq n^{1+\log_2 d}D(d,n)≤n1+log2​d. The statement covers every graph with the paper's labels and condition i), rather than a selected realization or a graph with extra geometric properties. The source prints the inequality for Δu(d,n)\Delta_u(d,n)Δu​(d,n), its polyhedron-diameter notation, and says in the same sentence that it continues to hold for the base abstraction (p. 3).

Significance

The bound shows that a common-facet path property is enough to prevent arbitrarily large diameter once ddd and nnn are fixed. It extends an upper bound first attached to graphs of polyhedra to a class specified entirely by finite sets and graph paths. This locates a concrete portion of the polyhedral diameter theory within combinatorics. The same paper gives an almost-quadratic lower bound for the abstraction, so the path property alone does not imply a linear bound in nnn (p. 3).

A machine-checked proof would connect a precise finite graph model with the real-exponent inequality, without relying on a polyhedron representation. Prove2Me already has a proved theorem called Hirsch.kalai_kleitman_bound for a bounded HHH-polytope, with a different exponent log⁡2d+2\log_2 d+2log2​d+2, and a proved Hirsch.kalai_kleitman_recursion for that polytope model. It also has Hirsch.clf_endpoint_peeling for connected layer families. Those results concern different objects or bounds and are useful prior work, but none is a formalization of the present claim for all G∈Bd,nG\in\mathcal B_{d,n}G∈Bd,n​.

Difficulty

The class permits edges and vertex families with no geometric realization, so a proof cannot use metric or convexity properties of a particular polyhedron. Condition i) guarantees one path preserving the common labels of two endpoints; it does not directly bound that path's length. Connectivity by itself would be too weak: a connected graph on the allowed labels can have a long shortest path unless its labels and the path condition constrain it. The formal task is to obtain the claimed uniform length bound from that condition across all finite families of ddd-sets. The Dagstuhl abstract calls the argument an adaptation of Kalai and Kleitman's proof but states no intermediate lemmas, so this proposal does not assign any in-source milestones (p. 3).

Formalization scope

Lean represents [n][n][n] by Fin n, with labels 0,…,n−10,\ldots,n-10,…,n−1 instead of 1,…,n1,\ldots,n1,…,n; only the number of labels matters. A vertex set is a Finset of Finset (Fin n), and a SimpleGraph has exactly those sets as vertices. This representation enforces distinct labels. InB d n V G requires VVV to be nonempty, every label to have cardinality ddd, and CondI for every ordered pair. The walk in CondI is equivalent to the source's path: repeated vertices can be erased, and the endpoints themselves contain u∩vu\cap vu∩v. Nonemptiness and condition i) together give the source's connectedness. A path through all graph vertices without the common-intersection condition would be a weaker and inadmissible model.

The theorem takes d≥1d\geq1d≥1, so the logarithm has its ordinary positive-dimension interpretation. The source does not print a logarithm base; the statement reads it as base two, following the cited Kalai–Kleitman convention. This interpretation is explicit in the theorem's prose and is an open source-reading question. The exponent is real: the path length is cast from a natural number to R\mathbb RR for comparison with (n:R)1+log⁡2d(n:\mathbb R)^{1+\log_2 d}(n:R)1+log2​d. No floor is applied to the exponent. The theorem asserts a bounded walk for each pair; that avoids defining D(d,n)D(d,n)D(d,n) through a supremum with an artificial value for an empty class. At n=0n=0n=0 with d≥1d\geq1d≥1, InB has no instance, as there are no ddd-element labels; at d=1d=1d=1, the bound is nnn, consistent with a connected graph on at most nnn labels.

The development needs reusable finite-set reasoning about common intersections, graph walks and their lengths, and real-logarithm bounds. Contributions establishing those facts for the paper's CondI model are within scope. A proof only for polyhedron graphs, for graphs with extra edge restrictions, or for a single chosen graph would leave the goal open.

Selected references

  • F. Eisenbrand, N. Hähnle, A. Razborov and T. Rothvoß, Diameter of Polyhedra: Limits of Abstraction, Dagstuhl Seminar Proceedings 10211, 2010. Extended abstract.
  • F. Eisenbrand, N. Hähnle, A. Razborov and T. Rothvoß, Diameter of Polyhedra: Limits of Abstraction, Mathematics of Operations Research 35, 786–794, 2010. DOI.
2 thms2 active usersReviewed
🏆Completed
Linear OptimizationOperations Research·Captain: mikedeng1

Solving Large-Scale Zero-One Linear Programming Problems: A Minimal Cover Inequality Cuts Off x̄ iff the Knapsack Problem (2.12) Has Optimal Value Less Than OneResearch Paper

Motivation

Large zero–one linear programs can contain constraints involving only a small fraction of their variables. Crowder, Johnson and Padberg study how to extract useful inequalities from one such row while solving the larger program. Their computational method identifies a violated inequality at a current linear programming solution, adds it, and resolves the relaxation. The mathematical question behind that step is whether a minimal cover inequality can be found by a separate optimization problem. The authors answer it in Section 2.3 of their 1983 paper, after developing cover and configuration inequalities in Section 2.2. Crowder, Johnson and Padberg (1983)

The target is a known result from that paper. It is a precise statement about a finite knapsack row and a point in the unit cube, independent of the paper's implementation and numerical experiments. The paper also discusses (1,k)(1,k)(1,k)-configurations and lifting, which turn the identified inequalities into valid cuts involving additional row variables. Those statements supply the mission's milestones and make the separation result useful in its original setting. Crowder, Johnson and Padberg (1983), Sections 2.2–2.4

Setting

Fix a finite index set KKK, positive rational coefficients aja_jaj​ for j∈Kj\in Kj∈K, and a rational right-hand side a0a_0a0​. A zero–one solution is a vector with each xj∈{0,1}x_j\in\{0,1\}xj​∈{0,1} satisfying the single row

∑j∈Kajxj≤a0.(2.5)\sum_{j\in K}a_jx_j\le a_0. \tag{2.5}j∈K∑​aj​xj​≤a0​.(2.5)

The model uses the support set x={j:xj=1}x=\{j:x_j=1\}x={j:xj​=1} for such a vector. A set S⊆KS\subseteq KS⊆K is a minimal cover when its total coefficient exceeds a0a_0a0​, but removing any one member makes the total at most a0a_0a0​:

∑j∈Saj>a0,∑j∈Saj−ak≤a0(k∈S).(2.6)\sum_{j\in S}a_j>a_0, \qquad \sum_{j\in S}a_j-a_k\le a_0\quad(k\in S). \tag{2.6}j∈S∑​aj​>a0​,j∈S∑​aj​−ak​≤a0​(k∈S).(2.6)

It yields the cover inequality ∑j∈Sxj≤∣S∣−1\sum_{j\in S}x_j\le |S|-1∑j∈S​xj​≤∣S∣−1 for every feasible zero–one vector. The right-hand side is interpreted as a real number, so the expression also has its usual meaning when SSS is empty.

A (1,k)(1,k)(1,k)-configuration consists of S∗⊆KS^*\subseteq KS∗⊆K, t∉S∗t\notin S^*t∈/S∗ and an integer 2≤k≤∣S∗∣2\le k\le |S^*|2≤k≤∣S∗∣. The set S∗S^*S∗ itself fits the row, while Q∪{t}Q\cup\{t\}Q∪{t} is a minimal cover for every kkk-element subset QQQ of S∗S^*S∗. For any k≤r≤∣S∗∣k\le r\le |S^*|k≤r≤∣S∗∣ and rrr-element T⊆S∗T\subseteq S^*T⊆S∗, its inequality is (r−k+1)xt+∑j∈Txj≤r(r-k+1)x_t+\sum_{j\in T}x_j\le r(r−k+1)xt​+∑j∈T​xj​≤r. Crowder, Johnson and Padberg (1983), pp. 810–811

For a point xˉ∈[0,1]K\bar x\in[0,1]^Kxˉ∈[0,1]K, the separation problem asks for a cover s⊆Ks\subseteq Ks⊆K minimizing ∑j∈s(1−xˉj)\sum_{j\in s}(1-\bar x_j)∑j∈s​(1−xˉj​), subject to the strict condition ∑j∈saj>a0\sum_{j\in s}a_j>a_0∑j∈s​aj​>a0​. This is problem (2.12). The set of attainable values is finite, but it is empty if the row has no cover. An optimal value zzz is therefore asserted only when a least attainable value exists.

Formalization targets

Cover separation

The goal is the paper's Section 2.3 equivalence:

(∃S⊆K minimal:∑j∈Sxˉj>∣S∣−1)⟺z<1,\left(\exists S\subseteq K\text{ minimal}: \sum_{j\in S}\bar x_j>|S|-1\right) \quad\Longleftrightarrow\quad z<1,​∃S⊆K minimal:j∈S∑​xˉj​>∣S∣−1​⟺z<1,

where zzz is the attained optimum of (2.12). A cover inequality cuts off xˉ\bar xxˉ precisely when its left-hand side exceeds ∣S∣−1|S|-1∣S∣−1. If there is no cover, (2.12) has no optimal value; the theorem does not assign it an artificial value. Crowder, Johnson and Padberg (1983), pp. 812–813

Valid inequalities and lifting

The milestones state the validity of (2.7) and every inequality in (2.9), the two assertions used in the separation equivalence, and the relaxed lifting claims of Section 2.4. For lifting, zkz_kzk​ is the maximum integer objective value in (2.10), zˉk\bar z_kzˉk​ is the maximum in its linear relaxation, and zk∗=⌊zˉk⌋z_k^*=\lfloor\bar z_k\rfloorzk∗​=⌊zˉk​⌋. The claims are zk≤zk∗z_k\le z_k^*zk​≤zk∗​ and validity after adding variable kkk with coefficient fk=f0−zk∗f_k=f_0-z_k^*fk​=f0​−zk∗​. They concern attained optima, as the paper's “maximum” language requires. Crowder, Johnson and Padberg (1983), pp. 811, 814

Significance

The equivalence turns a geometric question about which cover inequality excludes xˉ\bar xxˉ into a finite optimization test with a numerical threshold of one. It identifies when a minimal cover cut exists for a row, while the validity milestones certify that the inequalities can be added without removing zero–one feasible points. The lifting claims explain how an inequality first written on a subset of variables remains valid as further variables enter it. These are the mathematical guarantees used by the paper's cutting plane procedure. Crowder, Johnson and Padberg (1983), Sections 2.2–2.4

The paper proves the separation equivalence and states the surrounding validity claims. This mission records their exact statements in Lean; its theorem proofs remain open. A complete development would add machine-checked proofs for the finite cover argument, configuration inequalities, and lifting validity. The resulting definitions of feasible supports, attained optimization values and intermediate inequality validity can also be reused in other finite knapsack formalizations.

Difficulty

Checking all subsets of KKK directly grows rapidly with the row size. A separation result must relate a minimum over all covers to an inequality indexed by a minimal cover, while accounting for objective coefficients that may be zero when xˉj=1\bar x_j=1xˉj​=1. The same boundary matters for lifting: the zero–one maximum is compared with a continuous relaxation, and rounding is justified by the integer coefficients of the current inequality. If the lifting problem is infeasible, neither maximum exists. Crowder, Johnson and Padberg (1983), pp. 812–814

Formalization scope

Lean uses an arbitrary finite index type for KKK, replacing the paper's indices 1,…,n1,\ldots,n1,…,n. A zero–one vector is a finite support set. Row coefficients and a0a_0a0​ are rational, as on p. 810; the point xˉ\bar xxˉ and separation values are real. The target asks only that xˉ\bar xxˉ lie in the unit cube. Although the paper obtains it as an optimum of the full LP relaxation (2.11), that additional property is not used in the row-level equivalence. No sign condition is imposed on a0a_0a0​.

The strict knapsack condition in (2.12) is retained. “Chops off” means strict violation of (2.7). “Optimal value” means membership and leastness in the attainable value set; for lifting it means membership and greatestness. Thus rows with no cover or infeasible lifting subproblem do not acquire a default zero optimum. The size expressions ∣S∣−1|S|-1∣S∣−1 and r−k+1r-k+1r−k+1 are evaluated in the reals, so natural-number truncation cannot change them. Integer lifting coefficients and the floor of the relaxed optimum express the paper's “truncating to its integer part”; the feasible relaxed problem includes the zero vector, so this agrees with truncation toward zero there.

Validity during an intermediate lifting step ranges over feasible zero–one vectors supported in the current set SSS. After adding kkk, it ranges over support in S∪{k}S\cup\{k\}S∪{k}. The final inequality covers the whole row when that set is KKK. The development must preserve the full quantifier over every kkk-subset in (2.8) and every rrr and TTT in (2.9); restricting these would weaken the paper's claim. Contributions proving the named milestones, adding concrete examples, or developing general finite optimization lemmas are in scope. The facet assertions attributed to Padberg are outside the target. A related lifted-cover facet statement already appears on the platform as NemhauserWolsey.partitioned_cover_lifting_defines_facet; its conventions and claim differ from this mission's validity results.

Selected references

  • H. Crowder, E. L. Johnson and M. Padberg, Solving Large-Scale Zero-One Linear Programming Problems, Operations Research 31(5), 803–834, 1983. DOI: 10.1287/opre.31.5.803
8 thms2 active usersReviewed
🏆Completed
Operations ResearchTheoretical Computer Science·Captain: mikedeng1

On the Approximability of Single-Machine Scheduling with Precedence Constraints 2: Every Convex Bipartite Order Has a Realizer of Size 3Research Paper

Motivation

The single-machine scheduling problem 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_j C_j1∣prec∣∑wj​Cj​ asks for an order in which to process nnn jobs, each with a processing time and a weight, on one machine, respecting a partial order of precedence constraints, so that the weighted sum of completion times is as small as possible. It is NP-hard, and the best known approximation ratio for general precedence constraints is 222. Ambühl, Mastrolilli, Mutsanas and Svensson (Math. Oper. Res. 36(4), 2011) relate the problem to the dimension theory of partial orders: their Theorem 3.2 states that when the precedence constraints are given together with a realizer of size kkk, the problem has a (2−2/k)(2 - 2/k)(2−2/k)-approximation algorithm. Small realizers of structured precedence orders therefore translate directly into better approximation ratios.

This mission formalizes one such structural result, Lemma 4.1 of the paper: every convex bipartite order has a realizer of size 333. Convex bipartite orders form a class of precedence constraints that lies strictly between strong bipartite orders and general bipartite orders (Möhring, 1989). Combined with Theorem 3.2, the lemma gives a 4/34/34/3-approximation for 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_j C_j1∣prec∣∑wj​Cj​ on this class. According to the authors, the bound dim⁡≤3\dim \le 3dim≤3 for convex bipartite orders had not been known before.

Setting

A poset is a set NNN with a reflexive, antisymmetric, transitive relation PPP. Two elements x,yx, yx,y are incomparable, written x∥yx \parallel yx∥y, when neither (x,y)∈P(x,y) \in P(x,y)∈P nor (y,x)∈P(y,x) \in P(y,x)∈P; the set of incomparable pairs is inc(P)\mathrm{inc}(P)inc(P). A linear extension of PPP is a linear order LLL on NNN with P⊆LP \subseteq LP⊆L. A linear order LLL reverses the pair (x,y)(x,y)(x,y) when y<xy < xy<x in LLL. A realizer of size ttt is a family L1,…,LtL_1, \dots, L_tL1​,…,Lt​ of linear extensions of PPP such that every incomparable pair (x,y)(x,y)(x,y) is reversed by at least one LiL_iLi​. The dimension dim⁡(P)\dim(P)dim(P) is the least ttt for which a realizer of size ttt exists.

A convex bipartite order has jobs N=J−∪J+N = J^- \cup J^+N=J−∪J+ split into minus jobs J−={j1,…,ja}J^- = \{j_1, \dots, j_a\}J−={j1​,…,ja​} and plus jobs J+={ja+1,…,jn}J^+ = \{j_{a+1}, \dots, j_n\}J+={ja+1​,…,jn​}. Each plus job jkj_kjk​ carries two indices 1≤l(k)≤r(k)≤a1 \le l(k) \le r(k) \le a1≤l(k)≤r(k)≤a, and a minus job jij_iji​ precedes jkj_kjk​ exactly when l(k)≤i≤r(k)l(k) \le i \le r(k)l(k)≤i≤r(k). There are no other precedences: the predecessors of each plus job form a nonempty interval of consecutive minus jobs.

For the Appendix construction, the incomparable pairs are sorted into three sets E1,E2,E3E_1, E_2, E_3E1​,E2​,E3​ according to the kinds of the two jobs, the order of their indices and, for a pair (plus job jij_iji​, minus job jjj_jjj​), whether jjj_jjj​ precedes some plus job of larger index than iii. The sets Eˉm=Em∪P\bar E_m = E_m \cup PEˉm​=Em​∪P describe what the mmm-th linear order of the realizer must contain.

Formalization targets

Goal: Lemma 4.1

For every convex bipartite order (N,P):∃ L1,L2,L3 linear extensions of P with ∀(x,y)∈inc(P) ∃m: y<Lmx.\text{For every convex bipartite order } (N, P):\quad \exists\, L_1, L_2, L_3 \text{ linear extensions of } P \text{ with } \forall (x,y) \in \mathrm{inc}(P)\ \exists m:\ y <_{L_m} x .For every convex bipartite order (N,P):∃L1​,L2​,L3​ linear extensions of P with ∀(x,y)∈inc(P) ∃m: y<Lm​​x.

Equivalently, dim⁡(P)≤3\dim(P) \le 3dim(P)≤3. The goal holds for every numbering of the jobs and every choice of interval ends l,rl, rl,r.

Milestones

  • Lemma A.1. E1,E2,E3E_1, E_2, E_3E1​,E2​,E3​ partition inc(P)\mathrm{inc}(P)inc(P), and for each mmm, (x,y)∈Em(x,y) \in E_m(x,y)∈Em​ implies (y,x)∉Em(y,x) \notin E_m(y,x)∈/Em​.
  • Lemma A.2. Under the Appendix's numbering of the plus jobs (i<ji < ji<j implies l(i)≤l(j)l(i) \le l(j)l(i)≤l(j)), each Eˉm\bar E_mEˉm​ is contained in a linear order on NNN.

Significance

The result. Lemma 4.1 places convex bipartite orders among the classes of precedence constraints for which 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_j C_j1∣prec∣∑wj​Cj​ has an approximation ratio strictly below 222, namely 4/34/34/3 through Theorem 3.2 of the paper. The bound is tight: a bipartite order has dimension 222 exactly when it is a strong bipartite order (Möhring), so the bound 333 cannot be lowered for the class of convex bipartite orders. Since dim⁡(P)=3\dim(P) = 3dim(P)=3 also gives χ(GP)=3\chi(G_P) = 3χ(GP​)=3 for the graph of incomparable pairs, a 333-realizer yields an optimal colouring of that graph.

Formalizing it. The result is proved in the paper, with a short case analysis in the Appendix. No machine-checked proof is known to exist. Mathlib has no dimension theory of posets: realizers, reversals of incomparable pairs and the construction of linear extensions containing a given acyclic relation have to be set up here. The definitions of realizer and linear extension are stated for an arbitrary relation and can be reused for other dimension bounds (interval orders, semiorders, the other missions of this series).

Difficulty

The three relations Eˉm\bar E_mEˉm​ are not partial orders. Eˉ2\bar E_2Eˉ2​, for instance, is not transitive: with two minus jobs and one plus job j3j_3j3​ with l(3)=r(3)=2l(3) = r(3) = 2l(3)=r(3)=2, the pairs (j1,j2)∈E2(j_1, j_2) \in E_2(j1​,j2​)∈E2​ and (j2,j3)∈P(j_2, j_3) \in P(j2​,j3​)∈P are present but (j1,j3)(j_1, j_3)(j1​,j3​) lies in E1E_1E1​. The step that requires work is to show that each Eˉm\bar E_mEˉm​ has no cycle, so that it extends to a linear order. For Eˉ2\bar E_2Eˉ2​ this depends on the convexity of the predecessor intervals and on the numbering of the plus jobs by their left ends: it is not a consequence of bipartiteness alone. The goal itself carries no numbering assumption, so any proof must also account for renumbering the plus jobs.

Formalization scope

All declarations live in the namespace SingleMachinePrec.ConvexBipartite.

  • Relations. A relation on NNN is a predicate N → N → Prop. Linear orders use Mathlib's unbundled class IsLinearOrder. A linear extension of PPP is a linear order LLL with P⊆LP \subseteq LP⊆L; LLL reverses (x,y)(x,y)(x,y) when L y xL\,y\,xLyx and y≠xy \ne xy=x. A realizer of size ttt is a function Fin t → N → N → Prop, so repeated members are allowed, as in the paper's multisets.
  • Convex bipartite orders. The jobs are the type Fin a ⊕ Fin b: Sum.inl i is the minus job ji+1j_{i+1}ji+1​, Sum.inr k the plus job ja+k+1j_{a+k+1}ja+k+1​, with indices starting at 000. The interval ends are maps l r : Fin b → Fin a with l k ≤ r k, and PPP is equality together with the pairs (minus iii, plus kkk) with l(k)≤i≤r(k)l(k) \le i \le r(k)l(k)≤i≤r(k). A file-level instance records that PPP is a partial order. Any convex bipartite order in the paper's sense is one of these after naming its jobs. The cases a=0a = 0a=0 and b=0b = 0b=0 are allowed.
  • E1,E2,E3E_1, E_2, E_3E1​,E2​,E3​. Each set is defined as "incomparable and satisfies the paper's case condition". In the plus–minus clauses, "k>ik > ik>i" compares plus indices, because kkk exceeds the index of a plus job.
  • Lemma A.2. "Eˉm\bar E_mEˉm​ is an extension of PPP" is formalized as "there is a linear order containing Eˉm\bar E_mEˉm​", which is what the paper's proof establishes (absence of cycles) and what its final paragraph uses. Reading it as "Eˉm\bar E_mEˉm​ is a partial order" would make the lemma false (see Difficulty). The Appendix's numbering assumption is the hypothesis Monotone l of this milestone only.
  • Algorithmic wording not formalized. The paper states Lemma 4.1 as "a realizer of size 3 can be computed in polynomial time". What is formalized is the existence of the realizer; the running time is not.

A trivializing formalization would let the members of the realizer be arbitrary relations, which makes reversing every pair free. Here every member must be a linear order containing PPP.

Contributions welcome: proofs of the two milestones and of the goal; a general lemma that a relation whose transitive closure is antisymmetric extends to a linear order on a finite type; and the renumbering argument that removes the numbering assumption.

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
  • G. R. Brightwell, E. R. Scheinerman, Fractional dimension of partial orders, Order 9(2):139–158, 1992.
  • R. H. Möhring, Computationally tractable classes of ordered sets, in I. Rival (ed.), Algorithms and Order, NATO ASI Series 255, Kluwer, 1989, pp. 105–193.
  • W. T. Trotter, Combinatorics and Partially Ordered Sets: Dimension Theory, Johns Hopkins University Press, 1992.
7 thms2 active usersReviewed
🏆Completed
Graph TheoryOperations Research·Captain: mikedeng1

On the Graph Structure of Convex Polyhedra in n-Space II: Whitney's Theorem, a Graph Is n-Tuply Connected iff Any Two Points Are Joined by n Disjoint PathsResearch Paper

Motivation

Vertex connectivity measures how robust a network is against the failure of nodes. It can be measured in two ways that look different. One way counts the fewest nodes whose removal disconnects the network. The other counts the routes between two nodes that share no intermediate node. Whitney's theorem (1932) says that the two measures agree for every pair of nodes. It is the vertex form of Menger's theorem, and it underlies reliability analysis of communication and transportation networks, the design of fault-tolerant routing, and much of structural graph theory.

M. L. Balinski's 1961 paper On the graph structure of convex polyhedra in n-space proves that the graph of a bounded full-dimensional polyhedron in nnn-space is nnn-tuply connected (the subject of Mission I of this series). It then invokes Whitney's theorem to conclude that any two vertices of such a polyhedron are joined by nnn disjoint paths. Balinski gives a short new proof of Whitney's theorem through the max-flow min-cut theorem of Ford and Fulkerson and of Dantzig and Fulkerson. That makes the theorem a consequence of linear programming duality. This mission formalizes that part of the paper: the network vocabulary, the max-flow min-cut theorem with capacities on both points and lines, the integrality of maximum flows, and Whitney's theorem itself.

Timeline.

  • 1927: Menger states the disjoint-paths theorem for separating sets.
  • 1932: Whitney proves the characterization of nnn-connected graphs by nnn disjoint paths between every pair of points.
  • 1956: Ford and Fulkerson and Dantzig and Fulkerson prove the max-flow min-cut theorem.
  • 1961: Balinski derives Whitney's theorem from it with a unit-capacity network.

Setting

A graph GGG consists of a finite set VVV of points and a set of lines, each line being a pair of distinct points. A path from psp_sps​ to pkp_kpk​ is a sequence of lines (p1,p2),(p2,p3),…,(pm,pm+1)(p_1,p_2),(p_2,p_3),\dots,(p_m,p_{m+1})(p1​,p2​),(p2​,p3​),…,(pm​,pm+1​) with p1=psp_1 = p_sp1​=ps​, pm+1=pkp_{m+1} = p_kpm+1​=pk​ and m≥1m \ge 1m≥1. Paths are disjoint if they have no point in common except possibly their first and last points.

GGG is nnn-tuply connected if it has at least n+1n+1n+1 points and, for every set XXX of fewer than nnn points, the graph G−XG - XG−X remaining after deleting XXX is connected. GGG has nnn disjoint paths from psp_sps​ to pkp_kpk​ if there are nnn pairwise distinct paths from psp_sps​ to pkp_kpk​, none of which repeats a point, and no two of which share a point other than psp_sps​ and pkp_kpk​.

A network is a connected graph with a capacity c(x)≥0c(x) \ge 0c(x)≥0 on every point and c(e)≥0c(e) \ge 0c(e)≥0 on every line, and with a distinguished source psp_sps​ and sink pkp_kpk​. A flow assigns a number f(C)≥0f(C) \ge 0f(C)≥0 to every path CCC from psp_sps​ to pkp_kpk​, such that for every point xxx and every line eee

∑C∋xf(C)≤c(x),∑C∋ef(C)≤c(e).\sum_{C \ni x} f(C) \le c(x), \qquad \sum_{C \ni e} f(C) \le c(e).C∋x∑​f(C)≤c(x),C∋e∑​f(C)≤c(e).

Its value is val⁡(f)=∑Cf(C)\operatorname{val}(f) = \sum_C f(C)val(f)=∑C​f(C). A disconnecting set is a pair (X,F)(X,F)(X,F) of points and lines that meets every walk from psp_sps​ to pkp_kpk​. Its value is ∑x∈Xc(x)+∑e∈Fc(e)\sum_{x\in X} c(x) + \sum_{e \in F} c(e)∑x∈X​c(x)+∑e∈F​c(e).

The unit network of the proof has capacity 111 on every point except psp_sps​ and pkp_kpk​, and capacity n+1n+1n+1 on every line except the line pspkp_sp_kps​pk​ (if present), which has capacity 111. In Lean these are IsNTuplyConnected, HasNDisjointPaths, IsFlow, flowValue, IsDisconnecting, cutValue, unitCapV and unitCapE, all in the namespace Balinski61.Whitney.

Formalization targets

Goal: Whitney's theorem (p. 434)

For a finite graph GGG with at least two points and any n≥0n \ge 0n≥0:

G is n-tuply connected  ⟺  for all ps≠pk, G has n disjoint paths from ps to pk.G \text{ is } n\text{-tuply connected} \iff \text{for all } p_s \ne p_k,\ G \text{ has } n \text{ disjoint paths from } p_s \text{ to } p_k.G is n-tuply connected⟺for all ps​=pk​, G has n disjoint paths from ps​ to pk​.

Both directions are part of the goal.

Milestones, in the order of the proof

  1. Max-flow min-cut (p. 433). In every network there is a number MMM that is the value of some flow and of some disconnecting set, with every flow of value at most MMM and every disconnecting set of value at least MMM.
  2. Integrality (p. 434). If all capacities are integers, some maximum flow has only integer path flows.
  3. Min-cut in the unit network (p. 434). If GGG is nnn-tuply connected and ps≠pkp_s \ne p_kps​=pk​, every disconnecting set of the unit network has value at least nnn.
  4. Paths from unit flows (p. 434). An integral flow of value at least nnn in the unit network yields nnn disjoint paths from psp_sps​ to pkp_kpk​.
  5. Sufficiency (p. 434). If every pair of distinct points is joined by nnn disjoint paths, GGG is nnn-tuply connected.

Significance

Whitney's theorem turns a statement about all small deletion sets into the existence of explicit, verifiable path systems, and back again. In applications it certifies connectivity by exhibiting paths, and it certifies that connectivity is no larger by exhibiting a separating set. It is the base of the theory of kkk-connected graphs: ear decompositions, the fan lemma, and the structure of minimally kkk-connected graphs all use it. Inside this paper it supplies the COROLLARY that any two vertices of a bounded full-dimensional polyhedron in nnn-space are joined by nnn disjoint edge paths.

None of these results is formalized for vertex connectivity at this Mathlib revision. Mathlib has edge connectivity and connected components but no vertex Menger theorem. The Prove2Me library has max-flow min-cut statements for arc capacities only and integrality results for basic solutions of network LPs, but no flow model with capacities on points. A completed development gives a reusable vertex-capacitated max-flow min-cut theorem for undirected graphs and the first machine-checked Whitney theorem in this library. The results are classical and proved; the remaining work is the formalization.

Difficulty

The sufficiency direction is elementary. The necessity direction needs a global object (a flow, or a family of paths) to exist from purely local hypotheses about deletions. The obvious induction on nnn, which deletes a point and applies the hypothesis to a smaller graph, does not keep the path systems disjoint. In Balinski's route the weight falls on max-flow min-cut and integrality for path flows with capacities on points, neither of which exists in the library. A second difficulty sits in a case the paper's proof skips: a disconnecting set of the unit network may use the line pspkp_sp_kps​pk​ (capacity 111) together with up to n−2n-2n−2 points, and the deletion hypothesis of nnn-tuple connectedness speaks only about points.

Formalization scope

Graphs are Mathlib SimpleGraphs on a Fintype with decidable equality and adjacency. Paths are walks with IsPath. A flow is a real function on the finite type G.Path ps pk of simple paths; "through a point" and "through a line" mean membership in the walk's support and edge list. Capacities are functions V → ℝ and Sym2 V → ℝ. Nonnegativity, connectivity of GGG and ps≠pkp_s \ne p_kps​=pk​ are hypotheses of the network theorems.

The following readings of loose phrases are explicit in the statements:

  • "dropping out n−1n-1n−1 or fewer points" is ∣X∣<n|X| < n∣X∣<n;
  • "nnn disjoint paths" means nnn pairwise distinct simple paths. The printed path syntax permits repeated vertices; in a graph with lines psap_s aps​a, psbp_s bps​b, and pspkp_s p_kps​pk​, the distinct walks ps,a,ps,pkp_s,a,p_s,p_kps​,a,ps​,pk​ and ps,b,ps,pkp_s,b,p_s,p_kps​,b,ps​,pk​ share only their endpoints even though the graph is not 222-tuply connected. The theorem therefore uses its conventional simple-path reading;
  • path flows live on simple paths (merging and shortcutting changes no maximum value);
  • the paper leaves the capacities of psp_sps​ and pkp_kpk​ in the unit network unassigned, and here they are n+1n+1n+1;
  • "the condition is sufficient is obvious" is the full statement that nnn-tuple connectedness follows;
  • the hypothesis ∣V∣≥2|V| \ge 2∣V∣≥2 is added to the goal and to sufficiency, because the paper's "any pair of points" presupposes it and the equivalence fails for a one-point graph.

The max-flow min-cut milestone states that the maximum and the minimum are attained. A statement that only bounds some flow by every cut is satisfied by the zero flow. A connectivity notion without the n+1n+1n+1 point count would make every complete graph nnn-connected for all nnn. Both trivializations are excluded.

Contributions welcome: a vertex-capacitated augmenting-path or LP-duality proof of max-flow min-cut for path flows, integrality by an augmenting-path argument, the unit-network lemmas, and direct combinatorial proofs of Whitney's theorem that bypass flows.

Selected references

  • M. L. Balinski, On the graph structure of convex polyhedra in n-space, Pacific J. Math. 11 (1961), 431–434. https://doi.org/10.2140/pjm.1961.11.431
  • H. Whitney, Congruent graphs and the connectivity of graphs, Amer. J. Math. 54 (1932), 150–168. https://doi.org/10.2307/2371086
  • L. R. Ford, Jr. and D. R. Fulkerson, Maximal flow through a network, Canadian J. Math. 8 (1956), 399–404. https://doi.org/10.4153/CJM-1956-045-5
  • G. B. Dantzig and D. R. Fulkerson, On the max-flow min-cut theorem of networks, in Linear Inequalities and Related Systems, Ann. of Math. Stud. 38, Princeton Univ. Press, 1956, 215–221. https://doi.org/10.1515/9781400881987
  • K. Menger, Zur allgemeinen Kurventheorie, Fund. Math. 10 (1927), 96–115. https://doi.org/10.4064/fm-10-1-96-115
9 thms2 active usersReviewed
🏆Completed
Operations ResearchTheoretical Computer Science·Captain: mikedeng1

On the Approximability of Single-Machine Scheduling with Precedence Constraints 6: The Optimal Value of S_G Lies Between n² − an²(ln 1/a + 2) and n² − an²Research Paper

Motivation

The problem 1 ∣ prec ∣ ∑wjCj1\,|\,\mathrm{prec}\,|\,\sum w_jC_j1∣prec∣∑wj​Cj​ asks for a single-machine sequence of jobs, respecting precedence constraints, that minimizes the weighted sum of completion times. It has been known to be strongly NP-hard since Lawler (1978) and Lenstra and Rinnooy Kan (1978), several different 2-approximation algorithms are known, and closing the approximability gap is listed by Schuurman and Woeginger (1999) as one of ten outstanding open problems in scheduling theory. Ambühl, Mastrolilli, Mutsanas and Svensson (Math. Oper. Res. 36(4), 2011) give the first inapproximability result for this problem: under a widely believed complexity assumption it has no polynomial-time approximation scheme (PTAS). The bridge to that result is a quantitative link, Lemma 9.1, between the optimal value of a special bipartite scheduling instance and the maximum edge biclique of a bipartite graph, a problem whose hardness of approximation was established by Ambühl, Mastrolilli and Svensson (FOCS 2007). This mission formalizes that link.

Setting

A schedule of a finite job set is a sequence σ\sigmaσ listing every job once; the machine processes the jobs in that order from time 000 without idle time or pre-emption. Job jjj has a processing time pjp_jpj​ and a weight wjw_jwj​; its completion time CjC_jCj​ is the sum of the processing times of the jobs up to and including jjj, and the value of σ\sigmaσ is val(σ)=∑jwjCj\mathrm{val}(\sigma)=\sum_j w_jC_jval(σ)=∑j​wj​Cj​. Precedence constraints are a relation PPP on jobs: (i,j)∈P(i,j)\in P(i,j)∈P with i≠ji\ne ji=j means job iii must be completed before job jjj starts. A schedule respecting all of them is feasible, and a feasible schedule σ∗\sigma^*σ∗ of least value is optimal.

Let G=(U,V,E)G=(U,V,E)G=(U,V,E) be an nnn-by-nnn bipartite graph: ∣U∣=∣V∣=n|U|=|V|=n∣U∣=∣V∣=n and E⊆U×VE\subseteq U\times VE⊆U×V. An edge biclique is a pair A⊆UA\subseteq UA⊆U, B⊆VB\subseteq VB⊆V with A×B⊆EA\times B\subseteq EA×B⊆E, of value ∣A∣⋅∣B∣|A|\cdot|B|∣A∣⋅∣B∣; the maximum edge biclique problem (Definition 9.1) asks for one of largest value. The bipartite scheduling instance SGS_GSG​ has jobs U∪VU\cup VU∪V and precedence constraints

P=(U×V)∖E,P=(U\times V)\setminus E,P=(U×V)∖E,

so u∈Uu\in Uu∈U must precede v∈Vv\in Vv∈V exactly when (u,v)(u,v)(u,v) is not an edge. Jobs of UUU have p=1p=1p=1, w=0w=0w=0; jobs of VVV have p=0p=0p=0, w=1w=1w=1. Thus val(σ)=∑v∈VCv\mathrm{val}(\sigma)=\sum_{v\in V}C_vval(σ)=∑v∈V​Cv​, where CvC_vCv​ is the number of UUU-jobs scheduled before vvv. For i≥1i\ge1i≥1, σ(i)\sigma(i)σ(i) denotes the number of VVV-jobs scheduled before iii jobs of UUU have been scheduled.

In the Lean development these are weightedCompletion, IsOptimalSchedule, IsEdgeBiclique, maxBicliqueValue, precSG, procSG, weightSG, valSG, IsOptimalSG and vBefore in the namespace SingleMachinePrec.Biclique.

Formalization targets

Goal: Lemma 9.1 (p. 666)

If a maximum edge biclique of GGG has value an2an^2an2 with a∈(0,1]a\in(0,1]a∈(0,1], then SGS_GSG​ has an optimal schedule and every optimal schedule σ∗\sigma^*σ∗ satisfies

n2−an2(ln⁡1a+2)≤val(σ∗)≤n2−an2.n^2-an^2\Bigl(\ln\frac1a+2\Bigr)\le\mathrm{val}(\sigma^*)\le n^2-an^2 .n2−an2(lna1​+2)≤val(σ∗)≤n2−an2.

Milestones (proof of Lemma 9.1, §9, p. 666)

  1. For every edge biclique (A,B)(A,B)(A,B), a schedule in the block order U∖A→B→A→V∖BU\setminus A\to B\to A\to V\setminus BU∖A→B→A→V∖B exists, and every such schedule is feasible with
val(σ)=n2−∣A∣⋅∣B∣.\mathrm{val}(\sigma)=n^2-|A|\cdot|B| .val(σ)=n2−∣A∣⋅∣B∣.
  1. For every schedule, σ(n+1)=n\sigma(n+1)=nσ(n+1)=n and
val(σ)=∑i=1n(σ(i+1)−σ(i))i=n2−∑i=1nσ(i).\mathrm{val}(\sigma)=\sum_{i=1}^n\bigl(\sigma(i+1)-\sigma(i)\bigr)i=n^2-\sum_{i=1}^n\sigma(i).val(σ)=i=1∑n​(σ(i+1)−σ(i))i=n2−i=1∑n​σ(i).
  1. For every feasible schedule and i=1,…,ni=1,\dots,ni=1,…,n,
σ(i)(n−i+1)≤an2,σ(i)≤n.\sigma(i)(n-i+1)\le an^2,\qquad \sigma(i)\le n .σ(i)(n−i+1)≤an2,σ(i)≤n.

Significance

Lemma 9.1 shows that the optimal value of SGS_GSG​ determines the maximum edge biclique of GGG up to a factor of order ln⁡(1/a)\ln(1/a)ln(1/a) in the "area above the work line" n2−val(σ∗)n^2-\mathrm{val}(\sigma^*)n2−val(σ∗). Combined with the hardness of approximating maximum edge biclique (Theorem 9.1, cited from Ambühl, Mastrolilli and Svensson 2007) it yields Theorem 9.2: 1 ∣ prec ∣ ∑wjCj1\,|\,\mathrm{prec}\,|\,\sum w_jC_j1∣prec∣∑wj​Cj​ has no PTAS unless SAT can be decided by a probabilistic algorithm in time 2Nϵ2^{N^\epsilon}2Nϵ for every ϵ>0\epsilon>0ϵ>0. It also makes precise the two-dimensional Gantt chart picture of Eastman, Even and Isaacs (1964) and of Goemans and Williamson (2000), in which every point on the work line of a schedule defines an edge biclique.

The lemma is proved in the paper; it is not formalized anywhere to our knowledge. A formal proof certifies the combinatorial core of the no-PTAS result independently of the complexity-theoretic layer, and its definitions (the bipartite instance SGS_GSG​, edge bicliques, the profile σ(i)\sigma(i)σ(i)) are reusable for the gap inequality behind Theorem 9.2.

Difficulty

The upper bound is a direct computation on one explicit schedule. The lower bound is a statement about every feasible schedule, of which there are exponentially many, and it must hold with the explicit constant 222 and the factor ln⁡(1/a)\ln(1/a)ln(1/a) for every a∈(0,1]a\in(0,1]a∈(0,1]. The printed argument splits the sum at i=(1−a)ni=(1-a)ni=(1−a)n and uses ⌊an⌋\lfloor an\rfloor⌊an⌋, treating ananan as an integer; for general aaa (for example n=3n=3n=3, value 222, an=2/3an=2/3an=2/3) the split point is not an integer, so the printed estimate does not apply verbatim and the constant 222 has to be re-checked for non-integral ananan. On the formal side, the value identity requires relating completion times in a list to counting UUU-jobs before each VVV-job, with ties among zero-length jobs.

Formalization scope

Jobs are the disjoint union U ⊕ V of two finite types with Fintype.card U = Fintype.card V = n; EEE is a relation U → V → Prop. A schedule is a duplicate-free list containing every job; feasibility is the published LawlerPrec.MinMax.IsFeasible and completion times are the published MooreLateJobs.Shared.completionTime (time 000 start, no idle time). Processing times and weights are reals, here in {0,1}\{0,1\}{0,1}. The maximum edge biclique value is the maximum of ∣A∣⋅∣B∣|A|\cdot|B|∣A∣⋅∣B∣ over all edge bicliques, the empty ones included, so the hypothesis a>0a>0a>0 means E≠∅E\ne\emptysetE=∅. The logarithm is natural (Real.log).

Conventions and readings committed to:

  • The goal is stated for every optimal schedule, and the existence of an optimal schedule is a separate conclusion, so the bounds cannot hold vacuously. Proving the bounds for one particular schedule, or for an optimal value defined as an infimum that could be a junk default, would not be this lemma.
  • No integrality hypothesis on ananan is added.
  • Milestones 2 and 3 are stated for every schedule (respectively every feasible schedule), not only for σ∗\sigma^*σ∗; milestone 1 states the value of the block-order schedule as an equality, where the paper writes "≤⋯=\le\cdots=≤⋯=".
  • The paper's P=(U×V)∖EP=(U\times V)\setminus EP=(U×V)∖E is irreflexive; feasibility only constrains distinct jobs, so it agrees with the reflexive partial order of §1.

Not formalized: Theorem 9.1 (cited hardness of maximum edge biclique) and Theorem 9.2 (no PTAS under a complexity assumption); no polynomial-time or complexity-theoretic statement appears in the mission. Contributions welcome: proofs of the three milestones and of the goal; Mathlib's bounds on harmonic numbers (Mathlib/NumberTheory/Harmonic/Bounds.lean) are the relevant library.

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
  • C. Ambühl, M. Mastrolilli, O. Svensson, Inapproximability results for sparsest cut, optimal linear arrangement, and precedence constraint scheduling, Proc. 48th IEEE FOCS, 329–337, 2007 (reference [4] of the paper).
  • W. L. Eastman, S. Even, I. M. Isaacs, Bounds for the optimal scheduling of n jobs on m processors, Management Science 11(2):268–279, 1964 (reference [11]).
  • M. X. Goemans, D. P. Williamson, Two-dimensional Gantt charts and a scheduling algorithm of Lawler, SIAM J. Discrete Math. 13(3):281–294, 2000 (reference [15]).
  • P. Schuurman, G. J. Woeginger, Polynomial time approximation algorithms for machine scheduling: ten open problems, J. Scheduling 2(5):203–213, 1999 (reference [36]).
8 thms2 active usersReviewed
🏆Completed
Graph TheoryOperations ResearchTheoretical Computer Science·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 thms2 active usersReviewed
PreviousPage 10 of 15Next

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