Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Theoretical Computer Science

464 missions · 237 completed

The mathematical foundations of computation: which problems can be solved, by what algorithms, and at what cost in time, space, or communication. Distinguished by its emphasis on rigor and unconditional lower bounds, it spans computational complexity, algorithm design, automata and computability, cryptography, and the analysis of Boolean functions.

Missions

Open227Completed237All464
🏆Completed
Graph TheoryLinear Optimization·Captain: mikedeng1

Multicommodity Max-Flow Min-Cut Theorems and Their Use in Designing Approximation Algorithms I: Every Uniform Multicommodity Flow Problem Has Ω(𝒮/log n) ≤ f ≤ 𝒮Research Paper

Motivation

A multicommodity flow problem asks to route several commodities simultaneously through a capacitated network, each from its own source to its own sink. For a single commodity, the max-flow min-cut theorem of Ford and Fulkerson says that the largest routable flow equals the capacity of the smallest separating cut. For many commodities the natural analogue fails: the max-flow can be strictly smaller than the min-cut, already for four commodities (Okamura and Seymour, 1981). Leighton and Rao asked how large the gap can be, and showed that for uniform demands it is at most logarithmic in the number of nodes. The result turned the min-cut, an NP-hard quantity, into something approximable by a linear program, and became the basis of the first polylogarithmic approximation algorithms for sparsest cut, balanced cut, minimum feedback arc set, crossing number and VLSI layout (Leighton–Rao, J. ACM 1999).

Timeline:

  • 1956–1963. Ford and Fulkerson prove max-flow = min-cut for one commodity; Hu (1963) proves equality for two commodities.
  • 1981–1990. Okamura and Seymour give a 4-commodity example with max-flow 3/4 and min-cut 1. Schrijver (1990), combining Lomonosov and Papernov, characterizes demand graphs for which equality always holds. Shahrokhi and Matula (1990) prove equality for uniform demands under a condition on the dual.
  • 1988. Leighton and Rao (FOCS 1988) prove the Θ(log⁡n)\Theta(\log n)Θ(logn) gap for uniform demands; the journal version (1999) is the source of this mission.
  • 1995. Linial, London and Rabinovich, and independently Aumann and Rabani, extend the bound to arbitrary demands with O(log⁡k)O(\log k)O(logk) for kkk commodities, via metric embeddings (Combinatorica 1995).

Setting

A network is a finite set VVV of nnn nodes with a capacity C(u,v)=C(v,u)≥0C(u,v)=C(v,u)\ge0C(u,v)=C(v,u)≥0 for every pair; the pair is an edge when C(u,v)>0C(u,v)>0C(u,v)>0. For U⊆VU\subseteq VU⊆V with complement Uˉ\bar UUˉ, the cut capacity C(U,Uˉ)C(U,\bar U)C(U,Uˉ) is the total capacity of the edges joining UUU to Uˉ\bar UUˉ. As in the paper, the network is assumed connected: C(U,Uˉ)>0C(U,\bar U)>0C(U,Uˉ)>0 for every nonempty proper UUU.

In the uniform multicommodity flow problem (UMFP) every unordered pair of nodes is a commodity with demand one. A routing of value λ\lambdaλ sends λ\lambdaλ units between every pair simultaneously, so that the total flow of all commodities over each edge, in both directions together, is at most its capacity. The max-flow fff is the largest such λ\lambdaλ (the concurrent max-flow).

The min-cut (sparsest cut) of the UMFP is

S=min⁡∅≠U⊊VC(U,Uˉ)∣U∣ ∣Uˉ∣,\mathcal S=\min_{\emptyset\ne U\subsetneq V}\frac{C(U,\bar U)}{|U|\,|\bar U|},S=∅=U⊊Vmin​∣U∣∣Uˉ∣C(U,Uˉ)​,

the minimum over cuts of the cut capacity divided by the demand separated by the cut. The quantity C(U,Uˉ)/(∣U∣∣Uˉ∣)C(U,\bar U)/(|U||\bar U|)C(U,Uˉ)/(∣U∣∣Uˉ∣) is the ratio cost of the cut.

The milestones use the dual objects. A distance function ddd assigns a nonnegative length to each edge; d(u,v)d(u,v)d(u,v) is then the shortest-path distance, d(T,u)=min⁡t∈Td(t,u)d(T,u)=\min_{t\in T}d(t,u)d(T,u)=mint∈T​d(t,u), and the total weight is W=∑eC(e)d(e)W=\sum_{e}C(e)d(e)W=∑e​C(e)d(e). The distance constraint (5) is ∑{u,v}d(u,v)≥1\sum_{\{u,v\}}d(u,v)\ge1∑{u,v}​d(u,v)≥1, summed over unordered pairs. A set of nodes has radius at most Δ\DeltaΔ if some centre reaches all its nodes by walks inside the set of length at most Δ\DeltaΔ. Throughout, log⁡n=log⁡2n\log n=\log_2 nlogn=log2​n.

Formalization targets

Goal: Theorem 2 (p. 795)

Ω ⁣(Slog⁡n)≤f≤S\Omega\!\left(\frac{\mathcal S}{\log n}\right)\le f\le\mathcal SΩ(lognS​)≤f≤S

for every uniform multicommodity flow problem: there is one absolute constant c>0c>0c>0 such that c S/log⁡n≤f≤Sc\,\mathcal S/\log n\le f\le\mathcal ScS/logn≤f≤S for every connected network with n≥2n\ge2n≥2 nodes. The goal fixes no constant, so it is the form any improved argument would still prove.

Milestones, in the order the argument uses them

  1. Weak duality (§1.2, pp. 789–790): every concurrent flow of value λ\lambdaλ satisfies λ ∣U∣∣Uˉ∣≤C(U,Uˉ)\lambda\,|U||\bar U|\le C(U,\bar U)λ∣U∣∣Uˉ∣≤C(U,Uˉ) for every cut.
  2. LP duality (§2.2, p. 796): every distance function satisfying (5) has W≥fW\ge fW≥f, and some such function has W=fW=fW=f.
  3. Lemma 3 (p. 797): for every Δ>0\Delta>0Δ>0, VVV splits into components of radius at most Δ\DeltaΔ with at most 4Wlog⁡n/Δ4W\log n/\Delta4Wlogn/Δ capacity between components.
  4. Corollary 4 (pp. 798–799): either a set of at least 2n/32n/32n/3 nodes has radius at most 1/(2n2)1/(2n^2)1/(2n2), or some cut has ratio cost at most 36Wlog⁡n36W\log n36Wlogn.
  5. Lemma 5 (p. 799): if ∣T∣≥2n/3|T|\ge 2n/3∣T∣≥2n/3 and ∑u∉Td(T,u)≥1/(2n)\sum_{u\notin T}d(T,u)\ge 1/(2n)∑u∈/T​d(T,u)≥1/(2n), some cut has ratio cost at most 6W6W6W.
  6. Lemma 6 (p. 800): a distance function satisfying (5) yields a cut of ratio cost at most 36Wlog⁡n36W\log n36Wlogn.

Companion

S≤36 flog⁡n\mathcal S\le 36\,f\log nS≤36flogn and f≤Sf\le\mathcal Sf≤S: Theorem 2 with the constant its proof produces.

Significance

The theorem says that the LP relaxation of uniform sparsest cut has integrality gap O(log⁡n)O(\log n)O(logn), and the expander example of §2.1 (Theorem 1) shows this is tight. Combined with the cut-finding argument, it gives an O(log⁡n)O(\log n)O(logn)-approximation for sparsest cut and, through the paper's Section 3, approximation algorithms for balanced separators, minimum cut linear arrangement, crossing number, VLSI layout area, feedback arc set and the emulation of networks. Region-growing decompositions of the kind in Lemma 3 reappear in low-diameter decompositions, multicut (Garg–Vazirani–Yannakakis), and metric embedding.

The result is proved; no machine-checked proof of it is known. A formalization would provide: concurrent multicommodity flow and its LP dual on finite graphs, the uniform sparsest cut, a verified region-growing decomposition, and the max-flow min-cut gap theorem itself. The LP-duality milestone, which the paper cites rather than proves, is a reusable piece of flow theory in its own right.

Difficulty

The upper bound f≤Sf\le\mathcal Sf≤S is a counting argument. The lower bound is the content. The natural idea, rounding an optimal dual distance function by a threshold and cutting the long edges, does not work: the paper notes (p. 797) that no thresholding rule is known that always yields a good cut, because edges with large distance need not form a sparse cut. Any argument must convert the fractional dual into a single cut while losing only a log⁡n\log nlogn factor against the total weight WWW, which requires a global property of the distance function rather than the lengths of individual edges. The LP-duality step needs a correct dual for the concurrent flow problem with shortest-path distances, which Mathlib does not provide off the shelf.

Formalization scope

  • Nodes and capacities. V is a Fintype with decidable equality; capacities are a symmetric nonnegative real function vanishing on the diagonal, and the graph is the simple graph of pairs with positive capacity. Goal-level statements assume 2 ≤ Fintype.card V and connectivity in the form C(U,Uˉ)>0C(U,\bar U)>0C(U,Uˉ)>0 for every nonempty proper UUU (the paper's standing assumption, p. 789).
  • Demands and flows. Each unordered pair is encoded as two ordered commodities with demand 12\tfrac1221​ (the paper's footnote 2). Flows are arc flows with conservation per commodity, and both directions of an edge share its capacity. The max-flow is the real supremum of feasible values; with n≥2n\ge2n≥2 that set is nonempty and bounded.
  • Cuts. The min-cut ranges over nonempty proper subsets only. Including U=∅U=\emptysetU=∅ or U=VU=VU=V would give 0/00/00/0, which Lean evaluates to 000 and which would make the goal trivially true; this formalization excludes it.
  • Distances. d(u,v)d(u,v)d(u,v) is an infimum over walks in the graph, meaningful under connectivity; the radius of a component is measured by walks inside it.
  • Constants. The goal keeps the page's Ω\OmegaΩ with a single universal constant chosen before the network. The milestones carry the constants computed in the paper's proofs (363636 in Corollary 4 and Lemma 6, 666 in Lemma 5), and each statement says so.
  • Algorithms. "We can find" is formalized as existence. Running times and the remark that the proofs yield polynomial-time algorithms (p. 801) are out of scope.

Contributions welcome: proofs of any milestone; a general LP-duality theorem for concurrent multicommodity flow, reusable for the product and directed variants; and a general region-growing lemma.

Selected references

  • T. Leighton and S. Rao, Multicommodity max-flow min-cut theorems and their use in designing approximation algorithms, J. ACM 46(6):787–832, 1999. https://doi.org/10.1145/331524.331526
  • F. Shahrokhi and D. W. Matula, The maximum concurrent flow problem, J. ACM 37(2):318–334, 1990. https://doi.org/10.1145/77600.77620
  • H. Okamura and P. D. Seymour, Multicommodity flows in planar graphs, J. Combin. Theory Ser. B 31:75–81, 1981. https://doi.org/10.1016/S0095-8956(81)80012-3
  • N. Linial, E. London and Y. Rabinovich, The geometry of graphs and some of its algorithmic applications, Combinatorica 15(2):215–245, 1995. https://doi.org/10.1007/BF01200757
  • V. Chvátal, Linear Programming, W. H. Freeman, 1983.
10 thms3 active usersReviewed
🏆Completed
CombinatoricsLinear algebra·Captain: mikedeng1

Quick Approximation to Matrices and Applications I: Every Real Matrix Is a Sum of at Most 1/ε² Cut Matrices plus an Error W with |W(S,T)| ≤ ε√(|S||T|)‖A‖_FResearch Paper

Motivation

Many dense combinatorial optimization problems (Max-Cut, graph bisection, quadratic assignment on dense instances) have objective functions that are sums of matrix entries over rectangles S×TS\times TS×T. If a matrix can be replaced by a simple surrogate that agrees with it on every rectangle up to a small error, such problems can be solved approximately on the surrogate. Frieze and Kannan's paper Quick approximation to matrices and applications (Combinatorica 19 (1999) 175–220) supplies this surrogate: a sum of a bounded number of cut matrices. Their decomposition is a "weak" form of Szemerédi's regularity lemma whose number of pieces grows only like 1/ε21/\varepsilon^21/ε2, against a tower-type bound for Szemerédi's lemma. It underlies constant-time sampling algorithms for dense instances, the theory of graph limits (the cut distance of Lovász and Szegedy), and the cut-norm approximation algorithm of Alon and Naor.

Timeline:

  • 1978: Szemerédi's regularity lemma, with a number of parts of tower type in 1/ε1/\varepsilon1/ε.
  • 1996: Frieze and Kannan announce the weak regularity lemma (FOCS 1996).
  • 1999: the journal paper proves the matrix version, its algorithmic forms, and applications to Max-Cut, Max-SNP problems and pseudo-regular partitions.
  • 2004–2006: Alon and Naor approximate the cut norm via Grothendieck's inequality; Lovász and Szegedy build graph limits on the cut distance.

Setting

Let RRR and CCC be finite sets, ∣R∣=m|R|=m∣R∣=m, ∣C∣=n|C|=n∣C∣=n, and let AAA be a real R×CR\times CR×C matrix. For S⊆RS\subseteq RS⊆R, T⊆CT\subseteq CT⊆C write

A(S,T)=∑(i,j)∈S×TA(i,j).A(S,T)=\sum_{(i,j)\in S\times T}A(i,j).A(S,T)=(i,j)∈S×T∑​A(i,j).

The Frobenius norm is ∥A∥F=(∑i,jA(i,j)2)1/2\|A\|_F=\big(\sum_{i,j}A(i,j)^2\big)^{1/2}∥A∥F​=(∑i,j​A(i,j)2)1/2, and the cut norm is ∥A∥C=max⁡S⊆R, T⊆C∣A(S,T)∣\|A\|_C=\max_{S\subseteq R,\,T\subseteq C}|A(S,T)|∥A∥C​=maxS⊆R,T⊆C​∣A(S,T)∣.

For S⊆RS\subseteq RS⊆R, T⊆CT\subseteq CT⊆C and d∈Rd\in\mathbb Rd∈R, the cut matrix CUT(S,T,d)\mathrm{CUT}(S,T,d)CUT(S,T,d) has entry ddd on S×TS\times TS×T and 000 elsewhere; it has rank at most one. A cut decomposition of width sss writes

A=D(1)+⋯+D(s)+W,D(t)=CUT(Rt,Ct,dt),A=D^{(1)}+\dots+D^{(s)}+W,\qquad D^{(t)}=\mathrm{CUT}(R_t,C_t,d_t),A=D(1)+⋯+D(s)+W,D(t)=CUT(Rt​,Ct​,dt​),

with error matrix WWW. In the Lean development these are blockSum, frobNorm, cutNorm, cutMatrix and cutSum in the namespace FriezeKannan.CutDecomp.

Formalization targets

Goal: Theorem 7 (p. 188)

For every real R×CR\times CR×C matrix AAA and every ε>0\varepsilon>0ε>0 there are s≤1/ε2s\le1/\varepsilon^2s≤1/ε2 cut matrices D(t)=CUT(Rt,Ct,dt)D^{(t)}=\mathrm{CUT}(R_t,C_t,d_t)D(t)=CUT(Rt​,Ct​,dt​) such that W(s)=A−(D(1)+⋯+D(s))W^{(s)}=A-(D^{(1)}+\dots+D^{(s)})W(s)=A−(D(1)+⋯+D(s)) satisfies

∣W(s)(S,T)∣≤ε∣S∣ ∣T∣  ∥A∥F∀S⊆R, T⊆C.(17)|W^{(s)}(S,T)|\le\varepsilon\sqrt{|S|\,|T|}\;\|A\|_F\qquad\forall S\subseteq R,\ T\subseteq C. \tag{17}∣W(s)(S,T)∣≤ε∣S∣∣T∣​∥A∥F​∀S⊆R, T⊆C.(17)

Milestones (the steps of the proof, p. 188)

  1. With d=W(S,T)/(∣S∣∣T∣)d=W(S,T)/(|S||T|)d=W(S,T)/(∣S∣∣T∣) for nonempty S,TS,TS,T:
∥W−CUT(S,T,d)∥F2−∥W∥F2=−∣S∣ ∣T∣ d2=−W(S,T)2∣S∣ ∣T∣.\|W-\mathrm{CUT}(S,T,d)\|_F^2-\|W\|_F^2=-|S|\,|T|\,d^2=-\frac{W(S,T)^2}{|S|\,|T|}.∥W−CUT(S,T,d)∥F2​−∥W∥F2​=−∣S∣∣T∣d2=−∣S∣∣T∣W(S,T)2​.
  1. If ∣W(S,T)∣>ε∣S∣∣T∣ ∥A∥F|W(S,T)|>\varepsilon\sqrt{|S||T|}\,\|A\|_F∣W(S,T)∣>ε∣S∣∣T∣​∥A∥F​, this drop is at most −ε2∥A∥F2-\varepsilon^2\|A\|_F^2−ε2∥A∥F2​.
  2. If ∥W(t)∥F2≤(1−ε2t)∥A∥F2\|W^{(t)}\|_F^2\le(1-\varepsilon^2t)\|A\|_F^2∥W(t)∥F2​≤(1−ε2t)∥A∥F2​, then either W(t)W^{(t)}W(t) satisfies (17), or one more cut matrix gives ∥W(t+1)∥F2≤(1−ε2(t+1))∥A∥F2\|W^{(t+1)}\|_F^2\le(1-\varepsilon^2(t+1))\|A\|_F^2∥W(t+1)∥F2​≤(1−ε2(t+1))∥A∥F2​.

Companion statements

  • The Remark after Theorem 7 (the case m=1m=1m=1): a real nnn-vector aaa is a weighted sum of at most ε−2\varepsilon^{-2}ε−2 0-1 vectors plus an error www with ∣w(S)∣≤εn ∣a∣|w(S)|\le\varepsilon\sqrt n\,|a|∣w(S)∣≤εn​∣a∣ for all SSS.
  • Inequality (1), p. 178: ∥M∥C≤sup⁡x≠0∥Mx∥1/∥x∥∞≤4∥M∥C\|M\|_C\le\sup_{x\ne0}\|Mx\|_1/\|x\|_\infty\le4\|M\|_C∥M∥C​≤supx=0​∥Mx∥1​/∥x∥∞​≤4∥M∥C​.
  • §2.3, p. 179: the symmetrised error Wˉ=(W+WT)/2\bar W=(W+W^T)/2Wˉ=(W+WT)/2 satisfies ∥Wˉ∥C≤∥W∥C\|\bar W\|_C\le\|W\|_C∥Wˉ∥C​≤∥W∥C​.

Significance

Theorem 7 says that the width needed to approximate a matrix on all rectangles, with error measured against ∣S∣∣T∣∥A∥F\sqrt{|S||T|}\|A\|_F∣S∣∣T∣​∥A∥F​, depends on ε\varepsilonε only, not on mmm or nnn. Since ∣S∣∣T∣≤mn|S||T|\le mn∣S∣∣T∣≤mn, it gives a cut decomposition with ∥W∥C≤εmn∥A∥F\|W\|_C\le\varepsilon\sqrt{mn}\|A\|_F∥W∥C​≤εmn​∥A∥F​, the error of the paper's algorithmic Theorem 2. It is the existence statement behind the paper's constant-time approximation schemes for dense Max-Cut and Max-SNP problems and behind its pseudo-regular partitions of graphs (§5.1), and it is the matrix form of the weak regularity lemma used throughout the theory of dense graph limits.

The result is proved in the paper, by an elementary argument. Its formalization adds a machine-checked weak regularity lemma for real matrices with explicit constants, stated in a form (finite index types, block sums, cut matrices) that downstream developments in graph limits, property testing and approximation algorithms can import. To the best of the curators' knowledge it has no machine-checked proof; Mathlib contains Szemerédi's regularity lemma for graphs but not this matrix decomposition or the cut norm of matrices.

Difficulty

The obvious decomposition, one cut matrix per nonzero entry, achieves W=0W=0W=0 but has width mnmnmn; the content is that width 1/ε21/\varepsilon^21/ε2 suffices uniformly in the dimensions. Bounding the error on a single rectangle is easy; bounding it on all 2m+n2^{m+n}2m+n rectangles at once, with an error that scales as ∣S∣∣T∣\sqrt{|S||T|}∣S∣∣T∣​ rather than as mn\sqrt{mn}mn​, is the point. The bound is also against the Frobenius norm of the original matrix AAA, not of the current error, so the argument must track a single global quantity across all steps. In Lean, the main friction is bookkeeping: real casts of cardinalities inside square roots, the division by ∣S∣∣T∣|S||T|∣S∣∣T∣, and the identification of the squared Frobenius norm with a double sum.

Formalization scope

  • Index sets are arbitrary finite types R, C with decidable equality; subsets are Finsets, and ∣S∣|S|∣S∣ is the cardinality cast to R\mathbb RR. Empty RRR or CCC is allowed.
  • ∥A∥F\|A\|_F∥A∥F​ is the real square root of the sum of squares (frobNorm); Mathlib's scoped matrix norms are not used. The cut norm is a genuine maximum over the finite nonempty set of pairs of subsets.
  • A family of sss cut matrices is three functions on Fin s; s=0s=0s=0 is allowed (for ε≥1\varepsilon\ge1ε≥1 or A=0A=0A=0 no cut matrix is needed). No nonemptiness, disjointness or sign condition is imposed on RtR_tRt​, CtC_tCt​, dtd_tdt​.
  • Added hypotheses, each disclosed in the statement's Formalization Note: ε>0\varepsilon>0ε>0 (the paper's implicit range of the error parameter; ε<1\varepsilon<1ε<1 is not assumed); nonempty S,TS,TS,T in milestone 1 (they make ddd well defined); the strict inequality >>> in milestone 2 (the negation of (17); the page prints ≥\ge≥); a nonempty column set for inequality (1) (the supremum over Rn∖{0}\mathbb R^n\setminus\{0\}Rn∖{0} needs n≥1n\ge1n≥1). Inequality (1) is stated without division, as the existence of a nonzero xxx with ∥M∥C∥x∥∞≤∥Mx∥1\|M\|_C\|x\|_\infty\le\|Mx\|_1∥M∥C​∥x∥∞​≤∥Mx∥1​ together with ∥Mx∥1≤4∥M∥C∥x∥∞\|Mx\|_1\le4\|M\|_C\|x\|_\infty∥Mx∥1​≤4∥M∥C​∥x∥∞​ for all xxx.
  • The width bound s≤1/ε2s\le1/\varepsilon^2s≤1/ε2 is the theorem; without it the statement is trivial (one cut matrix per entry gives W=0W=0W=0). A formalization that drops it, bounds the error by ∥W∥F\|W\|_F∥W∥F​ instead of ∥A∥F\|A\|_F∥A∥F​, or uses mn\sqrt{mn}mn​ in place of ∣S∣∣T∣\sqrt{|S||T|}∣S∣∣T∣​ proves a different theorem.
  • Reusable infrastructure: the identity ∥M∥F2=∑i,jM(i,j)2\|M\|_F^2=\sum_{i,j}M(i,j)^2∥M∥F2​=∑i,j​M(i,j)2, block sums of cut matrices and their sums, and the cut norm with its basic properties (transpose invariance, comparison with the ∞→1\infty\to1∞→1 norm). Contributions of these as standalone lemmas, and of the milestones in the listed order, are welcome.

Selected references

  • A. Frieze and R. Kannan, Quick approximation to matrices and applications, Combinatorica 19(2) (1999) 175–220. https://doi.org/10.1007/s004930050052
  • A. Frieze and R. Kannan, The regularity lemma and approximation schemes for dense problems, Proc. 37th FOCS (1996) 12–20. https://doi.org/10.1109/SFCS.1996.548459
  • E. Szemerédi, Regular partitions of graphs, Problèmes combinatoires et théorie des graphes, Colloq. Internat. CNRS 260 (1978) 399–401.
  • N. Alon and A. Naor, Approximating the cut-norm via Grothendieck's inequality, SIAM J. Comput. 35(4) (2006) 787–803. https://doi.org/10.1137/S0097539704441629
  • L. Lovász and B. Szegedy, Limits of dense graph sequences, J. Combin. Theory Ser. B 96(6) (2006) 933–957. https://doi.org/10.1016/j.jctb.2006.05.002
5 thms3 active usersReviewed
🏆Completed
Algorithmic Game TheoryOperations ResearchOptimization·Captain: mikedeng1

Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments IV: Uniform Pricing Earns OPT/(1 + ln(v_m/v_1)) in Unit-Demand Min-PricingResearch Paper

Pricing with unit-demand consumers

A seller with nnn distinct items faces mmm consumers, each of whom wants at most one item. Choosing item prices to maximise revenue in this setting is the unit-demand envy-free pricing problem, introduced by Rusmevichientong and studied in algorithmic game theory and revenue management since the mid-2000s. Optimal prices are hard to compute in general, so the guarantees of simple pricing rules matter. The simplest such rule is uniform pricing: put one price qqq on every item and choose qqq to maximise revenue.

Aggarwal, Feder, Motwani and Zhu (ICALP 2004) and Guruswami, Hartline, Karlin, Kempe, Kenyon and McSherry (SODA 2005) analysed uniform pricing; in the min-buying model it recovers a 1/(1+ln⁡m)1/(1+\ln m)1/(1+lnm) fraction of the optimum, a bound Berbeglia and Joret attribute to Aggarwal et al. Berbeglia and Joret (arXiv:1606.01371, Algorithmica 2020) observed that the min-buying problem is a special case of assortment optimisation under a regular discrete choice model, and that uniform pricing then coincides with the revenue-ordered assortments heuristic. Their general guarantee for that heuristic yields a second bound for uniform pricing, in terms of the spread of the valuations rather than the number of consumers. This mission formalizes that bound (Corollary 4.7) together with the reduction it rests on (Theorem 4.6).

Setting

Unit-demand min-pricing (UDPmin⁡\mathrm{UDP}_{\min}UDPmin​). There are items [n][n][n] and consumers [m][m][m], m≥1m\ge 1m≥1. Consumer iii is interested in a set Bi⊆[n]B_i\subseteq[n]Bi​⊆[n] and has a valuation vi>0v_i>0vi​>0. Given a price assignment p:[n]→R>0p:[n]\to\mathbb R_{>0}p:[n]→R>0​, consumer iii buys a cheapest item of BiB_iBi​ if its price is at most viv_ivi​, and nothing otherwise. The seller's revenue is

revUDP(p)=∑i: Bi≠∅, min⁡x∈Bip(x)≤vi min⁡x∈Bip(x),\mathrm{rev}_{\mathrm{UDP}}(p)=\sum_{i:\ B_i\ne\emptyset,\ \min_{x\in B_i}p(x)\le v_i}\ \min_{x\in B_i}p(x),revUDP​(p)=i: Bi​=∅, minx∈Bi​​p(x)≤vi​∑​ x∈Bi​min​p(x),

the optimum is OPTUDP=sup⁡p>0revUDP(p)\mathrm{OPT}_{\mathrm{UDP}}=\sup_{p>0}\mathrm{rev}_{\mathrm{UDP}}(p)OPTUDP​=supp>0​revUDP​(p), and uniform pricing earns UP=sup⁡q>0revUDP(q,…,q)\mathrm{UP}=\sup_{q>0}\mathrm{rev}_{\mathrm{UDP}}(q,\dots,q)UP=supq>0​revUDP​(q,…,q). Write vmax⁡=vmv_{\max}=v_mvmax​=vm​ and vmin⁡=v1v_{\min}=v_1vmin​=v1​ for the extreme valuations.

Regular choice models. A finite product set C\mathcal CC carries choice probabilities P(x,S)\mathcal P(x,S)P(x,S) for x∈Cx\in\mathcal Cx∈C, S⊆CS\subseteq\mathcal CS⊆C, with no-purchase probability 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 (i) all probabilities are nonnegative, (ii) P(x,S)=0\mathcal P(x,S)=0P(x,S)=0 for x∉Sx\notin Sx∈/S, (iii) ∑x∈SP(x,S)≤1\sum_{x\in S}\mathcal P(x,S)\le1∑x∈S​P(x,S)≤1, and (iv) 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}. With revenues r:C→R>0r:\mathcal C\to\mathbb R_{>0}r:C→R>0​, an assortment 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⁡Srev(S)\mathrm{OPT}=\max_S\mathrm{rev}(S)OPT=maxS​rev(S). If r1<⋯<rkr_1<\dots<r_kr1​<⋯<rk​ are the distinct revenues and Si={x:r(x)≥ri}S_i=\{x: r(x)\ge r_i\}Si​={x:r(x)≥ri​}, revenue-ordered assortments earn RO=max⁡irev(Si)\mathrm{RO}=\max_i\mathrm{rev}(S_i)RO=maxi​rev(Si​).

The reduction. From a UDPmin⁡\mathrm{UDP}_{\min}UDPmin​ instance the paper builds C=[n]×{v1,…,vm}\mathcal C=[n]\times\{v_1,\dots,v_m\}C=[n]×{v1​,…,vm​}, r((x,v))=m vr((x,v))=m\,vr((x,v))=mv, and P=1m∑iPi\mathcal P=\frac1m\sum_i\mathcal P_iP=m1​∑i​Pi​, where Pi\mathcal P_iPi​ spreads probability one over the set Qi(S)Q_i(S)Qi​(S) of pairs (x,v)∈S(x,v)\in S(x,v)∈S with x∈Bix\in B_ix∈Bi​, v≤viv\le v_iv≤vi​ and vvv minimal among the pairs of SSS whose item lies in BiB_iBi​.

Formalization targets

Goal: Corollary 4.7

OPTUDP ≤ (1+ln⁡vmax⁡vmin⁡)⋅UP.\mathrm{OPT}_{\mathrm{UDP}}\ \le\ \Big(1+\ln\frac{v_{\max}}{v_{\min}}\Big)\cdot\mathrm{UP}.OPTUDP​ ≤ (1+lnvmin​vmax​​)⋅UP.

The bound is the paper's "approximates the optimum revenue to within a factor of 1/(1+ln⁡ρ)1/(1+\ln\rho)1/(1+lnρ), ρ:=vm/v1\rho:=v_m/v_1ρ:=vm​/v1​", stated in product form.

Milestones

  1. Lemma 4.5. Some optimal price assignment takes only valuations as prices; in particular OPTUDP\mathrm{OPT}_{\mathrm{UDP}}OPTUDP​ is attained.
  2. Regularity of the reduction's choice model (proof of Theorem 4.6, pp. 17–18).
  3. Assortments to prices: rev(S)=revUDP(pS)\mathrm{rev}(S)=\mathrm{rev}_{\mathrm{UDP}}(p_S)rev(S)=revUDP​(pS​) with pS(x)=min⁡{v:(x,v)∈S}p_S(x)=\min\{v:(x,v)\in S\}pS​(x)=min{v:(x,v)∈S} (p. 19).
  4. Prices to assortments: for valuation-valued ppp, rev(Sp)=revUDP(p)\mathrm{rev}(S_p)=\mathrm{rev}_{\mathrm{UDP}}(p)rev(Sp​)=revUDP​(p) with Sp={(x,v):v≥p(x)}S_p=\{(x,v): v\ge p(x)\}Sp​={(x,v):v≥p(x)} (p. 19).
  5. Theorem 4.6 for the explicit instance: regularity, positive revenues, equal optima, and the two-way correspondence between threshold sets and uniform valuation prices.
  6. Theorem 3.2 for an arbitrary regular model:
OPT≤(∑i=1kri−ri−1ri)RO,∑i=1kri−ri−1ri≤1+ln⁡rkr1,r0:=0.\mathrm{OPT}\le\Big(\sum_{i=1}^k\frac{r_i-r_{i-1}}{r_i}\Big)\mathrm{RO},\qquad \sum_{i=1}^k\frac{r_i-r_{i-1}}{r_i}\le 1+\ln\frac{r_k}{r_1},\qquad r_0:=0.OPT≤(i=1∑k​ri​ri​−ri−1​​)RO,i=1∑k​ri​ri​−ri−1​​≤1+lnr1​rk​​,r0​:=0.

Significance

The corollary gives a uniform-pricing guarantee that is independent of the number of consumers and items: when valuations lie within a factor ρ\rhoρ of each other, a single price recovers a 1/(1+ln⁡ρ)1/(1+\ln\rho)1/(1+lnρ) fraction of the optimum, which improves on 1/(1+ln⁡m)1/(1+\ln m)1/(1+lnm) whenever ρ<m\rho<mρ<m. More broadly, Theorem 4.6 places UDPmin⁡\mathrm{UDP}_{\min}UDPmin​ inside the class of regular choice models, so every guarantee proved for revenue-ordered assortments under regularity transfers to uniform pricing.

All results here are proved in the paper. Nothing in this mission is formalized elsewhere as far as a search of the platform shows: there is no regular choice model, no UDPmin⁡\mathrm{UDP}_{\min}UDPmin​ problem and no revenue-ordered guarantee on Prove2Me. The mission produces the reduction as reusable infrastructure, a machine-checked version of the paper's Lemma 4.5 (whose appendix proof is an informal repair argument), and the first formal guarantee for uniform pricing.

Difficulty

The reduction is where the work sits. Regularity hinges on the no-purchase case of axiom (iv): enlarging an assortment can only make QiQ_iQi​ nonempty, never empty, and this must be checked for the specific sets Qi(S)Q_i(S)Qi​(S) whose definition depends on the whole of SSS. The revenue identities require tracking the minimum over BiB_iBi​ of a price function obtained from an assortment, including items absent from SSS, which receive the paper's "+∞+\infty+∞" price. Lemma 4.5 is needed to make the optimum of the uncountable pricing problem coincide with that of a finite assortment problem, and the statement is about a supremum over a continuum of price vectors whose attainment is not obvious. Theorem 3.2 has to relate the optimum to threshold sets of an arbitrary regular model, where products may share revenues and the optimal assortment need not be a threshold set.

The naive route of comparing uniform pricing to the optimum directly, consumer by consumer, does not give a bound in terms of vm/v1v_m/v_1vm​/v1​; the reduction to a regular model is what supplies it.

Formalization scope

  • Items and consumers are finite types X and M, with Nonempty M (m≥1m\ge1m≥1). Valuations satisfy vi>0v_i>0vi​>0 as a field of the instance; the paper says "non-negative", but its reduction (r=m v>0r=m\,v>0r=mv>0) and its ratio vm/v1v_m/v_1vm​/v1​ need positivity.
  • A consumer with Bi=∅B_i=\emptysetBi​=∅ pays 000. Ties among cheapest items do not affect payments, so payments are defined directly.
  • OPTUDP\mathrm{OPT}_{\mathrm{UDP}}OPTUDP​ and UP\mathrm{UP}UP are real suprema over positive price assignments and positive uniform prices; both index sets are nonempty and the revenues are bounded by ∑ivi\sum_i v_i∑i​vi​. Uniform pricing ranges over all q>0q>0q>0, not only valuations.
  • The no-purchase option is not a product: it 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), and regularity includes axiom (iv) for it. OPT is a maximum over all Finset C; RO is a maximum over the kkk threshold sets only. Distinct revenues are sorted 0-based (Fin k), with r0=0r_0=0r0​=0 handled as a separate case.
  • The paper's +∞+\infty+∞ price for pSp_SpS​ is the real 1+∑ivi1+\sum_i v_i1+∑i​vi​, which exceeds every valuation, as its Remark on p. 18 allows. The valuation set collapses equal valuations.
  • Approximation factors are in product form, OPT≤D⋅RO\mathrm{OPT}\le D\cdot\mathrm{RO}OPT≤D⋅RO; ln⁡\lnln is Real.log applied to ratios ≥1\ge1≥1.
  • Trivializing readings are excluded: Theorem 4.6 is stated for the explicit instance of its proof, not as the bare existence of a regular instance with the same optimum (a single product of revenue OPTUDP\mathrm{OPT}_{\mathrm{UDP}}OPTUDP​ would satisfy that); OPTUDP\mathrm{OPT}_{\mathrm{UDP}}OPTUDP​ ranges over all positive price assignments, not uniform ones; valuations are positive, so no logarithm takes a junk value.
  • Theorem 3.2 is restated here for this mission's copy of the regular model; it is also the goal of mission I of this series. Reusable pieces: the regular choice model and revenue-ordered value, and the UDPmin⁡\mathrm{UDP}_{\min}UDPmin​ model. Contributions of any milestone, and alternative proofs of Lemma 4.5, are welcome.

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; Algorithmica 82 (2020). https://arxiv.org/abs/1606.01371
  • G. Aggarwal, T. Feder, R. Motwani and A. Zhu, Algorithms for multi-product pricing, ICALP 2004. https://doi.org/10.1007/978-3-540-27836-8_8
  • V. Guruswami, J. D. Hartline, A. R. Karlin, D. Kempe, C. Kenyon and F. McSherry, On profit-maximizing envy-free pricing, SODA 2005. https://dl.acm.org/doi/10.5555/1070432.1070598
  • P. Rusmevichientong, B. Van Roy and P. W. Glynn, A nonparametric approach to multiproduct pricing, Operations Research 54(1), 2006. https://doi.org/10.1287/opre.1050.0236
  • 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.1754
13 thms3 active usersReviewed
🏆Completed
Machine LearningOperations Research·Captain: mikedeng1

The Weighted Majority Algorithm 2: On Any Partition into Segments, WML Makes at Most Σ min(l_i, (log(n/(βγ)) + m_i log(1/β))/log(1/u)) MistakesResearch Paper

Motivation

Prediction with expert advice asks a master algorithm to combine the binary predictions of a pool of nnn prediction algorithms, trial after trial, so that it makes few more mistakes than the best member of the pool, without any probabilistic assumption on the data. Littlestone and Warmuth's Weighted Majority Algorithm WM (Littlestone and Warmuth 1994) keeps a weight per member, predicts with the weighted vote, and multiplies the weights of the wrong members by a factor β<1\beta<1β<1 after each of its own mistakes. Its mistake bound is of order log⁡n+m\log n + mlogn+m, where mmm is the number of mistakes of the best member on the whole sequence.

That bound compares the master with a single fixed member. In many applications the member that predicts well changes over time: one algorithm is accurate for a while, then a second one, then a third. Against a shifting target of this kind, WM can do badly, because the weight of a member that was wrong for a long time has become so small that it takes many mistakes to recover it. Section 3 of the paper introduces the variant WML, which never decreases a weight that is already small relative to the total, and proves that it pays only a logarithmic cost per change of target. This was the first result of what is now called tracking the best expert; it was later refined by Herbster and Warmuth's Fixed-Share algorithm (Herbster and Warmuth 1998).

Timeline. 1989: a conference version of the paper appears at FOCS. 1994: the journal version, the source of this mission, gives Lemma 3.1 and Theorem 3.1 in their present form. 1998: Herbster and Warmuth's Tracking the Best Expert generalizes the shifting bound to other loss functions.

Setting

A run takes place over TTT trials t=0,…,T−1t = 0, \dots, T-1t=0,…,T−1. In trial ttt each member i∈{1,…,n}i \in \{1,\dots,n\}i∈{1,…,n} of the pool predicts xi(t)∈{0,1}x_i^{(t)} \in \{0,1\}xi(t)​∈{0,1}, the master predicts λ(t)∈{0,1}\lambda^{(t)} \in \{0,1\}λ(t)∈{0,1}, and the label ρ(t)∈{0,1}\rho^{(t)} \in \{0,1\}ρ(t)∈{0,1} is revealed. The predictions of the pool and the labels are arbitrary.

The master keeps weights wi(t)>0w_i^{(t)} > 0wi(t)​>0, the weights at the beginning of trial ttt. Let W(t)=∑jwj(t)W^{(t)} = \sum_j w_j^{(t)}W(t)=∑j​wj(t)​ be the total weight and qb(t)q_b^{(t)}qb(t)​ the total weight of the members predicting bbb. WML has parameters 0<β<10 < \beta < 10<β<1 and 0≤γ<120 \le \gamma < \tfrac120≤γ<21​:

  1. it predicts 000 if q0(t)>q1(t)q_0^{(t)} > q_1^{(t)}q0(t)​>q1(t)​, 111 if q1(t)>q0(t)q_1^{(t)} > q_0^{(t)}q1(t)​>q0(t)​, and either value on a tie;
  2. if λ(t)≠ρ(t)\lambda^{(t)} \ne \rho^{(t)}λ(t)=ρ(t), each member iii with xi(t)≠ρ(t)x_i^{(t)} \ne \rho^{(t)}xi(t)​=ρ(t) and wi(t)>γnW(t)w_i^{(t)} > \frac{\gamma}{n} W^{(t)}wi(t)​>nγ​W(t) has its weight multiplied by β\betaβ; all other weights stay;
  3. if λ(t)=ρ(t)\lambda^{(t)} = \rho^{(t)}λ(t)=ρ(t), no weight changes.

With γ=0\gamma = 0γ=0 this is WM. Write mmm for the number of mistakes of the master and

u=1+β2+(1−β)γ,u = \frac{1+\beta}{2} + (1-\beta)\gamma ,u=21+β​+(1−β)γ,

which lies in [12,1)[\tfrac12, 1)[21​,1) for the stated parameter ranges.

Formalization targets

Goal: Theorem 3.1 (p. 224)

Split the sequence into consecutive segments S1,…,Sk\mathcal S_1, \dots, \mathcal S_kS1​,…,Sk​, with lil_ili​ trials in Si\mathcal S_iSi​, and let mim_imi​ be the smallest number of mistakes any single pool member makes on Si\mathcal S_iSi​. If γ>0\gamma > 0γ>0 and every initial weight is at least βγ/n\beta\gamma/nβγ/n times the total initial weight, then

m≤∑i=1kmin⁡(li, log⁡(n/(βγ))+milog⁡(1/β)log⁡(1/u)).m \le \sum_{i=1}^{k} \min\left( l_i,\ \frac{\log\big(n/(\beta\gamma)\big) + m_i \log(1/\beta)}{\log(1/u)} \right).m≤i=1∑k​min(li​, log(1/u)log(n/(βγ))+mi​log(1/β)​).

The partition is arbitrary and unknown to the algorithm; the bound holds for all partitions simultaneously.

Milestone: Lemma 3.1 (pp. 223–224)

The single-segment case, k=1k = 1k=1, together with the statement that the weight floor survives the run: if every initial weight is at least βγ/n\beta\gamma/nβγ/n times the total, then

m≤log⁡(n/(βγ))+m0log⁡(1/β)log⁡(1/u),m \le \frac{\log\big(n/(\beta\gamma)\big) + m_0 \log(1/\beta)}{\log(1/u)},m≤log(1/u)log(n/(βγ))+m0​log(1/β)​,

and every final weight is at least βγ/n\beta\gamma/nβγ/n times the total final weight.

Milestone: the contraction step (p. 224, proof of Lemma 3.1)

On every mistake trial, W(t+1)≤u W(t)W^{(t+1)} \le u\, W^{(t)}W(t+1)≤uW(t).

Significance

The result. Theorem 3.1 shows that an online algorithm with no knowledge of when or how often the best predictor changes can compete with the best sequence of predictors, at a price of O(log⁡(n/(βγ)))O(\log(n/(\beta\gamma)))O(log(n/(βγ))) mistakes per segment. The min⁡\minmin with lil_ili​ says that a segment never costs more than its length. The bound is the template for later tracking and adaptive-regret results in online learning and for algorithms whose weights are kept bounded away from zero.

Formalizing it. The result is proved in the paper; no machine-checked proof of a shifting-target bound exists on Prove2Me or, to our knowledge, in Mathlib. This mission produces a formal model of WML in which ties are arbitrary (a run is a relation, so the bounds cover every tie-breaking rule), the per-mistake contraction, the single-sequence bound with its preserved invariant, and the segment bound. The model of trials, majority votes and mistake counts is reusable for the other algorithms of the paper.

Difficulty

The standard argument for WM compares the decrease of the total weight with the final weight of the best member. For a shifting target this fails twice. First, the best member of a later segment may have been wrong during all earlier segments; under WM its weight has shrunk by β\betaβ per mistake, and the master has to make many mistakes before that member regains influence. Second, the floor that prevents this also prevents some wrong members from being penalized, so the total weight shrinks by less on each mistake than under WM. The statement has to control both effects at once, and the per-segment accounting has to restart from whatever weights the previous segments left behind. A natural first idea, applying the single-sequence bound to each part of an arbitrary (interleaved) partition, does not work: for interleaved subsequences the inequality is false, and the theorem is about consecutive segments.

Formalization scope

The pool is Fin n with n>0n > 0n>0; trials are Fin T; binary values are Bool; weights are real, w : ℕ → Fin n → ℝ, with w t the weights at the beginning of trial ttt, w 0 the initial and w T the final weights. A run is the relation IsWMLRun β γ x ρ w lam; it fixes positive initial weights, the majority rule with arbitrary ties, and the update, whose threshold is strict and uses the total weight at the beginning of the trial. Logarithms are natural logarithms; the base cancels in every ratio. The minimum over the pool is an infimum in ℕ over Fin n. Segments are given by monotone boundaries b : Fin (k+1) → ℕ with b 0 = 0 and b (Fin.last k) = T; empty segments are allowed and contribute 000.

The paper allows γ=0\gamma = 0γ=0 for the algorithm, but then the bounds read log⁡(n/0)=+∞\log(n/0) = +\inftylog(n/0)=+∞. In Lean, n / 0 = 0, so the bounds assume γ>0\gamma > 0γ>0; the contraction milestone holds for 0≤γ<120 \le \gamma < \tfrac120≤γ<21​ and is stated so. The hypothesis γ<12\gamma < \tfrac12γ<21​ is kept as stated in the paper, not replaced by u<1u < 1u<1. A trivializing formalization, such as fixing the tie-breaking rule, dropping the strict weight test, letting the floor hypothesis be unsatisfiable, or allowing interleaved subsequences (which makes the goal false), is ruled out: a run exists for every input from every positive initial weights, and equal initial weights satisfy the floor.

Contributions welcome: proofs of the three statements, and lemmas on restricting a WML run to a block of trials, which the goal needs and which apply to any segment-wise analysis.

Selected references

  • N. Littlestone, M. K. Warmuth, The Weighted Majority Algorithm, Information and Computation 108(2), 212–261, 1994. https://doi.org/10.1006/inco.1994.1009
  • N. Littlestone, M. K. Warmuth, The Weighted Majority Algorithm, Proc. 30th IEEE FOCS, 256–261, 1989. https://doi.org/10.1109/SFCS.1989.63487
  • M. Herbster, M. K. Warmuth, Tracking the Best Expert, Machine Learning 32, 151–178, 1998. https://doi.org/10.1023/A:1007424614876
4 thms3 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchProbability·Captain: mikedeng1

Online Stochastic Matching: Beating 1-1/e 3: No Online Algorithm Beats Expected Approximation Factor 26/27 on a 6-Cycle, and None Reaches 1 − o(1) on Disjoint 6-CyclesResearch Paper

Motivation

Online bipartite matching asks an algorithm to assign arriving requests to resources it cannot reassign later. In display advertising, the requests are page views (impressions) and the resources are advertisers who have bought a fixed number of impressions in advance; each impression must be served immediately or lost. In the adversarial model, Karp, Vazirani and Vazirani (STOC 1990) showed that the RANKING algorithm achieves a 1−1/e1 - 1/e1−1/e fraction of the optimum in expectation, and that no online algorithm does better. Advertising systems, however, have historical traffic data, which motivates the i.i.d. model: the impression types are drawn independently from a distribution that is known in advance.

Feldman, Mehta, Mirrokni and Muthukrishnan (arXiv:0905.4100, FOCS 2009) were the first to beat 1−1/e1 - 1/e1−1/e in this model, with an algorithm that achieves ≈0.67\approx 0.67≈0.67 with high probability. Their Section 3 asks the complementary question: how close to 111 can any online algorithm get when the distribution is known? Their Theorem 3 answers that the expected approximation factor of every online algorithm is bounded strictly away from 111, already on a graph with six vertices. Later work in the same model (Manshadi, Oveis Gharan and Saberi, SODA 2011, arXiv:1007.1673) sharpened both the algorithmic and the hardness side.

Setting

An instance consists of a bipartite graph G=(A,I,E)G = (A, I, E)G=(A,I,E) between a finite set AAA of advertisers and a finite set III of impression types, a distribution DDD on III, and a number nnn of arrivals. In this mission DDD is the uniform distribution on III. Online, nnn impressions arrive one at a time; their types ω(0),…,ω(n−1)\omega(0), \dots, \omega(n-1)ω(0),…,ω(n−1) are independent draws from DDD. When impression ttt arrives, the algorithm must immediately either assign it to an advertiser aaa with (a,ω(t))∈E(a, \omega(t)) \in E(a,ω(t))∈E that has not yet been used, or leave it unassigned. Each advertiser can be used at most once, and decisions are final.

A deterministic online algorithm decides what to do with arrival ttt from the types ω(0),…,ω(t)\omega(0), \dots, \omega(t)ω(0),…,ω(t) seen so far; it knows GGG, DDD and nnn, but not the future. A randomized online algorithm is a probability distribution over deterministic ones. ALG(ω)\mathrm{ALG}(\omega)ALG(ω) is the number of impressions the algorithm assigns on the arrival sequence ω\omegaω. OPT(ω)\mathrm{OPT}(\omega)OPT(ω) is the size of a maximum matching of the realization graph, which has one node per arrival ttt, joined to every advertiser adjacent to ω(t)\omega(t)ω(t): the most impressions that could have been assigned with hindsight. The expected approximation factor of an algorithm is

E ⁣[ALG(ω)OPT(ω)],\mathbb E\!\left[\frac{\mathrm{ALG}(\omega)}{\mathrm{OPT}(\omega)}\right],E[OPT(ω)ALG(ω)​],

the expectation taken over the algorithm's randomness and the arrivals.

The 6-cycle instance has A={a,b,c}A = \{a, b, c\}A={a,b,c}, I={x,y,z}I = \{x, y, z\}I={x,y,z} and E={(x,a),(y,a),(y,b),(z,b),(z,c),(x,c)}E = \{(x,a),(y,a),(y,b),(z,b),(z,c),(x,c)\}E={(x,a),(y,a),(y,b),(z,b),(z,c),(x,c)}, with the uniform distribution and n=3n = 3n=3. The family Γk\Gamma_kΓk​ consists of kkk disjoint copies of the 6-cycle, with the uniform distribution on its 3k3k3k impression types and n=3kn = 3kn=3k arrivals.

Formalization targets

Goal: Theorem 3

On the 6-cycle, every randomized online algorithm satisfies

E ⁣[ALGOPT]≤2627,\mathbb E\!\left[\frac{\mathrm{ALG}}{\mathrm{OPT}}\right] \le \frac{26}{27},E[OPTALG​]≤2726​,

and there are a constant c<1c < 1c<1 and a threshold K0K_0K0​ such that for every k≥K0k \ge K_0k≥K0​ and every randomized online algorithm on Γk\Gamma_kΓk​,

E ⁣[ALGOPT]≤c.\mathbb E\!\left[\frac{\mathrm{ALG}}{\mathrm{OPT}}\right] \le c .E[OPTALG​]≤c.

The second part is the paper's "there exists a family of instances with n→∞n \to \inftyn→∞ for which no algorithm can achieve an expected approximation of 1−o(1)1 - o(1)1−o(1)", stated on the paper's own family. It leaves ccc unspecified: Appendix B estimates c≈0.9898c \approx 0.9898c≈0.9898 through approximate counts, and a sharper constant would not invalidate the goal.

Milestones

  1. The (x,y,y)(x, y, y)(x,y,y) scenario (§3, p. 4). For every deterministic online algorithm on the 6-cycle and every type uuu of the first arrival there is a type vvv with ALG(u,v,v)≤2\mathrm{ALG}(u, v, v) \le 2ALG(u,v,v)≤2 and OPT(u,v,v)=3\mathrm{OPT}(u, v, v) = 3OPT(u,v,v)=3.
  2. Theorem 3, first sentence (p. 4). The bound 26/2726/2726/27 on the 6-cycle, for every randomized online algorithm.

Significance

Theorem 3 sets the ceiling against which algorithms in the i.i.d. model are measured. It rules out an online algorithm with expected factor arbitrarily close to 111, even though the distribution is known and the instance is tiny, so a constant-factor gap between online and offline matching is intrinsic to the model and not an artefact of adversarial arrivals. The paper's positive results (1−1/e1 - 1/e1−1/e for the Suggested Matching algorithm, ≈0.67\approx 0.67≈0.67 for Two Suggested Matchings) sit between 1−1/e1 - 1/e1−1/e and this ceiling.

The single-instance bound has a short counting proof. The family statement is argued in the paper only in outline: Appendix B approximates the fractions of copies receiving 111, 222, 333 or more impressions, assumes the most favourable outcome on each, and reports the resulting ratio as approximately 0.98980.98980.9898. A machine-checked proof of the family statement would turn that sketch into a theorem with an explicit constant. To our knowledge neither part has been formalized before.

Difficulty

The single-cycle bound reduces to a finite check, but the paper's "without loss of generality (from the symmetry of the 6-cycle)" hides the cases the algorithm can choose: assigning the first impression to either neighbour, or not assigning it at all. Each case needs its own bad continuation. The passage from deterministic to randomized algorithms is an averaging step.

The family bound is harder. On Γk\Gamma_kΓk​ the algorithm sees all arrivals in all copies and may coordinate its decisions across copies, so the per-copy loss of the single cycle does not transfer by independence. Moreover the target is the expectation of a ratio, E[ALG/OPT]\mathbb E[\mathrm{ALG}/\mathrm{OPT}]E[ALG/OPT], not a ratio of expectations: a bound on the expected loss must be combined with concentration of the number of copies that receive exactly three impressions, and with a lower bound on OPT\mathrm{OPT}OPT that holds with high probability. The approximations "≃3/e3\simeq 3/e^3≃3/e3, 9/(2e3)9/(2e^3)9/(2e3), 27/(6e3)27/(6e^3)27/(6e3)" of Appendix B have unquantified errors and cannot be used as they stand.

Formalization scope

The model is formalized for general finite AAA and III with uniform arrivals. Arrival sequences are functions Fin n → I. A deterministic online algorithm is a function that receives the time ttt, the prefix of the first ttt types (Fin t → I) and the current type, and returns an advertiser or none; it cannot read later arrivals. A proposal of a non-adjacent or already used advertiser leaves the impression unassigned, so the encoding covers all online algorithms, including ones that skip an impression while a neighbour is free. A randomized online algorithm is a probability distribution (PMF) on the finite type of deterministic algorithms; by Kuhn's theorem this is equivalent to fresh random choices at each step. OPT\mathrm{OPT}OPT is the maximum size of an injective, edge-respecting partial assignment of arrivals to advertisers. The expected approximation factor is a finite average over all ∣I∣n|I|^n∣I∣n arrival sequences, weighted by the algorithm distribution.

The explicit instantiations of the paper's asymptotic wording are as follows:

  • "no algorithm can achieve an expected approximation of 1−o(1)1 - o(1)1−o(1)" becomes ∃ c<1, ∃K0, ∀k≥K0, ∀P: E[ALG/OPT]≤c\exists\, c < 1,\ \exists K_0,\ \forall k \ge K_0,\ \forall P:\ \mathbb E[\mathrm{ALG}/\mathrm{OPT}] \le c∃c<1, ∃K0​, ∀k≥K0​, ∀P: E[ALG/OPT]≤c on Γk\Gamma_kΓk​, with n=3kn = 3kn=3k tied to kkk, and with ccc and K0K_0K0​ chosen before kkk and before the algorithm.
  • The paper's constant 0.98980.98980.9898 is not part of the statement.

Lean's division sets x/0=0x/0 = 0x/0=0. A formalization of the form "there is an instance on which every algorithm has expected factor at most 26/2726/2726/27" would be satisfied trivially by a graph without edges, where OPT=0\mathrm{OPT} = 0OPT=0; the statements here are on the paper's explicit instances, where every impression type has an adjacent advertiser and OPT≥1\mathrm{OPT} \ge 1OPT≥1 on every arrival sequence.

A complete development needs finite case analysis on runs of online algorithms, averaging over mixed strategies, and, for the family, concentration inequalities for occupancy counts (Azuma or McDiarmid, both available on the platform) together with bounds on maximum matchings of disjoint unions. The online-algorithm model and the averaging lemmas are reusable for other lower bounds in the i.i.d. model. Proofs of either milestone, and independent proofs of the family statement with any explicit c<1c < 1c<1, are welcome.

Selected references

  • J. Feldman, A. Mehta, V. Mirrokni, S. Muthukrishnan, Online Stochastic Matching: Beating 1-1/e, FOCS 2009; arXiv:0905.4100v1. https://arxiv.org/abs/0905.4100
  • R. M. Karp, U. V. Vazirani, V. V. Vazirani, An optimal algorithm for on-line bipartite matching, STOC 1990. https://doi.org/10.1145/100216.100262
  • V. H. Manshadi, S. Oveis Gharan, A. Saberi, Online Stochastic Matching: Online Actions Based on Offline Statistics, SODA 2011; arXiv:1007.1673. https://arxiv.org/abs/1007.1673
5 thms3 active usersReviewed
🏆Completed
CombinatoricsGraph Theory·Captain: mikedeng1

On the Approximability of Single-Machine Scheduling with Precedence Constraints 3: The Incomparable-Pairs Graphs of Canonical Interval Orders Have Unbounded Chromatic NumberResearch Paper

Motivation

Single-machine scheduling with precedence constraints asks for an order in which one machine processes jobs while respecting prescribed comparisons between them. For the objective of minimizing total weighted completion time, the structure of the precedence order can influence the quality of approximation algorithms. Ambühl, Mastrolilli, Mutsanas, and Svensson connect one such approach to coloring a graph built from the order's incomparable pairs. A bounded number of colors would make a certain class of coloring-based guarantees uniform over the class of precedence orders. Their Theorem 4.3 shows that canonical interval orders have no such uniform bound for this particular graph. Ambühl et al., §4.2, p. 658

An interval order is a partial order whose elements can be represented by closed real intervals, with one element strictly before another when the first interval ends before the second begins. Interval orders are a familiar restriction of general precedence systems. The paper distinguishes them from semiorders, represented by intervals of equal length, for which it reports a three-element realizer bound and a corresponding scheduling guarantee. The chromatic obstruction here explains why that particular bounded-color route cannot be extended uniformly from semiorders to all interval orders. It does not assert that every scheduling approach to interval orders fails. Ambühl et al., §§3–4.2, pp. 656–658

Setting

For an integer n≥2n\ge2n≥2, let [n][n][n] be a linearly ordered set with nnn points. The canonical interval order InI_nIn​ has one element for every pair of distinct points a<ba<ba<b in [n][n][n], viewed as the closed interval [a,b][a,b][a,b]. For two such intervals, write S≤InTS\le_{I_n}TS≤In​​T when S=TS=TS=T or the right endpoint of SSS is strictly less than the left endpoint of TTT. Thus intervals that touch at an endpoint remain incomparable. The formalization represents each interval by its two-element endpoint set, and its order relation is equivalent to this endpoint rule. Ambühl et al., §4.2, p. 658

For any partial order PPP on a ground set NNN, an ordered incomparable pair is (x,y)(x,y)(x,y) with neither x≤Pyx\le_P yx≤P​y nor y≤Pxy\le_P xy≤P​x. Both (x,y)(x,y)(x,y) and (y,x)(y,x)(y,x) are vertices when xxx and yyy are incomparable. A linear extension of PPP is a linear order on the same elements that contains all comparisons of PPP. It reverses (x,y)(x,y)(x,y) if it places yyy before xxx. The graph GPG_PGP​ has the ordered incomparable pairs as vertices. Two distinct vertices form an edge when no linear extension reverses both, while each singleton can be reversed by some linear extension. This singleton condition is the minimality clause in the paper's definition through its hypergraph of incomparable pairs. Ambühl et al., §2.1, p. 655, §3, p. 656

A proper kkk-coloring assigns one of kkk colors to every vertex of GPG_PGP​ so that endpoints of an edge have different colors. The graph's chromatic number χ(GP)\chi(G_P)χ(GP​) is the least positive number of colors that works, when the graph is nonempty. The mission uses the equivalent predicate that no proper map into a fixed kkk-element color set exists once nnn is large enough. This also gives a clear interpretation at k=0k=0k=0.

Formalization targets

The goal is Theorem 4.3: for each fixed integer kkk, all sufficiently large canonical interval orders have incomparable-pairs graphs with chromatic number greater than kkk. With k≥0k\ge0k≥0, the Lean target is

∀k∈N, ∃n0≥2, ∀n≥n0,GIn is not k-colorable.\forall k\in\mathbb N,\ \exists n_0\ge2,\ \forall n\ge n_0,\quad G_{I_n}\text{ is not $k$-colorable}.∀k∈N, ∃n0​≥2, ∀n≥n0​,GIn​​ is not k-colorable.

The paper states kkk as an integer. Its negative cases are immediate from nonnegativity of chromatic numbers; the target states the substantive range. The threshold n0n_0n0​ may depend on kkk, while InI_nIn​ and GInG_{I_n}GIn​​ are constructed from nnn rather than supplied as arbitrary objects. Ambühl et al., Theorem 4.3, p. 658

Two assertions from the theorem's proof serve as milestones. First, χ(GIn)\chi(G_{I_n})χ(GIn​​) is nondecreasing as nnn increases. Second, for four endpoints i<j<ℓ<mi<j<\ell<mi<j<ℓ<m, the vertices ({i,j},{j,ℓ})(\{i,j\},\{j,\ell\})({i,j},{j,ℓ}) and ({j,ℓ},{ℓ,m})(\{j,\ell\},\{\ell,m\})({j,ℓ},{ℓ,m}) are adjacent in GInG_{I_n}GIn​​. They are recorded with their wording and provenance from §4.2. A separately published theorem on hypergraph Ramsey numbers is included as a reference because the source uses that theorem in its large-nnn argument. Ambühl et al., proof of Theorem 4.3, p. 658

Significance

The result provides a precise structural limitation: across the canonical interval orders, the chromatic numbers of GInG_{I_n}GIn​​ cannot be bounded by one constant. The graph GPG_PGP​ is a specific graph of ordered incomparable pairs, distinct from other graphs that can be associated with a scheduling instance. Thus the conclusion addresses the bounded-color approach attached to this graph, without asserting an algorithmic impossibility for interval-order scheduling. The paper contrasts this limitation with the bounded-realizer situation for semiorders. Ambühl et al., §§3–4.2, pp. 656–658

The theorem is proved in the 2011 paper. The remaining work in this mission is a machine-checked development of its definition layer, its two stated supporting assertions, and the unbounded-color conclusion. The partial-order property of InI_nIn​ is already proved in the definition file; the three theorem statements are open proof targets in the proposal. A completed formalization would also give reusable components for finite interval orders, ordered incomparable pairs, and coloring arguments in dimension theory.

Difficulty

A pair of intervals overlapping or touching is easy to recognize as incomparable, but graph adjacency has a stronger meaning: it quantifies over every linear extension of the entire partial order. A local test based only on the four displayed intervals would silently change the graph. Likewise, the fact that two vertices have different endpoint descriptions does not by itself make them adjacent. The formal argument must connect the canonical interval representation to the global extension-based edge condition, then show that a bounded coloring cannot persist as the endpoint set grows. Ambühl et al., proof of Theorem 4.3, p. 658

Formalization scope

The endpoint set is Fin n, hence indexed 0,…,n−10,\ldots,n-10,…,n−1 rather than the paper's 1,…,n1,\ldots,n1,…,n; the shift preserves order. Elements of InI_nIn​ are exactly two-element subsets of this set. The order relation includes equality and strict separation of endpoint sets and is proved to be a partial order. The graph is defined for an arbitrary binary relation so it can be reused, but every target here applies it to the concrete partial order InI_nIn​. Vertices are ordered incomparable pairs, and adjacency retains both singleton reversibility and failure of simultaneous reversibility. Defining adjacency directly from the desired four-point pattern would trivialize the target and is excluded.

Colorable n k means existence of a proper function from all vertices of GInG_{I_n}GIn​​ to Fin k. For an empty graph, zero-colorability is allowed by this representation. The theorem therefore includes a threshold n0≥2n_0\ge2n0​≥2 and quantifies over every n≥n0n\ge n_0n≥n0​; the assertion itself forces the threshold past any empty-graph exceptions. The formal target uses natural kkk; the paper's negative integer values need no separate theorem. No asymptotic notation or unspecified constants occur. The source's Ramsey number R(3:4,…,4)R(3:4,\ldots,4)R(3:4,…,4) is not hard-coded into the goal, because the source asserts existence of a suitable threshold, not that the least threshold equals that value.

The printed adjacency justification names an alternating-cycle pair set different from the two adjacent vertices in the preceding sentence. The milestone formalizes the adjacency assertion itself and records the printed explanation verbatim for review; it does not turn the mismatched pair set into a hypothesis. The mission needs finite order embeddings, linear extensions, graph colorings, and a finite Ramsey theorem. The existing published Ramsey theorem is referenced as a reusable result. No scheduling approximation algorithm, polynomial-time statement, or numerical scheduling guarantee is formalized in this mission.

Selected references

  • Christoph Ambühl, Monaldo Mastrolilli, Nikolaos Mutsanas, and Ola Svensson, On the Approximability of Single-Machine Scheduling with Precedence Constraints, Mathematics of Operations Research 36(4), 653–669, 2011. DOI: 10.1287/moor.1110.0512
7 thms3 active usersReviewed
🏆Completed
CombinatoricsGraph TheoryOperations Research·Captain: mikedeng1

On the Approximability of Single-Machine Scheduling with Precedence Constraints 1: A k-Fold Realizer of Size t Yields a Vertex Cover of Expected Weight at Most (2 − 2/(t/k)) Times OptimalResearch Paper

Motivation

Single-machine scheduling with precedence constraints, written 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_jC_j1∣prec∣∑wj​Cj​ in the notation of Graham et al., asks for an order of nnn weighted jobs on one machine that respects a given partial order and minimizes the weighted sum of completion times. The problem is strongly NP-hard (Lawler 1978; Lenstra and Rinnooy Kan 1978), and closing its approximability gap is listed by Schuurman and Woeginger among ten outstanding open problems in scheduling theory. Several 2-approximation algorithms are known (Schulz 1996; Hall et al. 1997; Chudak and Hochbaum 1999; Chekuri and Motwani 1999; Margot et al. 2003).

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 built from the precedence order. Ambühl, Mastrolilli, Mutsanas and Svensson (Math. Oper. Res. 36(4), 2011) observed that this graph is the graph of incomparable pairs of dimension theory, and used that identification to obtain (2−2/f)(2-2/f)(2−2/f)-approximations for orders of fractional dimension at most fff. This mission formalizes that framework: the identification of the two graphs and the rounding guarantee of the paper's Theorem 5.1.

Setting

An instance SSS consists of a finite set NNN of jobs, a partial order PPP on NNN (reflexive, antisymmetric, transitive; (i,j)∈P(i,j)\in P(i,j)∈P with i≠ji\ne ji=j means job iii finishes before job jjj starts), processing times pj≥0p_j\ge 0pj​≥0 and weights wj≥0w_j\ge 0wj​≥0.

Two jobs x,yx,yx,y are incomparable, x∥yx\parallel yx∥y, when neither (x,y)(x,y)(x,y) nor (y,x)(y,x)(y,x) lies in PPP. The set inc⁡(P)\operatorname{inc}(P)inc(P) of incomparable pairs consists of ordered pairs and is closed under swapping. A linear extension of PPP is a linear order L⊇PL\supseteq PL⊇P on NNN; it reverses (x,y)∈inc⁡(P)(x,y)\in\operatorname{inc}(P)(x,y)∈inc(P) when y<xy<xy<x in LLL. A nonempty multiset L1,…,LtL_1,\dots,L_tL1​,…,Lt​ of linear extensions is a k:tk:tk:t-realizer if every incomparable pair is reversed by at least kkk of them. The fractional dimension fdim⁡(P)\operatorname{fdim}(P)fdim(P) is the least ratio t/kt/kt/k over all k:tk:tk:t-realizers.

The vertex cover graph GPSG^S_PGPS​ has the incomparable pairs as nodes; 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,ℓ),(k,j)∈P(i,\ell),(k,j)\in P(i,ℓ),(k,j)∈P (symmetrically closed). Node (i,j)(i,j)(i,j) has weight w(i,j)=piwjw_{(i,j)}=p_iw_jw(i,j)​=pi​wj​, and w(C)=∑u∈Cwuw(C)=\sum_{u\in C}w_uw(C)=∑u∈C​wu​. OPT\mathrm{OPT}OPT is the minimum weight of a vertex cover of GPSG^S_PGPS​. The LP relaxation [CS-LP] asks for x∈[0,1]inc⁡(P)x\in[0,1]^{\operatorname{inc}(P)}x∈[0,1]inc(P) with xu+xv≥1x_u+x_v\ge1xu​+xv​≥1 on every edge, minimizing ∑uwuxu\sum_u w_ux_u∑u​wu​xu​. For a solution xxx write Va={u:xu=a}V_a=\{u: x_u=a\}Va​={u:xu​=a}, and for a linear extension LLL let I1/2(L)I_{1/2}(L)I1/2​(L) be the pairs of V1/2V_{1/2}V1/2​ reversed in LLL.

The graph of incomparable pairs GPG_PGP​ (Felsner and Trotter 2000) also has the incomparable pairs as vertices; two of them are adjacent when the pair of them is a minimal set of incomparable pairs that no linear extension reverses entirely.

Formalization targets

Goal: Theorem 5.1 (p. 658)

For an instance SSS, a k:tk:tk:t-realizer L1,…,LtL_1,\dots,L_tL1​,…,Lt​ of PPP, and a half-integral optimal solution xxx of [CS-LP], put Ci=V1∪(V1/2∖I1/2(Li))C_i=V_1\cup\bigl(V_{1/2}\setminus I_{1/2}(L_i)\bigr)Ci​=V1​∪(V1/2​∖I1/2​(Li​)). Then every CiC_iCi​ is a vertex cover of GPSG^S_PGPS​, and

1t∑i=1tw(Ci)  ≤  (2−2t/k)OPT.\frac1t\sum_{i=1}^t w(C_i)\;\le\;\Bigl(2-\frac{2}{t/k}\Bigr)\mathrm{OPT}.t1​i=1∑t​w(Ci​)≤(2−t/k2​)OPT.

Milestones

  • Proposition 3.2 (p. 657): GPS=GPG^S_P=G_PGPS​=GP​.
  • Footnote 4 (p. 659): the pairs reversed by a linear extension are independent in GPSG^S_PGPS​.
  • Eq. (4): 1t∣{i:Li reverses u}∣≥k/t\frac1t|\{i: L_i\text{ reverses }u\}|\ge k/tt1​∣{i:Li​ reverses u}∣≥k/t for every incomparable pair uuu.
  • Eq. (5): 1t∑iw(I1/2(Li))≥kt w(V1/2)\frac1t\sum_i w(I_{1/2}(L_i))\ge \frac kt\,w(V_{1/2})t1​∑i​w(I1/2​(Li​))≥tk​w(V1/2​).
  • Hochbaum's observation (§5, p. 659): for half-integral feasible xxx, V1∪CV_1\cup CV1​∪C covers GPSG^S_PGPS​ whenever CCC covers GPS[V1/2]G^S_P[V_{1/2}]GPS​[V1/2​].
  • Eqs. (6)–(8): 1t∑iw(Ci)≤w(V1)+(1−kt)w(V1/2)≤2(1−kt)(w(V1)+12w(V1/2))≤(2−2t/k)OPT\frac1t\sum_iw(C_i)\le w(V_1)+(1-\frac kt)w(V_{1/2})\le 2(1-\frac kt)(w(V_1)+\frac12w(V_{1/2}))\le(2-\frac2{t/k})\mathrm{OPT}t1​∑i​w(Ci​)≤w(V1​)+(1−tk​)w(V1/2​)≤2(1−tk​)(w(V1​)+21​w(V1/2​))≤(2−t/k2​)OPT when PPP is not a linear order.

Significance

The result. Combined with the cited Theorem 2.1 (Ambühl–Mastrolilli 2009; Correa–Schulz 2005), which turns an α\alphaα-approximate vertex cover of GPSG^S_PGPS​ into an α\alphaα-approximate schedule, Theorem 5.1 gives a (2−2/f)(2-2/f)(2−2/f)-approximation for 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_jC_j1∣prec∣∑wj​Cj​ whenever the precedence order has an efficiently samplable realizer with t/k≤ft/k\le ft/k≤f. The paper applies it to interval orders (3/23/23/2), convex bipartite orders and semiorders (4/34/34/3), orders of bounded degree and orders of interval dimension two; for the earlier special classes it matched or improved the best known ratios, and for the last two it gave the first results. Proposition 3.2 makes the dimension theory of posets (realizers, critical pairs, fractional dimension) directly available to the vertex cover approach.

Formalizing it. The results are proved in the paper; to the best of current knowledge none of them has a machine-checked proof. The mission produces a Lean development of incomparable pairs, linear extensions, kkk-fold realizers and the hypergraph of incomparable pairs, which other dimension-theory missions can reuse, and a verified LP-rounding argument for half-integral vertex cover solutions under a distribution of independent sets.

Difficulty

Most of the rounding argument is arithmetic over finite sums. The central nontrivial step is the inclusion GP⊆GPSG_P\subseteq G^S_PGP​⊆GPS​ in Proposition 3.2: for two incomparable pairs that are not adjacent under the three-case rule, one must construct a single linear extension reversing both. This needs an extension of PPP by two new comparabilities whose transitive closure is still antisymmetric, followed by Szpilrajn's theorem; checking that every potential cycle is excluded by the three cases is the actual content. The opposite inclusion, and footnote 4, follow from transitivity and antisymmetry of linear orders. A second point of care is the inequality t≥2kt\ge 2kt≥2k used in step (7): it is not part of the definition of a realizer, and follows from each linear extension reversing exactly one of (x,y)(x,y)(x,y) and (y,x)(y,x)(y,x).

Formalization scope

  • Jobs form a finite type N; the precedence order is a relation P : N → N → Prop with IsPartialOrder, carried by the structure Instance. Processing times and weights are nonnegative reals.
  • inc⁡(P)\operatorname{inc}(P)inc(P) is the subtype IncPair P of N × N; a linear extension is a relation with IsLinearOrder containing P; a k:tk:tk:t-realizer is a family Fin t → LinearExtension P with t>0t>0t>0. Reversal of (x,y)(x,y)(x,y) means y<xy<xy<x in LLL throughout; the page's "y>xy>xy>x" in Eq. (4) and "Prob[j>i]\mathrm{Prob}[j>i]Prob[j>i]" in Eq. (5) are the same family of inequalities because inc⁡(P)\operatorname{inc}(P)inc(P) is symmetric.
  • GPSG^S_PGPS​ is the symmetric closure of the printed three-case rule on distinct nodes. GPG_PGP​ is defined through linear extensions and hyperedge minimality, never through the three-case rule, so Proposition 3.2 is a genuine statement and not a definitional equality.
  • [CS-LP] drops the constant term ∑jpjwj+∑(i,j)∈Ppiwj\sum_jp_jw_j+\sum_{(i,j)\in P}p_iw_j∑j​pj​wj​+∑(i,j)∈P​pi​wj​ of [CS-IP], which does not affect optimality. OPT\mathrm{OPT}OPT is the minimum weight of a vertex cover of GPSG^S_PGPS​, taken over a finite nonempty family.
  • Not formalized: "efficiently samplable", "polynomial time" and "randomized algorithm". The expectation over a uniformly sampled LiL_iLi​ is stated as the average 1t∑i=1t\frac1t\sum_{i=1}^tt1​∑i=1t​, which is equivalent to it and stronger than the existence of one good index. The existence of a half-integral optimal [CS-LP] solution (Nemhauser–Trotter, cited) is a hypothesis on xxx. The conversion of a vertex cover into a schedule (Theorem 2.1, cited) is not formalized; the goal is stated for vertex covers of GPSG^S_PGPS​.
  • Constants: 2−2/(t/k)2-2/(t/k)2−2/(t/k) in real arithmetic; it equals 2−2k/t2-2k/t2−2k/t, and equals 222 when k=0k=0k=0.
  • The paper assumes fdim⁡(P)≥2\operatorname{fdim}(P)\ge2fdim(P)≥2, i.e. PPP is not a linear order. Eqs. (6)–(8) carry that hypothesis as the paper does; the goal omits it because for a linear order both sides are 000.
  • Conclusion (a), that each CiC_iCi​ is a vertex cover, is part of the goal and is not assumed.

Contributions welcome: proofs of the milestones in any order, and a reusable Szpilrajn-style lemma producing a linear extension that reverses a prescribed set of compatible incomparable pairs.

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), 2009 (reference [2] of the paper).
  • 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
  • G. R. Brightwell, E. R. Scheinerman, Fractional dimension of partial orders, Order 9(2):139–158, 1992 (reference [7]).
  • S. Felsner, W. T. Trotter, Dimension, graph and hypergraph coloring, Order 17(2):167–177, 2000 (reference [13]).
  • D. S. Hochbaum, Efficient bounds for the stable set, vertex cover and set packing problems, Discrete Appl. Math. 6(3):243–254, 1983 (reference [20]).
  • G. L. Nemhauser, L. E. Trotter, Vertex packings: structural properties and algorithms, Math. Programming 8(1):232–248, 1975 (reference [29]).
10 thms3 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchOptimization+1·Captain: mikedeng1

A Tight Linear Time (1/2)-Approximation for Unconstrained Submodular Maximization 2: Randomized Double Greedy Achieves 1/2 of the Optimum in ExpectationResearch Paper

Motivation

Many selection problems assign a value to each subset of a finite collection: the coverage supplied by chosen facilities, the influence reached by chosen seeds, or the value of a coalition. A submodular set function has diminishing returns in the precise sense that the combined value of two sets, counting their overlap once, does not exceed the sum of their separate values. When the function is also monotone, taking more elements never hurts. The unconstrained problem studied here permits nonmonotone functions, so both accepting and rejecting an element can matter. The question is what a single pass through the elements can guarantee when the function is available through value queries. Buchbinder et al., FOCS 2012

The randomized algorithm in this mission attains an expected one-half approximation for every nonnegative submodular function. The paper presents this as tight in the value-oracle setting: it recalls the earlier result of Feige, Mirrokni and Vondrák that a fixed improvement beyond one-half requires exponentially many queries. The contribution here is therefore both the guarantee and a short adaptive rule that attains it in a linear number of iterations. The local proposal follows the FOCS 2012 version of the paper; its theorem numbering differs from the later SIAM Journal on Computing article. Buchbinder et al., §I.A and Theorem I.2

Setting

Let N\mathcal NN be a finite ground set, and let f:2N→R≥0f:2^{\mathcal N}\to\mathbb R_{\ge0}f:2N→R≥0​ assign a nonnegative real value to every subset. The unconstrained submodular maximization problem asks for the largest value f(S)f(S)f(S) among all S⊆NS\subseteq\mathcal NS⊆N. Write OPTOPTOPT for that value when no confusion arises, and OOO for a set attaining it. Submodularity means

f(A∪B)+f(A∩B)≤f(A)+f(B)(A,B⊆N).f(A\cup B)+f(A\cap B)\le f(A)+f(B)\qquad(A,B\subseteq\mathcal N).f(A∪B)+f(A∩B)≤f(A)+f(B)(A,B⊆N).

There is no monotonicity or normalization assumption: f(∅)f(\varnothing)f(∅) and f(N)f(\mathcal N)f(N) may both be positive. A value oracle returns f(S)f(S)f(S) for a requested subset SSS. The paper's complexity claim counts such queries, assuming a query takes constant time. Buchbinder et al., §I and footnotes 1–2

Algorithm 2 visits the elements once in an arbitrary order u1,…,unu_1,\ldots,u_nu1​,…,un​. It keeps two sets, starting at X0=∅X_0=\varnothingX0​=∅ and Y0=NY_0=\mathcal NY0​=N. At step iii, it measures the gain aia_iai​ from adding uiu_iui​ to Xi−1X_{i-1}Xi−1​ and the gain bib_ibi​ from removing uiu_iui​ from Yi−1Y_{i-1}Yi−1​. It clips each gain at zero, giving ai′=max⁡(ai,0)a'_i=\max(a_i,0)ai′​=max(ai​,0) and bi′=max⁡(bi,0)b'_i=\max(b_i,0)bi′​=max(bi​,0). It adds uiu_iui​ to XXX with probability ai′/(ai′+bi′)a'_i/(a'_i+b'_i)ai′​/(ai′​+bi′​) and otherwise removes it from YYY. When both clipped gains vanish, the paper defines the add probability as one. After all elements have been processed, the two sets coincide, and the algorithm returns their common value. The state law is adaptive: its probability at step iii depends on the actual pair of sets produced by earlier choices. Buchbinder et al., Algorithm 2

Formalization targets

The main target is Theorem I.2 for this exact algorithm and for every enumeration of the ground set:

max⁡S⊆Nf(S)≤2 E[f(Xn)].\max_{S\subseteq\mathcal N}f(S)\le 2\,\mathbb E[f(X_n)].S⊆Nmax​f(S)≤2E[f(Xn​)].

The milestone statements retain the paper's key local quantities. For a comparison optimum OOO, set OPTi=(O∪Xi)∩YiOPT_i=(O\cup X_i)\cap Y_iOPTi​=(O∪Xi​)∩Yi​. Lemma II.1 asserts ai+bi≥0a_i+b_i\ge0ai​+bi​≥0. The endpoint statement identifies OPT0=OOPT_0=OOPT0​=O and OPTn=Xn=YnOPT_n=X_n=Y_nOPTn​=Xn​=Yn​. Inequality (3) bounds the conditional loss in the positive-gain case; Lemma III.1 compares the expected change of OPTiOPT_iOPTi​ with the expected combined change of XiX_iXi​ and YiY_iYi​. The telescoped display keeps the initial endpoint values f(∅)f(\varnothing)f(∅) and f(N)f(\mathcal N)f(N) before using nonnegativity. Buchbinder et al., Lemmas II.1 and III.1, inequality (3), proof of Theorem I.2

A companion target is Theorem I.4 via its second proof. For two normalized monotone submodular utilities f1,f2f_1,f_2f1​,f2​, let g(S)=f1(S)+f2(N∖S)g(S)=f_1(S)+f_2(\mathcal N\setminus S)g(S)=f1​(S)+f2​(N∖S). The maximum of ggg is exactly the optimal welfare of a two-player partition. Algorithm 2 on ggg is asked to satisfy

3max⁡S⊆Ng(S)≤4 E[g(Xn)].3\max_{S\subseteq\mathcal N}g(S)\le4\,\mathbb E[g(X_n)].3S⊆Nmax​g(S)≤4E[g(Xn​)].

This is the paper's three-quarter guarantee in its welfare application. Buchbinder et al., Theorem I.4 and Proof (2)

Significance

The main theorem gives a specific randomized rule whose expected value is at least half the best subset value, even when accepting an element can lower the objective. It applies without restricting the cardinality or shape of the chosen subset. The welfare corollary shows that keeping the initial endpoint values in the analysis yields a stronger guarantee for the objective formed from two monotone players. Buchbinder et al., Theorems I.2 and I.4

This mission formalizes the statement of the algorithm, its intermediate state laws, its comparison set, and the paper's numbered proof targets. The algorithmic guarantee is proved in the source paper; the local Lean theorem files are open statements with sorry and do not yet give machine-checked proofs of these results. A completed development would supply a reusable formal model of an adaptive finite random process over pairs of subsets, as well as the specific submodular inequalities. The published Submodular and OPT definitions from the earlier Feige–Mirrokni–Vondrák formalization are reused here.

Difficulty

The two possible updates cannot be assessed independently. The probability of each choice depends on the current state, and the comparison set OPTiOPT_iOPTi​ can gain or lose the processed element in a way that differs from the two algorithm sets. A bound on the expected value of XiX_iXi​ alone does not control the movement of OPTiOPT_iOPTi​. The proof must handle the clipped gains, including the case when both are zero, while preserving the exact joint law of (Xi,Yi)(X_i,Y_i)(Xi​,Yi​). Buchbinder et al., proof of Lemma III.1

Formalization scope

The ground set is a finite Lean type; subsets are Finset X, and values are real numbers. An order is a list with no repeated elements that covers the type, including the empty type. The run is an explicit finite mass function on pairs of subsets after every prefix of the list. Expectation is a finite weighted sum, so it has no integrability exception. The transition clips the two real marginal gains and handles 0/00/00/0 by assigning probability one to the add branch, exactly as Algorithm 2 specifies. The optimum is the published maximum over all subsets. No ratio divides by a possibly zero optimum.

The theorem fixes Algorithm 2 itself; an arbitrary process with nested sets or a process defined by its desired approximation property does not satisfy this scope. The Lean goal states the value bound and leaves the paper's linear-time claim outside the formal theorem. The algorithm uses four value evaluations per processed element in its printed rule; the Lean development represents those evaluations, not an implementation cost model. The statement that its two final sets coincide is a separate milestone.

The source's main-text decreasing-returns definition has an overbroad quantifier on the added element. This development uses the equivalent lattice inequality given in the paper's footnote, which permits nonmonotone functions. The proof of Lemma II.1 also has a set-index slip, and the proof of Theorem I.2 prints FFF for fff in one display; neither slip is copied into a formal statement. The one-step inequality (3) is stated for any nested pair with the processed element in Y∖XY\setminus XY∖X, a generalization of the conditioned reachable states in the paper. Contributions proving the endpoint invariant, conditional inequality, one-step expected estimate, and final bound are all within scope.

Selected references

  • Niv Buchbinder, Moran Feldman, Joseph Naor and Roy Schwartz, A Tight Linear Time (1/2)-Approximation for Unconstrained Submodular Maximization, Proceedings of the 53rd IEEE Symposium on Foundations of Computer Science, 2012. FOCS version used here.
  • Niv Buchbinder, Moran Feldman, Joseph Naor and Roy Schwartz, A Tight Linear Time (1/2)-Approximation for Unconstrained Submodular Maximization, SIAM Journal on Computing 44(5), 2015. DOI: 10.1137/130929205. The cited statement indices above refer to the FOCS version.
11 thms3 active usersReviewed
🏆Completed
CombinatoricsGraph TheoryOperations Research+1·Captain: mikedeng1

Secretary Problems: Weights and Discounts 5: A 3e-Competitive Algorithm for the Graphic Matroid Secretary ProblemResearch Paper

Motivation

In the secretary problem, nnn items with nonnegative values arrive one at a time in a uniformly random order, and an online algorithm must decide on each arrival, irrevocably, whether to keep it. The classical version keeps one item; the rule that observes a 1/e1/e1/e fraction of the arrivals and then takes the first item better than everything seen picks the best item with probability at least 1/e1/e1/e (Ferguson 1989).

Babaioff, Immorlica and Kleinberg (SODA 2007; journal version J. ACM 2018) introduced the matroid secretary problem: the kept set must be independent in a known matroid. It models online auctions in which the feasible sets of winners have matroid structure, for example hiring along the edges of a network without closing a cycle. They gave a 161616-competitive algorithm when the matroid is graphic, i.e. the items are the edges of a graph and a set is feasible when it contains no cycle.

Timeline for graphic matroids:

  • 2007, Babaioff–Immorlica–Kleinberg: 161616-competitive.
  • 2009, Babaioff–Dinitz–Gupta–Immorlica–Talwar (SODA 2009, Theorem 1.5): 3e≈8.153e\approx 8.153e≈8.15-competitive, through a random reduction to partition matroids. This mission formalizes that result.
  • 2009, Korula–Pál (ICALP 2009): 2e2e2e-competitive, by a different reduction.

Setting

Let G=(V,E)G=(V,E)G=(V,E) be a finite simple graph. Each edge eee has a value v(e)≥0v(e)\ge 0v(e)≥0. A set S⊆ES\subseteq ES⊆E is independent in the graphic matroid of GGG if the graph (V,S)(V,S)(V,S) has no cycle. The offline optimum is

OPT(G,v)=max⁡{∑e∈Sv(e):S⊆E acyclic}.\mathrm{OPT}(G,v)=\max\Big\{\sum_{e\in S}v(e): S\subseteq E\ \text{acyclic}\Big\}.OPT(G,v)=max{e∈S∑​v(e):S⊆E acyclic}.

The edges arrive in a uniformly random order. An algorithm sees each edge and its value on arrival and decides at once whether to select it. The selected set must be acyclic. The algorithm is α\alphaα-competitive if OPT(G,v)≤α⋅E[value of the selected set]\mathrm{OPT}(G,v)\le\alpha\cdot\mathbb E[\text{value of the selected set}]OPT(G,v)≤α⋅E[value of the selected set] for every GGG and every v≥0v\ge 0v≥0.

A partition matroid on a subset U′⊆EU'\subseteq EU′⊆E is given by a family PPP of nonempty, pairwise disjoint parts with union U′U'U′: a set is independent when it lies in U′U'U′ and meets each part at most once. Its max-weight base has value val(P,v)=∑p∈Pmax⁡e∈pv(e)\mathrm{val}(P,v)=\sum_{p\in P}\max_{e\in p}v(e)val(P,v)=∑p∈P​maxe∈p​v(e).

Definition 5.1. A random partition μ\muμ (a probability distribution on such families, chosen from GGG alone) is an α\alphaα-partition scheme if every partition in its support has only acyclic independent sets, and for every v≥0v\ge 0v≥0,

OPT(G,v)≤α⋅EP∼μ[val(P,v)].\mathrm{OPT}(G,v)\le \alpha\cdot\mathbb E_{P\sim\mu}[\mathrm{val}(P,v)].OPT(G,v)≤α⋅EP∼μ​[val(P,v)].

The random partition of Lemma 5.3. Pick an edge {u,w}\{u,w\}{u,w} uniformly at random. With probability 12\tfrac1221​ colour uuu red and www blue, otherwise the reverse. Colour every other vertex red or blue independently with probability 12\tfrac1221​. Each red vertex xxx gets a part: the red-blue edges at xxx. Then repeat on the edges with both endpoints blue, with fresh randomness.

The algorithm. Draw the partition, let the edges arrive, and on each part run the classical secretary rule on that part's arrivals. Output all selected edges.

Formalization targets

Goal: Theorem 1.5

For every finite simple graph GGG and every v≥0v\ge 0v≥0:

  1. every possible output of the algorithm is an acyclic set of edges of GGG;
OPT(G,v)≤3e⋅E[ALG].\mathrm{OPT}(G,v)\le 3e\cdot\mathbb E[\mathrm{ALG}].OPT(G,v)≤3e⋅E[ALG].

Part 1 is needed for the statement to have content: an algorithm that selects every edge would otherwise satisfy part 2.

Milestones

  • Section 2, p. 4. On m≥1m\ge1m≥1 arrivals, the classical rule selects the maximum with probability at least 1/e1/e1/e.
  • Theorem 5.4, first clause. For a fixed partition PPP, the per-part rule outputs a set independent in the partition matroid, and val(P,v)≤e⋅Eπ[ALG]\mathrm{val}(P,v)\le e\cdot\mathbb E_\pi[\mathrm{ALG}]val(P,v)≤e⋅Eπ​[ALG].
  • Lemma 5.3, independence. Every partition the random construction can produce is a partition matroid on a subset of EEE, and each of its independent sets is a forest.
  • Lemma 5.3. The construction is a 333-partition scheme.
  • Section 5, p. 10. Any α\alphaα-partition scheme for a graphic matroid, combined with the per-part rule, gives a feasible, eαe\alphaeα-competitive algorithm.

Significance

The theorem shows that the graphic matroid secretary problem admits a constant-competitive algorithm with a small explicit constant. It does so through a reduction: a random partition matroid that is feasible for the original matroid and loses only a constant factor in expectation. The reduction separates the combinatorics (Lemma 5.3) from the online part (Theorem 5.4). The same framework gives algorithms for uniform and transversal matroids and for the weighted and discounted variants on any matroid with an α\alphaα-partition property.

The result is proved in the paper; it has not been formalized. The mission contributes a machine-checked version of the reduction, a formal treatment of a recursively defined random partition, and the classical secretary bound in a reusable finite form. The constant 3e3e3e is not the best known for graphic matroids (Korula–Pál improve it to 2e2e2e), so the formal goal is this algorithm's guarantee, not the best possible ratio.

Difficulty

The online half is routine once the classical bound is available: the relative order of the edges in each part is uniform, and the parts are disjoint. The difficulty is Lemma 5.3. The natural idea of using a fixed optimal forest to build the partition is ruled out because the partition must be chosen before the values are seen. The expectation bound must therefore hold for every valuation at once, for a law that depends on the graph only. The construction is recursive and random: its expected value is not a closed-form sum, and any bound has to be carried through the random sequence of blue-blue subgraphs. Feasibility needs an invariant across rounds: the parts created later live inside the blue-blue edges of every earlier round.

Formalization scope

  • Graph. A SimpleGraph on a Fintype vertex type with decidable adjacency. The edges are G.edgeFinset, and acyclicity of SSS is (SimpleGraph.fromEdgeSet S).IsAcyclic. Multigraphs are not covered.
  • Values. Values are a real function v : Sym2 V → ℝ with ∀ e, 0 ≤ v e; only the values on edges matter.
  • OPT is a Finset.sup' over acyclic subsets of the edge set. A partition is a finite family of nonempty, pairwise disjoint parts inside the edge set. Its max-weight base value is the sum of the part maxima.
  • Random partition. A PMF defined by well-founded recursion on the number of edges. Empty parts are dropped, and edges with two red endpoints are discarded.
  • Random order. The edges are numbered by a fixed enumeration. An arrival order is a permutation of the numbers, and expectation over the order is the average over all ∣E∣!|E|!∣E∣! permutations.
  • Classical rule. It samples ⌊m/e⌋\lfloor m/e\rfloor⌊m/e⌋ arrivals of a part with mmm edges. Ties are broken by preferring the smaller edge number among equal values.
  • Constants. Competitiveness is multiplicative (OPT≤3e⋅E[ALG]\mathrm{OPT}\le 3e\cdot\mathbb E[\mathrm{ALG}]OPT≤3e⋅E[ALG]), so a zero expectation is not a loophole.
  • Ruling out trivial formalizations. In Definition 5.1 the random partition is fixed before the valuation, and the independence requirement holds for every partition in its support. A partition allowed to depend on vvv would make every matroid 111-partitionable.

A complete development needs the classical secretary bound in finite form, the uniformity of induced sub-orders of a uniform permutation, expectations of PMF.bind along a well-founded recursion, and facts about forests in SimpleGraph. The first two, and a general graphic-matroid layer, are reusable beyond this mission. Proofs of any milestone, alternative proofs of Lemma 5.3, and extensions to the uniform and transversal cases of Theorem 5.2 are welcome.

Selected references

  • M. Babaioff, M. Dinitz, A. Gupta, N. Immorlica, K. Talwar, Secretary Problems: Weights and Discounts, Proc. 20th ACM-SIAM Symposium on Discrete Algorithms (SODA), 2009. https://doi.org/10.1137/1.9781611973068.135
  • M. Babaioff, N. Immorlica, R. Kleinberg, Matroids, secretary problems, and online mechanisms, SODA 2007, pp. 434–443. https://dl.acm.org/doi/10.5555/1283383.1283429
  • M. Babaioff, N. Immorlica, D. Kempe, R. Kleinberg, Matroid Secretary Problems, Journal of the ACM 65(6), 2018. https://doi.org/10.1145/3212512
  • N. Korula, M. Pál, Algorithms for Secretary Problems on Graphs and Hypergraphs, ICALP 2009, LNCS 5556. https://doi.org/10.1007/978-3-642-02930-1_42
  • T. S. Ferguson, Who solved the secretary problem?, Statistical Science 4(3), 1989. https://doi.org/10.1214/ss/1177012493
10 thms3 active usersReviewed
🏆Completed
Operations ResearchOptimization·Captain: mikedeng1

A Simple Forward Algorithm to Solve General Dynamic Lot Sizing Models with n Periods in O(n log n) or O(n) Time: Minimal Optimal Predecessor Lists Are Characterized by Strictly Increasing BreakpointsResearch Paper

Motivation

The dynamic lot size model asks when, and how much, to order of a single item over a planning horizon of nnn periods with known, time-varying demands, setup costs, unit order costs and holding costs. It is the textbook model of production planning and the building block of material requirements planning, multi-item scheduling and many decomposition schemes for larger supply-chain problems.

Wagner and Whitin (1958) showed that some optimal policy orders only when inventory is zero, which turns the problem into a shortest-path recursion with O(n2)O(n^2)O(n2) running time. For more than thirty years this was the standard algorithm. In 1991 three groups independently reduced the complexity: Federgruen and Tzur (Management Science 37(8), 1991), Wagelmans, van Hoesel and Kolen (Operations Research 40, 1992) and Aggarwal and Park (Operations Research 41, 1993). Each obtained O(nlog⁡n)O(n \log n)O(nlogn) in general and O(n)O(n)O(n) under special cost structures. The Federgruen–Tzur algorithm is a forward algorithm: at iteration jjj it keeps a short list of periods that could still be the best last setup period for some future horizon, and updates it by local tests on neighbouring entries. This mission formalizes the theorem that justifies those tests.

Setting

For periods i=1,2,…i = 1, 2, \dotsi=1,2,… let did_idi​ be the demand, KiK_iKi​ the setup cost, cic_ici​ the variable per unit order cost and hih_ihi​ the cost of carrying a unit of inventory at the end of period iii. Write D(i)=∑k=1idkD(i) = \sum_{k=1}^{i} d_kD(i)=∑k=1i​dk​ and H(i)=∑k=1ihkH(i) = \sum_{k=1}^{i} h_kH(i)=∑k=1i​hk​, so D(0)=H(0)=0D(0) = H(0) = 0D(0)=H(0)=0. For i<ji < ji<j let cij=ci+hi+⋯+hj−1c_{ij} = c_i + h_i + \dots + h_{j-1}cij​=ci​+hi​+⋯+hj−1​, let C~(i)=ci−H(i−1)\tilde C(i) = c_i - H(i-1)C~(i)=ci​−H(i−1), and let

S(i,j)=∑r=ij−1hr (D(j)−D(r))S(i, j) = \sum_{r=i}^{j-1} h_r\,\bigl(D(j) - D(r)\bigr)S(i,j)=r=i∑j−1​hr​(D(j)−D(r))

be the carrying cost of an order placed in period iii that covers the demands of periods i,…,ji, \dots, ji,…,j.

The costs are given by the zero-inventory recursion (2): F(0)=0F(0) = 0F(0)=0 and, for 1≤l≤t1 \le l \le t1≤l≤t,

F(l,t)=F(l−1)+Kl+S(l,t)+cl [D(t)−D(l−1)],F(t)=min⁡1≤l≤tF(l,t).F(l, t) = F(l-1) + K_l + S(l, t) + c_l\,[D(t) - D(l-1)], \qquad F(t) = \min_{1 \le l \le t} F(l, t).F(l,t)=F(l−1)+Kl​+S(l,t)+cl​[D(t)−D(l−1)],F(t)=1≤l≤tmin​F(l,t).

F(l,t)F(l, t)F(l,t) is the cost of the first ttt periods when the last setup is in period lll.

For two periods k<lk < lk<l the difference Δk,l(t)=F(k,t)−F(l,t)\Delta_{k,l}(t) = F(k,t) - F(l,t)Δk,l​(t)=F(k,t)−F(l,t) is affine in D(t)D(t)D(t), with intercept A(k,l)A(k,l)A(k,l) given by (4) and slope ck,l−cl=C~(k)−C~(l)c_{k,l} - c_l = \tilde C(k) - \tilde C(l)ck,l​−cl​=C~(k)−C~(l). Its root G(k,l)G(k,l)G(k,l) is defined by (5): A(k,l)/(C~(l)−C~(k))A(k,l)/(\tilde C(l) - \tilde C(k))A(k,l)/(C~(l)−C~(k)) when the slopes differ, and +∞+\infty+∞ or −∞-\infty−∞ according to the sign of A(k,l)A(k,l)A(k,l) when they agree. It is extended symmetrically, G(l,k)=G(k,l)G(l,k) = G(k,l)G(l,k)=G(k,l).

At iteration jjj the future demands are unknown, so a future horizon has a potential cumulative demand x≥D(j)x \ge D(j)x≥D(j). The jjjth Minimal Optimal Predecessors list Ω(j)\Omega(j)Ω(j) is the set of periods l≤jl \le jl≤j that are the lowest-index optimal last setup period, among {1,…,j}\{1, \dots, j\}{1,…,j}, for every potential cumulative demand in some open interval above D(j)D(j)D(j).

Formalization targets

Goal: Theorem 1(a)

Let j≥1j \ge 1j≥1 and let S={i1,…,ir}S = \{i_1, \dots, i_r\}S={i1​,…,ir​} with Ω(j)⊆S⊆{1,…,j}\Omega(j) \subseteq S \subseteq \{1, \dots, j\}Ω(j)⊆S⊆{1,…,j}, ranked so that C~(i1)≥⋯≥C~(ir)\tilde C(i_1) \ge \dots \ge \tilde C(i_r)C~(i1​)≥⋯≥C~(ir​), with equal C~\tilde CC~-values in ascending order of index. Put g(1)=D(j)g(1) = D(j)g(1)=D(j) and g(l)=G(il,il−1)g(l) = G(i_l, i_{l-1})g(l)=G(il​,il−1​) for l=2,…,rl = 2, \dots, rl=2,…,r. Then

S=Ω(j)  ⟺  g(1)<g(2)<⋯<g(r)<∞.(6)S = \Omega(j) \iff g(1) < g(2) < \dots < g(r) < \infty. \tag{6}S=Ω(j)⟺g(1)<g(2)<⋯<g(r)<∞.(6)

Milestones

In attack order:

  • identity (1a) for the carrying costs;
  • Lemma 2(a)–(d), the linearity of Δk,l\Delta_{k,l}Δk,l​ and the sign test against its root G(k,l)G(k,l)G(k,l);
  • the claim that Ω(j)\Omega(j)Ω(j) contains an optimal last setup period for the horizon jjj;
  • the strict chains (7)–(8) of the Appendix;
  • Theorem 1(b), that under (6) the first entry i1i_1i1​ is an optimal last setup period l(j)l(j)l(j);
  • Theorem 1(c)(i)–(iii), the three elimination rules: g(2)≤D(j)g(2) \le D(j)g(2)≤D(j) removes i1i_1i1​, g(k+1)≤g(k)g(k+1) \le g(k)g(k+1)≤g(k) removes iki_kik​, and g(r)=∞g(r) = \inftyg(r)=∞ removes iri_rir​.

A supporting item potCost_spec certifies that the potential costs used to define Ω(j)\Omega(j)Ω(j) agree with the paper's F(l,t)F(l,t)F(l,t), up to a term that does not depend on lll.

Significance

Theorem 1 is what makes the forward algorithm correct. Part (a) reduces the minimality of a candidate list to a condition on consecutive pairs of a sorted list. Part (c) says which entry to delete when the condition fails. Part (b) says where to read off the optimal last setup period. With these, Ω(j)\Omega(j)Ω(j) is maintained by deletions at the ends and in the interior of a list ordered by C~\tilde CC~, and each period is inserted and deleted at most once; the O(nlog⁡n)O(n \log n)O(nlogn) bound, and the O(n)O(n)O(n) bound under the paper's special cost structures, follow from this bookkeeping. The same lower-envelope reasoning appears in the other 1991–1993 algorithms and in later extensions to backlogging and capacitated variants.

The result has a complete published proof. To our knowledge there is no machine-checked development of the Wagner–Whitin recursion or of any of the fast lot-sizing algorithms. This mission produces the model, the breakpoints and the Minimal Optimal Predecessors lists as reusable definitions, and a checked proof of the characterization. It also records two small corrections that a formal reading forces on the printed text (see Formalization scope).

Difficulty

Each piece in isolation is elementary algebra on affine functions. The difficulty is in the combinatorics of the lower envelope with ties. The natural argument "consecutive breakpoints increase, so each line owns an interval" must handle three things:

  • equal slopes, where G=±∞G = \pm\inftyG=±∞;
  • several lines meeting at one point;
  • the lowest-index tie-breaking that makes Ω(j)\Omega(j)Ω(j) minimal.

The "only if" direction needs every failure of (6) to be traced to an element that is never the unique lowest-index optimum on an interval. Ties are exactly where the printed definition of Ω(j)\Omega(j)Ω(j), read literally at a single demand value, breaks the theorem. A proof that ignores ties proves a statement that is false.

Formalization scope

  • Data and costs. The data are four functions N→R\mathbb N \to \mathbb RN→R bundled in a structure; values at index 000 are unused, and no sign conditions are imposed. FFF is defined by the recursion (2) with F(0)=0F(0) = 0F(0)=0. Its identification with the minimum cost over all feasible policies is the paper's Lemma 1 (Wagner–Whitin), which is not part of this mission. The horizon nnn is not a parameter.
  • Breakpoints. GGG and the critical values g(⋅)g(\cdot)g(⋅) take values in EReal, so ±∞\pm\infty±∞ are kept distinct from every real number. The final "<∞< \infty<∞" of (6) is part of the condition.
  • Ranked lists. A ranked set is a duplicate-free List ℕ. Lean lists are 0-based, so the paper's im+1i_{m+1}im+1​ and g(m+1)g(m+1)g(m+1) are entry mmm and gval j L m.
  • Disclosed change 1, Ω(j)\Omega(j)Ω(j). The page asks for a single potential cumulative demand D≥D(j)D \ge D(j)D≥D(j) at which lll is the lowest-index optimum. With that reading, Theorem 1(a) "only if" and Theorem 1(c) fail when two lines tie exactly at a breakpoint (an explicit five-period instance is in the definition's note). The formalization requires lll to be the lowest-index optimum on a nondegenerate open interval of potential demands above D(j)D(j)D(j). This is the paper's own description of the list on p. 915: "the unique optimal last setup period for any horizon … with potential cumulative demand g(k)<D<g(k+1)g(k) < D < g(k+1)g(k)<D<g(k+1)".
  • Disclosed change 2, Lemma 2(d). The printed hypothesis "ck,l<clc_{k,l} < c_lck,l​<cl​" duplicates part (c) and is read as "ck,l=clc_{k,l} = c_lck,l​=cl​". The equivalence "Δk,l≥0\Delta_{k,l} \ge 0Δk,l​≥0 iff D(t)≥G(k,l)D(t) \ge G(k,l)D(t)≥G(k,l)" is stated under A(k,l)≠0A(k,l) \ne 0A(k,l)=0, since A(k,l)=0A(k,l) = 0A(k,l)=0 gives G=+∞G = +\inftyG=+∞ by (5).
  • Ruling out trivial formalizations. The hypotheses of the goal are satisfiable for every j≥1j \ge 1j≥1: rank {1,…,j}\{1, \dots, j\}{1,…,j} itself. Ω(j)\Omega(j)Ω(j) is nonempty (a milestone). F(t)F(t)F(t) for t≥1t \ge 1t≥1 is a minimum over the nonempty set {1,…,t}\{1, \dots, t\}{1,…,t}, never a default value. GGG is never replaced by a real-valued junk value at equal slopes.
  • Out of scope. Lemma 1, Lemma 3, Corollaries 1–5, Theorem 2, the Algorithm's pseudo-code and its complexity analysis, and the submodularity discussion of §5.
  • Reusable infrastructure. The model, the recursion (2), AAA, GGG and Ω(j)\Omega(j)Ω(j) can be reused for the paper's algorithmic results and for related lot-sizing papers. Proofs of the milestones, in any order, are welcome.

Selected references

  • A. Federgruen and M. Tzur, A Simple Forward Algorithm to Solve General Dynamic Lot Sizing Models with n Periods in O(n log n) or O(n) Time, Management Science 37(8):909–925, 1991. https://doi.org/10.1287/mnsc.37.8.909
  • H. M. Wagner and T. M. Whitin, Dynamic Version of the Economic Lot Size Model, Management Science 5(1):89–96, 1958. https://doi.org/10.1287/mnsc.5.1.89
  • A. Wagelmans, S. van Hoesel and A. Kolen, Economic Lot-Sizing: An O(n log n) Algorithm That Runs in Linear Time in the Wagner-Whitin Case, Operations Research 40(1-supplement-1):S145–S156, 1992. https://doi.org/10.1287/opre.40.1.S145
  • A. Aggarwal and J. K. Park, Improved Algorithms for Economic Lot Size Problems, Operations Research 41(3):549–571, 1993. https://doi.org/10.1287/opre.41.3.549
16 thms3 active usersReviewed
🏆Completed
Operations Research·Captain: mikedeng1

On the Power of Randomization in On-Line Algorithms 4: Restarting a Bounded-Cost Algorithm Is (1+ε)α-Competitive in Games of Finite DiameterResearch Paper

Motivation

Competitive analysis compares an online algorithm, which must answer each request before seeing the next, with the optimal off-line solution of the same request sequence. Ben-David, Borodin, Karp, Tardos and Wigderson (Algorithmica 11, 1994) set up a general framework, request-answer games, in which paging, the KKK-server problem and metrical task systems are all instances, and used it to compare the power of randomized algorithms against several kinds of adversaries.

One of their results (Theorem 2.1) says that if a randomized algorithm is α\alphaα-competitive against every adaptive off-line adversary, then some deterministic algorithm is already α\alphaα-competitive. The argument is a game-theoretic existence proof: it says nothing about how to compute the deterministic algorithm. Section 4 of the paper, "A Constructive Version of Theorem 2.1", answers the natural follow-up question for a large class of games: under monotonicity, locality and a finite diameter, a deterministic algorithm that loses only a factor 1+ϵ1+\epsilon1+ϵ can be assembled from a finite object.

Setting

A request-answer game FFF has a request set RRR, a finite answer set AAA, and cost functions fn:Rn×An→Rf_n : R^n \times A^n \to \mathbb{R}fn​:Rn×An→R, n≥0n \ge 0n≥0; f0f_0f0​ is the cost of the empty play. The off-line optimum of a request sequence rrr of length nnn is c(r)=min⁡a∈Anfn(r,a)c(r) = \min_{a \in A^n} f_n(r, a)c(r)=mina∈An​fn​(r,a). A deterministic online algorithm GGG answers the iii-th request by a function gi(r1,…,ri)g_i(r_1, \dots, r_i)gi​(r1​,…,ri​) of the requests so far; its cost on rrr is cG(r)=fn(r,G(r))c_G(r) = f_n(r, G(r))cG​(r)=fn​(r,G(r)). It is α\alphaα-competitive if cG(r)≤α(c(r))c_G(r) \le \alpha(c(r))cG​(r)≤α(c(r)) for every rrr.

The game is monotone if fn+1(rt,ab)≥fn(r,a)f_{n+1}(rt, ab) \ge f_n(r, a)fn+1​(rt,ab)≥fn​(r,a) always, and local if for every h>0h > 0h>0 only finitely many request sequences have c(r)≤hc(r) \le hc(r)≤h. The discrepancy of two request-answer sequences is

δ((r,a),(r′,a′))=f(rr′,aa′)−f(r,a)−f(r′,a′),\delta((r,a),(r',a')) = f(rr', aa') - f(r,a) - f(r',a'),δ((r,a),(r′,a′))=f(rr′,aa′)−f(r,a)−f(r′,a′),

and the diameter D(F)D(F)D(F) is the supremum of ∣δ∣|\delta|∣δ∣ over all pairs.

For a real HHH, the set RHR_HRH​ consists of the request sequences all of whose proper prefixes have off-line optimum at most HHH. Given a deterministic algorithm AHA_HAH​, the restart algorithm simulates AHA_HAH​ and, as soon as the request sequence leaves RHR_HRH​, starts over as if it had received no previous requests. This cuts every request sequence into segments r=r(1) r(2)⋯r(t)r = r(1)\, r(2) \cdots r(t)r=r(1)r(2)⋯r(t), each a longest prefix in RHR_HRH​ of what remains.

Formalization targets

Goal: Theorem 4.1 through its construction

Let FFF be monotone and local, with finite nonempty RRR and AAA, f0≥0f_0 \ge 0f0​≥0, and diameter at most DDD. Let α(x)=d x\alpha(x) = d\,xα(x)=dx with d≥1d \ge 1d≥1, ϵ>0\epsilon > 0ϵ>0, and

H=(2+ϵ)Dϵ.H = \frac{(2+\epsilon) D}{\epsilon}.H=ϵ(2+ϵ)D​.

Then RHR_HRH​ is finite; the restart algorithm depends on AHA_HAH​ only through its values on RHR_HRH​; and for every AHA_HAH​ with cAH(r)≤α(c(r))c_{A_H}(r) \le \alpha(c(r))cAH​​(r)≤α(c(r)) on RHR_HRH​,

cRestart(r)≤(1+ϵ) α(c(r))for all request sequences r.c_{\mathrm{Restart}}(r) \le (1+\epsilon)\,\alpha(c(r)) \qquad \text{for all request sequences } r.cRestart​(r)≤(1+ϵ)α(c(r))for all request sequences r.

Milestones (proof of Theorem 4.1, pp. 17–18)

  1. RHR_HRH​ is finite.
  2. The restart rule produces the greedy decomposition into longest prefixes in RHR_HRH​.
  3. The restart algorithm answers AH(r(1)),AH(r(2)),…,AH(r(t))A_H(r(1)), A_H(r(2)), \dots, A_H(r(t))AH​(r(1)),AH​(r(2)),…,AH​(r(t)).
  4. c(r(i))≥Hc(r(i)) \ge Hc(r(i))≥H for i=1,…,t−1i = 1, \dots, t-1i=1,…,t−1.
  5. c(r)≥c(r(1))+∑i=2t(c(r(i))−D(F))c(r) \ge c(r(1)) + \sum_{i=2}^{t} (c(r(i)) - D(F))c(r)≥c(r(1))+∑i=2t​(c(r(i))−D(F)).
  6. cRestart(r)≤α(c(1))+∑i=2t(α(c(i))+D(F))c_{\mathrm{Restart}}(r) \le \alpha(c(1)) + \sum_{i=2}^{t} (\alpha(c(i)) + D(F))cRestart​(r)≤α(c(1))+∑i=2t​(α(c(i))+D(F)).

Significance

Theorem 2.1 shows that, against adaptive off-line adversaries, randomization gives no advantage, but only as an existence statement. Theorem 4.1 turns it into a recipe: a deterministic algorithm need only be good on the finite set RHR_HRH​, which can be prepared in advance, and restarting extends it to all inputs at a loss of 1+ϵ1+\epsilon1+ϵ. The paper illustrates this with KKK-server problems on finite graphs, where AHA_HAH​ is a finite table and each step of the resulting algorithm costs one dynamic-programming evaluation of an off-line optimum.

The restart construction is of independent use: it is a general way to turn a guarantee on bounded-cost inputs into a guarantee on all inputs when the cost is nearly additive over concatenation.

To our knowledge neither the abstract request-answer game model nor this theorem has a machine-checked proof. The mission produces a formal model of request-answer games with the monotonicity, locality and diameter conditions, the restart algorithm as an explicit definition, and the full chain of inequalities of the proof.

Difficulty

The individual inequalities are elementary; the work lies in the bookkeeping of the construction. The algorithm is defined online, one request at a time, while the analysis is phrased through the decomposition of the whole sequence. Relating the two requires showing that the online rule produces exactly the decomposition into longest prefixes in RHR_HRH​, and that the answers produced on each segment are those of AHA_HAH​ run from scratch, so that the cost of the whole play can be compared with the costs on the segments. The constants must close exactly: with H=(2+ϵ)D/ϵH = (2+\epsilon)D/\epsilonH=(2+ϵ)D/ϵ the additive losses at the t−1t-1t−1 cuts, on both the algorithm's side and the optimum's side, must be absorbed by the factor 1+ϵ1+\epsilon1+ϵ, and they do so only for ratios d≥1d \ge 1d≥1.

A tempting shortcut is to conclude Theorem 4.1 from Theorem 2.1 directly: a deterministic α\alphaα-competitive algorithm is trivially (1+ϵ)α(1+\epsilon)\alpha(1+ϵ)α-competitive when costs are nonnegative. This proves the sentence of the theorem without its point, and is ruled out below.

Formalization scope

  • Request and answer sequences are Lean Lists, oldest first; fn(r,a)f_n(r,a)fn​(r,a) is F.cost r a with r.length = a.length. Costs are real-valued (the paper allows +∞+\infty+∞; real costs are a special case).
  • The answer set is a Fintype and nonempty; in the goal the request set is a Fintype and nonempty.
  • A deterministic algorithm is one function List R → A. Competitiveness of a deterministic algorithm against adaptive off-line adversaries is stated as competitiveness on every request sequence, which is equivalent for deterministic algorithms (p. 8).
  • Finite diameter is a real bound DDD with ∣δ∣≤D|\delta| \le D∣δ∣≤D for all pairs, and every statement holds for every such DDD, in particular D=D(F)D = D(F)D=D(F); no real supremum is taken.
  • f0≥0f_0 \ge 0f0​≥0 is assumed in the goal; with monotonicity it makes every cost nonnegative. It holds in the paper's examples.
  • α(x)=d x\alpha(x) = d\,xα(x)=dx with d≥1d \ge 1d≥1 in the goal; the cost bound for the restart algorithm (milestone 6) holds for an arbitrary α\alphaα.
  • HHH is fixed to the paper's value (2+ϵ)D/ϵ(2+\epsilon)D/\epsilon(2+ϵ)D/ϵ.
  • The restart algorithm is a definition (a left fold over the requests), not a hypothesis. The paper's assumption of a randomized algorithm competitive against adaptive off-line adversaries is used only, through Theorem 2.1, to obtain AHA_HAH​; the goal quantifies over every AHA_HAH​ that is α\alphaα-competitive on RHR_HRH​. "Computable" has no precise meaning for real costs; its content is the finiteness of RHR_HRH​ together with the fact that the restart algorithm reads AHA_HAH​ only on RHR_HRH​. A formalization that proves only the existence of a deterministic (1+ϵ)α(1+\epsilon)\alpha(1+ϵ)α-competitive algorithm, without the construction, does not meet the goal.
  • Milestones 5 and 6 are stated for nonempty request sequences, so that t≥1t \ge 1t≥1; milestone 5 is stated for every decomposition into consecutive pieces.

Contributions welcome: proofs of the milestones, and a connection to the other missions of this series, where the randomized model and Theorem 2.1 are formalized.

Selected references

  • S. Ben-David, A. Borodin, R. Karp, G. Tardos, A. Wigderson, On the power of randomization in on-line algorithms, Algorithmica 11 (1994). https://doi.org/10.1007/BF01294260
  • M. Chrobak, H. Karloff, T. Payne, S. Vishwanathan, New results on server problems, SIAM Journal on Discrete Mathematics 4 (1991), 172–181. https://doi.org/10.1137/0404017
10 thms3 active usersReviewed
🏆Completed
Operations ResearchProbability·Captain: mikedeng1

On the Power of Randomization in On-Line Algorithms 3: An Augmented Potential Function Yields an Explicit Deterministic α∘β-Competitive AlgorithmResearch Paper

Motivation

Competitive analysis measures an online algorithm, which must answer each request before seeing the next, against the best off-line answer to the whole request sequence. Randomized online algorithms are often much better than deterministic ones against an oblivious adversary, who fixes the requests in advance; the paging problem is the standard example. Against an adaptive adversary, who sees the algorithm's answers before choosing the next request, the advantage can disappear.

Ben-David, Borodin, Karp, Tardos and Wigderson (Algorithmica 11, 1994; preliminary version STOC 1990) made this precise in an abstract framework of request-answer games. Their Corollary 2.1 says: if a game has a randomized algorithm that is α\alphaα-competitive against adaptive on-line adversaries and one that is β\betaβ-competitive against oblivious adversaries, then it has a deterministic α∘β\alpha\circ\betaα∘β-competitive algorithm. That proof is non-constructive: it goes through a game-theoretic determinacy argument. Section 3 of the paper gives a constructive version. Most competitive analyses of randomized algorithms against adaptive adversaries are carried out with a potential function, in the style of Manasse, McGeoch and Sleator (J. Algorithms 11, 1990). The paper shows that such a potential function, together with any oblivious-competitive algorithm HHH, determines an explicit deterministic algorithm MMM, answer by answer. This mission formalizes that construction and its guarantee.

Setting

A request-answer game has a request set RRR, a finite answer set AAA and cost functions fn:Rn×An→Rf_n : R^n \times A^n \to \mathbb Rfn​:Rn×An→R. The off-line optimum of r∈Rnr \in R^nr∈Rn is c(r)=min⁡a∈Anfn(r,a)c(r) = \min_{a \in A^n} f_n(r, a)c(r)=mina∈An​fn​(r,a). A deterministic online algorithm MMM is a sequence of maps mi:Ri→Am_i : R^i \to Ami​:Ri→A; on r=(r1,…,rn)r = (r_1,\dots,r_n)r=(r1​,…,rn​) it answers M(r)=(m1(r1),m2(r1,r2),…,mn(r))M(r) = (m_1(r_1), m_2(r_1,r_2), \dots, m_n(r))M(r)=(m1​(r1​),m2​(r1​,r2​),…,mn​(r)), at cost cM(r)=fn(r,M(r))c_M(r) = f_n(r, M(r))cM​(r)=fn​(r,M(r)). It is α\alphaα-competitive if cM(r)≤α(c(r))c_M(r) \le \alpha(c(r))cM​(r)≤α(c(r)) for all rrr. Throughout, α\alphaα and β\betaβ are affine maps R→R\mathbb R \to \mathbb RR→R (the paper's "linear functions").

A randomized online algorithm HHH is a probability distribution over deterministic algorithms HyH_yHy​; it is β\betaβ-competitive against any oblivious adversary if Ey[fn(r,Hy(r))]≤β(c(r))\mathbb E_y[f_n(r, H_y(r))] \le \beta(c(r))Ey​[fn​(r,Hy​(r))]≤β(c(r)) for all rrr. The algorithm GGG analysed by the potential function is described by its next-answer laws gn+1(rrn+1,a)g_{n+1}(r r_{n+1}, a)gn+1​(rrn+1​,a) on AAA, given the requests so far, the new request and its own past answers. An adaptive on-line adversary SSS chooses each request from the algorithm's past answers and answers it itself, before the algorithm does, for at most dQd_QdQ​ rounds; a configuration after nnn rounds is (r,a,b)∈Rn×An×An(r, a, b) \in R^n \times A^n \times A^n(r,a,b)∈Rn×An×An: requests, algorithm's answers, adversary's answers.

An augmented potential function for α\alphaα and GGG (Definition 3.1) is a family Φn:Rn×An×An→R\Phi_n : R^n \times A^n \times A^n \to \mathbb RΦn​:Rn×An×An→R with (1) Φ0=0\Phi_0 = 0Φ0​=0; (2) Φn(r,a,b)≤α(fn(r,b))−fn(r,a)\Phi_n(r,a,b) \le \alpha(f_n(r,b)) - f_n(r,a)Φn​(r,a,b)≤α(fn​(r,b))−fn​(r,a) for every configuration; (3) Ean+1∼gn+1(rrn+1,a)[Φn+1(rrn+1,aan+1,bbn+1)]≥Φn(r,a,b)\mathbb E_{a_{n+1} \sim g_{n+1}(r r_{n+1}, a)}[\Phi_{n+1}(r r_{n+1}, a a_{n+1}, b b_{n+1})] \ge \Phi_n(r,a,b)Ean+1​∼gn+1​(rrn+1​,a)​[Φn+1​(rrn+1​,aan+1​,bbn+1​)]≥Φn​(r,a,b) for every configuration, every rn+1∈Rr_{n+1} \in Rrn+1​∈R and every bn+1∈Ab_{n+1} \in Abn+1​∈A.

Formalization targets

Goal: Theorem 3.1 (p. 15)

Let Φ\PhiΦ be an augmented potential function for α\alphaα and GGG, and HHH a β\betaβ-competitive algorithm against oblivious adversaries. Say that MMM obeys the potential rule if for every r∈Rnr \in R^nr∈Rn and r′=rtr' = rtr′=rt,

Ey[Φn+1(r′,M(r) mn+1(r′),Hy(r′))] ≥ Ey[Φn(r,M(r),Hy(r))].\mathbb E_y\big[\Phi_{n+1}(r', M(r)\,m_{n+1}(r'), H_y(r'))\big] \ \ge\ \mathbb E_y\big[\Phi_n(r, M(r), H_y(r))\big].Ey​[Φn+1​(r′,M(r)mn+1​(r′),Hy​(r′))] ≥ Ey​[Φn​(r,M(r),Hy​(r))].

Then such an MMM exists, and every such MMM satisfies

cM(r)≤α(β(c(r)))for all r.c_M(r) \le \alpha\big(\beta(c(r))\big) \quad \text{for all } r .cM​(r)≤α(β(c(r)))for all r.

Both parts are part of the goal: the rule can be followed, and following it guarantees α∘β\alpha\circ\betaα∘β-competitiveness.

Milestones

  1. In every play of GGG against an adaptive on-line adversary, the expected final potential is nonnegative (proof of Lemma 3.1).
  2. Lemma 3.1, "if" direction: an augmented potential function for α\alphaα and GGG makes GGG α\alphaα-competitive against any adaptive on-line adversary, E[cG(S)]≤E[α(cS(G))]\mathbb E[c_G(S)] \le \mathbb E[\alpha(c_S(G))]E[cG​(S)]≤E[α(cS​(G))].
  3. For every rrr, ttt and every a∈Ana \in A^na∈An, some a′∈Aa' \in Aa′∈A satisfies Ey[Φn+1(rt,aa′,Hy(rt))]≥Ey[Φn(r,a,Hy(r))]\mathbb E_y[\Phi_{n+1}(rt, aa', H_y(rt))] \ge \mathbb E_y[\Phi_n(r, a, H_y(r))]Ey​[Φn+1​(rt,aa′,Hy​(rt))]≥Ey​[Φn​(r,a,Hy​(r))].
  4. If MMM obeys the rule, Ey[Φn(r,M(r),Hy(r))]≥0\mathbb E_y[\Phi_n(r, M(r), H_y(r))] \ge 0Ey​[Φn​(r,M(r),Hy​(r))]≥0 for every rrr.
  5. If MMM obeys the rule, fn(r,M(r))≤Ey[α(fn(r,Hy(r)))]f_n(r, M(r)) \le \mathbb E_y[\alpha(f_n(r, H_y(r)))]fn​(r,M(r))≤Ey​[α(fn​(r,Hy​(r)))] for every rrr: MMM is α\alphaα-competitive against the randomized adaptive adversary that serves its requests with HHH.

Significance

The theorem turns two separate analyses into one deterministic algorithm with an explicit description. The potential function certifies GGG against the strongest on-line adversary; the oblivious algorithm HHH need not be related to GGG, and the paper remarks that HHH may be GGG itself. The next answer of MMM is computable whenever the expected potential under HHH is (Corollary 3.1, stated informally in the paper), and the paper notes that for the potential functions used in the KKK-server literature this expectation is computable in time polynomial in the number of nodes and KKK. Read in this light, a potential-function proof for a randomized algorithm doubles as a deterministic algorithm.

The result is proved in the paper. As far as is known, neither this theorem nor the abstract framework of request-answer games with adaptive adversaries has a machine-checked formalization. The mission produces that framework and a checked derandomization principle that applies to every request-answer game with real costs, not to one problem.

Difficulty

The obvious argument for the existence of mn+1(r′)m_{n+1}(r')mn+1​(r′) averages property (3) of Φ\PhiΦ; the work is in seeing which configuration to apply it to. The rule compares MMM's configuration against HyH_yHy​'s answers, not against an adversary playing GGG, and the paper argues through an auxiliary on-line adversary that asks r′r'r′ and serves it with HyH_yHy​. Making this rigorous requires interchanging the expectation over HHH's coins with the finite expectation over GGG's next answer, and checking that the needed expectations are finite.

The second difficulty is the two kinds of randomness. GGG enters only through its next-answer laws, while HHH must be a single distribution over deterministic algorithms: the rule evaluates Hy(r)H_y(r)Hy​(r) and Hy(r′)H_y(r')Hy​(r′) with the same coins yyy. Replacing HHH by a behavioural description breaks the coupling between consecutive rounds.

Formalization scope

Requests and answers are Lean lists, oldest first, and fn(r,a)f_n(r,a)fn​(r,a) is F.cost r a on lists of common length; values on lists of different lengths are never used. Costs are real: the paper allows fn=+∞f_n = +\inftyfn​=+∞, so every statement here is about the real-valued games. The answer type is finite and nonempty, so the minimum c(r)c(r)c(r) exists. Affine maps are written α(x)=cx+d\alpha(x) = c x + dα(x)=cx+d. The goal additionally assumes α\alphaα nondecreasing: the last step of the paper's proof applies α\alphaα to an inequality, which needs it, and the paper's examples are positive ratios. In Lemma 3.1 and milestone 5 linearity of α\alphaα is kept as the paper's standing convention, although with α\alphaα inside the expectation the argument does not use it.

GGG is a map from (requests, own answers) to a probability mass function on AAA (behavioural form); its play against an adaptive on-line adversary is a probability mass function on final configurations, with finite support, and its expectations are finite sums. HHH is a probability measure on a coin space with a deterministic algorithm per coin, each answer measurable in the coins; expectations over HHH are Bochner integrals of functions with finitely many values. α\alphaα stays inside expectations, as in the paper's definition of competitiveness against adaptive adversaries. Adversaries stop by returning none and have a uniform depth bound.

The goal cannot be satisfied vacuously: it states the existence of an algorithm obeying the rule alongside the guarantee for every such algorithm, and Definition 3.1 is required at every configuration, not only at reachable ones.

Not formalized: the "only if" direction of Lemma 3.1, which the paper only sketches, and Corollary 3.1, whose notion of computability the paper leaves unspecified. The definitions of request-answer games, online algorithms, adversaries and competitiveness are reusable for the other missions of this paper and for any problem-specific competitive analysis. Contributions welcome: proofs of the milestones, and a lemma relating the mixed and behavioural forms of a randomized algorithm.

Selected references

  • S. Ben-David, A. Borodin, R. Karp, G. Tardos, A. Wigderson, On the power of randomization in on-line algorithms, Algorithmica 11 (1994), 2–14. https://doi.org/10.1007/BF01294260
  • M. Manasse, L. McGeoch, D. Sleator, Competitive algorithms for server problems, Journal of Algorithms 11 (1990), 208–230. https://doi.org/10.1016/0196-6774(90)90003-W
  • D. Sleator, R. Tarjan, Amortized efficiency of list update and paging rules, Communications of the ACM 28 (1985), 202–208. https://doi.org/10.1145/2786.2793
  • A. Borodin, R. El-Yaniv, Online Computation and Competitive Analysis, Cambridge University Press, 1998. ISBN 0-521-56392-5
9 thms3 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchOptimization·Captain: mikedeng1

A Tight Linear Time (1/2)-Approximation for Unconstrained Submodular Maximization 3: Fractional Double Greedy on the Multilinear Extension Achieves 1/2 of the OptimumResearch Paper

Motivation

Unconstrained submodular maximization (USM) asks for a subset SSS of a finite ground set N\mathcal NN maximizing a nonnegative submodular function fff. It contains Max-Cut, Max-DiCut and maximum facility location as special cases, and it is the basic subproblem of many constrained submodular maximization algorithms. Because fff is given only through a value oracle, the question is how close to the optimum a polynomial number of queries can get.

Timeline of the approximation ratio for USM in the value oracle model:

  • Feige, Mirrokni and Vondrák (FOCS 2007; SIAM J. Comput. 2011) showed that a uniformly random set achieves 1/41/41/4, local search achieves 1/31/31/3 and 2/52/52/5, and that no algorithm making polynomially many queries achieves 1/2+ε1/2 + \varepsilon1/2+ε for any fixed ε>0\varepsilon > 0ε>0.
  • Oveis Gharan and Vondrák (SODA 2011) reached 0.410.410.41 by simulated annealing; Feldman, Naor and Schwartz (ICALP 2011) reached 0.420.420.42.
  • Buchbinder, Feldman, Naor and Schwartz (FOCS 2012) closed the gap with the double greedy algorithms: a deterministic 1/31/31/3-approximation and a randomized 1/21/21/2-approximation, both linear in the number of oracle calls. Their Appendix A gives a third, fractional variant, which is the subject of this mission.

This is the third mission on the FOCS 2012 paper; the first two treat the deterministic and the randomized double greedy on sets.

Setting

Let N\mathcal NN be a finite ground set with nnn elements and f:2N→R≥0f : 2^{\mathcal N} \to \mathbb R_{\ge 0}f:2N→R≥0​. The function fff is submodular if

f(A)+f(B)≥f(A∪B)+f(A∩B)for all A,B⊆N.f(A) + f(B) \ge f(A \cup B) + f(A \cap B) \qquad \text{for all } A, B \subseteq \mathcal N .f(A)+f(B)≥f(A∪B)+f(A∩B)for all A,B⊆N.

Write f(OPT)=max⁡S⊆Nf(S)f(OPT) = \max_{S \subseteq \mathcal N} f(S)f(OPT)=maxS⊆N​f(S) and let OPTOPTOPT be a maximizing set.

The multilinear extension of fff is the function on vectors x∈[0,1]Nx \in [0,1]^{\mathcal N}x∈[0,1]N

F(x)=∑S⊆Nf(S)∏u∈Sxu∏u∉S(1−xu)=E[f(R(x))],F(x) = \sum_{S \subseteq \mathcal N} f(S) \prod_{u \in S} x_u \prod_{u \notin S} (1 - x_u) = \mathbb E\bigl[f(R(x))\bigr],F(x)=S⊆N∑​f(S)u∈S∏​xu​u∈/S∏​(1−xu​)=E[f(R(x))],

where the random set R(x)R(x)R(x) contains each element uuu independently with probability xux_uxu​. A set is identified with its characteristic vector, so FFF agrees with fff on {0,1}N\{0,1\}^{\mathcal N}{0,1}N, and {u}\{u\}{u} also denotes the unit vector at uuu. For vectors, x∨yx \vee yx∨y and x∧yx \wedge yx∧y are the coordinate-wise maximum and minimum.

Algorithm 4 (MultilinearUSM). Fix an arbitrary order u1,…,unu_1, \dots, u_nu1​,…,un​ of N\mathcal NN and start from x0=∅x_0 = \emptysetx0​=∅ and y0=Ny_0 = \mathcal Ny0​=N (the vectors 0\mathbf 00 and 1\mathbf 11). In iteration i=1,…,ni = 1, \dots, ni=1,…,n compute

ai=F(xi−1+{ui})−F(xi−1),bi=F(yi−1−{ui})−F(yi−1),a_i = F(x_{i-1} + \{u_i\}) - F(x_{i-1}), \qquad b_i = F(y_{i-1} - \{u_i\}) - F(y_{i-1}),ai​=F(xi−1​+{ui​})−F(xi−1​),bi​=F(yi−1​−{ui​})−F(yi−1​),

set ai′=max⁡{ai,0}a_i' = \max\{a_i, 0\}ai′​=max{ai​,0}, bi′=max⁡{bi,0}b_i' = \max\{b_i, 0\}bi′​=max{bi​,0}, and update

xi=xi−1+ai′ai′+bi′{ui},yi=yi−1−bi′ai′+bi′{ui},x_i = x_{i-1} + \frac{a_i'}{a_i' + b_i'} \{u_i\}, \qquad y_i = y_{i-1} - \frac{b_i'}{a_i' + b_i'} \{u_i\},xi​=xi−1​+ai′​+bi′​ai′​​{ui​},yi​=yi−1​−ai′​+bi′​bi′​​{ui​},

with the convention that the two fractions are 111 and 000 when ai′=bi′=0a_i' = b_i' = 0ai′​=bi′​=0. The output is the random set R(xn)R(x_n)R(xn​). Every choice before the output is deterministic; the algorithm queries FFF at four points per element.

For the analysis, OPTi=(OPT∨xi)∧yiOPT_i = (OPT \vee x_i) \wedge y_iOPTi​=(OPT∨xi​)∧yi​.

Formalization targets

Goal: Theorem A.1, oracle-access clause

For every nonnegative submodular fff and every order of the ground set,

xn=ynandf(OPT)≤2 F(xn)=2 E[f(R(xn))].x_n = y_n \qquad\text{and}\qquad f(OPT) \le 2\,F(x_n) = 2\,\mathbb E\bigl[f(R(x_n))\bigr].xn​=yn​andf(OPT)≤2F(xn​)=2E[f(R(xn​))].

Milestones, in the order the proof uses them

  1. ai+bi≥0a_i + b_i \ge 0ai​+bi​≥0 at every iteration (proof of Lemma A.2; the page cites Lemma II.1).
  2. Endpoints: OPT0=OPTOPT_0 = OPTOPT0​=OPT with F(OPT)=f(OPT)F(OPT) = f(OPT)F(OPT)=f(OPT), and OPTn=xn=ynOPT_n = x_n = y_nOPTn​=xn​=yn​.
  3. (4) and (5): if ai≥0a_i \ge 0ai​≥0 and bi>0b_i > 0bi​>0, then F(xi)−F(xi−1)=ai2/(ai+bi)F(x_i) - F(x_{i-1}) = a_i^2/(a_i+b_i)F(xi​)−F(xi−1​)=ai2​/(ai​+bi​) and F(yi)−F(yi−1)=bi2/(ai+bi)F(y_i) - F(y_{i-1}) = b_i^2/(a_i+b_i)F(yi​)−F(yi−1​)=bi2​/(ai​+bi​).
  4. (6): in the same case, F(OPTi−1)−F(OPTi)≤aibi/(ai+bi)F(OPT_{i-1}) - F(OPT_i) \le a_i b_i/(a_i + b_i)F(OPTi−1​)−F(OPTi​)≤ai​bi​/(ai​+bi​), whether or not ui∈OPTu_i \in OPTui​∈OPT.
  5. Lemma A.2: for every 1≤i≤n1 \le i \le n1≤i≤n,
F(OPTi−1)−F(OPTi)≤12[F(xi)−F(xi−1)+F(yi)−F(yi−1)].F(OPT_{i-1}) - F(OPT_i) \le \tfrac12\bigl[F(x_i) - F(x_{i-1}) + F(y_i) - F(y_{i-1})\bigr].F(OPTi−1​)−F(OPTi​)≤21​[F(xi​)−F(xi−1​)+F(yi​)−F(yi−1​)].
  1. Telescoped display: F(OPT0)−F(OPTn)≤12[F(xn)−F(x0)]+12[F(yn)−F(y0)]≤12(F(xn)+F(yn))F(OPT_0) - F(OPT_n) \le \tfrac12[F(x_n) - F(x_0)] + \tfrac12[F(y_n) - F(y_0)] \le \tfrac12(F(x_n) + F(y_n))F(OPT0​)−F(OPTn​)≤21​[F(xn​)−F(x0​)]+21​[F(yn​)−F(y0​)]≤21​(F(xn​)+F(yn​)).

Significance

The result. Theorem A.1 shows that the double greedy analysis survives a change of domain: the factor 1/21/21/2 is obtained by a procedure that never flips a coin until the end, and whose state is a pair of fractional points. The ratio matches the Feige–Mirrokni–Vondrák hardness bound, so it cannot be improved in the value oracle model. Its output is a fractional point together with an independent rounding, which separates the optimization from the rounding step.

Formalizing it. The result is proved on paper; no machine-checked proof of a double greedy guarantee is known. A complete development yields reusable facts about the multilinear extension of a submodular function on a finite type: FFF is affine in each coordinate, its coordinate increments are antitone in the other coordinates on [0,1]N[0,1]^{\mathcal N}[0,1]N, and FFF restricted to characteristic vectors is fff. These are the standard tools of every continuous-relaxation argument for submodular maximization.

Difficulty

The proof on the page is short, but it relies on two facts it does not prove. First, the page justifies ai+bi≥0a_i + b_i \ge 0ai​+bi​≥0 "by Lemma II.1", which is a statement about sets; for vectors it requires that the increment of FFF along a coordinate decreases as the other coordinates increase, a property of the multilinear extension of a submodular function that must be derived from the sum defining FFF. Second, inequality (6) is written out only for ui∉OPTu_i \notin OPTui​∈/OPT, and Case 2 of Lemma A.2 is omitted as analogous; the formal statements cover all cases. The main technical work is the bookkeeping of the run: that each coordinate is touched once, that xi−1(ui)=0x_{i-1}(u_i) = 0xi−1​(ui​)=0 and yi−1(ui)=1y_{i-1}(u_i) = 1yi−1​(ui​)=1 when it is touched, that xi≤OPTi≤yix_i \le OPT_i \le y_ixi​≤OPTi​≤yi​, and that every state stays in [0,1]N[0,1]^{\mathcal N}[0,1]N, where the antitonicity applies.

Formalization scope

  • The ground set is a Fintype XXX with decidable equality; sets are Finset X; fff is real-valued, with nonnegativity a hypothesis ∀ S, 0 ≤ f S wherever the page uses it (the goal and the telescoped display). Submodularity is the published NonmonotoneSubmod.Shared.Submodular, the lattice form f(S∪T)+f(S∩T)≤f(S)+f(T)f(S \cup T) + f(S \cap T) \le f(S) + f(T)f(S∪T)+f(S∩T)≤f(S)+f(T); f(OPT)f(OPT)f(OPT) is the published NonmonotoneSubmod.Shared.OPT; FFF is the published NonmonotoneSubmod.Shared.F, the sum above, defined for every x:X→Rx : X \to \mathbb Rx:X→R.
  • The order u1,…,unu_1, \dots, u_nu1​,…,un​ is a duplicate-free list containing every element; uiu_iui​ is the entry at index i−1i-1i−1, and nnn is the list's length. The state after iii iterations is obtained by folding one step over the first iii entries from (0,1)(\mathbf 0, \mathbf 1)(0,1). Statements hold for every such order.
  • The footnote's convention ai′/(ai′+bi′)=1a_i'/(a_i'+b_i') = 1ai′​/(ai′​+bi′​)=1, bi′/(ai′+bi′)=0b_i'/(a_i'+b_i') = 0bi′​/(ai′​+bi′​)=0 when ai′=bi′=0a_i' = b_i' = 0ai′​=bi′​=0 is an explicit case split; with Lean's 0/0=00/0 = 00/0=0 it would otherwise be reversed and the run would no longer end with xn=ynx_n = y_nxn​=yn​.
  • Corrected slips of the page: lines 3–4 of Algorithm 4 assign ai′,bi′a_i', b_i'ai′​,bi′​ but define ai,bia_i, b_iai​,bi​; "f:N→R+f : \mathcal N \to \mathbb R^+f:N→R+" means f:2N→R+f : 2^{\mathcal N} \to \mathbb R^+f:2N→R+; "F(x)≜E[R(x)]F(x) \triangleq \mathbb E[R(x)]F(x)≜E[R(x)]" means E[f(R(x))]\mathbb E[f(R(x))]E[f(R(x))]; "NSM" in Theorem A.1 means USM. The main text's one-line definition of submodularity, read literally, forces monotonicity; the footnote's lattice form is used.
  • Not formalized: the sampling clause of Theorem A.1 (ratio (1/2)−o(1)(1/2) - o(1)(1/2)−o(1) without oracle access to FFF, whose proof the paper refers to Calinescu, Chekuri, Pál and Vondrák) and the running time. The guarantee is stated for the algorithm as printed, so the trivial existence of a 1/21/21/2-approximation by exhaustive search does not satisfy it. A statement in which xnx_nxn​ is an arbitrary point, or the state any process with xi≤yix_i \le y_ixi​≤yi​, would not be this theorem.
  • Contributions welcome: the multilinear-extension facts above as general lemmas, the run invariants, and proofs of the milestones in any order.

Selected references

  • N. Buchbinder, M. Feldman, J. Naor, R. Schwartz, A Tight Linear Time (1/2)-Approximation for Unconstrained Submodular Maximization, FOCS 2012, 649–658. https://doi.org/10.1109/FOCS.2012.73 (journal version: SIAM J. Comput. 44(5), 2015, https://doi.org/10.1137/130929205; its numbering differs and is not used here).
  • U. Feige, V. S. Mirrokni, J. Vondrák, Maximizing Non-monotone Submodular Functions, SIAM J. Comput. 40(4), 2011, 1133–1153. https://doi.org/10.1137/090779346
  • S. Oveis Gharan, J. Vondrák, Submodular Maximization by Simulated Annealing, SODA 2011, 1098–1117. https://doi.org/10.1137/1.9781611973082.83
  • M. Feldman, J. Naor, R. Schwartz, Nonmonotone Submodular Maximization via a Structural Continuous Greedy Algorithm, ICALP 2011, 342–353. https://doi.org/10.1007/978-3-642-22006-7_29
  • G. Calinescu, C. Chekuri, M. Pál, J. Vondrák, Maximizing a Monotone Submodular Function Subject to a Matroid Constraint, SIAM J. Comput. 40(6), 2011, 1740–1766. https://doi.org/10.1137/080733991
11 thms3 active usersReviewed
🏆Completed
Operations ResearchProbability·Captain: mikedeng1

On the Power of Randomization in On-Line Algorithms 1: α-Competitiveness Against Adaptive On-Line and β Against Oblivious Adversaries Give a Deterministic α∘β-Competitive AlgorithmResearch Paper

Why randomization matters in online algorithms

An online algorithm must answer each request before it sees the next one. Its performance is compared with an optimum that may choose all its answers after seeing the complete request string. Randomization can improve an online algorithm's guarantee when the request string is fixed in advance. The comparison changes when an adversary chooses later requests after seeing the algorithm's earlier answers. Ben-David, Borodin, Karp, Tardos and Wigderson studied these choices of adversary in a common request-answer model and proved a general relation between their competitive guarantees (Ben-David et al., 1994, manuscript §§2–3).

The paper distinguishes three adversaries. An oblivious adversary fixes the request string before the algorithm's random choices affect any answer. An adaptive off-line adversary chooses the next request from previous answers but serves the resulting request string optimally after the play. An adaptive on-line adversary also chooses its own answer as each request arrives. The ability to react to answers makes the latter two adversaries materially different from the oblivious one for randomized algorithms (Ben-David et al., 1994, manuscript pp. 7–9).

Request-answer games and competitive cost

A request-answer game has a request set RRR, a finite nonempty answer set AAA, and a real cost fn(r,a)f_n(r,a)fn​(r,a) for a request string r∈Rnr\in R^nr∈Rn and an answer string a∈Ana\in A^na∈An. The off-line optimum for rrr is c(r)=min⁡a∈Anfn(r,a)c(r)=\min_{a\in A^n}f_n(r,a)c(r)=mina∈An​fn​(r,a). A deterministic online algorithm DDD returns its iiith answer from the first iii requests alone; it has no access to the rest of rrr or to the eventual stopping time. Its cost on rrr is cD(r)=fn(r,D(r))c_D(r)=f_n(r,D(r))cD​(r)=fn​(r,D(r)).

A randomized online algorithm is a distribution over deterministic online algorithms. With coins ω\omegaω, write GωG_\omegaGω​ for the resulting deterministic algorithm. For a fixed request string rrr, GGG is β\betaβ-competitive against oblivious adversaries when Eω[cGω(r)]≤β(c(r))\mathbb E_\omega[c_{G_\omega}(r)]\leq\beta(c(r))Eω​[cGω​​(r)]≤β(c(r)). The paper calls a transformation “linear” when it has the affine form x↦ux+vx\mapsto ux+vx↦ux+v (Ben-David et al., 1994, manuscript p. 7).

An adaptive off-line adversary QQQ has a rule from prior answer strings to either the next request or a stop signal, together with a common finite upper bound on play length. Let r(Gω,Q)r(G_\omega,Q)r(Gω​,Q) denote its request string and cQ(Gω)=c(r(Gω,Q))c_Q(G_\omega)=c(r(G_\omega,Q))cQ​(Gω​)=c(r(Gω​,Q)). Its competitiveness condition places the transformation inside the expectation: Eω[cGω(Q)]≤Eω[α(cQ(Gω))]\mathbb E_\omega[c_{G_\omega}(Q)]\leq\mathbb E_\omega[\alpha(c_Q(G_\omega))]Eω​[cGω​​(Q)]≤Eω​[α(cQ​(Gω​))]. An adaptive on-line adversary SSS has the same request rule and an additional answer rule; its own cost is cS(Gω)c_S(G_\omega)cS​(Gω​), and the corresponding condition uses Eω[α(cS(Gω))]\mathbb E_\omega[\alpha(c_S(G_\omega))]Eω​[α(cS​(Gω​))] on the right (Ben-David et al., 1994, manuscript pp. 8–9).

Formalization targets

Randomization against adaptive off-line adversaries

The first target is Theorem 2.1: if some randomized algorithm is α\alphaα-competitive against every adaptive off-line adversary, a deterministic algorithm has that same guarantee on every request string:

∃G  ∀Q,E[cG(Q)]≤E[α(cQ(G))]⟹∃D  ∀r,cD(r)≤α(c(r)).\exists G\;\forall Q,\quad \mathbb E[c_G(Q)]\leq\mathbb E[\alpha(c_Q(G))]\quad\Longrightarrow\quad\exists D\;\forall r,\quad c_D(r)\leq\alpha(c(r)).∃G∀Q,E[cG​(Q)]≤E[α(cQ​(G))]⟹∃D∀r,cD​(r)≤α(c(r)).

Composition of two guarantees

Theorem 2.2 takes an α\alphaα guarantee for GGG against adaptive on-line adversaries and a β\betaβ guarantee for another randomized algorithm against oblivious adversaries. It concludes that GGG has the composed guarantee against adaptive off-line adversaries:

E[cG(Q)]≤E[(α∘β)(cQ(G))]for every Q.\mathbb E[c_G(Q)]\leq\mathbb E[(\alpha\circ\beta)(c_Q(G))]\qquad\text{for every }Q.E[cG​(Q)]≤E[(α∘β)(cQ​(G))]for every Q.

The mission goal is Corollary 2.1, the deterministic consequence of these two results:

∃D  ∀r,cD(r)≤(α∘β)(c(r)).\exists D\;\forall r,\qquad c_D(r)\leq(\alpha\circ\beta)(c(r)).∃D∀r,cD​(r)≤(α∘β)(c(r)).

The milestones follow the paper's two theorems and the stated claims in their proofs, including the finite-horizon winning-position formulation and the adversary that simulates a fixed online algorithm (Ben-David et al., 1994, manuscript pp. 9–13).

What the result supplies

The corollary turns the existence of two randomized guarantees under different information rules into the existence of a deterministic online strategy with an explicit composed cost transformation. It is an existence result: it does not say that the deterministic strategy can be computed efficiently from the randomized algorithms. The paper itself notes that such a construction is unavailable in full generality and then examines settings where constructive versions are possible (Ben-David et al., 1994, manuscript p. 13).

The mathematical results were proved in the 1994 paper; this mission asks for machine-checked Lean proofs of the abstract model, the intermediate claims, and Corollary 2.1. The local draft currently contains compiled statements with proof placeholders, so it does not yet provide checked proofs. A completed development would make the adversary distinctions and the exact placement of expectations available for reuse in later online-algorithm formalizations.

Why the proof is difficult

The apparent shortcut is to treat an adaptive request sequence as fixed and apply a guarantee against oblivious adversaries directly. That loses the dependence of later requests on the algorithm's earlier answers. For Theorem 2.1, a winning request strategy must have one finite horizon that works for every answer path; separate finite horizons for each branch do not suffice when the answer set is infinite. For Theorem 2.2, the simulated adversary must make its own answers before the algorithm answers the current request, while still matching a fixed online benchmark along every resulting play. The expectation inequalities must remain valid when the request string itself depends on the algorithm's coins (Ben-David et al., 1994, manuscript pp. 9–11).

Formalization scope

Lean represents requests and answers as oldest-first lists. List index zero is request one in the paper. The general game is a separate definition; the algorithm, adversary, and competitiveness definitions build on it. An off-line adversary's rule returns Option R, where none is the stop signal, and has a uniform finite depth bound. A randomized algorithm consists of a coin probability space and a deterministic prefix algorithm for each coin; its answer events are measurable. Finiteness of AAA and bounded play depth make the cost of each fixed adversarial play take finitely many values, so its real expectation is an ordinary integrable expectation.

The formal game uses real-valued costs, a deliberate restriction of the paper's R∪{∞}\mathbb R\cup\{\infty\}R∪{∞} costs. The answer set is finite and nonempty, while the request set may be infinite. The transformations α\alphaα and β\betaβ are affine. Theorem 2.2 and the goal assume α\alphaα is monotone: the paper applies α\alphaα to an inequality in its proof, and its competitive-ratio examples have positive slope. Theorem 2.1 does not need this added assumption. The two randomized algorithms may have different coin spaces, each an arbitrary Lean type at the declaration's universe level. The off-line and on-line adaptive comparisons retain α\alphaα inside the expectation.

The target ranges over every equal-length request and answer play generated by these rules, including an adversary that stops without a request. It does not allow the deterministic algorithm to see future requests or choose a different policy for each adversary. Reusable contributions include the game interface, bounded adaptive plays, measurable randomized algorithms, and finite-horizon winning positions. The statements of all three principal results, their intervening claims, and proofs of those statements are within scope.

Selected references

  • S. Ben-David, A. Borodin, R. Karp, G. Tardos and A. Wigderson, On the Power of Randomization in On-Line Algorithms, Algorithmica 11, 1994. DOI: 10.1007/BF01294260. The local source is the authors' 20-page manuscript; citations above use its page numbers.
11 thms3 active usersReviewed
🏆Completed
Operations ResearchProbability·Captain: mikedeng1

On the Power of Randomization in On-Line Algorithms 2: The Bound α∘β Against Adaptive Off-Line Adversaries Is TightResearch Paper

Motivation

An on-line algorithm must answer each request as it arrives, without knowing the requests to come; paging, caching, the kkk-server problem and metrical task systems are standard examples. Its quality is measured by competitive analysis: its cost is compared with the cost of an optimal off-line solution that knows the whole request sequence. For randomized on-line algorithms the comparison depends on how much the adversary producing the requests is allowed to see. Ben-David, Borodin, Karp, Tardos and Wigderson (Algorithmica 11, 1994; conference version STOC 1990) introduced the three standard adversaries — oblivious, adaptive on-line and adaptive off-line — and related the competitive ratios achievable against each.

Their Theorem 2.2 (manuscript p. 10) shows that if a randomized algorithm is α\alphaα-competitive against adaptive on-line adversaries and some randomized algorithm is β\betaβ-competitive against oblivious adversaries, then the first algorithm is αβ\alpha\betaαβ-competitive against adaptive off-line adversaries. This mission formalizes the paper's claim (manuscript p. 11) that this product bound cannot be improved in general, together with the explicit construction on pp. 12–13 that proves it.

Setting

A request-answer game consists of a request set RRR, a finite answer set AAA, and cost functions fn:Rn×An→Rf_n : R^n \times A^n \to \mathbb Rfn​:Rn×An→R. For a request sequence r‾∈Rn\underline r \in R^nr​∈Rn, the off-line optimum is c(r‾)=min⁡a‾∈Anfn(r‾,a‾)c(\underline r) = \min_{\underline a \in A^n} f_n(\underline r, \underline a)c(r​)=mina​∈An​fn​(r​,a​). A deterministic on-line algorithm GGG answers the iii-th request with ai=gi(r1,…,ri)a_i = g_i(r_1, \dots, r_i)ai​=gi​(r1​,…,ri​); a randomized one is a probability distribution over deterministic algorithms GxG_xGx​, xxx being the coin tosses.

An adaptive off-line adversary QQQ chooses each request ri+1=qi(a1,…,ai)r_{i+1} = q_i(a_1, \dots, a_i)ri+1​=qi​(a1​,…,ai​) from the answers given so far, stops after at most dQd_QdQ​ requests, and pays the off-line optimum cQ(G)=c(r‾)c_Q(G) = c(\underline r)cQ​(G)=c(r​) of the requests it made; the algorithm pays cG(Q)=fn(r‾,a‾)c_G(Q) = f_n(\underline r, \underline a)cG​(Q)=fn​(r​,a​). An adaptive on-line adversary SSS must in addition answer each request itself, before the algorithm does, with bi+1=pi(a1,…,ai)b_{i+1} = p_i(a_1, \dots, a_i)bi+1​=pi​(a1​,…,ai​), and pays cS(G)=fn(r‾,b‾)c_S(G) = f_n(\underline r, \underline b)cS​(G)=fn​(r​,b​). An oblivious adversary fixes r‾\underline rr​ in advance and pays c(r‾)c(\underline r)c(r​). A randomized GGG is α\alphaα-competitive against oblivious adversaries if Ex[cGx(r‾)]≤α c(r‾)\mathbb E_x[c_{G_x}(\underline r)] \le \alpha\, c(\underline r)Ex​[cGx​​(r​)]≤αc(r​) for all r‾\underline rr​, and against adaptive on-line adversaries if Ex[cGx(S)]≤Ex[α cS(Gx)]\mathbb E_x[c_{G_x}(S)] \le \mathbb E_x[\alpha\, c_S(G_x)]Ex​[cGx​​(S)]≤Ex​[αcS​(Gx​)] for all SSS.

The construction uses the mates game: R=AR = AR=A is a set of 2t2t2t elements split into ttt pairs of mates, and for n≥2n \ge 2n≥2 the cost depends only on the first answer a1a_1a1​ and the second request r2r_2r2​: it is 111 if a1=r2a_1 = r_2a1​=r2​, MMM if a1a_1a1​ is the mate of r2r_2r2​, and mmm otherwise. The algorithm GGG draws a1a_1a1​ uniformly at random. The parameters solve

β=(2t−2)m+M+12t,α=1+(2t−1)M2+(2t−2)m.\beta = \frac{(2t-2)m + M + 1}{2t}, \qquad \alpha = \frac{1 + (2t-1)M}{2 + (2t-2)m}.β=2t(2t−2)m+M+1​,α=2+(2t−2)m1+(2t−1)M​.

Formalization targets

Goal: tightness of Theorem 2.2

For 1<β≤α1 < \beta \le \alpha1<β≤α (or α=β=1\alpha = \beta = 1α=β=1) and every C<αβC < \alpha\betaC<αβ, there are a game and a randomized algorithm GGG such that

G is α-competitive against adaptive on-line adversaries,G is β-competitive against oblivious adversaries,G \text{ is } \alpha\text{-competitive against adaptive on-line adversaries}, \qquad G \text{ is } \beta\text{-competitive against oblivious adversaries},G is α-competitive against adaptive on-line adversaries,G is β-competitive against oblivious adversaries,

and for every randomized algorithm KKK some adaptive off-line adversary QQQ achieves

E[cQ(K)]>0,E[cK(Q)]≥C⋅E[cQ(K)].\mathbb E[c_Q(K)] > 0, \qquad \mathbb E[c_K(Q)] \ge C\cdot \mathbb E[c_Q(K)].E[cQ​(K)]>0,E[cK​(Q)]≥C⋅E[cQ​(K)].

Milestones (pp. 12–13)

  1. The closed forms m(t)m(t)m(t), M(t)M(t)M(t) are the unique solution of the two equations.
  2. m(t)→βm(t) \to \betam(t)→β and M(t)→αβM(t) \to \alpha\betaM(t)→αβ as t→∞t \to \inftyt→∞.
  3. For all large ttt: M(t)≥max⁡(m(t)2,C)M(t) \ge \max(m(t)^2, C)M(t)≥max(m(t)2,C), 1≤m(t)≤M(t)1 \le m(t) \le M(t)1≤m(t)≤M(t), α(m(t)−1)≤M(t)−m(t)\alpha(m(t)-1) \le M(t) - m(t)α(m(t)−1)≤M(t)−m(t).
  4. GGG is β\betaβ-competitive against oblivious adversaries in the mates game.
  5. GGG is α\alphaα-competitive against adaptive on-line adversaries in the mates game.
  6. An adaptive off-line adversary makes every algorithm pay MMM while paying 111.

Significance

The result. Together with Theorem 2.2, the claim pins down exactly how much the adaptive off-line adversary can gain over the other two: the product αβ\alpha\betaαβ is an upper bound for every game and is approached by a single game for every admissible pair (α,β)(\alpha, \beta)(α,β). It shows that no general argument relating the three adversary models can give a bound better than the product, so any improvement for a specific problem (paging, kkk-server) must use the structure of that problem. The paging example cited on p. 11 (RANDOM against the three adversaries) gives one instance of tightness; the mates game gives tightness for every admissible pair.

Formalizing it. The result is proved in the paper, in about one page, with two steps left to the reader ("by inspection of the equations", "a simple case analysis"). No machine-checked proof of this or of any statement about adaptive adversaries is known to us. The formalization makes the model of §2 precise (sequences, stopping, the order in which adversary and algorithm commit, expectations over coins), checks the asymptotics of the parameters, and verifies the case analysis, which on inspection needs an inequality the page does not state. Two printed formulas on p. 12 contain typos; the formal statements carry the correct values.

Difficulty

The construction is explicit, but each competitiveness claim quantifies over all adversaries, which may adapt their requests to the algorithm's random answers, stop at any time, and (for the on-line adversary) commit to their own answers in advance. The algebra of α\alphaα-competitiveness is tight: the adversary's best expected advantage is exactly zero, so every case of its best reply must be checked with no slack. The page's condition M≥m2M \ge m^2M≥m2 does not suffice for this: when a1a_1a1​ is neither the adversary's first answer nor its mate, the reply "mate of a1a_1a1​" beats the reply "the adversary's own answer" only when α(m−1)≤M−m\alpha(m-1) \le M - mα(m−1)≤M−m, which holds for the solved parameters but is not implied by M≥m2M \ge m^2M≥m2. At β=1<α\beta = 1 < \alphaβ=1<α the solved parameter mmm is below 111 for every ttt, and the oblivious bound fails.

Formalization scope

All declarations live in OnlineRandomization.Tightness. The conventions:

  • Costs are real-valued; the paper allows +∞+\infty+∞, so the game class is a special case.
  • Answer sets are nonempty finite types; request sets are arbitrary types.
  • Sequences are Lean lists, oldest first; cost r a is fnf_nfn​ on lists of equal length nnn.
  • Adversaries return none for "stop" and carry a depth bound dQd_QdQ​; an on-line adversary's answer bi+1b_{i+1}bi+1​ depends only on a1,…,aia_1, \dots, a_ia1​,…,ai​.
  • Randomized algorithms are a probability space of coins with a deterministic algorithm per coin and measurable answers; expectations are Bochner integrals, with α\alphaα applied inside the expectation. In the goal, coin spaces range over Type.
  • Competitiveness uses the ratio functions x↦αxx \mapsto \alpha xx↦αx and x↦βxx \mapsto \beta xx↦βx, with no additive constant.
  • The mates game is on Fin t × Bool, with mate (i,b)↦(i,¬b)(i, b) \mapsto (i, \lnot b)(i,b)↦(i,¬b). The paper leaves the costs of plays with fewer than two requests undefined; the formalization sets f0=0f_0 = 0f0​=0 and f1≡1f_1 \equiv 1f1​≡1 (with f1≡0f_1 \equiv 0f1​≡0 the algorithm would not be α\alphaα-competitive).
  • Range. The goal assumes 1<β≤α1 < \beta \le \alpha1<β≤α or α=β=1\alpha = \beta = 1α=β=1; the page's case β=1<α\beta = 1 < \alphaβ=1<α is not covered by its construction and is left out. In fact the claim is false there for 1<C<α1 < C < \alpha1<C<α: an algorithm that is 111-competitive against oblivious adversaries answers optimally, almost surely, on every request sequence (its cost is never below the optimum and its expected cost does not exceed it), and an adaptive off-line adversary reaches only finitely many request sequences, so against K=GK = GK=G every adversary has E[cG(Q)]=E[cQ(G)]\mathbb E[c_G(Q)] = \mathbb E[c_Q(G)]E[cG​(Q)]=E[cQ​(G)], a ratio of 1<C1 < C1<C.

The positivity requirement E[cQ(K)]>0\mathbb E[c_Q(K)] > 0E[cQ​(K)]>0 in the goal is essential: without it the adversary that asks nothing satisfies E[cK(Q)]≥C⋅0\mathbb E[c_K(Q)] \ge C \cdot 0E[cK​(Q)]≥C⋅0 for every KKK, and the third clause would hold vacuously.

A complete development needs: finite expectations over a uniform coin, the evaluation of the play of an adversary against a constant algorithm, and limit and eventual-inequality arguments for rational functions of ttt. The model of §2 is shared with the other missions of this series and is reusable for any request-answer formulation of an on-line problem. Proofs of individual milestones are welcome.

Selected references

  • S. Ben-David, A. Borodin, R. Karp, G. Tardos, A. Wigderson, On the power of randomization in on-line algorithms, Algorithmica 11 (1994) 2–14. https://doi.org/10.1007/BF01294260 (cited from the authors' manuscript, manuscript pp. 7–13).
  • A. Borodin, R. El-Yaniv, Online Computation and Competitive Analysis, Cambridge University Press, 1998. ISBN 0-521-56392-5.
  • P. Raghavan, M. Snir, Memory versus randomization in on-line algorithms, IBM Journal of Research and Development 38 (1994) 683–707. https://doi.org/10.1147/rd.386.0683
9 thms3 active usersReviewed
🏆Completed
Operations ResearchProbability·Captain: mikedeng1

Secretary Problems: Weights and Discounts 3: An O(log n)-Competitive Algorithm for the Discounted Secretary ProblemResearch Paper

Motivation

In the secretary problem, nnn candidates with arbitrary values arrive one at a time in a uniformly random order, and an online decision maker must accept or reject each candidate on arrival, irrevocably, keeping at most one. The rule that observes the first n/en/en/e candidates and then accepts the first one better than everything seen so far selects the best candidate with probability about 1/e1/e1/e (Dynkin, 1963). The problem is a basic model of online selection and, read economically, of posted-price mechanisms for agents who arrive in random order: a rule that compares each agent only against a threshold set by earlier agents is truthful.

Babaioff, Dinitz, Gupta, Immorlica and Talwar (SODA 2009) study a variant in which time costs value. Selecting the candidate who arrives at time ttt earns that candidate's value multiplied by a discount d(t)d(t)d(t), for an arbitrary non-negative discount function ddd known in advance. Earlier work treated only specific discount shapes, such as geometric discounting d(t)=βtd(t)=\beta^td(t)=βt (Rasmussen and Pliska, 1976). For a general ddd the classical rule can fail badly: if all the discount mass sits in the first few time steps, a rule that waits through a sample of size n/en/en/e earns nothing. The paper shows that the best competitive ratio for arbitrary discounts lies between Ω(log⁡n/log⁡log⁡n)\Omega(\log n/\log\log n)Ω(logn/loglogn) (its Theorem 4.3) and O(log⁡n)O(\log n)O(logn) (its Theorem 4.4). This mission formalizes the upper bound.

Setting

There are n≥1n\ge1n≥1 elements, indexed {0,…,n−1}\{0,\dots,n-1\}{0,…,n−1}, with values v(e)≥0v(e)\ge 0v(e)≥0, and nnn times with discounts d(t)≥0d(t)\ge0d(t)≥0. A uniformly random permutation π\piπ fixes the order of arrivals: element π(t)\pi(t)π(t) arrives at time ttt. An algorithm knows ddd but not vvv; it sees each value on arrival and may select the current element, irrevocably, earning d(t) v(π(t))d(t)\,v(\pi(t))d(t)v(π(t)). Expectations over π\piπ are exact averages over the n!n!n! orders.

The offline optimum on order π\piπ is OPT(π)=max⁡td(t) v(π(t))\mathsf{OPT}(\pi)=\max_t d(t)\,v(\pi(t))OPT(π)=maxt​d(t)v(π(t)); it is a random variable, and the benchmark is its expectation Eπ[OPT]\mathbb E_\pi[\mathsf{OPT}]Eπ​[OPT] (p. 4 of the paper).

Let dmax⁡=max⁡td(t)d_{\max}=\max_t d(t)dmax​=maxt​d(t) and vmax⁡=max⁡ev(e)v_{\max}=\max_e v(e)vmax​=maxe​v(e). For c≥1c\ge1c≥1 the ccc-th discount class is the set of times

Pc={ i:2−cdmax⁡<d(i)≤2−(c−1)dmax⁡ }.P_c=\{\,i : 2^{-c}d_{\max}<d(i)\le 2^{-(c-1)}d_{\max}\,\}.Pc​={i:2−cdmax​<d(i)≤2−(c−1)dmax​}.

The quantity OPTc\mathsf{OPT}_cOPTc​ is the part of Eπ[OPT]\mathbb E_\pi[\mathsf{OPT}]Eπ​[OPT] earned when the optimal time (the smallest time attaining the maximum) lies in PcP_cPc​.

The classical secretary rule on mmm arrivals observes the first ⌊m/e⌋\lfloor m/e\rfloor⌊m/e⌋ and then selects the first arrival that ranks above every earlier arrival. Ranks use a fixed tie-break order: larger value first, and smaller element index among equal values.

The algorithm A\mathcal AA sets M=3⌈log⁡2n⌉+2M=3\lceil\log_2 n\rceil+2M=3⌈log2​n⌉+2, draws c∈{1,…,M}c\in\{1,\dots,M\}c∈{1,…,M} uniformly, and runs the classical rule on the subsequence of arrivals at the times of PcP_cPc​, ignoring all other arrivals.

Formalization targets

Goal: Theorem 4.4 with its explicit constant

Eπ[OPT]  ≤  4e (3⌈log⁡2n⌉+2)  E[A](n≥1, d≥0, v≥0).\mathbb E_\pi[\mathsf{OPT}]\;\le\;4e\,\bigl(3\lceil\log_2 n\rceil+2\bigr)\;\mathbb E[\mathcal A]\qquad(n\ge1,\ d\ge0,\ v\ge0).Eπ​[OPT]≤4e(3⌈log2​n⌉+2)E[A](n≥1, d≥0, v≥0).

The paper states E[OPT]/E[A]≤O(log⁡n)\mathbb E[\mathsf{OPT}]/\mathbb E[\mathcal A]\le O(\log n)E[OPT]/E[A]≤O(logn); the constant 4e4e4e is the one its proof yields.

Milestones

  1. The classical secretary rule selects the top-ranked of m≥1m\ge1m≥1 elements with probability at least 1/e1/e1/e (§2, p. 4).
  2. OPT1≥vmax⁡dmax⁡/n\mathsf{OPT}_1\ge v_{\max}d_{\max}/nOPT1​≥vmax​dmax​/n (proof of Theorem 4.4, p. 7).
  3. OPTc≤2−c 2n2dmax⁡vmax⁡\mathsf{OPT}_c\le 2^{-c}\,2n^2d_{\max}v_{\max}OPTc​≤2−c2n2dmax​vmax​ for every c≥1c\ge1c≥1 (p. 7).
  4. ∑c=13⌈log⁡2n⌉+1OPTc≥12Eπ[OPT]\sum_{c=1}^{3\lceil\log_2 n\rceil+1}\mathsf{OPT}_c\ge\tfrac12\mathbb E_\pi[\mathsf{OPT}]∑c=13⌈log2​n⌉+1​OPTc​≥21​Eπ​[OPT] (p. 7).
  5. E[Ac]≥OPTc/2e\mathbb E[\mathcal A_c]\ge\mathsf{OPT}_c/2eE[Ac​]≥OPTc​/2e for every c≥1c\ge1c≥1, where Ac\mathcal A_cAc​ is the classical rule on PcP_cPc​ (p. 7).

Significance

The theorem shows that a general discount function costs only a logarithmic factor against the offline benchmark, and that one algorithm achieves this without any knowledge of the values. Together with the lower bound of Theorem 4.3 it pins the competitive ratio of the discounted secretary problem between log⁡n/log⁡log⁡n\log n/\log\log nlogn/loglogn and log⁡n\log nlogn. The same scale-splitting idea, stated in the paper as Theorem 4.5 without full proof, extends the bound to the weighted discounted problem.

The result is proved in the paper; to our knowledge it has not been formalized. A complete development would also produce a machine-checked proof of the classical secretary guarantee for the rule with sample size exactly ⌊m/e⌋\lfloor m/e\rfloor⌊m/e⌋ at every finite mmm, with an explicit tie-break, which is reusable by every secretary-type mission. Milestone 1 is that statement. Sharper constants or a smaller class range are welcome as additional statements but do not replace the goal, which is about this algorithm with this MMM.

Difficulty

The obvious argument, running the classical rule on all nnn arrivals, fails because the discounts can be concentrated at times the rule spends sampling. Splitting by discount scale fixes this but creates two problems. First, there are unboundedly many scales, and one has to show that the offline optimum's mass outside the top O(log⁡n)O(\log n)O(logn) of them is negligible against E[OPT]\mathbb E[\mathsf{OPT}]E[OPT], a random quantity rather than a fixed maximum. Second, the classical rule on a class sees only a random subset of the elements, in random order, and the guarantee must be transferred to this subsequence, conditioning on which elements land in PcP_cPc​. Neither step is deep, but both require careful bookkeeping of permutations, and the classical 1/e1/e1/e bound at finite mmm with a floor in the sample size is itself a nontrivial estimate.

Formalization scope

Elements and times are Fin n, an order is π : Equiv.Perm (Fin n) read as time ↦\mapsto↦ element, and the paper's time t=1,…,nt=1,\dots,nt=1,…,n is index t−1t-1t−1. Values and discounts are Fin n → ℝ with non-negativity hypotheses. Every expectation is the finite average 1n!∑π\frac1{n!}\sum_\pin!1​∑π​; the algorithm's random class is the explicit average 1M∑c=1M\frac1M\sum_{c=1}^MM1​∑c=1M​. Maxima are suprema over the finite index set. The logarithm is base 2, ⌈log⁡2n⌉\lceil\log_2 n\rceil⌈log2​n⌉ is Nat.clog 2 n, and the sample size is Nat.floor (m / Real.exp 1). Ties are broken by the order on Lex (ℝ × (Fin n)ᵒᵈ) (larger value, then smaller index); distinct values are not assumed. The optimal time is the smallest maximizing time, so that the OPTc\mathsf{OPT}_cOPTc​ add up to E[OPT]\mathbb E[\mathsf{OPT}]E[OPT]. Competitiveness is stated multiplicatively, never as a quotient, so E[A]=0\mathbb E[\mathcal A]=0E[A]=0 is not a loophole.

The goal is a statement about the specific algorithm A\mathcal AA, not "there exists an algorithm": an existential over unrestricted algorithms is witnessed by a clairvoyant rule that reads the values in advance. A\mathcal AA sees the values only through comparisons among arrivals that have already occurred, and E[OPT]\mathbb E[\mathsf{OPT}]E[OPT] is the expected offline maximum over the same random order, not dmax⁡vmax⁡d_{\max}v_{\max}dmax​vmax​.

Needed infrastructure: averages over permutations and the fact that the elements landing at a fixed set of times form a uniformly random subset in uniformly random order; the finite-mmm analysis of the classical rule; and elementary estimates on geometric sums. Contributions of general lemmas about uniform permutations are welcome and reusable.

Selected references

  • M. Babaioff, M. Dinitz, A. Gupta, N. Immorlica, K. Talwar, Secretary Problems: Weights and Discounts, Proc. 20th ACM-SIAM Symposium on Discrete Algorithms (SODA), 2009. https://doi.org/10.1137/1.9781611973068.139
  • E. B. Dynkin, The optimum choice of the instant for stopping a Markov process, Soviet Math. Doklady 4, 1963.
  • T. S. Ferguson, Who solved the secretary problem?, Statistical Science 4(3), 1989. https://doi.org/10.1214/ss/1177012493
  • L. T. Rasmussen, S. R. Pliska, Choosing the maximum from a sequence with a discount function, Applied Mathematics and Optimization 2, 1976. https://doi.org/10.1007/BF01458209
  • M. Babaioff, N. Immorlica, R. Kleinberg, Matroids, secretary problems, and online mechanisms, SODA 2007. https://dl.acm.org/doi/10.5555/1283383.1283429
9 thms3 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchOptimization·Captain: mikedeng1

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

Motivation

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

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

Timeline:

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

Setting

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

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

The proof uses three auxiliary problems:

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

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

Formalization targets

Goal: Theorem 8

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

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

Milestones

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

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

Significance

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

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

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

Difficulty

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

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

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

Formalization scope

Conventions committed to in Lean:

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

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

Out of scope:

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

Infrastructure that a complete development needs:

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

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

Selected references

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

Online Primal-Dual Algorithms for Maximizing Ad-Auctions Revenue: The Competitive Ratio of the Primal-Dual Allocation AlgorithmResearch Paper

Motivation

Search engines sell advertisement slots next to their results through ad-auctions. Advertisers bid on keywords, and each advertiser also sets a daily budget: the most it is willing to pay in a day. Queries arrive one at a time and each must be assigned to an advertiser at once, with no knowledge of the queries still to come. The seller's revenue from an advertiser is capped by its budget, so an allocation rule that ignores budgets can exhaust a high bidder early and forgo revenue that a more even allocation would have collected. The question is how much of the offline optimum an online rule can guarantee against every arrival sequence.

Mehta, Saberi, Vazirani and Vazirani (FOCS 2005 / J. ACM 2007) gave a deterministic algorithm whose competitive ratio tends to 1−1/e1 - 1/e1−1/e when bids are small compared with budgets, and showed that no deterministic algorithm does better. Their algorithm builds on online bipartite matching (Karp, Vazirani and Vazirani, STOC 1990) and online bbb-matching (Kalyanasundaram and Pruhs, 2000). Buchbinder, Jain and Naor (ESA 2007) rederived the 1−1/e1 - 1/e1−1/e bound with an online primal-dual algorithm, which gives the ratio in closed form for every value of the bid-to-budget ratio and extends to multiple slots, stochastic information, bounded degree and budget flexibility. This mission formalizes the basic algorithm of that paper and its Theorem 1.

Setting

There is a finite nonempty set III of buyers. Buyer iii has a known budget B(i)>0B(i) > 0B(i)>0. Products j=1,…,mj = 1, \dots, mj=1,…,m arrive one by one; when product jjj arrives, every buyer's bid b(i,j)≥0b(i,j) \ge 0b(i,j)≥0 on it is revealed. The bid-to-budget ratio is

Rmax⁡=max⁡i∈I, jb(i,j)B(i).R_{\max} = \max_{i \in I,\, j} \frac{b(i,j)}{B(i)} .Rmax​=i∈I,jmax​B(i)b(i,j)​.

A fractional allocation y(i,j)≥0y(i,j) \ge 0y(i,j)≥0 assigns fractions of products to buyers; the revenue from buyer iii is the minimum of ∑jb(i,j) y(i,j)\sum_j b(i,j)\,y(i,j)∑j​b(i,j)y(i,j) and B(i)B(i)B(i).

The offline fractional problem is the packing LP, which the paper calls the dual:

max⁡∑j∑ib(i,j) y(i,j)s.t.∑iy(i,j)≤1  ∀j,∑jb(i,j) y(i,j)≤B(i)  ∀i,y≥0.\max \sum_{j}\sum_{i} b(i,j)\,y(i,j) \quad\text{s.t.}\quad \sum_i y(i,j) \le 1 \ \ \forall j,\qquad \sum_j b(i,j)\,y(i,j) \le B(i)\ \ \forall i,\qquad y \ge 0 .maxj∑​i∑​b(i,j)y(i,j)s.t.i∑​y(i,j)≤1  ∀j,j∑​b(i,j)y(i,j)≤B(i)  ∀i,y≥0.

Its LP dual, the paper's primal, is the covering LP:

min⁡∑iB(i) x(i)+∑jz(j)s.t.b(i,j) x(i)+z(j)≥b(i,j)  ∀i,j,x,z≥0.\min \sum_i B(i)\,x(i) + \sum_j z(j) \quad\text{s.t.}\quad b(i,j)\,x(i) + z(j) \ge b(i,j)\ \ \forall i,j,\qquad x, z \ge 0 .mini∑​B(i)x(i)+j∑​z(j)s.t.b(i,j)x(i)+z(j)≥b(i,j)  ∀i,j,x,z≥0.

The Allocation Algorithm has a parameter c>1c > 1c>1 and starts from x≡0x \equiv 0x≡0. When product jjj arrives it takes a buyer iii maximizing b(i,j)(1−x(i))b(i,j)(1 - x(i))b(i,j)(1−x(i)). If x(i)≥1x(i) \ge 1x(i)≥1, the product is not sold. Otherwise it charges iii the minimum of b(i,j)b(i,j)b(i,j) and iii's remaining budget, sets y(i,j)←1y(i,j) \leftarrow 1y(i,j)←1 and z(j)←b(i,j)(1−x(i))z(j) \leftarrow b(i,j)(1 - x(i))z(j)←b(i,j)(1−x(i)), and updates

x(i)←x(i)(1+b(i,j)B(i))+b(i,j)(c−1) B(i).x(i) \leftarrow x(i)\Big(1 + \frac{b(i,j)}{B(i)}\Big) + \frac{b(i,j)}{(c-1)\,B(i)} .x(i)←x(i)(1+B(i)b(i,j)​)+(c−1)B(i)b(i,j)​.

Its revenue is the total amount charged.

Formalization targets

Goal: Theorem 1

For every instance and every bound R>0R > 0R>0 with b(i,j)≤R B(i)b(i,j) \le R\,B(i)b(i,j)≤RB(i) for all i,ji, ji,j, the Allocation Algorithm run with c=(1+R)1/Rc = (1+R)^{1/R}c=(1+R)1/R, under any tie-breaking of the maximum, satisfies for every feasible y′y'y′ of the packing LP

Revenue  ≥  (1−1c)(1−R)∑j∑ib(i,j) y′(i,j).\mathrm{Revenue} \;\ge\; \Big(1 - \frac1c\Big)(1 - R)\sum_{j}\sum_{i} b(i,j)\,y'(i,j).Revenue≥(1−c1​)(1−R)j∑​i∑​b(i,j)y′(i,j).

With R=Rmax⁡R = R_{\max}R=Rmax​ this is the paper's statement that the algorithm is (1−1/c)(1−Rmax⁡)(1 - 1/c)(1 - R_{\max})(1−1/c)(1−Rmax​)-competitive; the fractional optimum bounds every integral offline allocation.

Milestones

The proof of Theorem 1 rests on three claims and three auxiliary facts, each a milestone:

  1. the inequality ln⁡(1+x)/x≥ln⁡(1+y)/y\ln(1+x)/x \ge \ln(1+y)/yln(1+x)/x≥ln(1+y)/y for 0<x≤y≤10 < x \le y \le 10<x≤y≤1;
  2. Claim (1): the final (x,z)(x, z)(x,z) is feasible for the covering LP;
  3. Claim (2): the covering cost of the run equals (1+1/(c−1))(1 + 1/(c-1))(1+1/(c−1)) times the packing value of the run's own yyy;
  4. Inequality (1): x(i)≥1c−1(c∑jb(i,j)y(i,j)/B(i)−1)x(i) \ge \frac{1}{c-1}\big(c^{\sum_j b(i,j) y(i,j)/B(i)} - 1\big)x(i)≥c−11​(c∑j​b(i,j)y(i,j)/B(i)−1) at every stage of the run;
  5. Claim (3): ∑jb(i,j) y(i,j)≤B(i)+max⁡jb(i,j)\sum_j b(i,j)\,y(i,j) \le B(i) + \max_j b(i,j)∑j​b(i,j)y(i,j)≤B(i)+maxj​b(i,j), and the amount charged to iii is at least (1−R)∑jb(i,j) y(i,j)(1 - R)\sum_j b(i,j)\,y(i,j)(1−R)∑j​b(i,j)y(i,j);
  6. weak duality for the LP pair above;

and, separately, the second sentence of Theorem 1,

lim⁡R→0+(1−1(1+R)1/R)(1−R)=1−1e.\lim_{R\to 0^+}\Big(1 - \frac{1}{(1+R)^{1/R}}\Big)(1-R) = 1 - \frac1e .R→0+lim​(1−(1+R)1/R1​)(1−R)=1−e1​.

Significance

Theorem 1 gives an explicit ratio for every value of Rmax⁡R_{\max}Rmax​, not only in the limit. It tends to the optimal deterministic ratio 1−1/e1 - 1/e1−1/e as bids become small, and it quantifies how the guarantee degrades as single bids become a larger share of a budget. The primal-dual analysis is the template for the paper's later sections and for a line of work on online packing and covering problems, surveyed in Buchbinder and Naor's monograph The Design of Competitive Online Algorithms via a Primal-Dual Approach (Foundations and Trends in TCS, 2009).

The result is proved in the paper, and the proof is short. What this mission adds is a machine-checked proof about an algorithm that is defined, not described: the run is computed by recursion from the instance, and the guarantee is proved for that run and every tie-breaking. A related private mission on the platform, The Design of Competitive Online Algorithms via a Primal-Dual Approach VI: Maximizing Ad-Auctions Revenue, states the monograph's Theorem 10.1, which is this theorem, in a form that takes the analysis's intermediate inequalities as hypotheses over arbitrary lists of won bids; the present mission states it for the algorithm itself. No machine-checked proof of Theorem 1 is known to this mission.

Difficulty

Each step of the proof is elementary; the difficulty is the bookkeeping of an online process. Claims (1) and (2) are statements about a single iteration that must be lifted to the whole run: Claim (1) uses that xxx only increases, and Claim (2) that each product is processed once. Inequality (1) is an induction over the iterations that allocate to one buyer, interleaved with iterations that allocate to others and is the only place where the value of ccc matters. Claim (3) needs a further invariant: the amount charged equals the minimum of the allocated bids and the budget.

A tempting shortcut is to take Inequality (1) and the "at most one undercharge" fact as hypotheses about some list of bids. That does not describe the algorithm and is not the theorem; here the only hypotheses are on the instance and on the tie-breaking rule.

Formalization scope

Buyers are a type I with [Fintype I] and [Nonempty I]; products are Fin m, whose order is the arrival order. Bids and budgets are real, with B(i)>0B(i) > 0B(i)>0 and b(i,j)≥0b(i,j) \ge 0b(i,j)≥0. The state of the algorithm records xxx, the amounts charged, yyy and zzz; one iteration is step, the run after kkk products is runPrefix, and revenue sums the charges of the final state. The tie-breaking rule is a function sel of the current xxx and the product, required to return a maximizer of b(i,j)(1−x(i))b(i,j)(1-x(i))b(i,j)(1−x(i)); the theorem holds for every such rule. The constant c=(1+R)1/Rc = (1+R)^{1/R}c=(1+R)1/R is a real power and requires R>0R > 0R>0. The theorem is stated for any bound RRR on the ratios, of which the exact maximum is one instance. Claims (1) and (2) are stated for every c>1c > 1c>1, which covers the paper's choice. The paper's inequality for ln⁡(1+x)/x\ln(1+x)/xln(1+x)/x allows x=0x = 0x=0, read as a limit; the Lean statement requires x>0x > 0x>0.

A statement over an unconstrained allocation, or one conditioned on the proof's own intermediate inequalities, would be trivially true or false; the targets here concern only the run the definitions compute.

The development needs finite sums, real powers and logarithms from Mathlib and an induction principle for the run. The LP pair and weak duality are reusable for the paper's extensions, and the run invariants for any primal-dual online algorithm with multiplicative updates. Proofs of any milestone are welcome, as are sharper variants, such as the exact-Rmax⁡R_{\max}Rmax​ form or the bound against integral allocations.

Selected references

  • N. Buchbinder, K. Jain, J. Naor, Online Primal-Dual Algorithms for Maximizing Ad-Auctions Revenue, Algorithms – ESA 2007, LNCS 4698, 2007. https://doi.org/10.1007/978-3-540-75520-3_24
  • A. Mehta, A. Saberi, U. Vazirani, V. Vazirani, AdWords and Generalized Online Matching, Journal of the ACM 54(5), 2007. https://doi.org/10.1145/1284320.1284321
  • R. M. Karp, U. V. Vazirani, V. V. Vazirani, An Optimal Algorithm for On-line Bipartite Matching, STOC 1990. https://doi.org/10.1145/100216.100262
  • B. Kalyanasundaram, K. R. Pruhs, An Optimal Deterministic Algorithm for Online b-Matching, Theoretical Computer Science 233(1–2), 2000. https://doi.org/10.1016/S0304-3975(99)00140-1
  • N. Buchbinder, J. Naor, The Design of Competitive Online Algorithms via a Primal-Dual Approach, Foundations and Trends in Theoretical Computer Science 3(2–3), 2009. https://doi.org/10.1561/0400000024
10 thms3 active usersReviewed
🏆Completed
Number TheoryProbabilityQuantum Information·Captain: mikedeng1

Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer 3: The Success Probability of Quantum Order FindingResearch Paper

Motivation

The security of the RSA cryptosystem rests on the assumed difficulty of factoring large integers, and the best known classical algorithms for factoring run in super-polynomial time. In 1994 Peter Shor showed that a quantum computer can factor an nnn-digit integer in time polynomial in nnn (Shor, SIAM J. Comput. 1997; conference version FOCS 1994). The algorithm has two parts. A classical reduction, due to Miller (1976), turns factoring into order finding: given xxx coprime to nnn, find the least r≥1r \ge 1r≥1 with xr≡1(modn)x^r \equiv 1 \pmod nxr≡1(modn). The quantum part solves order finding.

This mission formalizes the quantum part as Shor analyzes it in §5 of the journal paper: the construction of the quantum state, the probability of each measurement outcome, and the classical post-processing that reads rrr off the measured value. The paper's claim is that one run of this procedure returns rrr with probability at least φ(r)/3r\varphi(r)/3rφ(r)/3r.

Timeline:

  • 1976: Miller reduces factoring to order finding (with randomization).
  • 1985–1994: Deutsch, Bernstein–Vazirani and Simon give the quantum Fourier sampling ideas the algorithm builds on.
  • 1994: Shor's FOCS paper introduces the factoring and discrete logarithm algorithms.
  • 1997: the SIAM J. Comput. version gives the analysis formalized here, with qqq the power of 222 in [n2,2n2)[n^2, 2n^2)[n2,2n2).

Setting

Fix an integer n≥2n \ge 2n≥2 and an integer xxx coprime to nnn. Its order rrr is the least r≥1r \ge 1r≥1 with xr≡1(modn)x^r \equiv 1 \pmod nxr≡1(modn); since xxx is a unit, r≤φ(n)<nr \le \varphi(n) < nr≤φ(n)<n. Let q=2lq = 2^lq=2l be the power of 222 with n2≤q<2n2n^2 \le q < 2n^2n2≤q<2n2.

A quantum state on two registers, the first holding 0≤a<q0 \le a < q0≤a<q and the second a residue y∈Z/ny \in \mathbb{Z}/ny∈Z/n, is a complex vector ψ(a,y)\psi(a, y)ψ(a,y) indexed by the basis states ∣a,y⟩|a, y\rangle∣a,y⟩. Measuring it returns ∣a,y⟩|a, y\rangle∣a,y⟩ with probability ∣ψ(a,y)∣2|\psi(a, y)|^2∣ψ(a,y)∣2.

The Fourier matrix AqA_qAq​ is the q×qq \times qq×q matrix with entries (Aq)a,c=q−1/2exp⁡(2πiac/q)(A_q)_{a,c} = q^{-1/2}\exp(2\pi i a c/q)(Aq​)a,c​=q−1/2exp(2πiac/q), with rows indexing inputs and columns outputs. The algorithm

  1. prepares 1q1/2∑a=0q−1∣a⟩∣xa mod n⟩\frac{1}{q^{1/2}}\sum_{a=0}^{q-1}|a\rangle|x^a \bmod n\rangleq1/21​∑a=0q−1​∣a⟩∣xamodn⟩ (eq. (5.2)),
  2. applies AqA_qAq​ to the first register, obtaining 1q∑a,cexp⁡(2πiac/q)∣c⟩∣xa mod n⟩\frac1q\sum_{a,c}\exp(2\pi iac/q)|c\rangle|x^a \bmod n\rangleq1​∑a,c​exp(2πiac/q)∣c⟩∣xamodn⟩ (eq. (5.4)),
  3. measures, obtaining some ∣c,y⟩|c, y\rangle∣c,y⟩,
  4. rounds c/qc/qc/q to the nearest fraction with denominator smaller than nnn.

The observed ccc gives us rrr if some fraction with lowest-terms denominator below nnn is within 1/2q1/2q1/2q of c/qc/qc/q, and every such fraction has lowest-terms denominator exactly rrr. In the Lean development these objects are preFourierState, finalState, outcomeProb and yieldsOrder, in the namespace ShorAlgorithms.OrderFinding, and the shared definition ShorAlgorithms.Shared.fourierMatrix.

Formalization targets

Goal: success probability at least φ(r)/3r\varphi(r)/3rφ(r)/3r

For all sufficiently large nnn, with xxx, rrr and qqq as above,

Pr⁡[the observed c gives us r]  =  ∑c gives r ∑y∈Z/n∣Ψ(c,y)∣2  ≥  φ(r)3r,\Pr\bigl[\text{the observed } c \text{ gives us } r\bigr] \;=\; \sum_{c\ \text{gives}\ r}\ \sum_{y \in \mathbb{Z}/n} |\Psi(c, y)|^2 \;\ge\; \frac{\varphi(r)}{3r},Pr[the observed c gives us r]=c gives r∑​ y∈Z/n∑​∣Ψ(c,y)∣2≥3rφ(r)​,

where Ψ\PsiΨ is the state (5.4). The threshold on nnn is uniform in xxx and qqq; it is the paper's "for sufficiently large nnn" from the per-state bound.

Milestones

  1. Eqs. (5.5)–(5.6). For 0≤k<r0 \le k < r0≤k<r, the probability of ∣c,xk⟩|c, x^k\rangle∣c,xk⟩ equals ∣1q∑b=0⌊(q−k−1)/r⌋exp⁡(2πi(br+k)c/q)∣2\left|\frac1q\sum_{b=0}^{\lfloor (q-k-1)/r\rfloor}\exp(2\pi i(br+k)c/q)\right|^2​q1​∑b=0⌊(q−k−1)/r⌋​exp(2πi(br+k)c/q)​2.
  2. Eq. (5.11). For nnn past a threshold, every ∣c,xk⟩|c, x^k\rangle∣c,xk⟩ with −r/2≤rc−dq≤r/2-r/2 \le rc - dq \le r/2−r/2≤rc−dq≤r/2 for some integer ddd has probability at least 1/3r21/3r^21/3r2.
  3. Eq. (5.13). If n2≤qn^2 \le qn2≤q, at most one fraction with denominator below nnn lies within 1/2q1/2q1/2q of c/qc/qc/q.
  4. p. 1500. Such a fraction is a convergent of the continued fraction of c/qc/qc/q.
  5. p. 1501. At least φ(r)\varphi(r)φ(r) values of ccc are within 1/2q1/2q1/2q of some d/rd/rd/r with gcd⁡(d,r)=1\gcd(d, r) = 1gcd(d,r)=1; with the rrr distinct values of xkx^kxk this gives at least rφ(r)r\varphi(r)rφ(r) states ∣c,xk⟩|c, x^k\rangle∣c,xk⟩, and each such ccc gives us rrr.

Significance

The goal is the quantitative statement behind "order finding is in bounded-error quantum polynomial time": since φ(r)/r≥δ/log⁡log⁡r\varphi(r)/r \ge \delta/\log\log rφ(r)/r≥δ/loglogr for a constant δ\deltaδ (Hardy and Wright, Thm. 328), O(log⁡log⁡r)O(\log\log r)O(loglogr) repetitions find rrr with high probability, and Miller's reduction then factors nnn. Without the bound, the algorithm is a procedure with no guarantee.

The result is proved, in the paper and in textbooks (Nielsen and Chuang, 2000, §5.3), usually with a phase-estimation analysis rather than Shor's direct count. What this mission adds is a machine-checked proof of Shor's own argument, with his choice of qqq and his constants, starting from the state built by applying AqA_qAq​ to (5.2). Formal proofs of idealized versions exist elsewhere, for instance in the exact-period model where rrr divides qqq and the output is uniform on rrr peaks, but that model removes the approximation that the 1/3r21/3r^21/3r2 bound is about. Legendre's theorem on continued fractions is already on the platform (FamousTheorems.legendre_continued_fraction_theorem) and is included as a reference item.

Difficulty

The obvious route is to compute the output distribution in closed form. That works only when rrr divides qqq; here qqq is a power of 222 and rrr is arbitrary, so the amplitudes are geometric sums of ⌊(q−k−1)/r⌋+1\lfloor (q-k-1)/r\rfloor + 1⌊(q−k−1)/r⌋+1 terms whose phases do not cancel exactly. The per-state bound 1/3r21/3r^21/3r2 requires a lower bound on such a sum that is uniform in rrr, ccc and kkk, with error terms of order 1/q1/q1/q controlled against a main term of order 1/r21/r^21/r2. The constant 1/31/31/3 leaves only a small margin below the limiting value 4/π2≈0.4054/\pi^2 \approx 0.4054/π2≈0.405, so the errors must be bounded explicitly, not merely shown to vanish.

The second difficulty is the counting: distinct coprime numerators ddd must give distinct outcomes ccc in [0,q)[0, q)[0,q), and each good ccc must determine rrr uniquely, which uses r<nr < nr<n and n2≤qn^2 \le qn2≤q.

Formalization scope

Conventions the statements commit to:

  • States are functions Fin q × ZMod n → ℂ; the matrix convention is row = input, so applying AqA_qAq​ to the first register gives the amplitude ∑aψ(a,y)(Aq)a,c\sum_a \psi(a, y)(A_q)_{a,c}∑a​ψ(a,y)(Aq​)a,c​ at (c,y)(c, y)(c,y).
  • The final state is built by applying AqA_qAq​ to the state (5.2); the closed forms (5.5) and (5.6) are theorems, not definitions. No normalization hypothesis is assumed.
  • Probabilities are squared moduli; the probability of the event "ccc gives us rrr" sums over all y∈Z/ny \in \mathbb{Z}/ny∈Z/n, which is exact because yyy that are not powers of xxx have probability zero.
  • xxx is a natural number with gcd⁡(x,n)=1\gcd(x, n) = 1gcd(x,n)=1; rrr is orderOf (x : ZMod n). qqq enters through the three hypotheses q=2lq = 2^lq=2l, n2≤qn^2 \le qn2≤q, q<2n2q < 2n^2q<2n2, not through a function of nnn.
  • Fractions are rationals, and "in lowest terms" is Rat.den.
  • Thresholds "for sufficiently large nnn" are ∃N, ∀n≥N\exists N,\ \forall n \ge N∃N, ∀n≥N, with NNN quantified before xxx, qqq, ccc and kkk.
  • Condition (5.11) is stated in its equivalent form (5.12), with an integer ddd.
  • Printed slip. Eq. (5.13)'s justification says "Because q>n2q > n^2q>n2", but qqq was chosen with n2≤qn^2 \le qn2≤q, and q=n2q = n^2q=n2 when nnn is a power of 222. The uniqueness claim holds under n2≤qn^2 \le qn2≤q, and that is what is stated.

Typing the closed form (5.4)–(5.6) in as the definition of the final state would make milestone 1 trivial and hide whether the probability model is the paper's; the definitions exclude this by construction.

Not stated: the polynomial running time of any step, the O(log⁡log⁡r)O(\log\log r)O(loglogr) repetition count (no explicit constant), the reversible modular exponentiation of §3, and the post-processing heuristics on p. 1501. Needed infrastructure: bounds on geometric exponential sums, Euler's totient, Diophantine approximation by fractions with bounded denominator, and Mathlib's continued fractions. Lemmas on geometric sums of roots of unity and on the order of units mod nnn are reusable in the companion discrete logarithm mission.

Selected references

  • P. W. Shor, Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer, SIAM J. Comput. 26(5):1484–1509, 1997. https://doi.org/10.1137/S0097539795293172
  • P. W. Shor, Algorithms for quantum computation: discrete logarithms and factoring, Proc. 35th FOCS, 1994. https://doi.org/10.1109/SFCS.1994.365700
  • G. L. Miller, Riemann's hypothesis and tests for primality, J. Comput. System Sci. 13(3):300–317, 1976. https://doi.org/10.1016/S0022-0000(76)80043-8
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford, 1979 (Ch. X, continued fractions; Thm. 328).
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge, 2000. https://doi.org/10.1017/CBO9780511976667
12 thms3 active usersReviewed
🏆Completed
CombinatoricsComplexity TheoryOperations Research+2·Captain: mikedeng1

Maximizing Non-Monotone Submodular Functions V: Beating 1/2 for Symmetric Functions Requires Exponentially Many Value QueriesResearch Paper

Motivation

Maximizing a nonnegative submodular set function without constraints contains Max Cut, Max Directed Cut and facility-location problems as special cases. In the value-oracle model an algorithm knows nothing about the function except the values f(S)f(S)f(S) of the sets SSS it queries, and it is judged by the number of queries it makes. Feige, Mirrokni and Vondrák (SIAM J. Comput. 40(4), 2011) gave constant-factor algorithms in this model and matching limits on what any algorithm can do. For symmetric functions, such as cut functions of undirected graphs, a uniformly random set already achieves 12\tfrac1221​ of the optimum in expectation (Theorem 2.1 of the paper). The question this mission formalizes is whether any algorithm can do better, and the answer given by Theorem 4.5 is: not without exponentially many value queries. The same factor 12\tfrac1221​ was later shown to be achievable for general (non-symmetric) nonnegative submodular functions by Buchbinder, Feldman, Naor and Schwartz (FOCS 2012 / SIAM J. Comput. 2015), so the bound of Theorem 4.5 is the tight limit of the whole problem in the value-oracle model.

Timeline:

  • 2007 (FOCS) / 2011 (SIAM J. Comput.): Feige, Mirrokni and Vondrák prove that no algorithm with subexponentially many value queries achieves (12+ϵ)(\tfrac12 + \epsilon)(21​+ϵ) of the optimum on symmetric nonnegative submodular functions, and give 25\tfrac2552​ for general functions.
  • 2011: Vondrák's symmetry-gap framework (SIAM J. Comput. 42(1), 2013) generalizes the construction to constrained problems.
  • 2012: Buchbinder, Feldman, Naor and Schwartz give a randomized 12\tfrac1221​-approximation for general nonnegative submodular functions, matching the bound.

Setting

Let [n]={0,…,n−1}[n] = \{0, \dots, n-1\}[n]={0,…,n−1} be the ground set, with nnn even. A set function f:2[n]→Rf : 2^{[n]} \to \mathbb{R}f:2[n]→R is submodular if f(S∪T)+f(S∩T)≤f(S)+f(T)f(S \cup T) + f(S \cap T) \le f(S) + f(T)f(S∪T)+f(S∩T)≤f(S)+f(T) for all S,TS, TS,T, symmetric if f([n]∖S)=f(S)f([n]\setminus S) = f(S)f([n]∖S)=f(S) for all SSS, and OPT(f)=max⁡Sf(S)\mathrm{OPT}(f) = \max_{S} f(S)OPT(f)=maxS​f(S).

Fix an integer mmm with 1≤m≤n/21 \le m \le n/21≤m≤n/2 and write ϵ=m/n\epsilon = m/nϵ=m/n, so that ϵn\epsilon nϵn is an integer. For integers k,ℓk, \ellk,ℓ put

f(k,ℓ)={(k+ℓ)(n−k−ℓ)∣k−ℓ∣≤m,k(n−2ℓ)+(n−2k)ℓ+m2−2m∣k−ℓ∣∣k−ℓ∣>m.f(k,\ell) = \begin{cases} (k+\ell)(n-k-\ell) & |k-\ell| \le m,\\ k(n-2\ell) + (n-2k)\ell + m^2 - 2m|k-\ell| & |k-\ell| > m. \end{cases}f(k,ℓ)={(k+ℓ)(n−k−ℓ)k(n−2ℓ)+(n−2k)ℓ+m2−2m∣k−ℓ∣​∣k−ℓ∣≤m,∣k−ℓ∣>m.​

For a set C⊆[n]C \subseteq [n]C⊆[n] with ∣C∣=n/2|C| = n/2∣C∣=n/2 and D=[n]∖CD = [n] \setminus CD=[n]∖C, the hard instance is fC(S)=f(∣S∩C∣,∣S∩D∣)f_C(S) = f(|S\cap C|, |S\cap D|)fC​(S)=f(∣S∩C∣,∣S∩D∣). The cut function of the complete graph is g(S)=∣S∣(n−∣S∣)g(S) = |S|(n-|S|)g(S)=∣S∣(n−∣S∣), with maximum 14n2\tfrac14 n^241​n2. A set QQQ is balanced for CCC if ∣∣Q∩C∣−∣Q∩D∣∣≤m\bigl||Q\cap C| - |Q\cap D|\bigr| \le m​∣Q∩C∣−∣Q∩D∣​≤m; on balanced sets fC=gf_C = gfC​=g.

A deterministic adaptive qqq-query algorithm AAA chooses each query from the answers received so far, and after qqq answers outputs a set A(h)A(h)A(h) when run against an oracle hhh. A randomized algorithm is a distribution μ\muμ over deterministic ones, with expected value EA∼μ[h(A(h))]\mathbb{E}_{A\sim\mu}[h(A(h))]EA∼μ​[h(A(h))].

Formalization targets

Goal: Theorem 4.5 with the constants of its proof

For every such n,mn, mn,m:

  1. every fCf_CfC​ with ∣C∣=n/2|C| = n/2∣C∣=n/2 is nonnegative, symmetric and submodular, with
OPT(fC)=12n2(1−2ϵ+2ϵ2);\mathrm{OPT}(f_C) = \tfrac12 n^2 (1 - 2\epsilon + 2\epsilon^2);OPT(fC​)=21​n2(1−2ϵ+2ϵ2);
  1. for every q<eϵ2n/8q < e^{\epsilon^2 n/8}q<eϵ2n/8 and every randomized qqq-query algorithm μ\muμ there is a CCC with ∣C∣=n/2|C| = n/2∣C∣=n/2 and
EA∼μ[fC(A(fC))]≤14n2+(2e−ϵ2n/8+2e−ϵ2n/4) OPT(fC).\mathbb{E}_{A\sim\mu}\bigl[f_C(A(f_C))\bigr] \le \tfrac14 n^2 + \bigl(2e^{-\epsilon^2 n/8} + 2e^{-\epsilon^2 n/4}\bigr)\,\mathrm{OPT}(f_C).EA∼μ​[fC​(A(fC​))]≤41​n2+(2e−ϵ2n/8+2e−ϵ2n/4)OPT(fC​).

Hence the ratio attained is at most 12(1−2ϵ+2ϵ2)+4e−ϵ2n/8=12+ϵ+O(ϵ2)+4e−ϵ2n/8\frac{1}{2(1-2\epsilon+2\epsilon^2)} + 4e^{-\epsilon^2 n/8} = \tfrac12 + \epsilon + O(\epsilon^2) + 4e^{-\epsilon^2 n/8}2(1−2ϵ+2ϵ2)1​+4e−ϵ2n/8=21​+ϵ+O(ϵ2)+4e−ϵ2n/8.

Milestones

  • Theorem 1.2, the Chernoff bound for independent variables in [−1,1][-1,1][−1,1].
  • Submodularity of fCf_CfC​.
  • The value OPT(fC)=12n2(1−2ϵ+2ϵ2)\mathrm{OPT}(f_C) = \tfrac12 n^2(1 - 2\epsilon + 2\epsilon^2)OPT(fC​)=21​n2(1−2ϵ+2ϵ2), attained at S=CS = CS=C.
  • A fixed query is unbalanced for at most a 2e−ϵ2n/42e^{-\epsilon^2 n/4}2e−ϵ2n/4 fraction of the half-size sets CCC.
  • If all queries are balanced, the algorithm cannot distinguish fCf_CfC​ from ggg.
  • The deterministic case of the bound, averaged over CCC.

Significance

The theorem shows that the factor 12\tfrac1221​ for symmetric submodular maximization, and hence for unconstrained submodular maximization in general, cannot be improved by any algorithm that uses a subexponential number of value queries, whatever its running time. It is an information-theoretic bound and needs no complexity assumption. Along with the later matching 12\tfrac1221​-approximation, it settles the value-oracle approximability of the problem. The construction, a function equal to a symmetric function on "balanced" sets and larger elsewhere, is the prototype of the symmetry-gap technique used for many later oracle lower bounds.

The result is proved in the paper. No machine-checked version is known to exist: the platform has no value-oracle or query-lower-bound statement. Formalizing it requires a precise model of adaptive randomized query algorithms, a concentration bound for the hypergeometric distribution, and a finite verification of submodularity of an explicit two-regime function, and it fixes the constants that the printed statement leaves as O(⋅)O(\cdot)O(⋅) terms.

Difficulty

Two steps of the printed argument do not go through as written. First, the proof bounds the probability that a fixed query is unbalanced by citing the Chernoff bound for independent variables, but for a uniformly random half-size set CCC the count ∣Q∩C∣|Q\cap C|∣Q∩C∣ is hypergeometric, and the summands are not independent. A bound for sampling without replacement is needed instead. Replacing the balanced partition by independent coin flips is not an option: then ∣C∣≠n/2|C| \ne n/2∣C∣=n/2 in general, and the function is no longer the paper's instance.

Second, the argument counts only the queries, but the value an algorithm receives is fCf_CfC​ of its output, which equals ggg of the output only if the output is balanced as well. That event has to be controlled too.

Finally, submodularity of fCf_CfC​ must be checked across the boundary ∣k−ℓ∣=ϵn|k-\ell| = \epsilon n∣k−ℓ∣=ϵn between the two regimes, where the formula changes.

Formalization scope

  • The ground set is Fin n with nnn even; ϵn\epsilon nϵn is an integer mmm with 1≤m1 \le m1≤m and 2m≤n2m \le n2m≤n, following the paper's "assume that ϵn\epsilon nϵn is an integer". Sets are Finset (Fin n), and all values are real.
  • OPT\mathrm{OPT}OPT is Finset.sup' over all subsets. There is no junk value.
  • The partition (C,D)(C, D)(C,D) is uniform over half-size sets; probabilities over it are counts of n/2-subsets divided by (nn/2)\binom{n}{n/2}(n/2n​), written multiplied out.
  • A deterministic algorithm is a pair of decision rules query, output : List ℝ → Finset (Fin n) making exactly qqq adaptive queries with arbitrary real answers. A randomized algorithm is a PMF over deterministic algorithms, which covers every randomization with countable support. The algorithm sees fff only through query answers; it never receives CCC.
  • Pinned-down constants. The printed theorem, "fewer than eϵ2n/8e^{\epsilon^2 n/8}eϵ2n/8 queries" and "expected value at least (12+ϵ)OPT(\tfrac12+\epsilon)\mathrm{OPT}(21​+ϵ)OPT", is not what the proof gives for one and the same ϵ\epsilonϵ. On the proof's instances OPT=12n2(1−2ϵ+2ϵ2)\mathrm{OPT} = \tfrac12 n^2(1-2\epsilon+2\epsilon^2)OPT=21​n2(1−2ϵ+2ϵ2), and the ratio held is 12(1−2ϵ+2ϵ2)>12+ϵ\frac{1}{2(1-2\epsilon+2\epsilon^2)} > \tfrac12 + \epsilon2(1−2ϵ+2ϵ2)1​>21​+ϵ. The formal goal states the explicit bound the proof establishes. The literal printed pair, stated for the proof's family with the same ϵ\epsilonϵ, is false: the zero-query algorithm that outputs a fixed half-size set gets at least 14n2>(12+ϵ)OPT\tfrac14 n^2 > (\tfrac12+\epsilon)\mathrm{OPT}41​n2>(21​+ϵ)OPT.
  • Added term. The error term 2e−ϵ2n/42e^{-\epsilon^2 n/4}2e−ϵ2n/4 for the output set is added to the paper's 2e−ϵ2n/82e^{-\epsilon^2 n/8}2e−ϵ2n/8.
  • Ruled-out trivializations. A restricted algorithm class (non-adaptive, deterministic, or one that must return a queried set) would give a different, weaker theorem. So would a bound that lets the algorithm read CCC, which would make the statement false. Both the instance's properties (nonnegativity, symmetry, submodularity, the value of OPT) and the bound are part of the goal, so an empty or degenerate family cannot satisfy it. The quantifier order is: for every algorithm there is an instance.
  • Needed infrastructure: a value-oracle algorithm model; tail bounds for the hypergeometric distribution (Hoeffding's inequality for sampling without replacement), which Mathlib lacks; averaging over a PMF of algorithms. The algorithm model and the hypergeometric bound are reusable for other oracle lower bounds. Proofs of any milestone, and alternative derivations of the balance bound, are welcome.

Selected references

  • U. Feige, V. S. Mirrokni, J. Vondrák, Maximizing Non-Monotone Submodular Functions, SIAM J. Comput. 40(4):1133–1153, 2011. https://doi.org/10.1137/090779346
  • N. Alon, J. H. Spencer, The Probabilistic Method, Wiley (source of Theorem 1.2).
  • W. Hoeffding, Probability Inequalities for Sums of Bounded Random Variables, J. Amer. Statist. Assoc. 58(301):13–30, 1963. https://doi.org/10.1080/01621459.1963.10500830
  • J. Vondrák, Symmetry and Approximability of Submodular Maximization Problems, SIAM J. Comput. 42(1):265–304, 2013. https://doi.org/10.1137/110832318
  • N. Buchbinder, M. Feldman, J. Naor, R. Schwartz, A Tight Linear Time (1/2)-Approximation for Unconstrained Submodular Maximization, SIAM J. Comput. 44(5):1384–1402, 2015. https://doi.org/10.1137/130929205
12 thms3 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchOptimization+1·Captain: mikedeng1

Maximizing Non-Monotone Submodular Functions IV: Smooth Local Search Achieves 2/5 of the OptimumResearch Paper

Motivation

Many optimization problems ask for a subset of a finite ground set that maximizes a submodular function, a set function with diminishing marginal returns. Max Cut and Max Directed Cut in graphs, facility location with fixed costs, and the maximization of mutual information or entropy of a subset of random variables are all of this form. Unlike the monotone case, where a greedy algorithm achieves 1−1/e1 - 1/e1−1/e, a general nonnegative submodular function may decrease when elements are added, and the empty set and the full set can both be poor. The question is how large a constant fraction of the optimum a polynomial-time algorithm can guarantee when the function is given only through an oracle that returns f(S)f(S)f(S) for a queried set SSS.

Feige, Mirrokni and Vondrák (SIAM J. Comput. 40(4), 2011) gave the first constant-factor algorithms for this problem. A uniformly random set achieves 1/41/41/4 of the optimum, a deterministic local search achieves 1/3−ϵ/n1/3 - \epsilon/n1/3−ϵ/n, and a randomized smooth local search achieves 2/5−o(1)2/5 - o(1)2/5−o(1). The last result is the paper's best approximation for general nonnegative submodular functions (Table 1, p. 1136), and it is the subject of this mission.

Timeline. Feige, Mirrokni and Vondrák: 1/41/41/4, 1/31/31/3 and 2/52/52/5 (FOCS 2007; journal version 2011). Gharan and Vondrák (SODA 2011): about 0.410.410.41 by simulated annealing. Buchbinder, Feldman, Naor and Schwartz (FOCS 2012; SIAM J. Comput. 2015): a randomized double greedy algorithm achieving 1/21/21/2, which matches the 1/21/21/2 hardness in the value oracle model proved in the same paper by Feige, Mirrokni and Vondrák.

Setting

Let XXX be a finite ground set with n=∣X∣≥1n = |X| \ge 1n=∣X∣≥1 elements and f:2X→Rf : 2^X \to \mathbb{R}f:2X→R a function with f(S)≥0f(S) \ge 0f(S)≥0 for all SSS and

f(S∪T)+f(S∩T)≤f(S)+f(T)(S,T⊆X).f(S \cup T) + f(S \cap T) \le f(S) + f(T) \quad (S, T \subseteq X).f(S∪T)+f(S∩T)≤f(S)+f(T)(S,T⊆X).

Write OPT=max⁡S⊆Xf(S)OPT = \max_{S \subseteq X} f(S)OPT=maxS⊆X​f(S). The multilinear extension of fff is

F(x)=∑S⊆Xf(S)∏i∈Sxi∏j∉S(1−xj),F(x) = \sum_{S \subseteq X} f(S) \prod_{i \in S} x_i \prod_{j \notin S} (1 - x_j),F(x)=S⊆X∑​f(S)i∈S∏​xi​j∈/S∏​(1−xj​),

the expected value of fff on a random set containing each iii independently with probability xix_ixi​.

For A⊆XA \subseteq XA⊆X and δ∈[−1,1]\delta \in [-1,1]δ∈[−1,1], the random set R(A,δ)\mathcal{R}(A,\delta)R(A,δ) is sampled with bias δ\deltaδ based on AAA: each element of AAA is included independently with probability p=(1+δ)/2p = (1+\delta)/2p=(1+δ)/2, each element of B=X∖AB = X \setminus AB=X∖A with probability q=(1−δ)/2q = (1-\delta)/2q=(1−δ)/2. The potential is Φ(A)=E[f(R(A,δ))]\Phi(A) = \mathbf{E}[f(\mathcal{R}(A,\delta))]Φ(A)=E[f(R(A,δ))] and the smoothed marginal value of xxx is

ωA,δ(x)=E[f(R(A,δ)∪{x})]−E[f(R(A,δ)∖{x})].\omega_{A,\delta}(x) = \mathbf{E}[f(\mathcal{R}(A,\delta) \cup \{x\})] - \mathbf{E}[f(\mathcal{R}(A,\delta) \setminus \{x\})].ωA,δ​(x)=E[f(R(A,δ)∪{x})]−E[f(R(A,δ)∖{x})].

Algorithm SLS starts from A=∅A = \emptysetA=∅. At each iteration it obtains estimates ω~A,δ(x)\tilde\omega_{A,\delta}(x)ω~A,δ​(x) within ±1n2OPT\pm\frac{1}{n^2}OPT±n21​OPT of ωA,δ(x)\omega_{A,\delta}(x)ωA,δ​(x). If some x∉Ax \notin Ax∈/A has ω~A,δ(x)>2n2OPT\tilde\omega_{A,\delta}(x) > \frac{2}{n^2}OPTω~A,δ​(x)>n22​OPT it adds xxx; otherwise, if some x∈Ax \in Ax∈A has ω~A,δ(x)<−2n2OPT\tilde\omega_{A,\delta}(x) < -\frac{2}{n^2}OPTω~A,δ​(x)<−n22​OPT it removes xxx; otherwise it stops and returns a random set R(A,δ′)\mathcal{R}(A, \delta')R(A,δ′).

Formalization targets

Goal: Theorem 3.6 in the explicit form of its proof

With δ=1/3\delta = 1/3δ=1/3, and δ′=1/3\delta' = 1/3δ′=1/3 with probability 0.90.90.9 or δ′=−1\delta' = -1δ′=−1 with probability 0.10.10.1: for every run from ∅\emptyset∅ whose estimates are all accurate and which has terminated at AAA,

910 E[f(R(A,13))]+110 f(X∖A)≥(25−95n)OPT,\tfrac{9}{10}\,\mathbf{E}[f(\mathcal{R}(A,\tfrac13))] + \tfrac{1}{10}\,f(X \setminus A) \ge \Big(\frac{2}{5} - \frac{9}{5n}\Big) OPT,109​E[f(R(A,31​))]+101​f(X∖A)≥(52​−5n9​)OPT,

and, for every δ∈(0,1]\delta \in (0,1]δ∈(0,1], every run of kkk iterations with accurate estimates has k<n2/δk < n^2/\deltak<n2/δ (fewer than 3n23n^23n2 for δ=1/3\delta = 1/3δ=1/3).

Milestones, in attack order

  • Lemma 2.2: E[g(A(p))]≥(1−p)g(∅)+p g(A)\mathbf{E}[g(A(p))] \ge (1-p)g(\emptyset) + p\,g(A)E[g(A(p))]≥(1−p)g(∅)+pg(A).
  • Display (∗): for three independently sampled sets, E[f(A1(p1)∪A2(p2)∪A3(p3))]≥∑I⊆{1,2,3}∏i∈Ipi∏i∉I(1−pi)f(⋃i∈IAi)\mathbf{E}[f(A_1(p_1) \cup A_2(p_2) \cup A_3(p_3))] \ge \sum_{I \subseteq \{1,2,3\}} \prod_{i\in I} p_i \prod_{i \notin I}(1-p_i) f(\bigcup_{i \in I} A_i)E[f(A1​(p1​)∪A2​(p2​)∪A3​(p3​))]≥∑I⊆{1,2,3}​∏i∈I​pi​∏i∈/I​(1−pi​)f(⋃i∈I​Ai​).
  • The increment identity Φ(A∪{x})−Φ(A)=δ ωA,δ(x)\Phi(A \cup \{x\}) - \Phi(A) = \delta\,\omega_{A,\delta}(x)Φ(A∪{x})−Φ(A)=δωA,δ​(x) for x∉Ax \notin Ax∈/A, and its removal counterpart.
  • 0≤Φ(A)≤OPT0 \le \Phi(A) \le OPT0≤Φ(A)≤OPT.
  • The terminal upper estimates E[f(R∪(B∩C))], E[f(R∩(B∪C))]≤E[f(R)]+2nOPT\mathbf{E}[f(R \cup (B\cap C))],\ \mathbf{E}[f(R \cap (B \cup C))] \le \mathbf{E}[f(R)] + \frac{2}{n}OPTE[f(R∪(B∩C))], E[f(R∩(B∪C))]≤E[f(R)]+n2​OPT.
  • The two lower bounds in 272727ths on the same two expectations.
  • The final chain E[f(R)]+19f(B)+2nOPT≥49OPT\mathbf{E}[f(R)] + \frac19 f(B) + \frac2n OPT \ge \frac49 OPTE[f(R)]+91​f(B)+n2​OPT≥94​OPT.

Significance

The 2/52/52/5 bound showed that local search on a smoothed objective, the multilinear extension restricted to points with two coordinate values, beats both uniform sampling and plain local search for non-monotone submodular maximization. The multilinear extension later became the standard tool for submodular maximization under constraints, through continuous greedy methods, contention resolution schemes and the analysis of randomized rounding. Display (∗) and Lemma 2.2 are the basic sampling inequalities for submodular functions and are reused throughout that literature.

The result is proved in the paper; it has been superseded in ratio by later algorithms reaching 1/21/21/2. To our knowledge none of it is machine-checked: Mathlib has no multilinear extension and no submodular maximization results. This mission produces a checked form of the analysis with every constant explicit: the o(1)o(1)o(1) as 95n\frac{9}{5n}5n9​, "polynomial time" as n2/δn^2/\deltan2/δ iterations, and the dependence on the accuracy of the sampled estimates as an explicit hypothesis.

Difficulty

The iteration bound and the increment identity are routine once the multilinear extension is set up. The substance lies in the lower bounds. The returned set RRR is random, so the comparison with the optimal set CCC cannot be made through a single local-optimality inequality as in deterministic local search; and approximate local optimality holds only for the smoothed marginals ωA,δ\omega_{A,\delta}ωA,δ​, which are averages over the random set, not for fff at any fixed set. Sampling inequalities such as (∗) are stated for independent samples of arbitrary, possibly overlapping sets, and their expectations are sums over products of subsets; the bookkeeping of such sums is the main formalization burden. The constants must also balance exactly: with δ=1/3\delta = 1/3δ=1/3 the 272727ths add up so that the 910/110\frac{9}{10}/\frac{1}{10}109​/101​ mixture yields 2/52/52/5. A different split or a different δ\deltaδ gives a different constant.

Formalization scope

The ground set is a Fintype X with decidable equality; sets are Finset X; fff is real valued, with nonnegativity and submodularity as hypotheses. The standing assumptions of the paper are made explicit: f≥0f \ge 0f≥0 (§3), value-oracle access (modelled by fff itself), n=∣X∣n = |X|n=∣X∣, and n≥1n \ge 1n≥1 (Nonempty X), so that 1n2\frac{1}{n^2}n21​ and 95n\frac{9}{5n}5n9​ are not Lean's junk value of division by zero. OPTOPTOPT is Finset.sup' over all subsets, which has no junk value. Every expectation over independently sampled sets is an exact finite sum: the multilinear extension for R(A,δ)\mathcal{R}(A,\delta)R(A,δ), and iterated sums over subsets for the several independent samples in Lemma 2.2 and (∗). Sampling probabilities carry 0≤p≤10 \le p \le 10≤p≤1.

The algorithm is a relation, not a choice: any element meeting the step-3 condition may be added, and removal is allowed only when no addition applies. The goal quantifies over every run from ∅\emptyset∅, every choice of accurate estimates, recomputed at every iteration, and every termination point. Accuracy is non-strict ("within ±\pm±"); the thresholds are strict. The thresholds and accuracy use OPTOPTOPT itself, as the proof does, although step 1 of the algorithm says an estimate of OPTOPTOPT is used. The sampling that produces the estimates and its "with high probability" are not modelled, and the goal is conditional on accurate estimates. The value of the δ′=−1\delta' = -1δ′=−1 branch is written as E[f(R(A,−1))]\mathbf{E}[f(\mathcal{R}(A,-1))]E[f(R(A,−1))], which equals f(X∖A)f(X \setminus A)f(X∖A).

A statement "every AAA satisfying the terminal conditions gives 2/5−9/(5n)2/5 - 9/(5n)2/5−9/(5n)" would be the final milestone plus arithmetic and is not the goal. The goal fixes the start at ∅\emptyset∅, the thresholds, the step order, the accuracy of every estimate, termination and the 0.9/0.10.9/0.10.9/0.1 mixture.

A complete development needs basic calculus of the multilinear extension: affinity in one coordinate, translation f(⋅∪D)f(\cdot \cup D)f(⋅∪D) and restriction f(⋅∩D)f(\cdot \cap D)f(⋅∩D), and splitting a sample into disjoint pieces. These lemmas are reusable for any work on submodular maximization, and contributions of them as separate theorems are welcome. The golden-ratio variant δ=δ′\delta = \delta'δ=δ′ (proof omitted in the paper), the tight example and the hardness results of §4 are out of scope.

Selected references

  • U. Feige, V. S. Mirrokni, J. Vondrák, Maximizing Non-Monotone Submodular Functions, SIAM J. Comput. 40(4):1133–1153, 2011. https://doi.org/10.1137/090779346
  • S. O. Gharan, J. Vondrák, Submodular Maximization by Simulated Annealing, SODA 2011. https://doi.org/10.1137/1.9781611973082.83
  • N. Buchbinder, M. Feldman, J. Naor, R. Schwartz, A Tight Linear Time (1/2)-Approximation for Unconstrained Submodular Maximization, SIAM J. Comput. 44(5):1384–1402, 2015. https://doi.org/10.1137/130929205
  • G. Calinescu, C. Chekuri, M. Pál, J. Vondrák, Maximizing a Monotone Submodular Function Subject to a Matroid Constraint, SIAM J. Comput. 40(6):1740–1766, 2011. https://doi.org/10.1137/080733991
14 thms3 active usersReviewed
🏆Completed
Combinatorics·Captain: mikedeng1

Three Partition Refinement Algorithms 1: Distinguishing-Prefix Refinement Finds Every Distinguishing Prefix Within m′ Refinement StepsResearch Paper

Motivation

Sorting a collection of strings in lexicographic order is a basic subroutine in compilers, string indexing, suffix sorting and the construction of tries. The classical method, due to Aho, Hopcroft and Ullman (The Design and Analysis of Computer Algorithms, 1974), is a multipass radix sort that processes the strings from the last position to the first and runs in time proportional to the total length mmm of the input plus the alphabet size kkk. Mehlhorn showed that a straightforward first-to-last scan of equal-length strings needs Ω(km)\Omega(km)Ω(km) time.

Paige and Tarjan (Three Partition Refinement Algorithms, SIAM J. Comput. 16(6):973–989, 1987, doi:10.1137/0216062) observed that most of the input is usually irrelevant to the sorted order: only the shortest prefix of each string that tells it apart from the others matters. Their algorithm separates the problem into two steps: (1) find this distinguishing prefix of every string, and (2) sort the distinguishing prefixes. The first step is carried out by partition refinement, the same technique that the paper then applies to the relational coarsest partition problem (the second mission of this series) and to double lexical ordering. This mission formalizes the correctness and the termination bound of the first step.

Setting

Fix k≥1k \ge 1k≥1 and the alphabet Σ={1,2,…,k}\Sigma = \{1, 2, \dots, k\}Σ={1,2,…,k}, together with an end marker 000 that is smaller than every symbol. A string is a finite sequence over Σ∪{0}\Sigma \cup \{0\}Σ∪{0}; λ\lambdaλ is the empty string, ∣x∣|x|∣x∣ is the length of xxx, and x(i)x(i)x(i) is its iii-th symbol (positions start at 111). A string α\alphaα is a prefix of yyy if y=αzy = \alpha zy=αz for some string zzz; every string is a prefix of itself. Σ∗0\Sigma^*0Σ∗0 is the set of strings over Σ\SigmaΣ followed by one 000.

The input is a multiset U={x1,…,xn}⊆Σ∗0U = \{x_1, \dots, x_n\} \subseteq \Sigma^*0U={x1​,…,xn​}⊆Σ∗0 with n≥1n \ge 1n≥1; repeated strings are allowed. The distinguishing prefix xi′x'_ixi′​ of xix_ixi​ in UUU is

  1. the shortest prefix of xix_ixi​ that is not a prefix of any other string of UUU, if xix_ixi​ occurs only once in UUU;
  2. xix_ixi​ itself, if xix_ixi​ occurs more than once.

Write m′=∑i=1n∣xi′∣m' = \sum_{i=1}^n |x'_i|m′=∑i=1n​∣xi′​∣.

For a string α\alphaα, the labeled block BαB_\alphaBα​ is the multiset of strings of UUU having α\alphaα as a prefix, and α\alphaα is its associated prefix. The block is finished if α=xi′\alpha = x'_iα=xi′​ for some iii, and unfinished otherwise. A state PPP is a set of labeled blocks. For an unfinished block Bα∈PB_\alpha \in PBα​∈P,

split(Bα,P)=(P∖{Bα})∪{ Bαu:u∈Σ∪{0}, ∃x∈Bα, x(∣α∣+1)=u }.\mathrm{split}(B_\alpha, P) = \bigl(P \setminus \{B_\alpha\}\bigr) \cup \{\, B_{\alpha u} : u \in \Sigma \cup \{0\},\ \exists x \in B_\alpha,\ x(|\alpha| + 1) = u \,\}.split(Bα​,P)=(P∖{Bα​})∪{Bαu​:u∈Σ∪{0}, ∃x∈Bα​, x(∣α∣+1)=u}.

The refinement algorithm starts from P0={Bλ}P_0 = \{B_\lambda\}P0​={Bλ​} and repeatedly applies the step Refine: pick any unfinished block Bα∈PB_\alpha \in PBα​∈P and replace PPP by split(Bα,P)\mathrm{split}(B_\alpha, P)split(Bα​,P). A run with KKK steps is a sequence P0,P1,…,PKP_0, P_1, \dots, P_KP0​,P1​,…,PK​ produced this way; any choice of unfinished block is allowed at every step.

Formalization targets

Goal: Theorem 1 with its explicit bound

For every run P0,…,PKP_0, \dots, P_KP0​,…,PK​ of the refinement algorithm,

K≤m′and(no block of PK is unfinished)  ⟹  PK={Bx1′,…,Bxn′}.K \le m' \qquad\text{and}\qquad \bigl(\text{no block of } P_K \text{ is unfinished}\bigr) \implies P_K = \{B_{x'_1}, \dots, B_{x'_n}\}.K≤m′and(no block of PK​ is unfinished)⟹PK​={Bx1′​​,…,Bxn′​​}.

The paper states Theorem 1 as "The algorithm terminates and is correct"; the bound K≤m′K \le m'K≤m′ is the explicit count proved in the last sentence of its proof (p. 975). The goal leaves the choice of unfinished block free, so it holds for every refinement order.

Milestones

  1. End markers (§2, p. 974). No string of UUU is a proper prefix of another string of UUU.
  2. Lemma 1 (p. 975). Along every run, every finished block BβB_\betaBβ​ is contained in a block of the current state, and there is Bα∈PjB_\alpha \in P_jBα​∈Pj​ with α\alphaα a prefix of β\betaβ.
  3. Proof of Theorem 1, second sentence. For every state of every run, ∑Bα∈Pj∣α∣≤m′\sum_{B_\alpha \in P_j} |\alpha| \le m'∑Bα​∈Pj​​∣α∣≤m′.
  4. Proof of Theorem 1, third sentence. Every Refine step strictly increases ∑Bα∈P∣α∣\sum_{B_\alpha \in P} |\alpha|∑Bα​∈P​∣α∣.

A supplementary item states the fact on p. 974 that motivates the two-step design: distinct strings have distinct distinguishing prefixes, and xi≤xjx_i \le x_jxi​≤xj​ iff xi′≤xj′x'_i \le x'_jxi′​≤xj′​ in lexicographic order.

Significance

Theorem 1 is the correctness half of the lexicographic sorting algorithm: once the distinguishing prefixes are known, sorting UUU reduces to sorting strings of total length m′m'm′, which is how the paper obtains its O(m′+k)O(m' + k)O(m′+k) time bound in place of O(m+k)O(m + k)O(m+k). The paper notes that under a natural probability model m′m'm′ is of order nlog⁡knn \log_k nnlogk​n, so the gain is large whenever the strings are long. The invariant of Lemma 1 is also the pattern on which the later partition refinement algorithms of the paper are built: a target partition is shown to refine every intermediate partition, and a potential function bounds the number of refinement steps.

The result has been proved since 1987. What this mission adds is a machine-checked account of the algorithm at the level of its abstract refinement steps, with the explicit bound m′m'm′ rather than an asymptotic statement, and with the standing assumptions of §2 (non-empty input, end markers) made explicit. No formal proof of this algorithm is present in Mathlib or on the platform.

Difficulty

The obvious argument, "each step makes the labels longer, and labels never grow past the distinguishing prefixes", needs two facts that are not immediate from the definitions. First, the labels of a reached state must form a partition of UUU into non-empty blocks whose labels are pairwise incomparable under the prefix order; this is an invariant of the algorithm, not part of the definition of a state, and it fails for arbitrary sets of labels. Second, the label of an unfinished block must be a proper prefix of every finished label below it, which rests on the end marker: without end markers a string can be a proper prefix of another, a block can be unfinished yet have no strictly longer children, and the algorithm can stall. Bounding the sum of label lengths by m′m'm′ is not a label-by-label comparison: a state can have fewer labels than strings, and the distinguishing prefixes of different strings can share a label as a common prefix.

Formalization scope

Strings are List (Fin (k + 1)), with 0 : Fin (k + 1) the end marker and prefix the Mathlib relation <+:. The multiset UUU is a family x : Fin n → List (Fin (k + 1)), so repetitions are distinct indices. The hypotheses 0 < n and EndMarked x (U⊆Σ∗0U \subseteq \Sigma^*0U⊆Σ∗0) are the paper's standing assumptions. A block is identified by its label, not by its set of members: two different labels can carry the same strings (for instance Bλ=B1B_\lambda = B_1Bλ​=B1​ when every string begins with 111), and they are different blocks. A state is a Finset of labels; a run is a map Fin (K + 1) → Finset (List (Fin (k + 1))) starting at {[]}. Positions are 1-based in the paper and 0-based in Lean, so x(∣α∣+1)x(|\alpha|+1)x(∣α∣+1) is (x i)[α.length]?. The distinguishing prefix follows the literal definition; for n=1n = 1n=1 it is the empty string.

The running times O(m′+k)O(m' + k)O(m′+k) and O(n+k)O(n + k)O(n+k) space, the implementation with a queue and a global index, and step two (sorting the prefixes via the refinement tree) are RAM-model statements and are not formalized. Only the explicit step count m′m'm′ is. A formalization in which split added all k+1k + 1k+1 children, added none, or allowed refining a finished block would change the theorem (the bound fails or the goal becomes vacuous); the definitions add exactly the children realized by some string of the block and refine only unfinished blocks.

Contributions welcome: proofs of the milestones, a general API for prefix-closed partition refinement on List, and a proof of the order-preservation item via Mathlib's List.Lex.

Selected references

  • R. Paige, R. E. Tarjan, Three Partition Refinement Algorithms, SIAM Journal on Computing 16(6):973–989, 1987. https://doi.org/10.1137/0216062
  • A. V. Aho, J. E. Hopcroft, J. D. Ullman, The Design and Analysis of Computer Algorithms, Addison-Wesley, 1974.
  • K. Mehlhorn, Data Structures and Algorithms 1: Sorting and Searching, Springer, 1984. https://doi.org/10.1007/978-3-642-69672-5
7 thms3 active usersReviewed
🏆Completed
CombinatoricsGraph TheoryOperations Research·Captain: mikedeng1

A Linear-Time Algorithm for Finding a Sparse k-Connected Spanning Subgraph of a k-Connected Graph 2: FOREST's Forests E_i Preserve Local Node-Connectivity up to i in a Simple GraphResearch Paper

Motivation

Given a kkk-connected graph, many connectivity algorithms run in time that grows with the number of edges ∣E∣|E|∣E∣. A sparse certificate is a spanning subgraph with only O(k∣V∣)O(k|V|)O(k∣V∣) edges that is still kkk-connected; computing one first and running the expensive algorithm on it replaces ∣E∣|E|∣E∣ by k∣V∣k|V|k∣V∣ in the bound. Finding a kkk-connected spanning subgraph with the minimum number of edges is NP-complete for every fixed k≥2k \ge 2k≥2 (Garey and Johnson, problem GT31), so the question is how cheaply a sparse, not necessarily minimum, certificate can be found.

Nagamochi and Ibaraki (Algorithmica 7 (1992) 583–596) answered this with a single linear-time scanning procedure, FOREST, which partitions the edges into classes E1,E2,…,E∣E∣E_1, E_2, \dots, E_{|E|}E1​,E2​,…,E∣E∣​. They showed that the prefix unions E1∪⋯∪EkE_1 \cup \dots \cup E_kE1​∪⋯∪Ek​ are certificates for edge-connectivity and, for simple graphs, for node-connectivity. The node-connectivity result is the subject of this mission; the edge-connectivity result is the preceding mission of this series.

Timeline:

  • 1980: Galil gives an algorithm testing κ(G)≥k\kappa(G) \ge kκ(G)≥k whose running time depends on ∣E∣|E|∣E∣ (SIAM J. Comput. 9).
  • Before 1992 (as cited on p. 583): Suzuki et al. give O(∣E∣)O(|E|)O(∣E∣)-time algorithms for sparse 2- and 3-node-connected spanning subgraphs; Nishizeki and Poljak find, for general kkk, a kkk-node-connected spanning subgraph with at most k(∣V∣−1)k(|V|-1)k(∣V∣−1) edges in O(∣V∣1/2∣E∣2)O(|V|^{1/2}|E|^2)O(∣V∣1/2∣E∣2) time.
  • 1992: Nagamochi and Ibaraki prove that FOREST, which runs in O(∣V∣+∣E∣)O(|V| + |E|)O(∣V∣+∣E∣) time, yields a kkk-node-connected spanning subgraph GkG_kGk​ of every simple kkk-node-connected graph, with ∣E(Gk)∣≤k∣V∣−k(k+1)/2|E(G_k)| \le k|V| - k(k+1)/2∣E(Gk​)∣≤k∣V∣−k(k+1)/2. Their Theorem 3.1 states a stronger, local form.
  • 1993: Cheriyan, Kao and Thurimella isolate "scan-first search" as the general principle behind such certificates (SIAM J. Comput. 22 (1993)).

Setting

A graph G=(V,E)G = (V, E)G=(V,E) has a finite node set VVV with ∣V∣≥2|V| \ge 2∣V∣≥2 and a finite edge set EEE; each edge has an unordered pair of two distinct end nodes. In this mission the graph is simple: no two edges have the same end nodes. For F⊆EF \subseteq EF⊆E, (V,F)(V, F)(V,F) is the spanning subgraph with edge set FFF.

The local node-connectivity κ(x,y;H)\kappa(x, y; H)κ(x,y;H) of nodes x,yx, yx,y in a graph HHH on VVV is ∣V∣−1|V| - 1∣V∣−1 if xxx and yyy are adjacent in HHH, and otherwise the minimum size of a node set W⊆V−{x,y}W \subseteq V - \{x, y\}W⊆V−{x,y} whose deletion leaves no xxx–yyy path. The node connectivity is κ(G)=min⁡x,yκ(x,y;G)\kappa(G) = \min_{x, y} \kappa(x, y; G)κ(G)=minx,y​κ(x,y;G).

Procedure FOREST keeps a label r(v)≥0r(v) \ge 0r(v)≥0 on each node, initially 000. While some node is unscanned, it chooses an unscanned node xxx of largest label; for each unscanned edge e=(x,y)e = (x, y)e=(x,y) it puts eee into the class Er(y)+1E_{r(y)+1}Er(y)+1​, increases r(x)r(x)r(x) by one if r(x)=r(y)r(x) = r(y)r(x)=r(y), and increases r(y)r(y)r(y) by one; then it marks xxx scanned. Ties are broken arbitrarily. The time instants are the states between these elementary operations; Ei∗E^*_iEi∗​ denotes the class iii at an instant, and EiE_iEi​ its final value. Put

Gi=(V, E1∪E2∪⋯∪Ei).G_i = (V,\ E_1 \cup E_2 \cup \dots \cup E_i).Gi​=(V, E1​∪E2​∪⋯∪Ei​).

Formalization targets

Goal: Theorem 3.1

For a simple graph GGG and the classes of any completed run of FOREST, for 1≤i≤∣E∣1 \le i \le |E|1≤i≤∣E∣,

κ(x,y;Gi) ≥ min⁡{κ(x,y;G), i}for any x,y∈V.(3.1)\kappa(x, y; G_i) \ \ge\ \min\{\kappa(x, y; G),\ i\} \qquad \text{for any } x, y \in V. \tag{3.1}κ(x,y;Gi​) ≥ min{κ(x,y;G), i}for any x,y∈V.(3.1)

The statement is local: it holds pair by pair, not only for the global minimum, and for every tie-breaking of the procedure.

Milestones, in the order the proof uses them

  1. Lemma 2.2: at every instant, a node vvv meets EiE_iEi​ exactly for i=1,…,r(v)i = 1, \dots, r(v)i=1,…,r(v).
  2. Lemma 2.4(b): at every instant, a uuu–vvv path in Ej∗E^*_jEj∗​ yields uuu–vvv paths in every Ei∗E^*_iEi∗​, i<ji < ji<j.
  3. In-degree at most one (§2, p. 588): orienting each edge from the earlier-scanned to the later-scanned end, every node has at most one entering arc in each class.
  4. Lemma 3.1: if an xxx–yyy path of Ej∗E^*_jEj∗​ has the form x,u1,…,uk=w,yx, u_1, \dots, u_k = w, yx,u1​,…,uk​=w,y with k=1k = 1k=1 or u1u_1u1​ scanned before www, then any www–xxx and www–yyy paths in Ei∗E^*_iEi∗​ (i<ji < ji<j) share a node other than www.
  5. Lemma 3.2: for a node cut set W={w1,…,wi}W = \{w_1, \dots, w_i\}W={w1​,…,wi​} of Gi+1G_{i+1}Gi+1​ (in scan order) separating a component XXX from the rest YYY, immediately after wtw_twt​ is scanned every XXX–YYY path of Et∗E^*_tEt∗​ passes through wtw_twt​, and Ej∗E^*_jEj∗​ has no XXX–YYY path for t+1≤j≤i+1t + 1 \le j \le i + 1t+1≤j≤i+1.

A companion item states the paper's announcement in §3: GkG_kGk​ is kkk-node-connected for every 1≤k≤κ(G)1 \le k \le \kappa(G)1≤k≤κ(G).

Significance

Theorem 3.1 at i=ki = ki=k shows that the first kkk classes of FOREST form a kkk-node-connected spanning subgraph whenever GGG is, and the edge-count analysis of the companion mission bounds its size by k∣V∣−k(k+1)/2k|V| - k(k+1)/2k∣V∣−k(k+1)/2. Since FOREST runs in linear time, any algorithm testing κ(G)≥k\kappa(G) \ge kκ(G)≥k can be run on GkG_kGk​ instead of GGG; the paper uses this to improve the bound for testing κ(G)≥k\kappa(G) \ge kκ(G)≥k from O(max⁡{k2∣V∣1/2,k∣V∣}∣E∣)O(\max\{k^2|V|^{1/2}, k|V|\}|E|)O(max{k2∣V∣1/2,k∣V∣}∣E∣) to O(max⁡{k3∣V∣3/2,k2∣V∣2})O(\max\{k^3|V|^{3/2}, k^2|V|^2\})O(max{k3∣V∣3/2,k2∣V∣2}), and similar gains for computing the number of node-disjoint paths between two nodes. The local form (3.1) is what makes the sss–ttt applications possible.

The result is proved on paper. As far as a search of the platform shows, neither FOREST nor local node-connectivity has a machine-checked treatment there; Mathlib has no notion of vertex connectivity of a pair of nodes. This mission produces a formal model of FOREST as a nondeterministic transition system and the statements needed to verify the paper's proof step by step.

Difficulty

For edge-connectivity the analogous statement follows from a general principle: any sequence of maximal spanning forests, each taken in what remains of the graph, preserves local edge-connectivity up to its length. The obvious attempt is to prove (3.1) the same way, from the fact that each EiE_iEi​ is a maximal spanning forest of what the earlier classes leave. The paper gives no such argument for node-connectivity: its proof uses the specific scan order of FOREST in an essential way, through the orientation of edges from earlier- to later-scanned nodes and the in-degree bound of that orientation (Lemmas 3.1 and 3.2). The argument tracks, for a hypothetical node cut WWW of size iii in Gi+1G_{i+1}Gi+1​, the classes at the moments the nodes of WWW are scanned, which requires reasoning about intermediate states of the algorithm and about paths in several classes at once. None of this reduces to a static property of the output partition.

Formalization scope

  • Graphs. A node type V and an edge type E, both finite, with ends : E → Sym2 V; loop-freeness is ∀ e, ¬ (ends e).IsDiag and simplicity is Function.Injective ends. The standing assumptions of p. 583 and p. 589 (∣V∣≥2|V| \ge 2∣V∣≥2, no self-loop, a simple graph when node-connectivity is discussed) appear as hypotheses; §2 items are stated for loopless graphs, as on the page.
  • Connectivity. κ(x,y;(V,F))\kappa(x, y; (V, F))κ(x,y;(V,F)) is valued in N∞\mathbb N_\inftyN∞​: ∣V∣−1|V| - 1∣V∣−1 on adjacent pairs, the minimum node cut otherwise, and ⊤\top⊤ when x=yx = yx=y, where (3.1) holds trivially.
  • FOREST. A nondeterministic step relation with three steps (select, scan, finish), a run of length KKK from the initial state, and completion when every node is scanned. Every tie-breaking is allowed, so the theorems quantify over all completed runs. A time instant is a state of a run; "scanned before" compares positions in the run's selection order; "immediately after wtw_twt​ has been scanned" is the state right after the finish step of wtw_twt​.
  • Excluded. The running time "O(∣V∣+∣E∣)O(|V| + |E|)O(∣V∣+∣E∣)" and everything in §4 (the connectivity-testing algorithms and their bounds) are not stated: the paper fixes no machine model. The edge bounds on ∣Ei∣|E_i|∣Ei​∣ belong to the edge-connectivity mission.
  • Ruled out. Stating (3.1) for an arbitrary partition into maximal spanning forests, or reading the classes Ej∗E^*_jEj∗​ of Lemmas 3.1–3.2 off the final state, would state a different theorem from the one the paper proves; the classes are those of a run of FOREST at the instant the page specifies.
  • Infrastructure. Reachability avoiding a node set, walks and paths in SimpleGraph, and invariants of the FOREST transition system. The definitions duplicate those of the edge-connectivity mission by design and are candidates for a shared layer. Contributions of general lemmas about the run (label invariants, monotonicity of classes along a run) are welcome.

Selected references

  • H. Nagamochi, T. Ibaraki, A linear-time algorithm for finding a sparse kkk-connected spanning subgraph of a kkk-connected graph, Algorithmica 7 (1992), 583–596. https://doi.org/10.1007/BF01758778
  • Z. Galil, Finding the vertex connectivity of graphs, SIAM J. Comput. 9 (1980), 197–199. https://doi.org/10.1137/0209016
  • J. Cheriyan, M.-Y. Kao, R. Thurimella, Scan-first search and sparse certificates: an improved parallel algorithm for kkk-vertex connectivity, SIAM J. Comput. 22 (1993), 157–174. https://doi.org/10.1137/0222013
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, Freeman, 1979.
9 thms3 active usersReviewed
🏆Completed
CombinatoricsGraph TheoryOperations Research·Captain: mikedeng1

A Linear-Time Algorithm for Finding a Sparse k-Connected Spanning Subgraph of a k-Connected Graph 1: The FOREST Decomposition Preserves Local Edge-ConnectivityResearch Paper

Motivation

Many graph algorithms for connectivity questions run in time proportional to the number of edges. When the question is only whether a graph is kkk-edge-connected, or what its local edge-connectivities are up to a threshold kkk, most edges are irrelevant: a spanning subgraph with O(k∣V∣)O(k|V|)O(k∣V∣) edges already carries the answer. Such a subgraph is called a sparse certificate. Computing one first and running the expensive algorithm on it replaces ∣E∣|E|∣E∣ by k∣V∣k|V|k∣V∣ in the running time of connectivity testing, of Matula-type edge-connectivity algorithms, and of sss–ttt flow computations used as connectivity oracles.

Nagamochi and Ibaraki (Algorithmica 7, 1992) gave a procedure, FOREST, that computes such a certificate for edge-connectivity and, on simple graphs, for node-connectivity, with a single graph search. The same partition of the edges into forests is the engine of their deterministic minimum-cut algorithm for multigraphs (SIAM J. Discrete Math. 5, 1992), which later became the maximum-adjacency ordering of the Stoer–Wagner minimum-cut algorithm (J. ACM 44, 1997).

Timeline:

  • 1927: Menger identifies the minimum number of edges separating two nodes with the maximum number of edge-disjoint paths between them.
  • 1992: Nagamochi and Ibaraki publish Procedure FOREST and prove that its iii-th prefix preserves local edge-connectivity up to iii in multigraphs, and local node-connectivity up to iii in simple graphs.
  • 1993: Cheriyan, Kao and Thurimella (SIAM J. Comput. 22) obtain sparse certificates by scan-first search; Frank, Ibaraki and Nagamochi (J. Graph Theory 17) give a shorter proof for the node-connectivity case.
  • 1994: Nishizeki and Poljak (Discrete Appl. Math. 55) publish the forest-decomposition lemma (Lemma 2.1 below), found independently.

Setting

A graph G=(V,E)G = (V, E)G=(V,E) is finite and undirected, has ∣V∣≥2|V| \ge 2∣V∣≥2 nodes, may have multiple edges (several edges with the same pair of end nodes), and has no self-loop. It is simple if no two edges have the same end nodes. For F⊆EF \subseteq EF⊆E, (V,F)(V, F)(V,F) is the spanning subgraph with edge set FFF. It is a forest if it has no cycle; two parallel edges form a cycle. It is a maximal spanning forest in (V,H)(V, H)(V,H), for F⊆HF \subseteq HF⊆H, if adding any edge of H∖FH \setminus FH∖F to FFF creates a cycle.

The local edge-connectivity λ(x,y;H)\lambda(x, y; H)λ(x,y;H) is the minimum number of edges of HHH whose removal leaves no path from xxx to yyy; parallel edges count separately. It is ∞\infty∞ when x=yx = yx=y.

Procedure FOREST keeps a label r(v)∈Nr(v) \in \mathbb{N}r(v)∈N on every node, initially 000, and classes E1,E2,…,E∣E∣E_1, E_2, \dots, E_{|E|}E1​,E2​,…,E∣E∣​, initially empty. While an unscanned node exists, it picks an unscanned node xxx of largest label; for every unscanned edge eee from xxx to a node yyy, it puts eee into Er(y)+1E_{r(y)+1}Er(y)+1​, increases r(x)r(x)r(x) by one if r(x)=r(y)r(x) = r(y)r(x)=r(y), increases r(y)r(y)r(y) by one, and marks eee scanned; then it marks xxx scanned. Ties among nodes and the order of edges are free. On termination, Gi=(V,E1∪⋯∪Ei)G_i = (V, E_1 \cup \cdots \cup E_i)Gi​=(V,E1​∪⋯∪Ei​).

Formalization targets

Goal: Theorem 2.1 (pp. 588–589), without the running time

For every graph GGG and every completed execution of FOREST on GGG:

  1. every edge lies in exactly one class EiE_iEi​, 1≤i≤∣E∣1 \le i \le |E|1≤i≤∣E∣;
  2. for i=1,…,∣E∣i = 1, \dots, |E|i=1,…,∣E∣,
λ(x,y;Gi)≥min⁡{λ(x,y;G), i}for all x,y∈V;(2.1)\lambda(x, y; G_i) \ge \min\{\lambda(x, y; G),\ i\} \qquad \text{for all } x, y \in V; \tag{2.1}λ(x,y;Gi​)≥min{λ(x,y;G), i}for all x,y∈V;(2.1)
  1. ∣Ei∣≤∣V∣−1|E_i| \le |V| - 1∣Ei​∣≤∣V∣−1 for all iii;
  2. if GGG is simple, ∣Ei∣≤∣V∣−i|E_i| \le |V| - i∣Ei​∣≤∣V∣−i for i≤∣V∣−1i \le |V| - 1i≤∣V∣−1 and Ei=∅E_i = \emptysetEi​=∅ for i≥∣V∣i \ge |V|i≥∣V∣.

Milestones

  • Lemma 2.2 (p. 587): during the execution, a node vvv has incident edges in exactly the classes E1,…,Er(v)E_1, \dots, E_{r(v)}E1​,…,Er(v)​.
  • Lemma 2.3 (p. 587): each (V,Ei)(V, E_i)(V,Ei​) is a forest at every instant.
  • Lemma 2.4 (p. 588): (a) an edge (u,v)(u, v)(u,v) added to EiE_iEi​ has its end nodes joined by a path in Ei−1E_{i-1}Ei−1​; (b) a path in EjE_jEj​ between uuu and vvv yields a path in every EiE_iEi​, i<ji < ji<j.
  • Lemma 2.5 (p. 588): each output (V,Ei)(V, E_i)(V,Ei​) is a maximal spanning forest in G−E1∪⋯∪Ei−1G - E_1 \cup \cdots \cup E_{i-1}G−E1​∪⋯∪Ei−1​.
  • Lemma 2.1 (p. 584): any sequence of successive maximal spanning forests satisfies (2.1).

Companions

  • the sparse certificate (p. 589): if λ(x,y;G)≥k\lambda(x, y; G) \ge kλ(x,y;G)≥k for all x,yx, yx,y, then GkG_kGk​ is kkk-edge-connected with ∣E(Gk)∣≤k(∣V∣−1)|E(G_k)| \le k(|V| - 1)∣E(Gk​)∣≤k(∣V∣−1), and ∣E(Gk)∣≤k∣V∣−k(k+1)/2|E(G_k)| \le k|V| - k(k+1)/2∣E(Gk​)∣≤k∣V∣−k(k+1)/2 for simple GGG;
  • Lemma 2.6 (p. 589): for k≤δ(G)k \le \delta(G)k≤δ(G), GkG_kGk​ has a node of degree exactly kkk.

Significance

(2.1) says that one search produces, for every threshold kkk at once, a subgraph with at most k(∣V∣−1)k(|V|-1)k(∣V∣−1) edges that keeps every local edge-connectivity up to kkk. Any algorithm whose running time grows with ∣E∣|E|∣E∣ can then be run on GkG_kGk​ in place of GGG; §4 of the paper uses this to speed up kkk-connectivity tests and the computation of local connectivities. The same forest partition is the structural fact behind the Nagamochi–Ibaraki and Stoer–Wagner minimum-cut algorithms. Lemma 2.6 shows that the certificate is tight: its edge-connectivity is exactly kkk when λ(G)≥k\lambda(G) \ge kλ(G)≥k.

The results are proved in the paper. No machine-checked proof of them is known. Formalizing them means formalizing a graph search with free tie-breaking as a transition system, reasoning about invariants of all its executions, and proving a cut-counting statement for multigraphs. Mathlib's connectivity notions, such as SimpleGraph.IsEdgeReachable, do not see parallel edges, so the multigraph cut theory here is new.

Difficulty

Lemma 2.1 is a short cut argument once maximality is available. The difficulty is showing that FOREST, which assigns each edge to a class by looking only at the label of one end node, produces maximal forests in the successive residual graphs (Lemma 2.5). The obvious invariant, that the class of an edge is the first forest it does not close a cycle in, is not what line 7 computes. The label r(y)r(y)r(y) records only which classes touch yyy, not which component of each class contains yyy. The paper's argument needs Lemma 2.4(a): at the moment an edge is added to EiE_iEi​ its ends already lie in one tree of Ei−1E_{i-1}Ei−1​. That relies on the choice of the unscanned node of largest label, on the order of lines 8 and 9, and on an argument about the scan order of tree roots.

Formalization scope

  • Graphs. A graph is a finite node type V with ∣V∣≥2|V| \ge 2∣V∣≥2, a finite edge type E, and ends : E → Sym2 V with no diagonal value (no self-loops). Parallel edges are distinct elements of E. Simplicity is injectivity of ends, and edge subsets are Finset E. A forest is an edge set in which every edge is a bridge, a condition that sees parallel edges.
  • Connectivity. λ\lambdaλ is an infimum in ℕ∞ over separating edge sets, ∞\infty∞ at x=yx = yx=y. (2.1) is kept "for all x,yx, yx,y", as printed.
  • FOREST. FOREST is a nondeterministic step relation (select, scan, finish) on explicit states: labels, class index per edge (000 = unscanned), scanned nodes, current node and selection order. The theorems quantify over every run from the initial state, so no tie-breaking rule is fixed. "At some time instant" is a state of the run; "upon completion" is a run whose last state has every node scanned. Lemma 2.2 is stated at every state, not only after a scan block (the other steps change neither labels nor classes). Lemma 2.4(a) assumes i≥2i \ge 2i≥2, since E0E_0E0​ does not exist.
  • Exclusions. Theorem 2.1's clause "is found in O(∣V∣+∣E∣)O(|V| + |E|)O(∣V∣+∣E∣) time", the bucket implementation, and the time bound of the certificate are not formalized: the paper fixes no machine model. The goal consists of the structural conclusions only. The bound printed "if GGG is multiple" is stated for every loopless graph.
  • Non-triviality. The goal is about the classes of a run of FOREST. An arbitrary partition of EEE into maximal spanning forests is Lemma 2.1's hypothesis, not a formalization of Theorem 2.1. A statement in which the classes are unconstrained variables, or in which the run hypotheses cannot be met, would be trivial. A separate sanity file checks that a complete run on the triangle K3K_3K3​ exists and attains ∣E1∣=∣V∣−1|E_1| = |V| - 1∣E1​∣=∣V∣−1, ∣E2∣=∣V∣−2|E_2| = |V| - 2∣E2​∣=∣V∣−2.
  • Welcome contributions. Useful reusable infrastructure includes:
    • a cut and Menger layer for finite multigraphs;
    • forest and bridge lemmas for edge-indexed graphs;
    • invariant-style reasoning over runs.

Selected references

  • H. Nagamochi, T. Ibaraki, A linear-time algorithm for finding a sparse kkk-connected spanning subgraph of a kkk-connected graph, Algorithmica 7 (1992) 583–596. https://doi.org/10.1007/BF01758778
  • H. Nagamochi, T. Ibaraki, Computing edge-connectivity in multigraphs and capacitated graphs, SIAM J. Discrete Math. 5 (1992) 54–66. https://doi.org/10.1137/0405004
  • T. Nishizeki, S. Poljak, kkk-connectivity and decomposition of graphs into forests, Discrete Appl. Math. 55 (1994) 295–301. https://doi.org/10.1016/0166-218X(94)90014-0
  • J. Cheriyan, M.-Y. Kao, R. Thurimella, Scan-first search and sparse certificates: an improved parallel algorithm for kkk-vertex connectivity, SIAM J. Comput. 22 (1993) 157–174. https://doi.org/10.1137/0222013
  • A. Frank, T. Ibaraki, H. Nagamochi, On sparse subgraphs preserving connectivity properties, J. Graph Theory 17 (1993) 275–281. https://doi.org/10.1002/jgt.3190170302
  • M. Stoer, F. Wagner, A simple min-cut algorithm, J. ACM 44 (1997) 585–591. https://doi.org/10.1145/263867.263872
10 thms3 active usersReviewed
🏆Completed
Control TheoryOperations Research·Captain: mikedeng1

Supervisory Control of a Class of Discrete Event Processes I: Minimally Restrictive Supervisors Exist iff the Supremal Controllable Legal Sublanguage Contains the Minimal Acceptable LanguageResearch Paper

Motivation

Manufacturing cells, communication protocols, traffic systems and database transaction managers are naturally described not by differential equations but by sequences of discrete events: a machine starts, a part arrives, a message is lost. Ramadge and Wonham's 1987 paper (SIAM J. Control Optim. 25(1)) set up a control theory for such systems in which the plant is an automaton, the controller is another automaton that may disable some events, and specifications are formal languages. The framework, now called supervisory control theory or the Ramadge–Wonham framework, is the standard model for the logical control of discrete event systems and underlies the textbook treatment in Cassandras and Lafortune (Introduction to Discrete Event Systems, 2008) and Wonham and Cai (Supervisory Control of Discrete-Event Systems, 2019).

The question the paper answers is the basic synthesis question of the theory: given a plant, a set of legal behaviours and a set of minimally acceptable behaviours, when does a controller exist that keeps the plant legal, achieves at least the acceptable behaviour, and never deadlocks, and what is the least restrictive such controller?

Setting

A generator is G=(Q,Σ,δ,q0,Qm)\mathcal G = (Q, \Sigma, \delta, q_0, Q_m)G=(Q,Σ,δ,q0​,Qm​) with a state set QQQ, a finite alphabet Σ\SigmaΣ of events, a partial transition function δ:Σ×Q→Q\delta : \Sigma \times Q \to Qδ:Σ×Q→Q, an initial state q0q_0q0​ and marker states Qm⊆QQ_m \subseteq QQm​⊆Q. Extending δ\deltaδ to strings, the generated language L(G)L(\mathcal G)L(G) is the set of strings www for which δ(w,q0)\delta(w, q_0)δ(w,q0​) is defined, and the marked language Lm(G)L_m(\mathcal G)Lm​(G) is the subset of those that end in QmQ_mQm​. The closure Kˉ\bar KKˉ of a language KKK is its set of prefixes; KKK is closed if K=KˉK = \bar KK=Kˉ. Throughout, G\mathcal GG is assumed trim in the sense L(G)=Lˉm(G)L(\mathcal G) = \bar L_m(\mathcal G)L(G)=Lˉm​(G): every generated string can be completed to a marked one.

The alphabet is split into controllable events Σc\Sigma_cΣc​ and uncontrollable events Σu=Σ−Σc\Sigma_u = \Sigma - \Sigma_cΣu​=Σ−Σc​. A supervisor S=(S,ϕ)\mathcal S = (S, \phi)S=(S,ϕ) consists of a deterministic, accessible automaton S=(X,Σ,ξ,x0,Xm)S = (X, \Sigma, \xi, x_0, X_m)S=(X,Σ,ξ,x0​,Xm​), whose state set XXX may be infinite, and a map ϕ\phiϕ assigning to each state xxx the set of controllable events it enables; uncontrollable events are always enabled. The closed loop S/G\mathcal S/\mathcal GS/G runs SSS and G\mathcal GG in lockstep: an event occurs when the plant can execute it, the supervisor enables it, and the supervisor's automaton can follow it. This defines the languages L(S/G)L(\mathcal S/\mathcal G)L(S/G), Lm(S/G)L_m(\mathcal S/\mathcal G)Lm​(S/G) and the controlled language Lc(S/G)=L(S/G)∩Lm(G)L_c(\mathcal S/\mathcal G) = L(\mathcal S/\mathcal G) \cap L_m(\mathcal G)Lc​(S/G)=L(S/G)∩Lm​(G).

S\mathcal SS is complete if its automaton never refuses an event that the plant can execute and ϕ\phiϕ enables; it is proper if it is complete and

Lˉm(S/G)=Lˉc(S/G)=L(S/G),\bar L_m(\mathcal S/\mathcal G) = \bar L_c(\mathcal S/\mathcal G) = L(\mathcal S/\mathcal G),Lˉm​(S/G)=Lˉc​(S/G)=L(S/G),

that is, every closed-loop string can be completed to a marked task. A language KKK is controllable if K⊆L(G)K \subseteq L(\mathcal G)K⊆L(G) and KˉΣu∩L(G)⊆Kˉ\bar K \Sigma_u \cap L(\mathcal G) \subseteq \bar KKˉΣu​∩L(G)⊆Kˉ. For L⊆L(G)L \subseteq L(\mathcal G)L⊆L(G), CG(L)\mathbf C_{\mathcal G}(L)CG​(L) is the class of controllable sublanguages of LLL and FG(L)\mathbf F_{\mathcal G}(L)FG​(L) the class of sublanguages KKK of LLL with K=Kˉ∩Lm(G)K = \bar K \cap L_m(\mathcal G)K=Kˉ∩Lm​(G).

Given ∅≠La⊆Lg⊆Lm(G)\emptyset \neq L_a \subseteq L_g \subseteq L_m(\mathcal G)∅=La​⊆Lg​⊆Lm​(G), the Supervisory Marking Problem (SMP) asks for a proper S\mathcal SS with La⊆Lm(S/G)⊆LgL_a \subseteq L_m(\mathcal S/\mathcal G) \subseteq L_gLa​⊆Lm​(S/G)⊆Lg​, and the Supervisory Control Problem (SCP) for a proper S\mathcal SS with La⊆Lc(S/G)⊆LgL_a \subseteq L_c(\mathcal S/\mathcal G) \subseteq L_gLa​⊆Lc​(S/G)⊆Lg​.

Formalization targets

Goal: Theorem 7.1 (pp. 218–219)

SMP solvable  ⟺  sup⁡CG(Lg)⊇La,SCP solvable  ⟺  sup⁡{CG(Lg)∩FG(Lg)}⊇La,\text{SMP solvable} \iff \sup \mathbf C_{\mathcal G}(L_g) \supseteq L_a, \qquad \text{SCP solvable} \iff \sup\{\mathbf C_{\mathcal G}(L_g) \cap \mathbf F_{\mathcal G}(L_g)\} \supseteq L_a,SMP solvable⟺supCG​(Lg​)⊇La​,SCP solvable⟺sup{CG​(Lg​)∩FG​(Lg​)}⊇La​,

and in each case the solving supervisor can be taken minimally restrictive: its marked (respectively controlled) language equals the supremal element and contains that of every proper supervisor whose language lies in LgL_gLg​.

Milestones

  1. Proposition 4.1 (i), (ii): marking is independent of control; every K⊆Lm(G)K \subseteq L_m(\mathcal G)K⊆Lm​(G), or K⊆L∩Lm(G)K \subseteq L \cap L_m(\mathcal G)K⊆L∩Lm​(G) for an achievable closed LLL, is the marked language of a complete supervisor.
  2. Proposition 5.1: a complete supervisor realizes (Lm,Lc,L)=(K1,K2,K3)(L_m, L_c, L) = (K_1, K_2, K_3)(Lm​,Lc​,L)=(K1​,K2​,K3​) iff K1⊆K2K_1 \subseteq K_2K1​⊆K2​, K2=K3∩Lm(G)K_2 = K_3 \cap L_m(\mathcal G)K2​=K3​∩Lm​(G) and K3K_3K3​ is closed and controllable.
  3. Theorem 6.1 (i), (ii): a proper supervisor with Lm(S/G)=KL_m(\mathcal S/\mathcal G) = KLm​(S/G)=K exists iff KKK is controllable; one with Lc(S/G)=KL_c(\mathcal S/\mathcal G) = KLc​(S/G)=K exists iff KKK is controllable and Lm(G)L_m(\mathcal G)Lm​(G)-closed.
  4. Proposition 7.1: CG(L)\mathbf C_{\mathcal G}(L)CG​(L) and FG(L)\mathbf F_{\mathcal G}(L)FG​(L) contain ∅\emptyset∅ and are closed under arbitrary unions.
  5. The supremal elements (p. 218): sup⁡CG(L)\sup \mathbf C_{\mathcal G}(L)supCG​(L), sup⁡FG(L)\sup \mathbf F_{\mathcal G}(L)supFG​(L) and sup⁡{CG(L)∩FG(L)}\sup\{\mathbf C_{\mathcal G}(L) \cap \mathbf F_{\mathcal G}(L)\}sup{CG​(L)∩FG​(L)} belong to their classes.

Significance

Theorem 7.1 reduces the existence of a correct, non-blocking controller to a single language inclusion, and identifies the supremal controllable sublanguage as the behaviour of the least restrictive solution. That object is the backbone of the later theory: modular and decentralized control, control under partial observation, and the computational results that sup⁡CG(L)\sup \mathbf C_{\mathcal G}(L)supCG​(L) is regular and computable when G\mathcal GG is finite and LLL regular all start from it.

The result is proved in the paper and has been taught for decades; it is not open. Mathlib has no model of generators with partial transitions, supervisors or controllability (its DFA has a total transition function and a single accepted language), and no formalization of this theory exists on the platform. This mission provides one: a reusable Lean model of generators with partial transitions, supervisors with possibly infinite state, closed loops and controllability, together with the paper's existence theorems stated against it. A second mission on the same paper (quotients of supervisors, Theorem 10.1) builds on the same objects.

Difficulty

The combinatorial core of Proposition 7.1 is a short prefix-closure computation. The work lies in the constructions: to show existence, a supervisor must be built for an arbitrary controllable language, which in general is not regular, so the supervisor needs an infinite state set (for example strings of the target language) together with a proof that the closed loop generates exactly the intended language, is complete, and is proper. The converse directions require relating the closed-loop run to separate runs of the plant and of the supervisor's automaton. The naive shortcut of reading the "sup" as an arbitrary member of CG(Lg)\mathbf C_{\mathcal G}(L_g)CG​(Lg​) containing LaL_aLa​ skips the content of the supremal-element milestone; the goal is stated with the actual union.

Formalization scope

  • The alphabet is a type α with [Fintype α]; strings are List α, the empty string is [], and sσs\sigmasσ is s ++ [σ]. Languages are Set (List α); the closure is pre K = {s | ∃ t, s ++ t ∈ K}.
  • A generator is a structure with a state type Q : Type, a partial transition δ : α → Q → Option Q, an initial state and a marker set. No finiteness is assumed on states, of the plant or of the supervisor. Supervisors and generators live in Type 1.
  • ϕ\phiϕ maps states to Ec → Bool; an event outside Σc\Sigma_cΣc​ is enabled by definition, so uncontrollable events cannot be disabled.
  • The closed loop is the product run from (x0,q0)(x_0, q_0)(x0​,q0​); the accessible part in the paper's display (2.1) changes no language and is not built.
  • Standing assumptions carried by every theorem: Σ\SigmaΣ finite, L(G)=Lˉm(G)L(\mathcal G) = \bar L_m(\mathcal G)L(G)=Lˉm​(G), supervisor automata accessible (as a hypothesis on every input supervisor and a conjunct of every "there exists a supervisor").
  • sup⁡\supsup is sSup in the complete lattice Set (List α).
  • A formalization in which the closed loop ignores ϕ\phiϕ or the plant, in which completeness is dropped from properness, or in which the supervisor is restricted to finitely many states, is not the paper's theorem and is ruled out by the definitions above.

Welcome contributions: the basic run lemmas (closed-loop runs project to plant and supervisor runs), the string-state supervisor construction, and proofs of the milestones in the given order.

Selected references

  • P. J. Ramadge and W. M. Wonham, Supervisory Control of a Class of Discrete Event Processes, SIAM J. Control Optim. 25(1):206–230, 1987. https://doi.org/10.1137/0325013
  • W. M. Wonham and P. J. Ramadge, On the Supremal Controllable Sublanguage of a Given Language, SIAM J. Control Optim. 25(3):637–659, 1987. https://doi.org/10.1137/0325036
  • C. G. Cassandras and S. Lafortune, Introduction to Discrete Event Systems, 2nd ed., Springer, 2008. https://doi.org/10.1007/978-0-387-68612-7
  • W. M. Wonham and K. Cai, Supervisory Control of Discrete-Event Systems, Springer, 2019. https://doi.org/10.1007/978-3-319-77452-7
12 thms3 active usersReviewed
PreviousPage 3 of 10Next

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