Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Collections

Online Algorithms

Competitive analysis: paging, k-server, metrical task systems, online primal-dual, secretary problems, and online matching.

18 open missions

Missions

1–18 of 18
OpenCompletedAll
Operations ResearchTheoretical Computer Science·Captain: Shuze Chen

The k-Server ConjectureOpen Problem

Motivation

The kkk-server problem was introduced by Manasse, McGeoch, and Sleator (STOC 1988 / J. Algorithms 1990) as a common generalization of paging, weighted caching, and related sequential decision problems, and their kkk-server conjecture has since become the central open question of competitive analysis. The conjecture asserts that a single ratio — exactly kkk — governs deterministic online server management on every metric space.

Timeline

  • 1985. Sleator and Tarjan introduce competitive analysis — an online algorithm judged against the offline optimum on every input — for list update and paging, and ask for a theory of such guarantees.
  • 1988–1990. Manasse, McGeoch, and Sleator introduce the kkk-server problem (STOC 1988; J. Algorithms 1990) and settle its extremes: no deterministic algorithm beats ratio kkk on any space with more than kkk points (Corollary 7), two servers admit a 222-competitive algorithm (Theorem 5, algorithm RES), and kkk servers on k+1k+1k+1 points admit a kkk-competitive one (Theorem 4, algorithm BAL). Section 8 poses the kkk-server conjecture, in the symmetric finite setting of the paper.
  • 1990. Fiat, Rabani, and Ravid (FOCS 1990) give the first competitive ratio depending on kkk alone — exponential in kkk, but finite on every metric space.
  • 1991. Chrobak, Karloff, Payne, and Vishwanathan (SIAM J. Discrete Math.) prove the conjecture on the real line via Double Coverage; Chrobak and Larmore (SIAM J. Comput.) extend it to all tree metrics.
  • 1995. Koutsoupias and Papadimitriou (J. ACM) prove the Work Function Algorithm is (2k−1)(2k-1)(2k−1)-competitive on every metric space — the breakthrough, and still the best general bound. Their Conjecture 1.1 fixes the conjecture's modern form: for every metric space there is an online algorithm with competitive ratio kkk.
  • 1996. The same authors verify the conjecture on spaces of k+2k+2k+2 points via the dual 2-evader problem (Inf. Process. Lett. 57).
  • 2004. Bartal and Koutsoupias prove the WFA itself is kkk-competitive on the line, weighted stars, and all spaces of k+2k+2k+2 points.
  • 2021. Coester and Koutsoupias (ICALP) give a unifying potential for all known WFA analyses and push the frontier to the circle.
  • 2023. Bubeck, Coester, and Rabani (STOC) refute the randomized analogue: no o(log⁡2k)o(\log^2 k)o(log2k)-competitive randomized algorithm exists in general. The deterministic conjecture — this mission's goal — survives as the central open question, with the gap between kkk and 2k−12k-12k−1 unmoved since 1995.
  • 2026. Coester, Koutsoupias, and Zbysiński post The kkk-server conjecture is true (arXiv:2609.15979), a claimed proof of the full conjecture: the Work Function Algorithm itself is kkk-competitive on every metric space, via a matrix representation of work functions and a potential function built on it. The preprint is not yet peer-reviewed; this mission's goal stays open until a machine-checked proof exists.

Setting

Fix a metric space MMM with distance function ddd, and a number of servers k≥1k \ge 1k≥1. A configuration records where the kkk servers stand: it is a function CCC assigning to each server i∈{1,…,k}i \in \{1, \dots, k\}i∈{1,…,k} a point C(i)∈MC(i) \in MC(i)∈M. Moving the servers from configuration CCC to configuration C′C'C′ means server iii travels from C(i)C(i)C(i) to C′(i)C'(i)C′(i); the movement cost is the total distance traveled,

moveCost(C,C′)  =  ∑i=1kd(C(i), C′(i)).\mathrm{moveCost}(C, C') \;=\; \sum_{i=1}^{k} d\bigl(C(i),\, C'(i)\bigr).moveCost(C,C′)=i=1∑k​d(C(i),C′(i)).

A request sequence is a finite list σ=(r1,…,rn)\sigma = (r_1, \dots, r_n)σ=(r1​,…,rn​) of points of MMM, presented one at a time; write σ≤j=(r1,…,rj)\sigma_{\le j} = (r_1, \dots, r_j)σ≤j​=(r1​,…,rj​) for the list of the first jjj requests (so σ≤0\sigma_{\le 0}σ≤0​ is the empty list).

A deterministic online algorithm AAA is a rule that, for every finite request sequence ℓ\ellℓ, specifies a configuration A(ℓ)A(\ell)A(ℓ) — where the servers stand after serving the requests of ℓ\ellℓ in order. In particular A(empty list)A(\text{empty list})A(empty list) is the initial configuration, before any request arrives. Two points about this way of modeling an algorithm:

  • Online and deterministic, by construction. The configuration after jjj requests is A(σ≤j)A(\sigma_{\le j})A(σ≤j​), a function of those first jjj requests only — the algorithm cannot see the future, and makes no random choices.
  • The service constraint. Whenever a request sequence ends with a request rrr, some server must stand at rrr immediately after: for every list ℓ\ellℓ and every point rrr, the configuration reached after serving ℓ\ellℓ followed by rrr places at least one server at the point rrr.

Running AAA on σ=(r1,…,rn)\sigma = (r_1, \dots, r_n)σ=(r1​,…,rn​) produces the configurations A(σ≤0), A(σ≤1), …, A(σ≤n)A(\sigma_{\le 0}),\, A(\sigma_{\le 1}),\, \dots,\, A(\sigma_{\le n})A(σ≤0​),A(σ≤1​),…,A(σ≤n​), and its cost is the total movement along this trajectory:

costA(σ)  =  ∑j=1nmoveCost(A(σ≤j−1), A(σ≤j)).\mathrm{cost}_A(\sigma) \;=\; \sum_{j=1}^{n} \mathrm{moveCost}\bigl(A(\sigma_{\le j-1}),\, A(\sigma_{\le j})\bigr).costA​(σ)=j=1∑n​moveCost(A(σ≤j−1​),A(σ≤j​)).

For comparison, an offline schedule for σ\sigmaσ starting at a configuration C0C_0C0​ is any sequence of configurations S0=C0,S1,…,SnS_0 = C_0, S_1, \dots, S_nS0​=C0​,S1​,…,Sn​ in which SjS_jSj​ places a server at the request rjr_jrj​, for each jjj — chosen with the whole of σ\sigmaσ known in advance. The optimal offline cost OPT(C0,σ)\mathrm{OPT}(C_0, \sigma)OPT(C0​,σ) is the infimum, over all such schedules, of the total movement ∑j=1nmoveCost(Sj−1,Sj)\sum_{j=1}^{n} \mathrm{moveCost}(S_{j-1}, S_j)∑j=1n​moveCost(Sj−1​,Sj​).

Finally, AAA is ccc-competitive if there is a constant aaa — depending on the algorithm, hence possibly on the metric space and the initial configuration, but never on the request sequence — with

costA(σ)  ≤  c⋅OPT(A(empty list), σ)+afor every request sequence σ.\mathrm{cost}_A(\sigma) \;\le\; c \cdot \mathrm{OPT}\bigl(A(\text{empty list}),\, \sigma\bigr) + a \qquad \text{for every request sequence } \sigma.costA​(σ)≤c⋅OPT(A(empty list),σ)+afor every request sequence σ.

Formalization targets

Goal — the kkk-server conjecture

For every k≥1, every metric space M, and every initial configuration C0: ∃ A starting at C0 that is k-competitive.\text{For every } k \ge 1,\ \text{every metric space } M,\ \text{and every initial configuration } C_0:\ \exists\, A \text{ starting at } C_0 \text{ that is } k\text{-competitive.}For every k≥1, every metric space M, and every initial configuration C0​: ∃A starting at C0​ that is k-competitive.

The goal fixes no algorithm: any kkk-competitive construction settles it. This is the weakest stable form of the conjecture — it survives every improvement in constants or techniques short of a disproof.

Milestones — the known ladder

The milestones are the classical results between the trivial and the conjectured, each an existence or impossibility statement over the same definitions: the lower bound c≥kc \ge kc≥k on any space with at least k+1k+1k+1 points; the conjecture for k=2k = 2k=2; for spaces of exactly k+1k+1k+1 points; for the real line; the (2k−1)(2k-1)(2k−1) upper bound of the Work Function Algorithm on every space; the conjecture for spaces of exactly k+2k+2k+2 points; the conjecture for three servers in the Manhattan plane (R2,ℓ1)(\mathbb{R}^2, \ell^1)(R2,ℓ1) — the one settled case over a genuinely two-dimensional continuum (Bein–Chrobak–Larmore 2002; reproved by the unifying potential of Coester–Koutsoupias 2021); Coester–Koutsoupias's 2021 result that the Work Function Algorithm itself — not just some algorithm — is 333-competitive for three servers on trees, stated over an explicit formalization of the WFA; and the 2023 Bubeck–Coester–Rabani refutation of the randomized analogue: there are (k+1)(k+1)(k+1)-point spaces on which every randomized algorithm is Ω(log⁡2k)\Omega(\log^2 k)Ω(log2k)-competitive, stated over a mixed-strategy model of randomized online algorithms.

Significance

A proof of the conjecture would close the founding problem of competitive analysis and pin down the exact power of determinism in online optimization over arbitrary metrics; a disproof would separate general metric spaces from every special class where the ratio kkk is known tight. Either outcome recalibrates the field's standard model of adversarial request sequences.

None of these results — not even the lower bound — has a machine-checked proof, and online algorithms as a subject are absent from Mathlib. This mission builds the base layer: a faithful model of online service systems (configurations, online algorithms as prefix functions, offline schedules, competitiveness), the classical possibility and impossibility results over it, and, at the top, the Koutsoupias–Papadimitriou bound, whose potential-function argument is self-contained but delicate. The model is reusable for paging, weighted caching, metrical task systems, and the randomized kkk-server problem.

Difficulty

The obvious first idea — the greedy algorithm, moving the nearest server to each request — is not competitive for any constant, already on three points of the line: two nearby points can ping-pong one server forever while a server parked slightly farther away never moves. Every known competitive algorithm must sometimes move a server other than the nearest one, and the whole difficulty of the conjecture is quantifying exactly how much such foresight-free hedging can achieve. The Work Function Algorithm's analysis via a potential over offline work functions loses a factor of two for reasons nobody has been able to remove; on the lower-bound side, no metric space is known where the deterministic ratio exceeds kkk.

Formalization scope

The Lean model commits to: configurations as functions Fin k → M (labeled servers — equivalent in cost to the unlabeled multiset model, since offline can permute labels for free); algorithms as total functions List M → (Fin k → M) with the service constraint, so a step may move several servers (the standard laziness reduction makes this equivalent to one-move-per-request); costs in ℝ via Metric.dist; the offline optimum as an sInf over schedules, which agrees with the attained minimum on finite spaces; and the additive-constant form of competitiveness, quantified as ∃ a, ∀ σ.

Two conventions guard against trivialization. The additive constant is quantified before the request sequence — allowing it to depend on σ\sigmaσ would make every algorithm 111-competitive. And the lower-bound milestone requires k+1k+1k+1 distinct points (Finset.card = k + 1); on spaces with at most kkk points the conjecture is trivially true and the lower bound false.

Three further definitional layers extend the model. The work function workFunction C₀ σ C is the sInf of (schedule cost + final move to C) over schedules serving σ from C₀, and the Work Function Algorithm WFA is defined on finite spaces with k ≥ 1 servers: after each request it moves to a configuration containing the request minimizing (movement cost) + (work function of the history including the request), a minimizer existing by finiteness and ties broken by a fixed arbitrary choice — matching the standard definition with its "ties broken arbitrarily" (our fixed choice is one admissible instance). A tree is a finite metric space carrying a tree graph whose weighted path lengths realize the metric — exactly "the set of vertices of a tree" of the sources. A randomized algorithm is a mixed strategy: a probability measure over an index type together with a deterministic algorithm per outcome and measurable per-sequence cost; its expected cost is a lower Lebesgue integral in [0,∞][0,\infty][0,∞], and ccc-competitiveness from C0C_0C0​ demands every outcome start at C0C_0C0​ and one additive constant work for all request sequences.

Welcome contributions: proofs of any milestone in any order (the lower bound and the (k+1)(k+1)(k+1)-point case are the natural entry points); alternative algorithms for milestones already closed; and infrastructure lemmas about moveCost, schedules, and work functions published as reusable platform theorems.

Selected references

  • M. Manasse, L. McGeoch, D. Sleator, Competitive algorithms for server problems, J. Algorithms 11 (1990). doi:10.1016/0196-6774(90)90003-W
  • A. Fiat, Y. Rabani, Y. Ravid, Competitive k-server algorithms, FOCS 1990. doi:10.1109/FSCS.1990.89566
  • M. Chrobak, H. Karloff, T. Payne, S. Vishwanathan, New results on server problems, SIAM J. Discrete Math. 4 (1991). doi:10.1137/0404017
  • M. Chrobak, L. Larmore, An optimal on-line algorithm for k servers on trees, SIAM J. Comput. 20 (1991). doi:10.1137/0220008
  • E. Koutsoupias, C. Papadimitriou, On the k-server conjecture, J. ACM 42 (1995). doi:10.1145/210118.210128
  • E. Koutsoupias, C. Papadimitriou, The 2-evader problem, Inf. Process. Lett. 57(5) (1996), 249–252.
  • C. Coester, E. Koutsoupias, Towards the k-server conjecture: a unifying potential, pushing the frontier to the circle, ICALP 2021. arXiv:2102.10474
  • S. Bubeck, C. Coester, Y. Rabani, The randomized k-server conjecture is false!, STOC 2023. arXiv:2211.05753
  • E. Koutsoupias, The k-server problem (survey), Computer Science Review 3 (2009). doi:10.1016/j.cosrev.2009.04.002
122 thms11 active usersReviewed
Linear OptimizationOperations ResearchTheoretical Computer Science·Captain: mikedeng1

The Design of Competitive Online Algorithms via a Primal-Dual Approach I: The Online Packing-Covering FrameworkTextbook

Motivation

Many online problems — renting vs. buying equipment, routing traffic without knowing future demand, allocating advertising budget as bids arrive — share a common linear-programming shape: a covering (minimization) problem whose constraints appear one at a time, or its dual packing (maximization) problem whose variables appear one at a time, with no ability to revisit past decisions. Buchbinder and Naor's survey [1] recasts the online primal-dual method, originally developed for offline approximation in Section 2 of the same survey (already covered on this platform via the PrimalDualOnline namespace, citing Buchbinder's thesis [2]), as a general recipe for such online problems, unifying earlier ad hoc analyses of the ski-rental problem (Chapter 3) and paving the way for the online set-cover, routing, caching, and ad-auction algorithms formalized elsewhere in this series. This mission covers Chapter 4, "The Basic Approach": the framework itself and its three founding algorithms.

Setting

Fix a finite index set III of primal variables with non-negative cost coefficients cic_ici​, and a finite index set JJJ of covering constraints, revealed one at a time in the order enumerated by JJJ. Each constraint jjj is given by a set S(j)⊆IS(j) \subseteq IS(j)⊆I (the book's simplified setting, in which every non-zero coefficient equals 111 and every right-hand side equals 111; Chapter 14 removes this restriction) and asserts ∑i∈S(j)xi≥1\sum_{i \in S(j)} x_i \ge 1∑i∈S(j)​xi​≥1. An online covering algorithm may only increase the xix_ixi​, never decrease them, and upon a constraint's arrival must eventually make it hold. The online covering problem is to minimize ∑icixi\sum_i c_i x_i∑i​ci​xi​ subject to every revealed constraint, online. Its Lagrangian dual is the online packing problem: a dual variable yjy_jyj​ arrives together with constraint jjj, may only be increased while jjj is being processed, and the objective is to maximize ∑jyj\sum_j y_j∑j​yj​ subject to ∑j∣i∈S(j)yj≤ci\sum_{j \mid i \in S(j)} y_j \le c_i∑j∣i∈S(j)​yj​≤ci​ for every iii — the packing constraint on iii becomes fully known only once every jjj with i∈S(j)i \in S(j)i∈S(j) has arrived, so it, too, is revealed gradually. d:=max⁡j∣S(j)∣d := \max_j |S(j)|d:=maxj​∣S(j)∣, the largest constraint size, is carried as an explicit parameter throughout.

Section 4.2 gives three algorithms solving both problems simultaneously — the same run produces a covering solution xxx and a packing solution yyy — with the same worst-case guarantee but different flavors: Algorithm 1 is a discrete process (each processing round performs a whole number of identical multiplicative-plus-additive updates until its constraint is satisfied); Algorithm 2 is the continuous limit of Algorithm 1 (dual variables increase continuously and xix_ixi​ follows an explicit exponential of the accumulated dual sum); Algorithm 3 replaces the continuous update with one triggered by an approximate complementary-slackness condition, at the cost of a mild dual infeasibility. All three make essential use of the online order: an algorithm that saw the whole instance up front would trivially solve the offline LP.

Formalization targets

Theorem 4.3 (Algorithm 3 — the goal, p. 124):

(∀i, ∑j∣i∈S(j)yj≤ci(1+ln⁡d)) ∧ (∀x′′ feasible, ∑icixi≤2(1+ln⁡d)∑icixi′′) ∧ (∀y′′ feasible, ∑jyj′′≤2∑jyj).\Big(\forall i,\ \textstyle\sum_{j \mid i \in S(j)} y_j \le c_i(1+\ln d)\Big) \ \wedge\ \Big(\forall x''\text{ feasible},\ \textstyle\sum_i c_i x_i \le 2(1+\ln d)\sum_i c_i x''_i\Big) \ \wedge\ \Big(\forall y''\text{ feasible},\ \textstyle\sum_j y''_j \le 2\sum_j y_j\Big).(∀i, ∑j∣i∈S(j)​yj​≤ci​(1+lnd)) ∧ (∀x′′ feasible, ∑i​ci​xi​≤2(1+lnd)∑i​ci​xi′′​) ∧ (∀y′′ feasible, ∑j​yj′′​≤2∑j​yj​).

Theorem 4.1 (Algorithm 1) and Theorem 4.2 (Algorithm 2, p. 118 and p. 121) are the same three-part guarantee for the other two algorithms, with log⁡2(3d+1)\log_2(3d+1)log2​(3d+1) in place of 1+ln⁡d1+\ln d1+lnd for Algorithm 1 (and its packing solution genuinely integral), and with 2ln⁡(1+d)2\ln(1+d)2ln(1+d) in place of both 2(1+ln⁡d)2(1+\ln d)2(1+lnd) and 222 for Algorithm 2 (whose packing solution is exactly feasible, not merely approximately so). The competitive ratios are stated against an arbitrary offline-feasible comparison solution on each side (covering and packing) rather than against an unconstructed LP optimum — the standard weak-duality reformulation of "ccc-competitive", and the one the book's own proofs (which invoke weak duality directly, never LP optimality) actually establish.

Significance

The framework converts three qualitatively different design intuitions — the ski-rental-style discrete doubling, the continuous primal-dual differential equation, and complementary slackness — into three algorithms with an identical asymptotic guarantee, Θ(log⁡d)\Theta(\log d)Θ(logd), matching the Ω(log⁡d)\Omega(\log d)Ω(logd) (packing) and Ω(log⁡n)\Omega(\log n)Ω(logn) (covering) lower bounds the book proves in Section 4.3 (Lemmas 4.5-4.6, not part of this mission). This is the load-bearing substrate for the rest of the survey: Chapter 5's online set-cover algorithm, Chapter 13's bounded-allocation problem, and Chapter 14's general packing-covering constraints all restate this chapter's framework locally rather than re-deriving it, and are formalized as separate missions in this series. Formalizing it here, once, with the exact constants each proof establishes, is what lets those missions cite a single faithful statement instead of three independently-drifting restatements. No formal development of this framework was found on the platform as of 2026-09-20 (searches below); this mission is the first.

Difficulty

The central obstacle is not the algebra of any single algorithm's proof — each is a short, self-contained argument — but stating the guarantee for an online process using only its final output. An algorithm is characterized here by the closed-form relation its update rule establishes between the accumulated dual sum and the primal value (e.g., Algorithm 3's xi=min⁡(1,d−1exp⁡(di/ci−1))x_i = \min(1, d^{-1}\exp(d_i/c_i - 1))xi​=min(1,d−1exp(di​/ci​−1)) once activated, 000 before), together with primal feasibility as the hypothesis that a run has completed; a formalization that instead handed the algorithm the whole instance in advance, or dropped primal feasibility as a hypothesis, would either trivialize the online promise or make the stated bound simply false. A second difficulty specific to this formalization: Theorem 4.3's own proof bounds the primal cost by splitting it into a piece controlled by the final-state complementary-slackness conditions (immediate from the closed form and the ddd-bound) and a piece controlled by a derivative/telescoping argument over the continuous accumulation process itself — the latter is a genuinely dynamic fact about the trajectory, not just its endpoint, and is recorded as a documented simplification below rather than folded into the hypotheses, since the goal is a faithful statement, not a proof.

Formalization scope

CoveringInstance I J bundles S : J → Finset I, c : I → ℝ, and d : ℝ with 0 < d and ∀ j, (S j).card ≤ d as explicit hypotheses (never derived as d := ⨆ j, (S j).card, matching the book's own presentation and avoiding a vacuous formalization in which d is chosen after the fact to make the bound trivial). dualSum inst y i := ∑_{j \mid i \in S(j)} y_j. Each algorithm's output is a noncomputable def from the accumulated dual data to ℝ (alg1X, alg2X, alg3X), so that "the algorithm's output" is genuinely a function of its dual trajectory rather than an independently-constrained free variable — ruling out the trivializing formalization in which xxx and yyy are unrelated variables merely required to satisfy the conclusion's own inequalities. Reals throughout (ℝ, not ℝ≥0 or ENNReal); Finset.filter realizes "jjj such that i∈S(j)i \in S(j)i∈S(j)"; Real.logb 2 and Real.log (natural log) match the book's own log₂ and ln. Reusable beyond this mission: CoveringInstance, dualSum, and the weak-duality-style "competitive against any feasible comparison solution" pattern, which 05-online-set-cover, 13-bounded-allocation, and 14-general-packing-covering are expected to restate locally (per this series' own rule against cross-draft imports between concurrent missions) rather than import directly. Welcome contributions: completing any of the three sorrys, and formalizing the Section 4.3 lower bounds (Lemmas 4.5-4.6) as a follow-on mission.

Selected references

  • 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):93-263, 2009. https://doi.org/10.1561/0400000024
  • N. Buchbinder. Designing Competitive Online Algorithms via a Primal-Dual Approach. PhD thesis, Tel Aviv University, 2008. https://www.tau.ac.il/~nivb/download/phd-thsis.pdf
12 thms3 active users
Algorithmic Game TheoryOperations ResearchTheoretical Computer Science·Captain: mikedeng1

The Design of Competitive Online Algorithms via a Primal-Dual Approach VI: Maximizing Ad-Auctions RevenueTextbook

Motivation

Search-engine ad-auctions sell items (ad slots) to buyers who arrive with per-item bids and a fixed daily budget, online: each item must be allocated the moment it appears, with no ability to revisit past allocations, and no buyer may ever be charged more than its budget. Buchbinder and Naor's survey [1] models this as a generalization of online bipartite matching and derives an allocation algorithm through the same online primal-dual recipe formalized elsewhere in this series (Chapter 4's framework), applied to a genuinely different LP shape — a maximization (revenue) problem whose "packing" role and "covering" role are the reverse of Chapters 4, 5 and 7 — and to a different constraint structure (budget caps on a per-buyer accumulated sum, not a per-constraint 0/1 covering requirement). This mission covers Section 10.1, "The Basic Algorithm": the single-slot allocation algorithm and its (1−1/c)(1−Rmax⁡)(1-1/c)(1-R_{\max})(1−1/c)(1−Rmax​)-competitive analysis (Theorem 10.1). Sections 10.2-10.3 (multiple ad-slots via strong duality; stochastic per-buyer spending guarantees) are out of scope — see Formalization scope.

Setting

Fix a finite set III of buyers, each with budget Bi>0B_i > 0Bi​>0, and a finite set MMM of items; buyer iii bids bij≥0b_{ij} \ge 0bij​≥0 for item jjj, revealed one item at a time in the order enumerated by MMM. Let Rmax⁡:=max⁡i,jbij/BiR_{\max} := \max_{i,j} b_{ij}/B_iRmax​:=maxi,j​bij​/Bi​, carried as an explicit positive parameter. The Allocation algorithm (p. 212), upon each item jjj's arrival, allocates it to the buyer iii maximizing bij(1−xi)b_{ij}(1-x_i)bij​(1−xi​) (where xi∈[0,1]x_i \in [0,1]xi​∈[0,1] is buyer iii's current primal value); if xi≥1x_i \ge 1xi​≥1 already, nothing happens (the buyer is "full"). Otherwise it charges buyer iii the minimum of bijb_{ij}bij​ and its remaining budget, sets the dual allocation variable yij←1y_{ij} \leftarrow 1yij​←1, sets zj←bij(1−xi)z_j \leftarrow b_{ij}(1-x_i)zj​←bij​(1−xi​), and updates xi←xi(1+bij/Bi)+bij/((c−1)Bi)x_i \leftarrow x_i(1+b_{ij}/B_i) + b_{ij}/((c-1)B_i)xi​←xi​(1+bij​/Bi​)+bij​/((c−1)Bi​) for a constant ccc fixed by the analysis. The revenue actually collected from buyer iii is min⁡ ⁣(∑jbijyij, Bi)\min\!\big(\sum_j b_{ij}y_{ij},\,B_i\big)min(∑j​bij​yij​,Bi​) — the buyer is never charged more than its budget. Unlike Chapters 4/7, this LP's dual (the ad-auctions revenue objective, Fig. 10.1) is the maximization problem being solved online; the "primal" covering LP (xix_ixi​, zjz_jzj​ variables) exists only as a duality certificate.

Formalization targets

Theorem 10.1 (the goal, p. 212), given the milestone's dual near-feasibility bound and the fact that each buyer's total accrued bids exceed its budget by at most a factor of Rmax⁡R_{\max}Rmax​ (Claim (3)'s consequence):

∀ (x′′,z′′) feasible for Fig. 10.1’s covering LP,  ∑iactualCharge(i) ≥ (1−1c)(1−Rmax⁡)(∑iBixi′′+∑jzj′′),\forall\, (x'', z'')\text{ feasible for Fig. 10.1's covering LP},\ \ \textstyle\sum_i \mathrm{actualCharge}(i) \ \ge\ (1-\tfrac1c)(1-R_{\max}) \Big(\textstyle\sum_i B_i x''_i + \sum_j z''_j\Big),∀(x′′,z′′) feasible for Fig. 10.1’s covering LP,  ∑i​actualCharge(i) ≥ (1−c1​)(1−Rmax​)(∑i​Bi​xi′′​+∑j​zj′′​),

with c=(1+Rmax⁡)1/Rmax⁡c = (1+R_{\max})^{1/R_{\max}}c=(1+Rmax​)1/Rmax​ taken verbatim from the theorem's own statement — the exact formula, not an O(⋅)O(\cdot)O(⋅) instantiation. Inequality (10.1) (p. 213), the milestone, is the book's own induction-proved lower bound on a buyer's primal value in terms of its accrued bids: xi≥1c−1(c(∑jbijyij)/Bi−1)x_i \ge \frac{1}{c-1}\big(c^{(\sum_j b_{ij}y_{ij})/B_i} - 1\big)xi​≥c−11​(c(∑j​bij​yij​)/Bi​−1).

Significance

Theorem 10.1 is the entry point to a short but influential sub-line of the online primal-dual method — Section 10.2 extends it to multiple ad-slots via strong duality for maximum-weight matching (rather than the weak duality this framework otherwise relies on throughout), and Section 10.3 incorporates stochastic per-buyer spending guarantees, both reusing this section's constant c=(1+Rmax⁡)1/Rmax⁡c=(1+R_{\max})^{1/R_{\max}}c=(1+Rmax​)1/Rmax​ and its limit c→ec\to ec→e as Rmax⁡→0R_{\max}\to0Rmax​→0 (recovering the classic (1−1/e)(1-1/e)(1−1/e)-competitive ratio for the unweighted, unbudgeted case). It is also the first mission in this series to apply the online primal-dual method to a genuine revenue-maximization problem rather than a covering/packing pair with matching roles. No formal development of ad-auctions, budgeted online matching, or this constant was found on the platform as of 2026-09-20 (searches below); this mission is the first.

Difficulty

As with 04-framework's Algorithm 3 and 07-generalized-caching's Fractional Caching algorithm, the central obstacle is characterizing an online process by its own final output. Unlike those two chapters, however, step (3)'s update increment varies per allocation (bd, the specific bid of the item just won), so no closed-form solution of the recurrence exists in general; buyerX is instead defined as an explicit List.foldl realizing the exact per-step update, over the (temporally ordered) list of bids a buyer actually won — a faithful, if less immediately readable, transcription of "the algorithm's output as a function of its own trajectory," in the same spirit as this series' other closed-form definitions. A second difficulty specific to this chapter: the book's derivation of inequality (10.1) is itself an induction on iterations (not a single algebraic step, unlike Eq. (7.2) in 07-generalized-caching), and the theorem's final bound further combines it with a separate "at most one undercharged iteration" argument (Claim (3)'s conclusion, p. 214-215) turning the raw accrued-bid bound into one about the actually collected (budget-capped) revenue — both are stated here as explicit hypotheses (the milestone, and h_at_most_one_undercharge) rather than derived, since the goal is a faithful statement, not a proof; both sorrys are documented, not silently discharged.

Formalization scope

AdAuctionsInstance I M bundles b : I → M → ℝ, B : I → ℝ (hB_pos), and Rmax : ℝ with hRmax_pos : 0 < Rmax and hRmax_bound : ∀ i j, b i j ≤ Rmax * B i — the last two as explicit hypotheses, never derived via Finset.sup, matching 04-framework's d and 07-generalized-caching's k. cParam inst := (1+Rmax)^(1/Rmax) uses Real.rpow (ℝ^ℝ). buyerX inst i bids folds step (3)'s update over a list of won bids; revenue/actualCharge realize ∑jbijyij\sum_j b_{ij}y_{ij}∑j​bij​yij​ and its budget-capped charge. Reals throughout. Explicitly out of scope: Section 10.2's multiple-slot generalization (Theorem 10.2), which requires strong duality for maximum-weight bipartite matching as an explicit premise (the book: "our analysis... crucially relies on strong duality") — a substantially different LP structure (an integral matching LP, not this section's per-buyer budget LP) that this mission's AdAuctionsInstance does not model, and whose applicability of PrimalDualOnline.LP.strong_duality_adapter (this book's own Theorem 2.2) was not verified in the time available. Section 10.3's stochastic guarantee (Theorem 10.3) is likewise out of scope, and its own BRIEF.md-flagged ambiguity (whether its scalar g is min⁡igi\min_i g_imini​gi​ or another aggregate of the per-buyer vector gig_igi​) was not resolved. Both are natural follow-on missions, not attempted here. Welcome contributions: completing the two sorrys, and the Section 10.2-10.3 follow-on mission(s).

Selected references

  • 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):93-263, 2009. https://doi.org/10.1561/0400000024
8 thms3 active users
Operations ResearchOptimizationProbability+1·Captain: mikedeng1

An Optimal On-Line Algorithm for Metrical Task System 2: The Randomized Competitive Ratio of the Uniform Task System Lies Between H(n) and 2H(n)Research Paper

Motivation

Metrical task systems, introduced by Borodin, Linial and Saks (J. ACM 39(4), 1992), are a common abstraction of on-line problems in which a server occupies one of finitely many states, pays a cost for each task depending on its current state, and may pay a transition cost to change state first. Paging, list update and the kkk-server problem all fit into this framework. The paper's first main result is that every deterministic on-line algorithm on an nnn-state metrical task system has competitive ratio at least 2n−12n-12n−1, and that 2n−12n-12n−1 is attained. The lower bound comes from an adversary that always charges the state the algorithm currently occupies. That adversary needs to know the algorithm's state, which suggests that randomization can help.

Section 7 of the paper makes this precise for the simplest system, the uniform task system, in which all transitions cost 111. There the randomized competitive ratio against an oblivious adversary is between H(n)H(n)H(n) and 2H(n)2H(n)2H(n), where H(n)=1+12+⋯+1nH(n)=1+\tfrac12+\cdots+\tfrac1nH(n)=1+21​+⋯+n1​ is between ln⁡n\ln nlnn and 1+ln⁡n1+\ln n1+lnn. This was the first logarithmic bound for a task system.

Timeline:

  • 1985: Sleator and Tarjan introduce competitive analysis for list update and paging (CACM 28(2)).
  • 1987/1992: Borodin, Linial and Saks define metrical task systems, prove the deterministic ratio 2n−12n-12n−1, and prove H(n)≤wˉ≤2H(n)H(n)\le\bar w\le 2H(n)H(n)≤wˉ≤2H(n) for the uniform system (conference version STOC 1987; journal version cited above).
  • 1991: Fiat, Karp, Luby, McGeoch, Sleator and Young prove the analogous 2Hk2H_k2Hk​ upper bound for randomized paging (J. Algorithms 12(4)).

Setting

A task system (S,d)(S,d)(S,d) is a finite set SSS of nnn states with a transition-cost matrix ddd: d(i,i)=0d(i,i)=0d(i,i)=0, d(i,j)>0d(i,j)>0d(i,j)>0 for i≠ji\ne ji=j, and d(i,k)≤d(i,j)+d(j,k)d(i,k)\le d(i,j)+d(j,k)d(i,k)≤d(i,j)+d(j,k). In the uniform task system, d(i,j)=1d(i,j)=1d(i,j)=1 for all i≠ji\neq ji=j. A task is a vector T∈R≥0ST\in\mathbb R_{\ge0}^ST∈R≥0S​ of processing costs. Given an initial state s0s_0s0​ and tasks T=T1⋯Tm\mathbf T=T^1\cdots T^mT=T1⋯Tm, a schedule is σ:{0,…,m}→S\sigma:\{0,\dots,m\}\to Sσ:{0,…,m}→S with σ(0)=s0\sigma(0)=s_0σ(0)=s0​, of cost

c(T;σ)=∑i=1md(σ(i−1),σ(i))+∑i=1mTi(σ(i)).c(\mathbf T;\sigma)=\sum_{i=1}^m d(\sigma(i-1),\sigma(i))+\sum_{i=1}^m T^i(\sigma(i)).c(T;σ)=i=1∑m​d(σ(i−1),σ(i))+i=1∑m​Ti(σ(i)).

The off-line optimum c0(T)c_0(\mathbf T)c0​(T) is the least cost over all schedules.

A deterministic on-line algorithm chooses σ(i)\sigma(i)σ(i) from s0s_0s0​ and T1,…,TiT^1,\dots,T^iT1,…,Ti. A randomized on-line algorithm RRR chooses σ(i)\sigma(i)σ(i) at random, with a distribution that depends on s0s_0s0​, on T1,…,TiT^1,\dots,T^iT1,…,Ti and on the states σ(0),…,σ(i−1)\sigma(0),\dots,\sigma(i-1)σ(0),…,σ(i−1) already visited. The task sequence is fixed before any random choice is made (an oblivious adversary). With pr(σ∣T)\mathrm{pr}(\sigma\mid\mathbf T)pr(σ∣T) the probability that RRR follows σ\sigmaσ, the expected cost is cˉR(T)=∑σc(T;σ) pr(σ∣T)\bar c_R(\mathbf T)=\sum_\sigma c(\mathbf T;\sigma)\,\mathrm{pr}(\sigma\mid\mathbf T)cˉR​(T)=∑σ​c(T;σ)pr(σ∣T). For w>0w>0w>0, RRR is expected www-competitive if there is a constant KKK with

cˉR(T)≤w c0(T)+K\bar c_R(\mathbf T)\le w\,c_0(\mathbf T)+KcˉR​(T)≤wc0​(T)+K

for every finite task sequence and every initial state. The randomized competitive ratio wˉ(S,d)\bar w(S,d)wˉ(S,d) is the infimum of all such www over all RRR.

Formalization targets

Goal: Theorem 7.1

For the uniform task system on n≥1n\ge1n≥1 states,

H(n)  ≤  wˉ(S,d)  ≤  2H(n).H(n)\;\le\;\bar w(S,d)\;\le\;2H(n).H(n)≤wˉ(S,d)≤2H(n).

Milestones

  1. Upper bound (p. 759). Some randomized on-line algorithm is expected 2H(n)2H(n)2H(n)-competitive on the uniform task system.
  2. Lemma 7.2 (p. 759). Let DDD be a distribution on infinite task sequences over a finite task alphabet, with E(c0(Tj))→∞E(c_0(\mathbf T^j))\to\inftyE(c0​(Tj))→∞, and let mj=inf⁡AE(cA(Tj))m_j=\inf_A E(c_A(\mathbf T^j))mj​=infA​E(cA​(Tj)) over deterministic on-line algorithms. Then every achievable www satisfies
lim sup⁡j→∞mjE(c0(Tj))≤w.\limsup_{j\to\infty}\frac{m_j}{E(c_0(\mathbf T^j))}\le w .j→∞limsup​E(c0​(Tj))mj​​≤w.
  1. mj≥j/nm_j\ge j/nmj​≥j/n (p. 760) when the tasks are independent uniformly random unit elementary tasks UsU_sUs​ (cost 111 in sss, 000 elsewhere).
  2. Coupon collector (p. 760). For i.i.d. uniform states on SSS, the expected number of draws until every state has appeared is nH(n)nH(n)nH(n).
  3. Off-line cost (p. 760). Under the same distribution, E(c0(Tj))≤j/(nH(n))+CE(c_0(\mathbf T^j))\le j/(nH(n))+CE(c0​(Tj))≤j/(nH(n))+C for a constant CCC independent of jjj.

Significance

The theorem shows that randomization reduces the competitive ratio of the uniform task system from 2n−12n-12n−1 to Θ(log⁡n)\Theta(\log n)Θ(logn). That is an exponential improvement, and it identifies the adversary's knowledge of the algorithm's state as the source of the deterministic lower bound. Lemma 7.2 is a form of Yao's principle adapted to the additive-constant definition of competitiveness. It is the standard tool for randomized lower bounds in on-line computation, and the same argument shape reappears for paging and kkk-server lower bounds.

On the formalization side, the result is proved but, as far as is known, has not been machine-checked. A complete development yields a reusable model of randomized on-line algorithms with oblivious adversaries, a Yao-type lemma usable for other on-line problems, and a coupon-collector expectation in the product-measure setting. The upper half additionally needs the continuous-time reduction of the paper's Lemma 3.1 in randomized form, or a direct discrete-time algorithm.

Difficulty

For the upper bound, the natural algorithm is continuous-time. It proceeds in phases, and inside a phase it stays in a state until that state has accumulated cost 111. A discrete task can saturate several states at once and straddle a phase boundary. So a discrete algorithm must either simulate the continuous one or be analyzed directly, and the expected transition count per phase must be controlled with the first phase starting in a deterministic state.

For the lower bound, the first obstacle is that the natural statement "wˉ≥lim sup⁡mj/E(c0)\bar w\ge\limsup m_j/E(c_0)wˉ≥limsupmj​/E(c0​)" silently assumes that a randomized algorithm's expected cost, averaged over random inputs, is at least that of the best deterministic algorithm. With the behavioural (kernel) definition used here, this requires converting a kernel into a mixture of deterministic algorithms, which is Kuhn's theorem on each finite horizon. The second obstacle is that the paper's claim E(c0(Tj))≤j/(nH(n))+O(1)E(c_0(\mathbf T^j))\le j/(nH(n))+O(1)E(c0​(Tj))≤j/(nH(n))+O(1) is supported only by the elementary renewal theorem, which gives a limit of ratios; the additive bound needs a sharper renewal estimate. Mathlib has no renewal theory. The hypothesis E(c0(Tj))→∞E(c_0(\mathbf T^j))\to\inftyE(c0​(Tj))→∞ of Lemma 7.2 must also be established for the uniform distribution; the paper does not prove it separately.

Formalization scope

  • States form a finite nonempty type S, and nnn = Fintype.card S; no n≥2n\ge2n≥2 assumption is made (at n=1n=1n=1 the goal reads 1≤wˉ≤21\le\bar w\le21≤wˉ≤2, and wˉ=1\bar w=1wˉ=1). H(n)H(n)H(n) is Mathlib's harmonic n, cast to R\mathbb RR. The uniform system has unit transition cost.
  • Tasks are finite and nonnegative. The paper also allows +∞+\infty+∞ entries; these are excluded. Task sequences are Fin m → S → ℝ and schedules are Fin (m+1) → S with σ 0 = s₀. c0c_0c0​ is a finite minimum.
  • A randomized algorithm is a kernel S → List (S → ℝ) → List S → PMF S. This is the paper's scheduler–taskmaster description (p. 758), equivalent to a distribution over deterministic algorithms on every finite task sequence. pr(σ∣T)\mathrm{pr}(\sigma\mid\mathbf T)pr(σ∣T) is the product of kernel probabilities, and cˉR\bar c_RcˉR​ is a finite sum. That pr(⋅∣T)\mathrm{pr}(\cdot\mid\mathbf T)pr(⋅∣T) sums to 111 has been checked locally.
  • wˉ(S,d)\bar w(S,d)wˉ(S,d) is the real sInf of {w:∃R, R expected w-competitive}\{w : \exists R,\ R \text{ expected } w\text{-competitive}\}{w:∃R, R expected w-competitive}. On the empty set this would be 000, so the upper bound is stated as the existence of an expected 2H(n)2H(n)2H(n)-competitive algorithm, and Lemma 7.2 is stated for every achievable www. The goal's lower half forces the set to be nonempty. Statements of the form "wˉ≤c\bar w\le cwˉ≤c" alone are therefore not acceptable substitutes for milestones 1 and 2.
  • Lemma 7.2 is restricted to task sequences over a finite alphabet, with the product σ\sigmaσ-algebra and measurable singletons. This makes every E(cA(Tj))E(c_A(\mathbf T^j))E(cA​(Tj)) a genuine integral for every deterministic AAA, and it covers the paper's application. The lim sup⁡\limsuplimsup of Lemma 7.2 is taken in EReal.
  • The coupon-collector time takes values in [0,∞][0,\infty][0,∞] and its expectation is a lower Lebesgue integral. Milestone 5 renders the paper's O(1)O(1)O(1) as an explicit constant CCC chosen before jjj.

Contributions welcome: proofs of any milestone; a discrete-time randomized phase algorithm; a general Kuhn-type conversion from kernels to mixtures of deterministic algorithms; renewal-theoretic lemmas.

Selected references

  • A. Borodin, N. Linial, M. Saks, An Optimal On-Line Algorithm for Metrical Task System, J. ACM 39(4):745–763, 1992. https://doi.org/10.1145/146585.146588
  • D. D. Sleator, R. E. Tarjan, Amortized Efficiency of List Update and Paging Rules, Commun. ACM 28(2):202–208, 1985. https://doi.org/10.1145/2786.2793
  • A. Fiat, R. M. Karp, M. Luby, L. A. McGeoch, D. D. Sleator, N. E. Young, Competitive Paging Algorithms, J. Algorithms 12(4):685–699, 1991. https://doi.org/10.1016/0196-6774(91)90041-V
  • A. C.-C. Yao, Probabilistic Computations: Toward a Unified Measure of Complexity, FOCS 1977, 222–227. https://doi.org/10.1109/SFCS.1977.24
  • S. M. Ross, Applied Probability Models with Optimization Applications, Holden-Day, 1970 (the elementary renewal theorem cited as [20] in the paper).
9 thms1 active userReviewed
Machine LearningOperations ResearchProbability+1·Captain: mikedeng1

Competitive Caching with Machine Learned Advice: The Competitive Ratio of Predictive MarkerResearch Paper

Motivation

Caching (online paging) is one of the oldest problems in online algorithms: a fast memory of kkk slots serves a sequence of requests, and every request for an element not in the fast memory is a cache miss that forces the element to be loaded, possibly evicting another one. With the whole request sequence known in advance, evicting the element whose next request is furthest in the future is optimal (Bélády, 1966). Without that knowledge, no deterministic algorithm is better than kkk-competitive, and the best randomized algorithms are Θ(log⁡k)\Theta(\log k)Θ(logk)-competitive (Fiat, Karp, Luby, McGeoch, Sleator and Young, 1991).

Lykouris and Vassilvitskii asked what happens in between: an online algorithm receives, with every request, a machine-learned prediction of the element's next arrival time. A good predictor should make the algorithm nearly as good as Bélády's rule (consistency), and a bad predictor should never make it worse than a classical algorithm (robustness). Their paper (arXiv:1802.05399v4; J. ACM 2021) is one of the founding papers of learning-augmented algorithms, and its algorithm, Predictive Marker, is the reference point for the later literature on caching with predictions.

Timeline.

  • 1966: Bélády's furthest-in-future rule is optimal offline.
  • 1985: Sleator and Tarjan show that deterministic online paging is at best kkk-competitive.
  • 1991: Fiat et al. introduce the Marker algorithm, 2Hk2H_k2Hk​-competitive, and the clean-element lower bound on the optimum.
  • 2018: Lykouris and Vassilvitskii (arXiv:1802.05399) introduce Predictive Marker, with ratio 2min⁡(1+2Sℓ(ϵ),2Hk)2\min(1+2S_\ell(\epsilon), 2H_k)2min(1+2Sℓ​(ϵ),2Hk​) for an ϵ\epsilonϵ-accurate predictor.
  • 2020: Rohatgi (arXiv:1910.12172, SODA 2020) and Wei (arXiv:2005.13716, APPROX/RANDOM 2020) improve the dependence on the error.

Setting

A request sequence σ=(z1,…,zn)\sigma = (z_1, \dots, z_n)σ=(z1​,…,zn​) lists elements of a set ZZZ. A cache of size k≥1k \ge 1k≥1 starts empty. A request for a cached element is a hit; otherwise it is a miss, the element is loaded, and if the cache is full some element is evicted first. The offline optimum Opt(σ)\mathrm{Opt}(\sigma)Opt(σ) is the least number of misses over all eviction schedules chosen with knowledge of σ\sigmaσ.

With each request ziz_izi​ the algorithm receives a real prediction hih_ihi​. The true label yiy_iyi​ is the position of the next request of ziz_izi​, or n+1n+1n+1 if there is none. For a loss function ℓ≥0\ell \ge 0ℓ≥0, the error of the predictions is ηℓ(h,σ)=∑iℓ(yi,hi)\eta_\ell(h,\sigma) = \sum_i \ell(y_i, h_i)ηℓ​(h,σ)=∑i​ℓ(yi​,hi​), and the predictions are ϵ\epsilonϵ-accurate when ηℓ(h,σ)≤ϵ⋅Opt(σ)\eta_\ell(h,\sigma) \le \epsilon \cdot \mathrm{Opt}(\sigma)ηℓ​(h,σ)≤ϵ⋅Opt(σ).

The spread of ℓ\ellℓ measures how cheaply a predictor can get the order of arrivals completely wrong: Sℓ(m)S_\ell(m)Sℓ​(m) is the least length T≥1T \ge 1T≥1 such that every strictly increasing integer sequence a1<⋯<aTa_1 < \dots < a_Ta1​<⋯<aT​ and every non-increasing real sequence b1≥⋯≥bTb_1 \ge \dots \ge b_Tb1​≥⋯≥bT​ have total loss ∑iℓ(ai,bi)≥m\sum_i \ell(a_i, b_i) \ge m∑i​ℓ(ai​,bi​)≥m.

Predictive Marker (Algorithm 1) works in the phases of the Marker algorithm. Requested elements are marked. A phase ends when the cache is full, every cached element is marked, and a miss occurs; then all marks are removed. An element requested in a phase but not in the previous one is clean, and Q(σ)Q(\sigma)Q(σ) is the total number of clean elements. Each clean miss starts a chain. An element evicted in the current phase that is requested again (a stale miss) extends the chain in which it was evicted. Evictions are among unmarked elements. As long as the chain's length n(r,c)n(r,c)n(r,c) is at most Hk=1+12+⋯+1kH_k = 1 + \tfrac12 + \dots + \tfrac1kHk​=1+21​+⋯+k1​, the evicted element is one with the largest prediction. After that it is chosen uniformly at random. The expected number of misses of Predictive Marker is costPM(σ)\mathrm{cost}_{PM}(\sigma)costPM​(σ).

Formalization targets

Goal: Theorem 3.3

If SSS is concave on [0,∞)[0,\infty)[0,∞) and majorizes the spread, then for every ϵ≥0\epsilon \ge 0ϵ≥0, every tie-breaking rule, and every sequence with ϵ\epsilonϵ-accurate predictions,

E[costPM(σ)]≤2⋅min⁡(1+2S(ϵ), 2Hk)⋅Opt(σ).\mathbb E\bigl[\mathrm{cost}_{PM}(\sigma)\bigr] \le 2\cdot\min\bigl(1 + 2S(\epsilon),\ 2H_k\bigr)\cdot \mathrm{Opt}(\sigma).E[costPM​(σ)]≤2⋅min(1+2S(ϵ), 2Hk​)⋅Opt(σ).

Milestones

  • Claim 1 (Fiat et al.): Q(σ)≤2 Opt(σ)Q(\sigma) \le 2\,\mathrm{Opt}(\sigma)Q(σ)≤2Opt(σ).
  • Proof of Theorem 3.3, last sentence: Opt(σ)≤Q(σ)\mathrm{Opt}(\sigma) \le Q(\sigma)Opt(σ)≤Q(σ).
  • Lemma 3.3: a chain that evicts by the predictions only has length n(r,c)≤1+S(ηr,c)n(r,c) \le 1 + S(\eta_{r,c})n(r,c)≤1+S(ηr,c​), where ηr,c\eta_{r,c}ηr,c​ is the error of the predictions on the elements evicted into it.
  • Lemma 3.4: E[n(r,c)]≤E[min⁡(1+2S(ηr,c),2Hk)]\mathbb E[n(r,c)] \le \mathbb E[\min(1 + 2S(\eta_{r,c}), 2H_k)]E[n(r,c)]≤E[min(1+2S(ηr,c​),2Hk​)].

Significance

The result. Theorem 3.3 gives both guarantees at once. For an exact predictor (ϵ=0\epsilon = 0ϵ=0) the ratio is a constant, 2(1+2S(0))2(1 + 2S(0))2(1+2S(0)), independent of kkk; for an arbitrary predictor it is 4Hk4H_k4Hk​, within a constant factor of the optimal randomized ratio. In between, the ratio degrades with the error at the rate of the spread: for the absolute loss the spread grows like m\sqrt mm​, so the ratio grows like ϵ\sqrt\epsilonϵ​. The spread and the chain decomposition are the tools later papers build on to trade consistency against robustness.

Formalizing it. The theorem is proved on paper; no machine-checked proof of it, of the Marker analysis, or of the clean-element bound of Fiat et al. is known. A formalization supplies a precise model of a randomized online algorithm with predictions. It also settles the details the paper leaves implicit: the eviction missing from the clean branch of Algorithm 1 as printed, the cap 2Hk2H_k2Hk​ printed as 2log⁡k2\log k2logk in Lemma 3.4, and the behaviour of the spread at 000.

Difficulty

The obvious argument charges every miss to a chain and bounds each chain separately. That works for chains that follow the predictions, but a chain that switches to random evictions interacts with every other chain of the phase, because all of them evict from the same pool of unmarked elements. A bound on its expected length must hold whatever the other chains evict, including evictions that depend on earlier coin flips. A second difficulty is summing. The chain errors ηr,c\eta_{r,c}ηr,c​ and the chain lengths are both random, while the hypothesis controls only the total error ηℓ(h,σ)\eta_\ell(h,\sigma)ηℓ​(h,σ) against Opt(σ)\mathrm{Opt}(\sigma)Opt(σ), not the number of chains Q(σ)Q(\sigma)Q(σ) in which the error is spread.

Formalization scope

Elements form a type with decidable equality. A request sequence is a list; predictions are one real per request, and every real sequence is allowed. Labels are 1-based next-arrival positions, with n+1n+1n+1 for elements never requested again. The paper prints the label with equal features; the element is meant. Opt\mathrm{Opt}Opt is computed as the minimum over all demand-paging schedules from the empty cache, which loses no generality. HkH_kHk​ is harmonic k as a real number, never log⁡k\log klogk.

Predictive Marker is a PMF over final states. The random eviction of line 21 is uniform over the unmarked cached elements, and ties in the arg max are a parameter quantified universally. The eviction of lines 23–24 is also performed after a clean miss; as printed, it sits only in the stale branch. The expected cost lies in [0,∞][0,\infty][0,∞].

The spread takes real arguments and lengths T≥1T \ge 1T≥1. SSS must be concave on [0,∞)[0,\infty)[0,∞), finite, and at least the spread. It must also be continuous at 000, which the paper does not say: without it the chain lemma fails for losses whose minimal reversed-order loss stays 000 over several lengths. ϵ\epsilonϵ-accuracy is the pointwise condition on the given pair (σ,h)(\sigma, h)(σ,h). The competitive ratio is written as a product, so Opt(σ)=0\mathrm{Opt}(\sigma) = 0Opt(σ)=0 needs no special case. Lemma 3.3 is stated pointwise for chains without random evictions, as its proof shows. Lemma 3.4 has 2Hk2H_k2Hk​ in place of the printed 2log⁡k2\log k2logk, with the minimum inside the expectation because ηr,c\eta_{r,c}ηr,c​ is random.

The statement must not be trivialized. Opt\mathrm{Opt}Opt is the true offline optimum, not Bélády's rule applied to the predictions. The expectation is taken over Predictive Marker's own run, never compared with itself. The spread hypothesis is satisfiable; for example, the constant loss 111 has spread max⁡(1,⌈m⌉)≤m+1\max(1,\lceil m\rceil) \le m + 1max(1,⌈m⌉)≤m+1.

Out of scope: Lemma 3.2 (the special-marking algorithm SM, which enters only through Lemma 3.4's proof); Lemma 3.1 and Corollaries 1–2, whose printed constants are false for small mmm or disagree with Theorem 3.3; the lower bounds of §3.1 and §3.4; the extensions of §4; the experiments of §5; running time and learnability.

Welcome contributions: the Marker phase structure and its equivalence with the combinatorial phases, the clean-element bounds Q/2≤Opt≤QQ/2 \le \mathrm{Opt} \le QQ/2≤Opt≤Q (reusable for any marking algorithm), and a bound on the expected number of misses caused by elements evicted uniformly at random within a phase.

Selected references

  • T. Lykouris, S. Vassilvitskii, Competitive Caching with Machine Learned Advice, arXiv:1802.05399v4, 2020; J. ACM 68(4), 2021. https://arxiv.org/abs/1802.05399v4
  • A. Fiat, R. M. Karp, M. Luby, L. A. McGeoch, D. D. Sleator, N. E. Young, Competitive paging algorithms, J. Algorithms 12(4), 1991. https://doi.org/10.1016/0196-6774(91)90041-V
  • L. A. Bélády, A study of replacement algorithms for a virtual-storage computer, IBM Systems Journal 5(2), 1966. https://doi.org/10.1147/sj.52.0078
  • D. D. Sleator, R. E. Tarjan, Amortized efficiency of list update and paging rules, Comm. ACM 28(2), 1985. https://doi.org/10.1145/2786.2793
  • D. Rohatgi, Near-optimal bounds for online caching with machine learned advice, SODA 2020. https://arxiv.org/abs/1910.12172
  • A. Wei, Better and simpler learning-augmented online caching, APPROX/RANDOM 2020. https://arxiv.org/abs/2005.13716
10 thms1 active userReviewed
CombinatoricsOperations ResearchOptimization+1·Captain: mikedeng1

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

Motivation

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

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

Timeline:

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

Setting

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

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

The proof uses three auxiliary problems:

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

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

Formalization targets

Goal: Theorem 8

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

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

Milestones

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

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

Significance

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

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

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

Difficulty

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

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

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

Formalization scope

Conventions committed to in Lean:

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

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

Out of scope:

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

Infrastructure that a complete development needs:

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

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

Selected references

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

A Polylogarithmic-Competitive Algorithm for the k-Server Problem: Randomized k-Server Is O(log² k · log³ n · log log n)-Competitive on Every n-Point MetricResearch Paper

Motivation

The k-server problem (Manasse, McGeoch and Sleator, 1990) is the central problem of online computation: kkk servers sit on points of a metric space, requests arrive one at a time at points of the space, and each request must be served by moving a server to it, at a cost equal to the distance travelled. An online algorithm decides without knowing future requests; its quality is its competitive ratio, the worst-case ratio between its cost and the cost of an optimal offline schedule. Paging (caching) is the special case of a uniform metric, and weighted paging the case of a weighted star.

Timeline of the upper bounds for general metrics:

  • 1990: Manasse, McGeoch and Sleator prove that every deterministic algorithm has ratio at least kkk and conjecture that kkk is achievable.
  • 1991: Fiat, Rabani and Ravid give the first ratio depending on kkk only (exponential in kkk).
  • 1995: Koutsoupias and Papadimitriou prove that the work function algorithm is (2k−1)(2k-1)(2k−1)-competitive.
  • For randomized algorithms against an oblivious adversary, the conjectured answer is O(log⁡k)O(\log k)O(logk), achieved for paging (Fiat et al., 1991), but until 2011 nothing better than the deterministic 2k−12k-12k−1 was known for general metrics, even when the ratio may depend on the number of points nnn.
  • 2011: Bansal, Buchbinder, Mądry and Naor give the first polylogarithmic bound, O(log⁡2klog⁡3nlog⁡log⁡n)O(\log^2 k\log^3 n\log\log n)O(log2klog3nloglogn) (arXiv:1110.1580; J. ACM 62(5), 2015, DOI 10.1145/2783434), the result of this mission.

Setting

Let (M,dist)(M,\mathrm{dist})(M,dist) be a finite metric space with nnn points and kkk a number of servers. A configuration C:{1,…,k}→MC:\{1,\dots,k\}\to MC:{1,…,k}→M places server iii at C(i)C(i)C(i). A deterministic online algorithm maps each prefix of the request sequence to a configuration that has a server at the last request; its cost on a sequence ρ\rhoρ is the total distance travelled. OPT(C0,ρ)\mathrm{OPT}(C_0,\rho)OPT(C0​,ρ) is the least cost of any schedule serving ρ\rhoρ from the initial configuration C0C_0C0​. A randomized algorithm is a probability distribution over deterministic online algorithms, all starting at C0C_0C0​; it is ccc-competitive if there is a constant aaa such that its expected cost on every request sequence ρ\rhoρ is at most c⋅OPT(C0,ρ)+ac\cdot\mathrm{OPT}(C_0,\rho)+ac⋅OPT(C0​,ρ)+a.

The paper works with three auxiliary objects. A σ-HST is a rooted tree whose leaves are the points, in which all edges from a node to its children have one common length, equal to 1/σ1/\sigma1/σ times the length of the edge above that node; the distance between two leaves is the length of the tree path. A weighted σ-HST only requires that the edge above a non-root internal node be at least σ\sigmaσ times each edge below it. In the fractional k-server problem on a tree, the state is a vector xxx of server probabilities on the leaves with 0≤xi≤10\le x_i\le10≤xi​≤1 and ∑ixi=k\sum_i x_i=k∑i​xi​=k, a request at leaf iii forces xi=1x_i=1xi​=1, and moving from xxx to x′x'x′ costs ∑vW(v) ∣xv′−xv∣\sum_v W(v)\,|x'_v-x_v|∑v​W(v)∣xv′​−xv​∣, where xvx_vxv​ is the mass below node vvv and W(v)W(v)W(v) the length of the edge above vvv. In the allocation problem on a weighted star with weights wiw_iwi​, requests carry a location iti^tit, a monotone cost vector ht(0)≥⋯≥ht(k)≥0h^t(0)\ge\dots\ge h^t(k)\ge0ht(0)≥⋯≥ht(k)≥0 (the cost of serving with jjj servers there) and a server quota κ(t)≤k\kappa(t)\le kκ(t)≤k.

Formalization targets

Goal: Theorem 1

There is a universal constant C>0C>0C>0 such that for all k≥2k\ge2k≥2, every metric space MMM with n≥3n\ge3n≥3 points and every initial configuration C0C_0C0​, some randomized online algorithm starting at C0C_0C0​ is

C log⁡2k log⁡3n log⁡log⁡n-competitive.C\,\log^2 k\,\log^3 n\,\log\log n\text{-competitive.}Clog2klog3nloglogn-competitive.

Milestones

In the order the proof uses them:

  1. Claim 15: the fix-stage inequality behind the allocation algorithm's analysis.
  2. Theorem 5: for every 0<ε≤10<\varepsilon\le10<ε≤1, a fractional allocation algorithm whose hit cost is at most (1+ε)(Opt+wmax⁡g(κ))+a(1+\varepsilon)(\mathrm{Opt}+w_{\max}g(\kappa))+a(1+ε)(Opt+wmax​g(κ))+a and whose movement cost is at most O(log⁡(k/ε))(Opt+wmax⁡g(κ))+aO(\log(k/\varepsilon))(\mathrm{Opt}+w_{\max}g(\kappa))+aO(log(k/ε))(Opt+wmax​g(κ))+a, where g(κ)=∑t∣κ(t)−κ(t−1)∣g(\kappa)=\sum_t|\kappa(t)-\kappa(t-1)|g(κ)=∑t​∣κ(t)−κ(t−1)∣.
  3. Theorem 6: given such allocation algorithms, an O(ℓlog⁡(kℓ))O(\ell\log(k\ell))O(ℓlog(kℓ))-competitive fractional k-server algorithm on every weighted σ-HST of depth ℓ\ellℓ with σ=Ω(ℓlog⁡(kℓ))\sigma=\Omega(\ell\log(k\ell))σ=Ω(ℓlog(kℓ)).
  4. Theorem 8: every σ-HST with nnn leaves becomes a weighted σ-HST of depth O(log⁡n)O(\log n)O(logn) on the same leaves, with distances distorted by at most 2σ/(σ−1)2\sigma/(\sigma-1)2σ/(σ−1).
  5. Lemma 25 and Theorem 24: on a σ-HST with σ>5\sigma>5σ>5, randomized states consistent with a changing fractional state can be maintained online at cost O(ct)O(c_t)O(ct​) per step.
  6. Theorem 7: on a σ-HST with σ>5\sigma>5σ>5, a ccc-competitive fractional algorithm yields an O(c)O(c)O(c)-competitive randomized one.

Significance

The theorem broke the exponential gap between the Ω(log⁡k)\Omega(\log k)Ω(logk) lower bound and the 2k−12k-12k−1 upper bound for randomized k-server, and showed that randomization helps on every finite metric, not only on uniform or specially structured ones. Its two-level method (a fractional algorithm on trees driven by per-node allocation problems, followed by an online rounding) became the template for later work, including the O(log⁡2k)O(\log^2 k)O(log2k) bound on HSTs of Bubeck, Cohen, Lee, Lee and Mądry (STOC 2018) and Lee's O(log⁡6k)O(\log^6 k)O(log6k) bound on general metrics (FOCS 2018).

The result is proved, in this paper. As far as is known it has no machine-checked proof. Formalizing it means formalizing the analysis of an online algorithm driven by a continuous-time process, a potential-function argument with exact constants, a tree contraction with a distortion bound, and an online randomized rounding against a transportation cost. The allocation, HST and rounding statements are reusable for other online problems on trees (metrical task systems, weighted paging).

Difficulty

For a deterministic or randomized algorithm on a tree, the natural recursion splits the servers of each node among its children. Coté, Meyerson and Poplawski showed that this works if each node solves an allocation problem with a strong guarantee: hit cost within a factor 1+ε1+\varepsilon1+ε of optimal. Integral allocation algorithms cannot achieve this; the integrality gap example of the paper (p. 8) gives a factor Ω(k)\Omega(k)Ω(k). The fractional relaxation avoids the gap, but then the rounding step must keep a randomized state consistent with a fractional state at constant-factor cost, and the HSTs obtained from general metrics have depth growing with the aspect ratio, which a depth-dependent ratio cannot afford. Each of the three reductions (allocation to fractional k-server, deep HST to shallow weighted HST, fractional to randomized) loses only polylogarithmic or constant factors, and the main theorem needs all three at once.

Formalization scope

The k-server model, randomized algorithms and competitiveness are the published definitions KServer_model and KServer_randomized; competitiveness carries an additive constant fixed before the request sequence. Trees are finite rooted trees with a parent map, a depth function and positive edge lengths; points of the k-server problem are the leaves, and the theorems take an arbitrary finite metric space together with a bijection to the leaves and the hypothesis that the distance equals the tree distance. Fractional k-server states have exactly kkk units of mass, each leaf at most 111, and fractional algorithms are measured against the integral offline optimum. The allocation optimum is the integral optimum; cost vectors are finite, non-negative and non-increasing; the diameter of the star is wmax⁡=max⁡iwiw_{\max}=\max_i w_iwmax​=maxi​wi​. The cost of changing a randomized state is the transportation cost over couplings, with minimum-matching cost between configurations. Every O(⋅)O(\cdot)O(⋅) is an explicit constant quantified before the instance, except that in Theorems 7 and 24 and Lemma 25 it may depend on σ\sigmaσ.

Formalizations that make the targets trivial are excluded: the fractional state must place a full server on every request and stay in [0,1][0,1][0,1], the benchmark is the integral optimum (not the algorithm's own or the fractional cost), and no constant may depend on kkk, nnn, the metric or the tree, since otherwise Theorem 1 would follow from the 2k−12k-12k−1 bound.

The proof of Theorem 1 also uses the embedding of Fakcharoenphol, Rao and Talwar [18] of a finite metric into a distribution over σ-HSTs with expected distortion O(σlog⁡σn)O(\sigma\log_\sigma n)O(σlogσ​n). It is an external ingredient, not a result of this paper, and is not a milestone; contributions formalizing it (or Bartal's earlier embedding) are welcome, as are formalizations of the integral optimum's properties on trees (Lemmas 21–22 of the paper), which are not stated here.

Selected references

  • N. Bansal, N. Buchbinder, A. Mądry, J. Naor, A Polylogarithmic-Competitive Algorithm for the k-Server Problem, arXiv:1110.1580v1, 2011; J. ACM 62(5), 2015. https://arxiv.org/abs/1110.1580, https://doi.org/10.1145/2783434
  • M. Manasse, L. McGeoch, D. Sleator, Competitive algorithms for server problems, J. Algorithms 11, 1990. https://doi.org/10.1016/0196-6774(90)90003-W
  • E. Koutsoupias, C. Papadimitriou, On the k-server conjecture, J. ACM 42(5), 1995. https://doi.org/10.1145/210118.210128
  • A. Fiat, R. Karp, M. Luby, L. McGeoch, D. Sleator, N. Young, Competitive paging algorithms, J. Algorithms 12, 1991. https://doi.org/10.1016/0196-6774(91)90041-V
  • J. Fakcharoenphol, S. Rao, K. Talwar, A tight bound on approximating arbitrary metrics by tree metrics, J. Comput. Syst. Sci. 69(3), 2004. https://doi.org/10.1016/j.jcss.2004.04.011
  • A. Coté, A. Meyerson, L. Poplawski, Randomized k-server on hierarchical binary trees, STOC 2008. https://doi.org/10.1145/1374376.1374474
14 thms3 active usersReviewed
Dynamic ProgrammingOperations ResearchProbability·Captain: mikedeng1

Uniformly Bounded Regret in the Multi-Secretary Problem 1: The Budget-Ratio Policy Has Regret at Most a₁M(ε), Uniformly in the Number of Candidates n and the Budget kResearch Paper

Motivation

The multi-secretary problem is the simplest model of capacity allocation under uncertainty: a decision maker sees nnn candidates one at a time and may hire at most kkk of them, with every decision final. The same structure underlies single-resource revenue management (accepting or rejecting booking requests against a fixed inventory; see Talluri and van Ryzin, The Theory and Practice of Revenue Management, 2004), online knapsack and packing problems, and dynamic assortment of limited stock.

The performance of an online policy is measured against the offline benchmark, the value of the best kkk candidates chosen with full hindsight. The gap between the two is the regret.

  • In the version where the values arrive as a uniform random permutation, Kleinberg (2005) proved that the minimal regret is of order k\sqrt kk​ and gave an algorithm attaining it (as summarized in Remark 1 of the paper below).
  • Arlotto and Gurvich (arXiv:1710.07719, 2017; Stochastic Systems 2019) showed that when the values have a finite support, the optimal online policy, and an explicit simple policy, have regret bounded by a constant that does not depend on nnn or kkk. The constant depends only on the smallest probability mass.

This mission formalizes that upper bound.

Setting

Abilities take values in a finite set A={am<am−1<⋯<a1}\mathcal A=\{a_m<a_{m-1}<\dots<a_1\}A={am​<am−1​<⋯<a1​} of distinct positive reals, with probabilities fj=P(X=aj)>0f_j=\mathbb P(X=a_j)>0fj​=P(X=aj​)>0, ∑jfj=1\sum_j f_j=1∑j​fj​=1. Write Fˉ(aj)=f1+⋯+fj−1\bar F(a_j)=f_1+\dots+f_{j-1}Fˉ(aj​)=f1​+⋯+fj−1​ for the mass strictly above aja_jaj​, and

ϵ=12min⁡{fm,…,f1}.\epsilon=\tfrac12\min\{f_m,\dots,f_1\}.ϵ=21​min{fm​,…,f1​}.

The abilities X1,…,XnX_1,\dots,X_nX1​,…,Xn​ are independent with this distribution. Budget pairs range over the triangle T={(n,k):0≤k≤n}\mathcal T=\{(n,k):0\le k\le n\}T={(n,k):0≤k≤n}.

  • Offline value. Voff∗(n,k)=E[max⁡{∑tXtσt:σ∈{0,1}n, ∑tσt≤k}]V^*_{\mathrm{off}}(n,k)=\mathbb E\big[\max\{\sum_t X_t\sigma_t:\sigma\in\{0,1\}^n,\ \sum_t\sigma_t\le k\}\big]Voff∗​(n,k)=E[max{∑t​Xt​σt​:σ∈{0,1}n, ∑t​σt​≤k}].
  • Online policies. A policy decides σt∈{0,1}\sigma_t\in\{0,1\}σt​∈{0,1} using only X1,…,XtX_1,\dots,X_tX1​,…,Xt​ and must select at most kkk candidates on every realization. Π(n,k)\Pi(n,k)Π(n,k) is the set of such policies, Vonπ(n,k)=E[∑tXtσtπ]V^\pi_{\mathrm{on}}(n,k)=\mathbb E[\sum_t X_t\sigma^\pi_t]Vonπ​(n,k)=E[∑t​Xt​σtπ​], and Von∗(n,k)=max⁡π∈Π(n,k)Vonπ(n,k)V^*_{\mathrm{on}}(n,k)=\max_{\pi\in\Pi(n,k)}V^\pi_{\mathrm{on}}(n,k)Von∗​(n,k)=maxπ∈Π(n,k)​Vonπ​(n,k).
  • Counts. ZjrZ^r_jZjr​ is the number of aja_jaj​-candidates among the first rrr. The offline sort selects Sjr=min⁡{Zjr,(k−∑i<jZir)+}\mathfrak S^r_j=\min\{Z^r_j,(k-\sum_{i<j}Z^r_i)_+\}Sjr​=min{Zjr​,(k−∑i<j​Zir​)+​} of them. Sjπ,rS^{\pi,r}_jSjπ,r​ counts those selected by π\piπ.
  • Action index. j0(n,k)j_0(n,k)j0​(n,k) is the largest jjj with Fˉ(aj)+12fj≤k/n\bar F(a_j)+\tfrac12f_j\le k/nFˉ(aj​)+21​fj​≤k/n, or 111 if there is none.
  • Thresholds. T1=0T_1=0T1​=0, Tj=Fˉ(aj)+12fjT_j=\bar F(a_j)+\tfrac12 f_jTj​=Fˉ(aj​)+21​fj​ for 2≤j≤m2\le j\le m2≤j≤m, and Tm+1=+∞T_{m+1}=+\inftyTm+1​=+∞.
  • Budget-Ratio (BR) policy. With remaining budget KtK_tKt​ (K0=kK_0=kK0​=k), at time t+1t+1t+1 the policy finds jjj with Tj≤Kt/(n−t)<Tj+1T_j\le K_t/(n-t)<T_{j+1}Tj​≤Kt​/(n−t)<Tj+1​. It selects Xt+1X_{t+1}Xt+1​ if and only if Kt>0K_t>0Kt​>0 and Xt+1≥ajX_{t+1}\ge a_jXt+1​≥aj​.
  • Stopping times. For 0<δ<ϵ0<\delta<\epsilon0<δ<ϵ, τ0\tau_0τ0​ is the first time the budget ratio comes within δ/2\delta/2δ/2 of a threshold, or the cut-off n−2δ−1−1n-2\delta^{-1}-1n−2δ−1−1. The time τ\tauτ of (20) is the first later time the ratio leaves the δ\deltaδ-band around that threshold, or the cut-off.

Formalization targets

Goal: Theorem 1 (first display)

For every ϵ>0\epsilon>0ϵ>0 there is a constant MMM such that for every instance with 12min⁡jfj=ϵ\tfrac12\min_jf_j=\epsilon21​minj​fj​=ϵ and all (n,k)∈T(n,k)\in\mathcal T(n,k)∈T, br∈Π(n,k)\mathrm{br}\in\Pi(n,k)br∈Π(n,k) and

Voff∗(n,k)−Von∗(n,k)≤Voff∗(n,k)−Vonbr(n,k)≤a1M.V^*_{\mathrm{off}}(n,k)-V^*_{\mathrm{on}}(n,k)\le V^*_{\mathrm{off}}(n,k)-V^{\mathrm{br}}_{\mathrm{on}}(n,k)\le a_1M.Voff∗​(n,k)−Von∗​(n,k)≤Voff∗​(n,k)−Vonbr​(n,k)≤a1​M.

No constant is fixed. Only the shape is asserted: a bound uniform in nnn, kkk, the support size and the distribution, given ϵ\epsilonϵ.

Milestones, in the order the proof uses them

  • The benchmark inequality Vonπ≤Voff∗V^\pi_{\mathrm{on}}\le V^*_{\mathrm{off}}Vonπ​≤Voff∗​ (p. 5).
  • The sort identity Voff∗=∑jajE[Sjn]V^*_{\mathrm{off}}=\sum_ja_j\mathbb E[\mathfrak S^n_j]Voff∗​=∑j​aj​E[Sjn​] (4).
  • The binomial overshoot bound E[(B−k)+]≤1/(4ε)\mathbb E[(B-k)_+]\le1/(4\varepsilon)E[(B−k)+​]≤1/(4ε) (Lemma 2).
  • The offline decomposition Voff∗=∑i<jaiE[Zin]+ajE[Sjn]+aj+1E[Sj+1n]±a1/(4ϵ)V^*_{\mathrm{off}}=\sum_{i<j}a_i\mathbb E[Z^n_i]+a_j\mathbb E[\mathfrak S^n_j]+a_{j+1}\mathbb E[\mathfrak S^n_{j+1}]\pm a_1/(4\epsilon)Voff∗​=∑i<j​ai​E[Zin​]+aj​E[Sjn​]+aj+1​E[Sj+1n​]±a1​/(4ϵ) (Proposition 1).
  • The sufficient condition: four properties (i)–(iv) of a policy up to a stopping time imply regret at most 3a1M+a1/(4ϵ)3a_1M+a_1/(4\epsilon)3a1​M+a1​/(4ϵ) (Proposition 2).
  • The identification j0(n,k)=jj_0(n,k)=jj0​(n,k)=j on k/n∈[Tj,Tj+1)k/n\in[T_j,T_{j+1})k/n∈[Tj​,Tj+1​) (p. 17).
  • The BR selection probability and the jump bound ∣Kt/(n−t)−Kt+1/(n−t−1)∣≤δ/2|K_t/(n-t)-K_{t+1}/(n-t-1)|\le\delta/2∣Kt​/(n−t)−Kt+1​/(n−t−1)∣≤δ/2 (p. 13).
  • E[τ]≥n−M\mathbb E[\tau]\ge n-ME[τ]≥n−M (Theorem 2).
  • BR and τ\tauτ satisfy (i)–(iv) (Corollary 1).
  • The state-space reduction vℓ(w,κ)=w+gℓ(κ)v_\ell(w,\kappa)=w+g_\ell(\kappa)vℓ​(w,κ)=w+gℓ​(κ) of the Bellman recursion (Proposition 5).

Significance

The result. Bounded regret means that the loss from not knowing the future is a fixed number of candidates' worth of value, however long the horizon and however large the budget. The bound holds uniformly over all distributions with the same ϵ\epsilonϵ. It is attained by an explicit, adaptive, non-randomized rule that compares one ratio with mmm fixed thresholds. The companion result of the same paper shows that every non-adaptive policy suffers regret of order n\sqrt nn​ in the interior regime. Together they quantify the value of adapting to the remaining budget. Lemma 1 of the paper shows the dependence on ϵ\epsilonϵ cannot be removed.

Formalizing it. The result is proved in the paper, but no part of it is machine-checked; there is no multi-secretary or bounded-regret development on the platform. The mission produces several pieces of machinery: a reusable finite model of sequential selection with online policies and the offline benchmark; an explicit online policy with its stopping-time analysis; and a binomial overshoot bound usable elsewhere. The constant MMM is not made explicit in the paper. A formal proof would give one, and sharper constants are welcome.

Difficulty

The offline decomposition and the sufficient condition are bookkeeping with counts and one concentration bound. The hard step is Theorem 2: showing that the budget ratio Kt/(n−t)K_t/(n-t)Kt​/(n−t) stays within δ\deltaδ of its attracting threshold until a bounded expected number of periods before the end. Near the horizon a single selection moves the ratio by about 1/(n−t)1/(n-t)1/(n−t), so the band becomes easy to leave. Equivalently, the target δ(n−τ0−u)\delta(n-\tau_0-u)δ(n−τ0​−u) that the deviation process must exceed shrinks to zero. A standard martingale or drift argument with a fixed band therefore does not give a bound uniform in nnn. The paper combines the mean-reverting drift of the deviation process with an exponential tail bound (its Proposition 4) and a Lyapunov argument. A second subtlety is uniformity: every constant must depend on ϵ\epsilonϵ (and δ\deltaδ) only, never on mmm, the aja_jaj​, nnn or kkk.

Formalization scope

The source is arXiv:1710.07719v2; its printed page numbers equal the PDF page numbers.

Representation.

  • Ability levels are Fin m, with index 0 the largest value a1a_1a1​; Lean index iii is the paper's i+1i+1i+1.
  • Each instance carries aaa strictly decreasing and positive, fff positive with ∑f=1\sum f=1∑f=1.
  • Expectations are finite sums over sequences x:Fin n→Fin mx:\mathrm{Fin}\,n\to\mathrm{Fin}\,mx:Finn→Finm weighted by ∏tf(xt)\prod_tf(x_t)∏t​f(xt​), so no measure theory is needed.
  • Policies are deterministic selection rules σ(x,t)\sigma(x,t)σ(x,t) that are non-anticipating and feasible. Von∗V^*_{\mathrm{on}}Von∗​ is a maximum over this finite set. The paper allows randomized policies; for this finite problem the optimal values coincide (p. 39). In any case, restricting to deterministic policies can only lower Von∗V^*_{\mathrm{on}}Von∗​ and so does not weaken the goal.
  • Voff∗V^*_{\mathrm{off}}Voff∗​ is defined as an expected maximum over selection vectors, not by the sort formula. The sort formula is a milestone.

Quantifiers. The constant MMM in the goal is chosen after ϵ\epsilonϵ and before mmm, the instance, nnn and kkk. A statement with MMM chosen after the instance, or after nnn, is trivial (regret ≤a1n\le a_1n≤a1​n) and is excluded.

Corrections to the printed text, disclosed in the items.

  1. In Theorem 2 and Corollary 1, MMM depends on the auxiliary δ∈(0,ϵ)\delta\in(0,\epsilon)δ∈(0,ϵ) as well, because τ\tauτ does. δ\deltaδ is quantified before MMM. The goal itself is δ\deltaδ-free.
  2. Lemma 2's conditions p+ε≤k/np+\varepsilon\le k/np+ε≤k/n, k/n≤p−εk/n\le p-\varepsilonk/n≤p−ε are stated as (p+ε)n≤k(p+\varepsilon)n\le k(p+ε)n≤k, k≤(p−ε)nk\le(p-\varepsilon)nk≤(p−ε)n, the form used in its proof. This avoids a false case at n=0n=0n=0.
  3. The BR rule is applied at every time t+1∈{1,…,n}t+1\in\{1,\dots,n\}t+1∈{1,…,n}; p. 11 writes {1,…,n−1}\{1,\dots,n-1\}{1,…,n−1}.
  4. τ\tauτ is capped at nnn, which matters only when n=0n=0n=0.
  5. In Proposition 5 the recursions are imposed for κ≥1\kappa\ge1κ≥1 (boundary conditions at κ=0\kappa=0κ=0), and only identity (49) is stated.

Infrastructure. The model definitions (instance, offline value, online policies, counts, thresholds, action index) and the binomial overshoot lemma are reusable for other finite-support online selection and revenue-management results. All of the following are welcome:

  • proofs of individual milestones;
  • an explicit constant;
  • a formal derivation of Von∗(n,k)=vn(0,k)V^*_{\mathrm{on}}(n,k)=v_n(0,k)Von∗​(n,k)=vn​(0,k) connecting Proposition 5 to Von∗V^*_{\mathrm{on}}Von∗​.

Selected references

  • A. Arlotto, I. Gurvich, Uniformly Bounded Regret in the Multi-Secretary Problem, arXiv:1710.07719v2, 2018; Stochastic Systems 9(3), 2019. https://arxiv.org/abs/1710.07719
  • R. Kleinberg, A multiple-choice secretary algorithm with applications to online auctions, SODA 2005. https://dl.acm.org/doi/10.5555/1070432.1070519
  • K. T. Talluri, G. J. van Ryzin, The Theory and Practice of Revenue Management, Springer, 2004. https://doi.org/10.1007/b139000
  • D. P. Bertsekas, S. E. Shreve, Stochastic Optimal Control: The Discrete-Time Case, Academic Press, 1978.
  • S. Boucheron, G. Lugosi, P. Massart, Concentration Inequalities, Oxford University Press, 2013. https://doi.org/10.1093/acprof:oso/9780199535255.001.0001
14 thms2 active usersReviewed
Operations ResearchProbability·Captain: mikedeng1

Uniformly Bounded Regret in the Multi-Secretary Problem 2: When (f₁+ε)n ≤ k ≤ (1−fₘ−ε)n, Every Non-Adaptive Policy Has Regret at Least M√nResearch Paper

Motivation

The multi-secretary problem is the basic model of selecting under a budget from a stream of offers. Hiring a fixed number of candidates, accepting a fixed number of requests for a perishable resource, and admitting customers into a capacity-limited service all have this structure. Each item must be accepted or rejected on arrival, and the comparison point is the offline decision maker, who sees the whole sequence and keeps the best kkk items. The gap between the two expected values is the regret.

A common class of heuristics in revenue management and online resource allocation does not react to the realised history. These policies fix in advance, period by period, a probability of accepting each type of item, and then follow it until the budget runs out; static bid-price and randomised-acceptance rules are of this kind (Talluri and van Ryzin 2004). Arlotto and Gurvich (arXiv:1710.07719v2, Theorem 1) show that when abilities take finitely many values, an adaptive policy has regret bounded uniformly in the horizon nnn and the budget kkk. Their Theorem 3 shows that the restriction to non-adaptive policies costs order n\sqrt nn​ over a wide range of budgets. Read together, these two results separate adaptive from non-adaptive control by an unbounded factor. This mission formalizes the non-adaptive half.

Setting

There are nnn candidates with abilities X1,…,XnX_1,\dots,X_nX1​,…,Xn​, independent and identically distributed on mmm values 0<am<am−1<⋯<a10<a_m<a_{m-1}<\dots<a_10<am​<am−1​<⋯<a1​, with masses fj=P(X1=aj)>0f_j=\mathbb P(X_1=a_j)>0fj​=P(X1​=aj​)>0 and ∑jfj=1\sum_jf_j=1∑j​fj​=1. Write ϵ=12min⁡jfj\epsilon=\tfrac12\min_jf_jϵ=21​minj​fj​ and Fˉ(aj)=f1+⋯+fj−1\bar F(a_j)=f_1+\dots+f_{j-1}Fˉ(aj​)=f1​+⋯+fj−1​. The budget is kkk, with 0≤k≤n0\le k\le n0≤k≤n.

The offline value is

Voff∗(n,k)=E[max⁡{∑tXtσt:σ∈{0,1}n, ∑tσt≤k}].V^*_{\mathrm{off}}(n,k)=\mathbb E\Big[\max\Big\{\textstyle\sum_tX_t\sigma_t:\sigma\in\{0,1\}^n,\ \sum_t\sigma_t\le k\Big\}\Big].Voff∗​(n,k)=E[max{∑t​Xt​σt​:σ∈{0,1}n, ∑t​σt​≤k}].

A non-adaptive policy is a matrix π={pj,t∈[0,1]}\pi=\{p_{j,t}\in[0,1]\}π={pj,t​∈[0,1]}. At time ttt, if budget remains and Xt=ajX_t=a_jXt​=aj​, the candidate is selected with probability pj,tp_{j,t}pj,t​, independently of everything else. The selection coins BtB_tBt​ are then independent Bernoulli variables with qt=E[Bt]=∑jpj,tfjq_t=\mathbb E[B_t]=\sum_jp_{j,t}f_jqt​=E[Bt​]=∑j​pj,t​fj​. The policy selects until kkk coins have come up. Its value Vonπ(n,k)V^\pi_{\mathrm{on}}(n,k)Vonπ​(n,k) is the expected total ability selected, and

Vna∗(n,k)=sup⁡πVonπ(n,k).V^*_{\mathrm{na}}(n,k)=\sup_{\pi}V^\pi_{\mathrm{on}}(n,k).Vna∗​(n,k)=πsup​Vonπ​(n,k).

The deterministic relaxation replaces the random counts Zjn=#{t:Xt=aj}Z^n_j=\#\{t:X_t=a_j\}Zjn​=#{t:Xt​=aj​} by their means. Its value is

DR(n,k)=max⁡{∑jajsj:0≤sj≤nfj, ∑jsj≤k},DR(n,k)=\max\Big\{\textstyle\sum_ja_js_j:0\le s_j\le nf_j,\ \sum_js_j\le k\Big\},DR(n,k)=max{∑j​aj​sj​:0≤sj​≤nfj​, ∑j​sj​≤k},

with solution sj∗=min⁡{nfj,(k−nFˉ(aj))+}s^*_j=\min\{nf_j,(k-n\bar F(a_j))_+\}sj∗​=min{nfj​,(k−nFˉ(aj​))+​}. The index policy takes its probabilities from s∗s^*s∗: pj,t=sj∗/(nfj)p_{j,t}=s^*_j/(nf_j)pj,t​=sj∗​/(nfj​).

Formalization targets

Goal: Theorem 3 (p. 25)

For every ϵ>0\epsilon>0ϵ>0, mmm and aaa there is M=M(ϵ,m,a)>0M=M(\epsilon,m,a)>0M=M(ϵ,m,a)>0 such that, for all masses with 12min⁡jfj=ϵ\tfrac12\min_jf_j=\epsilon21​minj​fj​=ϵ and all (n,k)(n,k)(n,k) with (f1+ϵ)n≤k≤(1−fm−ϵ)n(f_1+\epsilon)n\le k\le(1-f_m-\epsilon)n(f1​+ϵ)n≤k≤(1−fm​−ϵ)n,

Mn≤Voff∗(n,k)−Vna∗(n,k).M\sqrt n\le V^*_{\mathrm{off}}(n,k)-V^*_{\mathrm{na}}(n,k).Mn​≤Voff∗​(n,k)−Vna∗​(n,k).

The constant does not depend on the masses beyond ϵ\epsilonϵ, nor on nnn or kkk.

Milestones

  • Lemma 2 (p. 8): binomial overshoot, E[(B−k)+]≤1/(4ε)\mathbb E[(B-k)_+]\le1/(4\varepsilon)E[(B−k)+​]≤1/(4ε) when kkk exceeds the mean by εn\varepsilon nεn, and the symmetric bound.
  • Remark 2 (pp. 10–11): s∗s^*s∗ solves the relaxation, and Voff∗≤DRV^*_{\mathrm{off}}\le DRVoff∗​≤DR.
  • Lemma 3 (p. 25): the index policy satisfies DR−Vnaid≤ε−1a1nDR-V^{\mathrm{id}}_{\mathrm{na}}\le\varepsilon^{-1}a_1\sqrt nDR−Vnaid​≤ε−1a1​n​ when k/n≥εk/n\ge\varepsilonk/n≥ε, so the order n\sqrt nn​ is attained.
  • Lemma 5 (p. 26): for a centred Bernoulli sum with variance ς2\varsigma^2ς2, E[(±N−Υς)+]≥β1ς−(2+32)\mathbb E[(\pm N-\Upsilon\varsigma)_+]\ge\beta_1\varsigma-(2+3\sqrt2)E[(±N−Υς)+​]≥β1​ς−(2+32​) with β1(Υ)>0\beta_1(\Upsilon)>0β1​(Υ)>0, and E[(N+Υς)+2]≤β2ς2\mathbb E[(N+\Upsilon\varsigma)_+^2]\le\beta_2\varsigma^2E[(N+Υς)+2​]≤β2​ς2.
  • Lemma 7 (p. 27): an optimal non-adaptive policy exists, and any optimal one has f1/2≤qt≤1−fm/2f_1/2\le q_t\le1-f_m/2f1​/2≤qt​≤1−fm​/2 outside 2Mn2M\sqrt n2Mn​ periods, so ∑tqt(1−qt)≥f1fm4(n−2Mn)\sum_tq_t(1-q_t)\ge\tfrac{f_1f_m}4(n-2M\sqrt n)∑t​qt​(1−qt​)≥4f1​fm​​(n−2Mn​).
  • Lemma 4 (p. 25): for k≤n(f1−ϵ)k\le n(f_1-\epsilon)k≤n(f1​−ϵ) the non-adaptive regret is at most a2/(4ϵ)a_2/(4\epsilon)a2​/(4ϵ).
  • Lemma 8 and Proposition 6 (p. 40): E[Sjn]=sj∗±Mn\mathbb E[\mathfrak S^n_j]=s^*_j\pm M\sqrt nE[Sjn​]=sj∗​±Mn​, and 0≤DR−Voff∗≤Mn0\le DR-V^*_{\mathrm{off}}\le M\sqrt n0≤DR−Voff∗​≤Mn​ in general and ≤a1m/(4ϵ′)\le a_1m/(4\epsilon')≤a1​m/(4ϵ′) when k/nk/nk/n is ϵ′\epsilon'ϵ′ away from the jump points of Fˉ\bar FFˉ.

Significance

Theorem 3 is the lower half of the separation in Theorem 1 of the paper. The Budget-Ratio policy and the dynamic-programming policy have regret O(1)O(1)O(1), uniformly in (n,k)(n,k)(n,k), while every non-adaptive policy has regret Ω(n)\Omega(\sqrt n)Ω(n​) when k/nk/nk/n lies strictly between f1f_1f1​ and 1−fm1-f_m1−fm​. The order n\sqrt nn​ of fluid and static randomised policies is therefore a property of the whole class, not of a poor choice inside it. Lemma 4 shows that the budget range cannot be removed: with a small budget a non-adaptive policy is as good as any.

The result is proved in the source but has not been machine-checked. A complete development would formalize, inside one finite probabilistic model: the binomial overshoot bound, a uniform anti-concentration estimate for Bernoulli sums, the structure of optimal non-adaptive policies, and the comparison with the offline sort. The source's proof of Theorem 3 also relies on a lemma that fails as printed (see Formalization scope), so a formal proof would close a real gap in the published argument.

Difficulty

The upper bound of order n\sqrt nn​ (Lemma 3) follows from a variance computation. The lower bound must hold for every non-adaptive policy, including time-varying ones, and the obvious argument does not cover them. That argument compares a policy with the index policy and shows the index policy loses n\sqrt nn​. A policy can, however, differ from the index policy by order n\sqrt nn​ in its expected selection counts and still have regret of the same order. The step "small regret forces sj(π)≈sj∗s_j(\pi)\approx s^*_jsj​(π)≈sj∗​", which the source uses, is exactly the step that fails.

What has to be shown is that the selection count ∑tBt\sum_tB_t∑t​Bt​ of an optimal policy fluctuates by order n\sqrt nn​, uniformly in the policy. A policy that runs out of budget early then misses top-value candidates late in the horizon, and one that keeps budget wastes slots. Both effects must be bounded below by a multiple of n\sqrt nn​ that is uniform over all masses with the same ϵ\epsilonϵ. Lemma 5 needs a normal approximation with an explicit, qqq-independent error. Lemma 7 needs the existence of an optimal policy, which is a maximisation over a continuum of matrices.

Formalization scope

The source is the arXiv preprint arXiv:1710.07719v2 (1 June 2018). Its printed page numbers equal the PDF page numbers.

  • Indices. The value and mass vectors are a f : Fin m → ℝ. Lean index jjj is the paper's index j+1j+1j+1, so a 0 =a1=a_1=a1​ is the largest value and f (Fin.rev 0) =fm=f_m=fm​ is the mass of the smallest. The standing assumptions of Sec. 2 are IsValues a (strictly decreasing, positive) and IsMasses f (positive, summing to one).
  • Expectations. All expectations are finite sums over outcome sequences. For the offline problem these are x:Fin n→Fin mx:\mathrm{Fin}\,n\to\mathrm{Fin}\,mx:Finn→Finm with weight ∏tfxt\prod_tf_{x_t}∏t​fxt​​. For a non-adaptive policy they are pairs (Xt,Bt)(X_t,B_t)(Xt​,Bt​) with weight ∏tfxt pxt,tbt(1−pxt,t)1−bt\prod_tf_{x_t}\,p_{x_t,t}^{b_t}(1-p_{x_t,t})^{1-b_t}∏t​fxt​​pxt​,tbt​​(1−pxt​,t​)1−bt​. No measure theory is used.
  • Selection rule. A candidate is selected iff its coin is 111 and fewer than kkk earlier coins were 111. This equals the paper's "up to the stopping time ν\nuν" for k≥1k\ge1k≥1. At k=0k=0k=0 the printed ν=1\nu=1ν=1 would allow a selection without budget, and the feasible rule is used.
  • Suprema. Vna∗V^*_{\mathrm{na}}Vna∗​ is a supremum over all matrices with entries in [0,1][0,1][0,1], not over 0/10/10/1 matrices or the index policy alone. DRDRDR is the supremum of its linear program; it is not defined by the formula ∑jajsj∗\sum_ja_js^*_j∑j​aj​sj∗​, which is a milestone.
  • Index policy. jidj_{\mathrm{id}}jid​ is the largest index with Fˉ(ajid)≤k/n\bar F(a_{j_{\mathrm{id}}})\le k/nFˉ(ajid​​)≤k/n. As printed the defining inequality has no solution at k=nk=nk=n.
  • Constants. Each constant is quantified after (ϵ,m,a)(\epsilon,m,a)(ϵ,m,a) and before (f,n,k)(f,n,k)(f,n,k). The goal's MMM and Lemma 5's β1\beta_1β1​ are strictly positive; with M=0M=0M=0 the goal would reduce to Vna∗≤Voff∗V^*_{\mathrm{na}}\le V^*_{\mathrm{off}}Vna∗​≤Voff∗​. Theorem 3 is posed for all nnn in the range, as printed, without a threshold on nnn. Lemma 2 is stated in the multiplied form (p+ε)n≤k(p+\varepsilon)n\le k(p+ε)n≤k of its proof. Lemma 4 adds m≥2m\ge2m≥2, so that a2a_2a2​ exists.
  • Disclosed gaps in the source. The source's proof of Theorem 3 relies on a lemma that fails as printed (Lemma 6, p. 26), so Lemma 6 is not part of this mission. The statement of Theorem 3 is posed as in the source. The printed argument for the second inequality of Lemma 7's (36) does not go through, and a corrected one also uses am−1a_{m-1}am−1​. Lemma 7's constant is therefore quantified after all of aaa.

Contributions of any kind are welcome. Reusable pieces include binomial overshoot bounds, anti-concentration for sums of independent Bernoulli variables (for example via a Wasserstein normal approximation, which Mathlib lacks), and compactness arguments for optimal randomised policies.

Selected references

  • A. Arlotto, I. Gurvich, Uniformly Bounded Regret in the Multi-Secretary Problem, arXiv:1710.07719v2, 2018; Stochastic Systems 9(3), 2019. https://arxiv.org/abs/1710.07719v2
  • K. T. Talluri, G. J. van Ryzin, The Theory and Practice of Revenue Management, Springer, 2004. https://doi.org/10.1007/b139000
  • N. Ross, Fundamentals of Stein's method, Probability Surveys 8, 2011. https://doi.org/10.1214/11-PS182
  • S. Boucheron, G. Lugosi, P. Massart, Concentration Inequalities, Oxford University Press, 2013. https://doi.org/10.1093/acprof:oso/9780199535255.001.0001
11 thms1 active userReviewed
Linear OptimizationOperations ResearchTheoretical Computer Science·Captain: mikedeng1

Online Primal-Dual Algorithms for Covering and Packing 1: The Online Fractional Packing Scheme Is B-Competitive and Violates Each Packing Constraint by at Most 2 log(1 + n·a_i(max)/a_i(min))/BResearch Paper

Motivation

Many resource-allocation problems arrive one request at a time and must be answered immediately: a bandwidth request is admitted or refused when it appears, an advertiser's budget is charged when a query arrives, a job is accepted before later jobs are seen. Their linear-programming relaxations are packing problems: maximize a total profit subject to capacity constraints, where the variables are revealed online and each must be set irrevocably when it is revealed. A standard yardstick for such an online algorithm is its competitive ratio, the worst-case ratio between the offline optimum and the algorithm's value.

Buchbinder and Naor (Math. Oper. Res. 2009) gave a single online primal-dual scheme for the general online fractional packing problem, together with a matching scheme for covering. The scheme raises the newly revealed packing variable while increasing the dual covering variables along an exponential curve, and its analysis is a short primal-dual argument. Earlier online algorithms for throughput-competitive routing and for set cover (Alon et al. 2009) can be read as instances of it, and the same template later became the basis of a monograph on the primal-dual approach to online algorithms (Buchbinder, Naor 2009).

This mission formalizes the paper's headline result, Theorem 3.1, for the scheme exactly as the paper defines it.

Setting

Fix a finite set III of n≥1n\ge1n≥1 packing constraints (equivalently, primal covering variables), with known capacities c(i)>0c(i)>0c(i)>0. Packing variables y(1),…,y(m)y(1),\dots,y(m)y(1),…,y(m) arrive one per round; in round jjj the variable y(j)y(j)y(j) is revealed together with its non-negative column a(i,j)a(i,j)a(i,j), i∈Ii\in Ii∈I. The offline problems form the primal-dual pair of Figure 1 of the paper:

(P) min⁡∑ic(i)x(i)  s.t. ∑ia(i,j)x(i)≥1 ∀j, x≥0;(D) max⁡∑jy(j)  s.t. ∑ja(i,j)y(j)≤c(i) ∀i, y≥0.\text{(P)}\ \min\sum_i c(i)x(i)\ \text{ s.t. } \sum_i a(i,j)x(i)\ge1\ \forall j,\ x\ge0;\qquad \text{(D)}\ \max\sum_j y(j)\ \text{ s.t. } \sum_j a(i,j)y(j)\le c(i)\ \forall i,\ y\ge0.(P) mini∑​c(i)x(i)  s.t. i∑​a(i,j)x(i)≥1 ∀j, x≥0;(D) maxj∑​y(j)  s.t. j∑​a(i,j)y(j)≤c(i) ∀i, y≥0.

The profit of every y(j)y(j)y(j) is normalized to 111. Every column is assumed to have a positive entry; otherwise the packing problem is unbounded. An online algorithm may set y(j)y(j)y(j) only in round jjj and never changes it later.

The scheme with parameter B>0B>0B>0 keeps a primal vector xxx (initially 000) and the dual vector yyy. In round jjj it computes the prefix maximum ai(max⁡)=max⁡k≤ja(i,k)a_i(\max)=\max_{k\le j}a(i,k)ai​(max)=maxk≤j​a(i,k). If the new covering constraint ∑ia(i,j)x(i)≥1\sum_i a(i,j)x(i)\ge1∑i​a(i,j)x(i)≥1 already holds, it sets y(j)=0y(j)=0y(j)=0. Otherwise it sets y(j)y(j)y(j) to the least t≥0t\ge0t≥0 at which the constraint holds after every x(i)x(i)x(i) is replaced by

max⁡{x(i), 1n ai(max⁡)[exp⁡(B2c(i)∑k=1ja(i,k)y(k))−1]},y(j)=t.\max\Big\{x(i),\ \frac{1}{n\,a_i(\max)}\Big[\exp\Big(\frac{B}{2c(i)}\sum_{k=1}^{j}a(i,k)y(k)\Big)-1\Big]\Big\},\qquad y(j)=t.max{x(i), nai​(max)1​[exp(2c(i)B​k=1∑j​a(i,k)y(k))−1]},y(j)=t.

After rrr rounds, X(r)=∑ic(i)x(i)X(r)=\sum_i c(i)x(i)X(r)=∑i​c(i)x(i) is the primal value and Y(r)=∑k≤ry(k)Y(r)=\sum_{k\le r}y(k)Y(r)=∑k≤r​y(k) the dual value. For the analysis, ai(max⁡)a_i(\max)ai​(max) and ai(min⁡)a_i(\min)ai​(min) also denote the largest and the smallest non-zero coefficient of row iii over all mmm columns.

Formalization targets

Goal: Theorem 3.1

For every B>0B>0B>0, the scheme's dual solution is non-negative; after every round rrr, every non-negative y′y'y′ that satisfies the packing constraints restricted to the first rrr columns has ∑k≤ry′(k)≤B Y(r)\sum_{k\le r}y'(k)\le B\,Y(r)∑k≤r​y′(k)≤BY(r) (the scheme is BBB-competitive); and after all mmm rounds, for every iii,

∑k=1ma(i,k) y(k) ≤ c(i)⋅2log⁡(1+n ai(max⁡)/ai(min⁡))B.\sum_{k=1}^{m}a(i,k)\,y(k)\ \le\ c(i)\cdot\frac{2\log\big(1+n\,a_i(\max)/a_i(\min)\big)}{B}.k=1∑m​a(i,k)y(k) ≤ c(i)⋅B2log(1+nai​(max)/ai​(min))​.

The paper states the second part as c(i)⋅O((log⁡n+log⁡(ai(max⁡)/ai(min⁡)))/B)c(i)\cdot O\big((\log n+\log(a_i(\max)/a_i(\min)))/B\big)c(i)⋅O((logn+log(ai​(max)/ai​(min)))/B); the bound above is the one its proof establishes.

Milestones: the three claims of the proof

  1. Claim (i): X(r)≤B⋅Y(r)X(r)\le B\cdot Y(r)X(r)≤B⋅Y(r) after every round rrr.
  2. Claim (ii): after every round, x≥0x\ge0x≥0 and xxx satisfies every covering constraint revealed so far; no x(i)x(i)x(i) ever decreases.
  3. Claim (iii): the violation bound of the goal.

Significance

Theorem 3.1 says that a solution within factor BBB of the optimum can be maintained online at the price of overloading each packing constraint by a factor of order (log⁡n+log⁡(amax⁡/amin⁡))/B(\log n+\log(a_{\max}/a_{\min}))/B(logn+log(amax​/amin​))/B. Scaling the output down by the overload gives a feasible online packing solution with competitive ratio O(log⁡n+log⁡(amax⁡/amin⁡))O(\log n+\log(a_{\max}/a_{\min}))O(logn+log(amax​/amin​)), and Lemma 3.1 of the paper shows that no online algorithm does better up to constant factors. The same trade-off underlies the paper's online rounding results for routing (§5.2) and, through the covering counterpart, for set cover (§5.1).

The result is proved in the paper. As far as the platform's record shows, it is not formalized: the published OnlinePrimalDual.GeneralPacking.theorem14_1 (a restatement of the monograph's version) takes the inequality X≤BYX\le BYX≤BY and primal feasibility as hypotheses on arbitrary vectors x,yx,yx,y and derives competitiveness by weak duality; it does not mention the scheme and has no violation bound. The published OnlinePrimalDual.GeneralPacking.lemma14_2 is the matching lower bound (Lemma 3.1) and is not part of this mission. A formalization here would give the first machine-checked guarantee for the scheme itself and a reusable analysis pattern (a potential bound integrated along a monotone path) for the other schemes of the paper.

Difficulty

The paper's argument is a derivative comparison along a continuous process: while y(j)y(j)y(j) rises, ∂X/∂y(j)≤B\partial X/\partial y(j)\le B∂X/∂y(j)≤B. In the discrete formulation each round jumps directly to the least admissible y(j)y(j)y(j), and each x(i)x(i)x(i) is a maximum of its old value and an exponential, so it is continuous but not differentiable where the maximum switches; the comparison must be turned into an integral inequality over [0,y(j)][0,y(j)][0,y(j)] for such functions. The prefix maximum ai(max⁡)a_i(\max)ai​(max) changes between rounds, and the claim that this never lowers or raises the primal value needs an invariant (x(i)x(i)x(i) is always at least the current increment value). The violation bound rests on a second invariant, x(i)≤1/ai(min⁡)x(i)\le1/a_i(\min)x(i)≤1/ai​(min), which holds because y(j)y(j)y(j) is the least admissible value; any formalization that loses minimality (for example by taking an arbitrary admissible ttt) loses claim (iii).

Formalization scope

  • The instance is the published OnlinePrimalDual.GeneralPacking.GeneralInstance I (Fin m) (costs c>0c>0c>0, coefficients a≥0a\ge0a≥0); n=∣I∣n=|I|n=∣I∣ with [Nonempty I]; columns are Fin m, arrive in index order, and are 0-based in Lean, so "the first rrr columns" is (k : ℕ) < r. m≥1m\ge1m≥1 is [NeZero m].
  • ai(max⁡)a_i(\max)ai​(max), ai(min⁡)a_i(\min)ai​(min) over all columns are the published aMax and aMin. For a row with no non-zero coefficient aMin is 000 and both sides of the violation bound are 000.
  • The scheme is a function, stateAfter inst B r. The continuous loop of the paper is replaced by its discrete implementation, which the paper itself prescribes (p. 4): y(j)y(j)y(j) is the least t≥0t\ge0t≥0 restoring the new covering constraint, written as sInf. Every theorem assumes that every column has a positive entry (the paper's standing assumption, p. 4); this makes the infimum attained. If the prefix maximum is 000, Lean's 1/0=01/0=01/0=0 gives increment 000, which agrees with the paper's bracket being 000.
  • Explicit constant. The paper's c(i)⋅O((log⁡n+log⁡(ai(max⁡)/ai(min⁡)))/B)c(i)\cdot O((\log n+\log(a_i(\max)/a_i(\min)))/B)c(i)⋅O((logn+log(ai​(max)/ai​(min)))/B) is instantiated as c(i)⋅2log⁡(1+n ai(max⁡)/ai(min⁡))/Bc(i)\cdot 2\log(1+n\,a_i(\max)/a_i(\min))/Bc(i)⋅2log(1+nai​(max)/ai​(min))/B, from claim (iii) of the proof. Logarithms are natural (Real.log), since they invert Real.exp.
  • BBB-competitiveness is stated against every non-negative feasible packing solution of every prefix of the input, not only at the end.
  • A formalization in which the inequality X≤BYX\le BYX≤BY, primal feasibility, or the bound x(i)≤1/ai(min⁡)x(i)\le1/a_i(\min)x(i)≤1/ai​(min) is assumed rather than derived from the scheme is ruled out: every statement here is about the vectors the scheme computes from the instance.
  • Needed infrastructure: monotonicity and continuity of the per-round primal path, attainment of the infimum, an integral (or mean-value) form of the derivative comparison for maxima of exponentials, and weak duality for finite LPs. The weak-duality step and the per-round integration lemma are reusable for missions 2–4 of this series. Proofs of the milestones, and of helper lemmas such as the invariant x(i)≤1/ai(min⁡)x(i)\le1/a_i(\min)x(i)≤1/ai​(min), are welcome.

Selected references

  • N. Buchbinder, J. Naor, Online Primal-Dual Algorithms for Covering and Packing, Mathematics of Operations Research, 2009. https://doi.org/10.1287/moor.1080.0363
  • 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
  • N. Alon, B. Awerbuch, Y. Azar, N. Buchbinder, J. Naor, The Online Set Cover Problem, SIAM Journal on Computing 39(2), 2009. https://doi.org/10.1137/060661946
8 thms1 active userReviewed
Linear OptimizationOperations ResearchTheoretical Computer Science·Captain: mikedeng1

Online Primal-Dual Algorithms for Covering and Packing 4: Online Rounding of the Fractional Routing Scheme Respects Capacities and Is O(log P(max)·[exp(1 + 2 ln m/u(min)) − 1])-CompetitiveResearch Paper

Motivation

In online routing of virtual circuits, connection requests between pairs of nodes of a capacitated network arrive one at a time. Each request must be accepted and routed on a single path with bandwidth 111, or rejected, immediately and irrevocably, and no edge may carry more than its capacity. The goal is to maximize the number of accepted requests (the throughput). The model goes back to Awerbuch, Azar and Plotkin (FOCS 1993), whose deterministic algorithm has a logarithmic competitive ratio when edge capacities are at least logarithmic in the size of the network, and it underlies the analysis of admission control in circuit-switched and bandwidth-reserved networks.

Buchbinder and Naor (Math. Oper. Res. 2009) recover an algorithm with the same competitive factor from a general recipe: an online primal–dual scheme first produces a feasible fractional routing online, and an online version of Raghavan's pessimistic estimator (J. Comput. Syst. Sci. 1988) then rounds it, also online. This mission formalizes that construction and its guarantee (Section 5.2 of the paper, with the Section 3 scheme it uses).

Timeline:

  • 1987–1988: Raghavan and Thompson introduce randomized rounding for multicommodity flow; Raghavan derandomizes it with pessimistic estimators.
  • 1993: Awerbuch, Azar and Plotkin give the deterministic throughput-competitive online routing algorithm.
  • 2005–2009: Buchbinder and Naor's primal–dual framework (ESA 2005; MOR 2009) derives an algorithm with the same factor systematically.

Setting

Let EEE be a finite set of mmm edges with capacities u(e)>0u(e) > 0u(e)>0, and u(min⁡)=min⁡eu(e)u(\min) = \min_e u(e)u(min)=mine​u(e). Requests r1,r2,…r_1, r_2, \dotsr1​,r2​,… arrive online; request rir_iri​ comes with a finite list P(ri)\mathcal P(r_i)P(ri​) of admissible paths, each a set of edges, all of size at most P(max⁡)P(\max)P(max).

A fractional routing assigns flows f(ri,P)≥0f(r_i, P) \ge 0f(ri​,P)≥0; it is feasible when ∑P∈P(ri)f(ri,P)≤1\sum_{P \in \mathcal P(r_i)} f(r_i, P) \le 1∑P∈P(ri​)​f(ri​,P)≤1 for every request and the load ∑ri∑P∋ef(ri,P)\sum_{r_i}\sum_{P \ni e} f(r_i, P)∑ri​​∑P∋e​f(ri​,P) of every edge is at most u(e)u(e)u(e). Its value is val(f)=∑ri∑Pf(ri,P)\mathrm{val}(f) = \sum_{r_i}\sum_P f(r_i, P)val(f)=∑ri​​∑P​f(ri​,P); OPT\mathrm{OPT}OPT is the largest value of a feasible routing of the arrived requests, an upper bound on the integral optimum.

The fractional scheme. The covering LP paired with the routing LP (the paper's primal, Fig. 3) has variables x(e)x(e)x(e) (cost u(e)u(e)u(e)) and Z(ri)Z(r_i)Z(ri​) (cost 111) with constraints ∑e∈Px(e)+Z(ri)≥1\sum_{e \in P} x(e) + Z(r_i) \ge 1∑e∈P​x(e)+Z(ri​)≥1. When rir_iri​ arrives, its paths are visited in order; for each path whose constraint fails, f(ri,P)f(r_i,P)f(ri​,P) is raised from 000 to the least value restoring it, while x(e)=max⁡(x(e),1ℓ(eB′Fe/(2u(e))−1))x(e) = \max\big(x(e), \tfrac1\ell(e^{B' F_e/(2u(e))}-1)\big)x(e)=max(x(e),ℓ1​(eB′Fe​/(2u(e))−1)) for e∈Pe \in Pe∈P and Z(ri)=max⁡(Z(ri),1ℓ(eB′f(ri)/2−1))Z(r_i) = \max\big(Z(r_i), \tfrac1\ell(e^{B' f(r_i)/2}-1)\big)Z(ri​)=max(Z(ri​),ℓ1​(eB′f(ri​)/2−1)) follow the flow (FeF_eFe​ the load of eee, f(ri)f(r_i)f(ri​) the flow of rir_iri​). The parameters are ℓ=P(max⁡)+1\ell = P(\max)+1ℓ=P(max)+1 and B′=2ln⁡(1+ℓ)B' = 2\ln(1+\ell)B′=2ln(1+ℓ).

The rounding. With the rounding scale B=exp⁡(1+ln⁡(2m)/u(min⁡))−1B = \exp(1 + \ln(2m)/u(\min)) - 1B=exp(1+ln(2m)/u(min))−1, the integral edge usage χ(e)\chi(e)χ(e) and the number sss of served requests, the potential is Φ=Φ1+Φ2\Phi = \Phi_1 + \Phi_2Φ=Φ1​+Φ2​,

Φ1=12exp⁡(val(f)2B−sln⁡2),Φ2=12m∑eexp⁡((1+ln⁡2mu(e))χ(e)−Fe).\Phi_1 = \tfrac12\exp\Big(\frac{\mathrm{val}(f)}{2B} - s\ln 2\Big), \qquad \Phi_2 = \frac1{2m}\sum_{e}\exp\Big(\Big(1+\frac{\ln 2m}{u(e)}\Big)\chi(e) - F_e\Big).Φ1​=21​exp(2Bval(f)​−sln2),Φ2​=2m1​e∑​exp((1+u(e)ln2m​)χ(e)−Fe​).

After the fractional round of rir_iri​, the algorithm serves rir_iri​ on a path P∈P(ri)P \in \mathcal P(r_i)P∈P(ri​) (adding 111 to χ(e)\chi(e)χ(e) for e∈Pe \in Pe∈P) if this gives potential at most the potential Φstart\Phi^{\mathrm{start}}Φstart before the round; otherwise it rejects rir_iri​.

Formalization targets

Goal: Lemma 5.4

For every request sequence, the algorithm never exceeds a capacity, and for every feasible fractional routing fff,

χ(e)≤u(e)  ∀e,∑riχ(ri) ≥ val(f)4Bln⁡2⋅ln⁡(P(max⁡)+2)−1.\chi(e) \le u(e)\ \ \forall e, \qquad \sum_{r_i}\chi(r_i) \ \ge\ \frac{\mathrm{val}(f)}{4B\ln 2\cdot\ln(P(\max)+2)} - 1 .χ(e)≤u(e)  ∀e,ri​∑​χ(ri​) ≥ 4Bln2⋅ln(P(max)+2)val(f)​−1.

Milestone: Theorem 3.2 (packing half, on routing)

The fractional scheme's flows falgf^{\mathrm{alg}}falg are feasible and val(f)≤2ln⁡(P(max⁡)+2) val(falg)\mathrm{val}(f) \le 2\ln(P(\max)+2)\,\mathrm{val}(f^{\mathrm{alg}})val(f)≤2ln(P(max)+2)val(falg) for every feasible fff.

Milestone: Lemma 5.3

Φ≤1\Phi \le 1Φ≤1 initially, Φ>0\Phi > 0Φ>0 always, and whenever the flows of a request are raised by a non-negative amount of total at most 111, serving the request on some path or rejecting it does not increase Φ\PhiΦ.

Significance

The result shows that a deterministic online algorithm for throughput-competitive routing, previously designed by hand, falls out of two generic components: an online fractional packing scheme and an online pessimistic estimator. A side product is that the fractional phase alone produces, online, a near-optimal routing that respects all capacities exactly, independently of their size. When u(min⁡)≥log⁡nu(\min) \ge \log nu(min)≥logn the rounding loses only a constant factor and the algorithm is O(log⁡P(max⁡))O(\log P(\max))O(logP(max))-competitive, as in Awerbuch–Azar–Plotkin.

The results are proved in the paper; none is machine-checked. A formalization supplies a checked instance of the online primal–dual method together with derandomized online rounding, and fixes the constants the paper leaves inside O(⋅)O(\cdot)O(⋅). Related platform content: the monograph's OnlinePrimalDual.Routing.per_copy_guarantee and routing_competitive concern the Buchbinder–Naor (1,O(log⁡n))(1, O(\log n))(1,O(logn))-competitive algorithm with copies of the graph, a different scheme.

Difficulty

The fractional guarantee is argued in the paper continuously (rates of change of the primal and dual values), while the scheme as formalized is discrete: each flow is the least value restoring a constraint, and every primal variable is a maximum whose branch may switch during the increase. The continuous argument does not transfer verbatim, and feasibility depends on the least value restoring the constraint with equality.

The rounding is a derandomization. The existence of a good path or a good rejection is established in the paper as an expectation over a random trial; a deterministic statement about finitely many alternatives is what the mission asks for, with all m+1m+1m+1 exponential terms of Φ\PhiΦ under control at once. A frequent first attempt compares with the potential after the fractional round; the rule compares with Φstart\Phi^{\mathrm{start}}Φstart, before the flow increase, and the guarantee is stated for that comparison.

Formalization scope

  • Edges are a non-empty Fintype E; capacities are real and positive. A request is a List (Finset E) of paths; the request sequence is a list. Simple paths of a graph are a special case; nothing in the argument uses graph structure. P(max⁡)P(\max)P(max) is a parameter with every path of size at most P(max⁡)P(\max)P(max).
  • A routing is a List (List ℝ) of the shape of the request sequence. OPT\mathrm{OPT}OPT is quantified as "every feasible fractional routing".
  • The continuous increase is its discrete equivalent (an attained sInf). Ties among good paths are broken by list order; only paths that received flow in the current round are candidates for serving, so a request whose flow was not increased is rejected.
  • Explicit constants replacing O(⋅)O(\cdot)O(⋅): Theorem 3.2's O(log⁡ℓ)O(\log \ell)O(logℓ) becomes 2ln⁡(1+ℓ)=2ln⁡(P(max⁡)+2)2\ln(1+\ell) = 2\ln(P(\max)+2)2ln(1+ℓ)=2ln(P(max)+2); Lemma 5.4's O(log⁡P(max⁡)⋅[exp⁡(1+2ln⁡m/u(min⁡))−1])O(\log P(\max)\cdot[\exp(1+2\ln m/u(\min))-1])O(logP(max)⋅[exp(1+2lnm/u(min))−1]) becomes 4Bln⁡2⋅ln⁡(P(max⁡)+2)4B\ln 2\cdot\ln(P(\max)+2)4Bln2⋅ln(P(max)+2) with additive −1-1−1, where B=exp⁡(1+ln⁡(2m)/u(min⁡))−1B = \exp(1+\ln(2m)/u(\min))-1B=exp(1+ln(2m)/u(min))−1 is the scale chosen on p. 15. The printed "2ln⁡m2\ln m2lnm" differs from the proof's "ln⁡2m\ln 2mln2m"; the proof's constant is used (it is at least as strong for m≥2m \ge 2m≥2). All logarithms are natural.
  • The algorithm is a fully specified function: a formalization in which requests are never served, or OPT\mathrm{OPT}OPT is a free variable pinned by hypotheses, would trivialize the goal and is ruled out.
  • Not included: the covering half of Theorem 3.2, and the remark on u(min⁡)≥log⁡nu(\min) \ge \log nu(min)≥logn.

Contributions welcome: proofs of the milestones, a reusable lemma "convex combination ≤\le≤ value ⇒\Rightarrow⇒ some outcome ≤\le≤ value" for derandomization, and the discrete-to-continuous bridge for the Section 3 scheme.

Selected references

  • N. Buchbinder, J. Naor, Online Primal-Dual Algorithms for Covering and Packing, Mathematics of Operations Research, 2009. https://doi.org/10.1287/moor.1080.0363
  • B. Awerbuch, Y. Azar, S. Plotkin, Throughput-Competitive On-Line Routing, Proc. 34th FOCS, pp. 32–40, 1993. https://doi.org/10.1109/SFCS.1993.366884
  • P. Raghavan, Probabilistic construction of deterministic algorithms: approximating packing integer programs, J. Comput. Syst. Sci. 37(2), 1988. https://doi.org/10.1016/0022-0000(88)90003-7
  • P. Raghavan, C. D. Thompson, Randomized rounding: a technique for provably good algorithms and algorithmic proofs, Combinatorica 7(4), 1987. https://doi.org/10.1007/BF02579324
7 thms1 active userReviewed
CombinatoricsOperations ResearchProbability+1·Captain: mikedeng1

Online Stochastic Matching: Beating 1-1/e 1: When OPT = Ω(n), the Two Suggested Matchings Algorithm Achieves ALG/OPT ≥ (1 − 2/e²)/(4/3 − 2/(3e)) − ε ≈ 0.670 with Probability 1 − e^(−Ω(n))Research Paper

Motivation

Online bipartite matching models a platform that must commit each arriving request to a resource immediately. The motivating application of Feldman, Mehta, Mirrokni and Muthukrishnan is display advertising: an ad server knows from past traffic how many impressions of each type (web page, audience segment) to expect, sells them to advertisers in advance, and must assign each impression to an interested advertiser the moment a user loads the page. The goal is to fill as many contracted impressions as possible.

When arrivals are chosen by an adversary, the best ratio an online algorithm can guarantee is 1−1/e≈0.6321 - 1/e \approx 0.6321−1/e≈0.632, achieved by the RANKING algorithm of Karp, Vazirani and Vazirani (STOC 1990). The ad server, however, is not facing an adversary: it has a forecast. The i.i.d. model captures this: the graph and the distribution of impression types are known in advance, and the impressions are independent draws. The paper asks whether this knowledge allows an online algorithm to beat 1−1/e1 - 1/e1−1/e, and answers yes.

Timeline.

  • 1990: Karp, Vazirani and Vazirani give RANKING, with ratio 1−1/e1 - 1/e1−1/e for adversarial arrivals, and show this is optimal in that model.
  • 2005: Mehta, Saberi, Vazirani and Vazirani obtain 1−1/e1 - 1/e1−1/e for the budgeted AdWords generalization.
  • 2009: Feldman, Mehta, Mirrokni and Muthukrishnan (arXiv:0905.4100, FOCS 2009) show that in the i.i.d. model the two suggested matchings algorithm achieves about 0.6700.6700.670 with high probability when OPT is linear in nnn, the first ratio above 1−1/e1 - 1/e1−1/e for this model, and that no online algorithm reaches 26/2726/2726/27 in expectation.

Setting

An instance is a bipartite graph G=(A,I,E)G = (A, I, E)G=(A,I,E) with a finite set AAA of advertisers, a finite set III of impression types, and edges E⊆A×IE \subseteq A \times IE⊆A×I recording which advertisers want which types. The mission treats the case analysed throughout §4.2 of the paper, in which one impression of each type is expected (ei=1e_i = 1ei​=1). So n=∣I∣n = |I|n=∣I∣ impressions arrive one at a time, with types ω(0),…,ω(n−1)\omega(0), \dots, \omega(n-1)ω(0),…,ω(n−1) drawn independently and uniformly from III. On arrival an impression must be assigned at once and irrevocably to a still unassigned advertiser adjacent to its type, or discarded. ALG(ω)\mathrm{ALG}(\omega)ALG(ω) is the number of impressions an algorithm assigns. 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 aaa with (a,ω(t))∈E(a, \omega(t)) \in E(a,ω(t))∈E.

The two suggested matchings (TSM) algorithm works offline first. Its boosted flow graph GfG_fGf​ has a source arc of capacity 222 into every advertiser, a unit-capacity arc along every edge of EEE, and an arc of capacity 222 from every type to a sink. The algorithm takes the edge set EfE_fEf​ of an integral maximum flow. Every vertex then has at most two edges of EfE_fEf​, so EfE_fEf​ splits into vertex-disjoint paths and cycles. The algorithm colours each component blue and red:

  • on cycles, the colours alternate;
  • on odd paths, the colours alternate, with more blue than red;
  • on even paths between advertisers, the colours alternate;
  • on even paths between types, the first two edges are blue, then the colours alternate, ending in blue.

Online, the first arrival of type iii tries the advertiser along iii's blue edge, the second tries the one along its red edge, and later arrivals are discarded. A tried advertiser that is already taken is not reassigned. The advertisers fall into four classes by their coloured edges: ABRA_{BR}ABR​ (one blue, one red), ABBA_{BB}ABB​ (two blue), ABA_BAB​ (one blue only) and ARA_RAR​ (one red only).

Formalization targets

Goal: Theorem 5, first sentence, ei=1e_i = 1ei​=1

Let

α=1−2/e24/3−2/(3e)≈0.67029.\alpha = \frac{1 - 2/e^2}{4/3 - 2/(3e)} \approx 0.67029 .α=4/3−2/(3e)1−2/e2​≈0.67029.

The goal has three parts. First, every maximum flow edge set admits a colouring that follows the rules. Second, for every ε>0\varepsilon > 0ε>0 and c>0c > 0c>0 there are δ>0\delta > 0δ>0 and NNN such that, for every instance with n≥Nn \ge Nn≥N, every maximum flow edge set and every rule-following colouring,

Pr⁡ω[ OPT≥c n  ⟹  ALG≥(α−ε) OPT ]  ≥  1−e−δn.\Pr_\omega\big[\ \mathrm{OPT} \ge c\,n \implies \mathrm{ALG} \ge (\alpha - \varepsilon)\,\mathrm{OPT}\ \big] \;\ge\; 1 - e^{-\delta n}.ωPr​[ OPT≥cn⟹ALG≥(α−ε)OPT ]≥1−e−δn.

Third, α>1−1/e\alpha > 1 - 1/eα>1−1/e.

Milestones

  1. Facts 1 and 2: concentration for two balls-in-bins statistics.
  2. The note of §4.2.1: each type has no coloured edge, one blue edge, or one blue and one red edge.
  3. Equation (1): ∣Ef∣=2∣ABR∣+2∣ABB∣+∣AB∣+∣AR∣|E_f| = 2|A_{BR}| + 2|A_{BB}| + |A_B| + |A_R|∣Ef​∣=2∣ABR​∣+2∣ABB​∣+∣AB​∣+∣AR​∣.
  4. Equation (2): with high probability, ALG≥(1−1/e2)∣ABB∣+(1−2/e2)∣ABR∣+(1−3/(2e))(∣AB∣+∣AR∣)−4εn\mathrm{ALG} \ge (1 - 1/e^2)|A_{BB}| + (1 - 2/e^2)|A_{BR}| + (1 - 3/(2e))(|A_B| + |A_R|) - 4\varepsilon nALG≥(1−1/e2)∣ABB​∣+(1−2/e2)∣ABR​∣+(1−3/(2e))(∣AB​∣+∣AR​∣)−4εn.
  5. Equation (3): ∣Ef∣=2(∣AT∣+∣IS∣)+∣Eδ∣|E_f| = 2(|A_T| + |I_S|) + |E_\delta|∣Ef​∣=2(∣AT​∣+∣IS​∣)+∣Eδ​∣ for the surgered residual cut (S,T)(S,T)(S,T) of GfG_fGf​.
  6. Equation (4): with high probability, OPT≤∣ABR∣+∣ABB∣+12(∣AB∣+∣AR∣)+(12−1e)∣Eδ∣+εn\mathrm{OPT} \le |A_{BR}| + |A_{BB}| + \tfrac12(|A_B| + |A_R|) + (\tfrac12 - \tfrac1e)|E_\delta| + \varepsilon nOPT≤∣ABR​∣+∣ABB​∣+21​(∣AB​∣+∣AR​∣)+(21​−e1​)∣Eδ​∣+εn.
  7. Lemma 1: ∣Eδ∣≤23∣ABR∣+43∣ABB∣+∣AB∣+13∣AR∣|E_\delta| \le \tfrac23|A_{BR}| + \tfrac43|A_{BB}| + |A_B| + \tfrac13|A_R|∣Eδ​∣≤32​∣ABR​∣+34​∣ABB​∣+∣AB​∣+31​∣AR​∣.

Significance

The theorem separates the i.i.d. model from the adversarial one: knowing the distribution is worth a constant factor above 1−1/e1 - 1/e1−1/e. The suggested matching algorithm of the same paper (Theorem 4) shows that following a single offline matching gets exactly 1−1/e1 - 1/e1−1/e, so the second, red matching is what crosses the barrier. The paper's question started a line of work on the i.i.d. and random-order models, with later improvements to the constant by other authors under further assumptions.

The result has a written proof but, as far as the platform record shows, no machine-checked one. The mission formalizes the paper's own argument: the flow-and-colouring construction, the balls-in-bins concentration facts, the cut-based bound on OPT and the combinatorial Lemma 1. It also fixes two slips in the printed statements (see Formalization scope). The pieces are reusable beyond this paper. The occupancy concentration (Fact 1) and the satisfied-sequences bound (Fact 2) recur in analyses of online algorithms with stochastic input. The degree-capped flow encoding and its path/cycle decomposition are standard tools for 2-matchings.

Difficulty

The upper bound on OPT is the delicate part. A cut of the flow graph bounds the maximum matching of the realization graph only after a second surgery that depends on the random arrivals. Its size must then be compared with the colour classes, which are defined by a different structure (the components of EfE_fEf​). Lemma 1 bridges the two, and it depends on the exact colouring rules: a colouring that only satisfies local degree conditions can put red edges at both ends of an even advertiser path, which breaks the inequality ∣AB∣≥∣AR∣|A_B| \ge |A_R|∣AB​∣≥∣AR​∣ behind (2). On the probabilistic side, the advertisers of ABRA_{BR}ABR​ share impression types with each other, so the success events are dependent, and Fact 2 needs a bounded-differences argument in which one ball affects up to ddd sequences.

Formalization scope

All declarations live in the namespace OnlineStochMatching.TSM. Advertisers and types are finite types A I : Type, and EEE is a Finset (A × I). Probabilities are counting ratios #{ω:Fin n→I∣P ω}/∣I∣n\#\{\omega : \mathrm{Fin}\ n \to I \mid P\,\omega\}/|I|^n#{ω:Fin n→I∣Pω}/∣I∣n, so there are no measurability side conditions. OPT is a maximum over the finite, nonempty set of partial injective assignments. An integral flow of GfG_fGf​ is its set of saturated middle edges, i.e. a subset of EEE with at most two edges per vertex; EfE_fEf​ is such a set of maximum cardinality. A colouring is given by a listing of the components of EfE_fEf​ as vertex sequences. Every theorem quantifies over every maximum EfE_fEf​ and every colouring the rules allow, since the paper fixes neither.

The paper's asymptotic phrases are replaced by explicit quantifiers that come from its own proofs:

  • "with probability 1−e−Ω(n)1 - e^{-\Omega(n)}1−e−Ω(n)" and "with high probability" (Theorem 5, (2), (4)) become: ∃ δ>0, ∃ N\exists\, \delta > 0,\ \exists\, N∃δ>0, ∃N, chosen before the instance, with probability at least 1−e−δn1 - e^{-\delta n}1−e−δn for all n≥Nn \ge Nn≥N;
  • "as long as OPT =Ω(n)= \Omega(n)=Ω(n)" becomes the event OPT≥c n\mathrm{OPT} \ge c\,nOPT≥cn for an arbitrary c>0c > 0c>0 fixed before δ\deltaδ and NNN;
  • the O(1)O(1)O(1) term in the bound on ∣Aδ∗∣|A^*_\delta|∣Aδ∗​∣ (p. 8) is absorbed into εn\varepsilon nεn for n≥Nn \ge Nn≥N.

Corrections to the printed statements:

  • Theorem 5 prints ALG/OPT−ϵ≥α\mathrm{ALG}/\mathrm{OPT} - \epsilon \ge \alphaALG/OPT−ϵ≥α; the proof concludes ALG/OPT+ϵ≥α\mathrm{ALG}/\mathrm{OPT} + \epsilon \ge \alphaALG/OPT+ϵ≥α, so the goal states ALG≥(α−ε)OPT\mathrm{ALG} \ge (\alpha - \varepsilon)\mathrm{OPT}ALG≥(α−ε)OPT;
  • Fact 1 prints the failure probability 2e−ϵn/22e^{-\epsilon n/2}2e−ϵn/2; its proof gives 2e−ϵ2n/22e^{-\epsilon^2 n/2}2e−ϵ2n/2, which is used;
  • Fact 2 states two hypotheses its proof uses: the bins of a sequence are distinct, and c2<nc^2 < nc2<n.

The ratio is multiplied out, so no division by OPT occurs. A colouring condition that is unsatisfiable, or a flow set that is not maximum, would make the goal vacuous or false; part (a) of the goal rules out the first, and every statement requires maximality. Not included: the reduction to general integer eie_iei​ (§4.2.4), the tightness sentence of Theorem 5 (§4.2.5), and footnote 7's variant of the algorithm.

Useful infrastructure: bounded-differences (McDiarmid/Azuma) inequalities for functions of i.i.d. uniform variables, which exist on the platform as separate theorems; path/cycle decomposition of graphs of maximum degree two; and max-flow min-cut for unit-capacity bipartite networks. Contributions are welcome on Facts 1 and 2 independently of the combinatorics, and on Lemma 1 and equations (1) and (3), which are deterministic.

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
  • A. Mehta, A. Saberi, U. Vazirani, V. Vazirani, AdWords and generalized online matching, FOCS 2005; J. ACM 54(5), 2007. https://doi.org/10.1145/1284320.1284321
13 thms1 active userReviewed
CombinatoricsLinear OptimizationOperations Research+1·Captain: mikedeng1

Online Primal-Dual Algorithms for Covering and Packing 3: A Deterministic O(log d log(n/OPT))-Competitive Algorithm for Online Unweighted Set CoverResearch Paper

Motivation

Online set cover is the basic covering problem in which the requests arrive over time. A ground set of elements and a family of sets are known in advance, but which elements must be covered is revealed one element at a time, and each arriving element has to be covered at once by a set chosen irrevocably. The problem models resource placement under unknown demand (facilities, servers, sensors that must serve clients as they appear) and is the prototype for a family of online covering problems.

Alon, Awerbuch, Azar, Buchbinder and Naor (SIAM J. Comput. 2009) gave the first deterministic algorithm, with competitive ratio O(log⁡mlog⁡n)O(\log m\log n)O(logmlogn) for nnn elements and mmm sets, and showed that no deterministic algorithm does better than Ω(log⁡mlog⁡n/(log⁡log⁡m+log⁡log⁡n))\Omega(\log m\log n/(\log\log m+\log\log n))Ω(logmlogn/(loglogm+loglogn)) on some instances. Buchbinder and Naor (Math. Oper. Res. 2009) recast the fractional part of that algorithm as an instance of a general online primal-dual scheme for covering and packing linear programs, and turned an offline pessimistic estimator of Srinivasan into an online potential function. The result, in their Section 5.1, is a deterministic algorithm whose ratio O(log⁡dlog⁡(n/OPT))O(\log d\log(n/OPT))O(logdlog(n/OPT)) depends on the maximum element frequency ddd instead of the number of sets mmm, and on the ratio n/OPTn/OPTn/OPT instead of nnn.

Timeline:

  • 2003 (conference), 2009 (journal): Alon et al., deterministic O(log⁡mlog⁡n)O(\log m\log n)O(logmlogn) for unweighted online set cover, with a potential ∑j∉Cn2wj\sum_{j\notin C} n^{2w_j}∑j∈/C​n2wj​.
  • 2005 (conference), 2009 (journal): Buchbinder and Naor, the general online fractional covering/packing scheme, and the derandomized rounding of this mission.

Setting

A set-cover instance consists of a finite ground set XXX of nnn elements and a finite family S\mathcal SS of mmm sets. For an element eee, Se\mathcal S_eSe​ is the collection of sets containing eee, and ddd bounds its size: ∣Se∣≤d|\mathcal S_e|\le d∣Se​∣≤d for every eee (the frequency). In the unweighted problem every set costs 111.

Elements arrive in a list σ\sigmaσ. The algorithm maintains:

  • fractional weights w(s)≥0w(s)\ge 0w(s)≥0 for the sets, produced by the paper's Section 3 scheme with {0,1}\{0,1\}{0,1} coefficients: when an element eee arrives that is not yet fractionally covered (∑s∈Sew(s)<1\sum_{s\in\mathcal S_e}w(s)<1∑s∈Se​​w(s)<1), its dual variable y(e)y(e)y(e) is raised to the least value at which
w(s)=max⁡{w(s), 1d(exp⁡(B2c(s)∑k: s∋eky(ek))−1)}(s∋e)w(s)=\max\Big\{w(s),\ \tfrac1d\Big(\exp\Big(\tfrac{B}{2c(s)}\textstyle\sum_{k:\ s\ni e_k}y(e_k)\Big)-1\Big)\Big\}\qquad(s\ni e)w(s)=max{w(s), d1​(exp(2c(s)B​∑k: s∋ek​​y(ek​))−1)}(s∋e)

gives ∑s∈Sew(s)≥1\sum_{s\in\mathcal S_e}w(s)\ge 1∑s∈Se​​w(s)≥1; here B>0B>0B>0 is a parameter;

  • a cover C⊆S\mathcal C\subseteq\mathcal SC⊆S that only grows; CCC is the set of elements covered by C\mathcal CC.

With f(e)=min⁡{1,exp⁡(−α+α∑s∋ew(s))}f(e)=\min\{1,\exp(-\alpha+\alpha\sum_{s\ni e}w(s))\}f(e)=min{1,exp(−α+α∑s∋e​w(s))}, the potential is Φ=Φ1+Φ2\Phi=\Phi_1+\Phi_2Φ=Φ1​+Φ2​ with

Φ1=1−∏e∈X∖C(1−f(e)),Φ2=exp⁡(∑s∈S((ln⁡2)χC(s)−αw(s))−OPT),\Phi_1=1-\prod_{e\in X\setminus C}\big(1-f(e)\big),\qquad \Phi_2=\exp\Big(\sum_{s\in\mathcal S}\big((\ln 2)\chi_{\mathcal C}(s)-\alpha w(s)\big)-OPT\Big),Φ1​=1−e∈X∖C∏​(1−f(e)),Φ2​=exp(s∈S∑​((ln2)χC​(s)−αw(s))−OPT),

where OPTOPTOPT is the optimum number of sets covering the arrived elements, assumed known, r=eln⁡(e/(e−1))r=e\ln(e/(e-1))r=eln(e/(e−1)) and α=max⁡{1,ln⁡(rn/OPT)}\alpha=\max\{1,\ln(rn/OPT)\}α=max{1,ln(rn/OPT)}. The rounding rule: each time the weight of a set sss is augmented, sss is added to C\mathcal CC if this does not increase Φ\PhiΦ.

Formalization targets

Goal: Lemma 5.2

Every arriving element is covered by C\mathcal CC, and at every time

∣C∣ ≤ (2αln⁡(1+d)+1) OPTln⁡2,α=max⁡{1,ln⁡rnOPT}.|\mathcal C|\ \le\ \frac{\big(2\alpha\ln(1+d)+1\big)\,OPT}{\ln 2},\qquad \alpha=\max\Big\{1,\ln\frac{rn}{OPT}\Big\}.∣C∣ ≤ ln2(2αln(1+d)+1)OPT​,α=max{1,lnOPTrn​}.

This is the paper's OPT⋅O(log⁡dlog⁡(n/OPT))OPT\cdot O(\log d\log(n/OPT))OPT⋅O(logdlog(n/OPT)) with the constant its proof gives.

Milestones

  1. Theorem 3.2 (covering half): for every B>0B>0B>0, the {0,1}\{0,1\}{0,1} scheme with frequency bound ℓ\ellℓ yields a fractional cover of cost at most 2ln⁡(1+ℓ)2\ln(1+\ell)2ln(1+ℓ) times that of any fractional cover.
  2. Lemma 5.1 (i): initially Φ≤1\Phi\le 1Φ≤1; Φ>0\Phi>0Φ>0 in every state.
  3. Lemma 5.1 (ii): after the weight of a set is augmented by δ≥0\delta\ge0δ≥0, taking the set or excluding it leaves Φ\PhiΦ no larger than before.

Significance

The bound improves Alon et al.'s O(log⁡mlog⁡n)O(\log m\log n)O(logmlogn) whenever sets are many but each element lies in few of them (d≪md\ll md≪m), and whenever the optimum is large compared with nnn. It also shows that the fractional part and the rounding part of an online covering algorithm can be designed separately: any online fractional solution with a competitive guarantee can be rounded deterministically by an online potential function. The same method gives the routing result of Section 5.2 of the paper.

As far as is known, none of these results is machine-checked. Related formal work on the platform covers Alon et al.'s algorithm and its O(log⁡mlog⁡n)O(\log m\log n)O(logmlogn) bound (a different potential and a different fractional update), and the general covering scheme of Buchbinder and Naor's monograph; neither states Lemma 5.1 or Lemma 5.2, nor Theorem 3.2 for the {0,1}\{0,1\}{0,1} scheme with ℓ\ellℓ in place of nnn. A complete development here would give a verified derandomized rounding argument, reusable for other online covering problems.

Difficulty

The difficulty is not in the final inequality, which follows from Φ2≤1\Phi_2\le1Φ2​≤1 in one line, but in keeping Φ≤1\Phi\le1Φ≤1 throughout. The decision to take a set must be made with no knowledge of future elements, and the potential has to account for elements that may never arrive: Φ1\Phi_1Φ1​ ranges over the whole ground set. Lemma 5.1 (ii) asks that, for every current state, one of the two decisions does not increase a non-linear function of all uncovered elements at once and this must hold for every state, not only for the states a particular run reaches. On the fractional side, Theorem 3.2 is only asserted, "along the same lines" as Theorem 3.1, so its constant has to be re-derived with ℓ\ellℓ in place of nnn and with the scheme's continuous increase made discrete.

A tempting shortcut is to bound ∣C∣|\mathcal C|∣C∣ by the number of rounds or by ∑sw(s)\sum_s w(s)∑s​w(s) directly; neither gives a logarithmic factor in n/OPTn/OPTn/OPT, which comes only from the choice of α\alphaα in Φ\PhiΦ.

Formalization scope

The instance is the published OnlinePrimalDual.OnlineSetCover.SetCoverInstance (elements E, set indices T, incidence elemSets, positive costs c), with elementWeight and coveredBy. Unit costs are the hypothesis ∀ s, inst.c s = 1 in Lemma 5.2; Theorem 3.2 is stated for general positive costs. Logarithms are natural (Real.log), because they invert Real.exp.

Committed conventions:

  • The fractional scheme is the discrete form of the continuous increase: in each round y(e)y(e)y(e) is the least t≥0t\ge0t≥0 (an sInf) at which the new constraint holds. Arrival lists may repeat elements; an element already covered changes nothing.
  • The rounding treats each set's increase within a round as one augmentation; the augmented sets are processed one at a time in the order of a list ord containing every set, and a set is added when Φ(w+δs1s,C∪{s})≤Φ(w,C)\Phi(w+\delta_s\mathbf 1_s,\mathcal C\cup\{s\})\le\Phi(w,\mathcal C)Φ(w+δs​1s​,C∪{s})≤Φ(w,C). Lemma 5.2 holds for every order.
  • The algorithm is a function, so every run exists.
  • OPTOPTOPT is a natural number ≥1\ge1≥1 bounding the size of some cover of the arrived elements; d≥1d\ge1d≥1 bounds the frequency of every element. At the true optimum and the maximum frequency this is the paper's statement.
  • Explicit constants replacing O(⋅)O(\cdot)O(⋅): Theorem 3.2's O(log⁡ℓ)O(\log\ell)O(logℓ) is 2ln⁡(1+ℓ)2\ln(1+\ell)2ln(1+ℓ); Lemma 5.2's OPT⋅O(log⁡dlog⁡(n/OPT))OPT\cdot O(\log d\log(n/OPT))OPT⋅O(logdlog(n/OPT)) is (2αln⁡(1+d)+1) OPT/ln⁡2(2\alpha\ln(1+d)+1)\,OPT/\ln2(2αln(1+d)+1)OPT/ln2 with α=max⁡{1,ln⁡(rn/OPT)}\alpha=\max\{1,\ln(rn/OPT)\}α=max{1,ln(rn/OPT)}.
  • The proof of Lemma 5.1 (i) on p. 13 prints ddd where rrr is meant ("exp⁡(−dne−α)\exp(-dne^{-\alpha})exp(−dne−α)", "α≥ln⁡(dn/OPT)\alpha\ge\ln(dn/OPT)α≥ln(dn/OPT)"); the statement uses rrr, and rrr is printed "eln⁡(e/e−1)e\ln(e/e-1)eln(e/e−1)" for eln⁡(e/(e−1))e\ln(e/(e-1))eln(e/(e−1)).

Not in scope: the packing half of Theorem 3.2, Theorem 3.1 for general coefficients, and the doubling wrapper of p. 12 that removes the assumption that OPTOPTOPT is known. The guarantee is for the algorithm run with the stated α\alphaα; a formalization in which Φ\PhiΦ's product ranges only over arrived elements, or in which the weights or the chosen family are free variables constrained by hypotheses instead of being produced by the algorithm (OPTOPTOPT is an input of the algorithm, as the page assumes it known), would be a different and weaker statement and is ruled out.

Contributions welcome: proofs of the three milestones and of the goal, and general lemmas about the fractional round (attainment of the least ttt, monotonicity of the weights) that other online covering missions can reuse.

Selected references

  • N. Buchbinder, J. Naor, Online Primal-Dual Algorithms for Covering and Packing, Mathematics of Operations Research 34(2), 2009. https://doi.org/10.1287/moor.1080.0363
  • N. Alon, B. Awerbuch, Y. Azar, N. Buchbinder, J. Naor, The Online Set Cover Problem, SIAM Journal on Computing 39(2), 2009. https://doi.org/10.1137/060661946
  • A. Srinivasan, Improved approximation guarantees for packing and covering integer programs, SIAM Journal on Computing 29(2), 1999. https://doi.org/10.1137/S0097539796314240
  • 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 thms1 active userReviewed
CombinatoricsOperations ResearchProbability+1·Captain: mikedeng1

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

Motivation

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

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

Setting

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

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

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

Formalization targets

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

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

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

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

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

Significance

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

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

Difficulty

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

Formalization scope

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

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

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

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

Selected references

  • Jon Feldman, Aranyak Mehta, Vahab Mirrokni, and S. Muthukrishnan, Online Stochastic Matching: Beating 1-1/e, arXiv:0905.4100v1, 2009; FOCS 2009. Preprint
10 thms1 active userReviewed
CombinatoricsOperations ResearchProbability+1·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 thms1 active userReviewed
CombinatoricsOperations ResearchProbability+1·Captain: mikedeng1

An Optimal Algorithm for On-line Bipartite Matching 1: RANKING Finds a Matching of Expected Size at Least n(1 − 1/e) − o(n) on Every Graph with a Perfect MatchingResearch Paper

Motivation

Online matching describes allocation when requests must be answered as they arrive. A matching decision uses only the edges revealed so far and cannot be revised when later requests appear. In the bipartite setting studied by Karp, U. Vazirani, and V. Vazirani, the arriving vertices are girls and the possible partners are boys. Even when the full graph has a perfect matching, a fixed greedy priority can leave many girls unmatched. The paper introduced RANKING, which randomly chooses the boys' priority order once and then uses that order for every arrival.

The question is quantitative: how many pairs does RANKING guarantee in expectation against a graph and an arrival order chosen before its random ranking? The paper's target is a fraction approaching 1−1/e1-1/e1−1/e of the nnn pairs in a perfect matching. Its printed Theorem 1 concerns an auxiliary algorithm called EARLY; the RANKING statement follows the intended chain through Lemmas 3 and 5. The original EARLY analysis has a gap for general upper-triangular matrices, so this mission states the RANKING target and the earlier, unaffected lemmas separately. This distinction matters because an assertion about EARLY would be a different formalization target.

Setting

A bipartite graph has a boy side UUU and a girl side VVV, each with nnn vertices. An edge (u,v)(u,v)(u,v) means that boy uuu may be paired with girl vvv. A matching is a collection of edges in which no boy or girl occurs twice. The standing hypothesis for the performance guarantee is that the graph has a perfect matching: some bijection from boys to girls selects an edge for every boy. The graph is otherwise arbitrary.

Girls arrive in a predetermined order. When girl vvv arrives, only her incident edges are revealed. RANKING first chooses a uniformly random permutation π\piπ of the boys. For each arriving girl, it selects the highest-ranked adjacent boy who is still unmatched, if one exists. Write MR(G,π)M_{\mathrm R}(G,\pi)MR​(G,π) for the final matching. The expected size is the finite average over all n!n!n! rankings. The guarantee must hold for every graph and every predetermined arrival order; relabeling girls lets the formal statement fix their order to n,n−1,…,1n,n-1,\ldots,1n,n−1,…,1.

The paper also uses a dual rows-arrive view. Boys arrive in an order and choose among eligible girls, whose priority order is fixed. This is the same greedy rule with the sides exchanged. For its triangular reduction, the columns are numbered 1,…,n1,\ldots,n1,…,n and column nnn has highest priority. An upper-triangular matrix with unit diagonal is a graph with every edge (i,i)(i,i)(i,i) and with an edge (i,j)(i,j)(i,j) only when i≤ji\le ji≤j. On such a graph, the auxiliary algorithm EARLY declines to match row iii if column iii is already covered by EARLY's own matching. For any matching MMM, D(M)D(M)D(M) denotes the indices for which both row iii and column iii are covered.

Formalization targets

RANKING guarantee

For every ε>0\varepsilon>0ε>0, one threshold NNN must work for every size n≥Nn\ge Nn≥N and every graph GGG with a perfect matching:

Eπ∼Unif(Sn)∣MR(G,π)∣≥(1−e−1−ε)n.\mathbb E_{\pi\sim\mathrm{Unif}(S_n)}|M_{\mathrm R}(G,\pi)| \ge (1-e^{-1}-\varepsilon)n.Eπ∼Unif(Sn​)​∣MR​(G,π)∣≥(1−e−1−ε)n.

This is the uniform lower-bound reading of the paper's n(1−1/e)−o(n)n(1-1/e)-o(n)n(1−1/e)−o(n) target. It does not prescribe a finite-nnn additive constant. The mission goal is the RANKING assertion drawn from the paper's Theorem 1 and Lemmas 3 and 5. The six milestones formalize the paper's Lemmas 1–5 and the corollary to Lemma 4, in their source order. They cover the side-exchange identity, arbitrary refusal algorithms, the triangular reduction, the matching-count identity, its expected form, and the pointwise comparison with EARLY.

Significance

The guarantee gives a concrete worst-case floor for a simple randomized allocation rule: as the graph size grows, RANKING matches at least an asymptotic 1−1/e1-1/e1−1/e fraction of the pairs available in a perfect matching, in expectation. The order is selected before the random permutation, matching the paper's performance measure. The lower bound remains meaningful on sparse graphs and does not rely on a density assumption. The original paper also studies the limit on what any randomized online algorithm can guarantee; that upper-bound result is treated in a separate mission.

A formal development here would provide reusable finite definitions for online greedy matching, arbitrary state-dependent refusal, fixed-order duality, and uniform expectation over permutations. The goal is a theorem statement awaiting a machine-checked proof; compiling the draft declarations verifies their Lean syntax and types, not their truth. The local mission separates the valid early structural statements from later statements whose published EARLY argument does not justify them on general upper-triangular graphs.

Difficulty

The random permutation does not make the fate of different vertices independent. Matching one girl removes a boy who might be essential to a later girl, so a per-arrival probability estimate cannot simply be added across all arrivals. The dual and triangular views capture useful structure, but turning that structure into a uniform bound for every graph is the main obstacle. In particular, reasoning about EARLY as if its matched-column set had the same monotonicity as unrestricted RANKING fails on some upper-triangular matrices. A proof of the goal must establish the RANKING guarantee without treating those later EARLY statements as available facts.

Formalization scope

Boys and girls are both Fin n; an adjacency matrix is a relation between them. A published predicate represents a perfect matching as an edge-preserving bijection. The generic greedy run processes arrival times in increasing order and interprets a smaller priority index as higher rank. It makes a new matching decision only from the current matching and the arriving vertex. RANKING's returned edges are consistently ordered as (boy, girl), even though girls arrive. The dual run has rows arriving in an arbitrary permutation and takes column n−1n-1n−1 as highest priority in Lean's zero-based numbering. EARLY's refusal test consults the matching constructed by EARLY itself.

All expectations are finite averages over Equiv.Perm (Fin n), scaled by 1/n!1/n!1/n!. The formal o(n)o(n)o(n) claim is ∀ε>0, ∃N, ∀n≥N, ∀G\forall\varepsilon>0,\ \exists N,\ \forall n\ge N,\ \forall G∀ε>0, ∃N, ∀n≥N, ∀G with a perfect matching, the displayed lower bound. Placing NNN before GGG preserves the worst-case meaning. The perfect-matching condition is essential: the empty graph cannot satisfy a positive linear guarantee. The triangular matrix in Lemma 3 retains every diagonal edge, so a zero matrix cannot witness the reduction. These conventions rule out vacuous versions of the goal and reduction.

The definitions of partial matching, covered vertices, greedy run, RANKING, EARLY, and uniform average are part of the mission. A complete proof may build further finite counting and permutation machinery; such lemmas can be shared beyond this paper. Contributions toward the six milestone statements, the goal, and faithful supporting results are in scope. Lemmas 6–12 and the printed EARLY Theorem 1 are outside this mission because their analysis depends on claims that fail for some matrices allowed by their surrounding assumptions.

Selected references

  • Richard M. Karp, Umesh V. Vazirani, and Vijay V. Vazirani, An Optimal Algorithm for On-line Bipartite Matching, Proceedings of the 22nd Annual ACM Symposium on Theory of Computing, 1990, pp. 352–358. DOI: 10.1145/100216.100262.
10 thms2 active usersReviewed
CombinatoricsOperations ResearchProbability+1·Captain: mikedeng1

An Optimal Algorithm for On-line Bipartite Matching 2: No Randomized On-line Algorithm Guarantees an Expected Matching Larger Than n(1 − 1/e) + o(n)Research Paper

Motivation

An online matching algorithm must commit to each assignment when a request arrives, before it sees later requests. This limitation arises whenever waiting for the full instance is impossible: an available resource can be assigned to a current request or saved for an unknown future request. The quality of the assignment is measured by how many requests can be matched. The central question is how much the lack of future information costs, even when an algorithm uses randomness. Karp, U. Vazirani, and V. Vazirani studied this question for bipartite graphs and proved an asymptotic ceiling of 1−1/e1-1/e1−1/e for the expected fraction matched by any randomized online rule on graphs with a perfect matching (Karp, Vazirani, and Vazirani, 1990).

The paper also analyzes RANKING and gives a matching asymptotic guarantee from below. This mission concerns the separate upper-bound result: the existence of hard instances for every algorithm. The upper bound matters independently of any particular proposed rule. It says that improving an algorithm's decisions cannot remove the worst-case loss due to decisions made before all columns are revealed. Its quantifiers make the adversarial model precise: the graph is selected with knowledge of the algorithm, before that algorithm's random choices are made.

Setting

There are nnn boys, represented by rows, and nnn girls, represented by columns. A bipartite graph G⊆[n]×[n]G\subseteq[n]\times[n]G⊆[n]×[n] records which boy and girl pairs may be matched. The graph is assumed to contain a perfect matching: some bijection between the two sides uses only edges of GGG. This assumption supplies an offline benchmark of nnn matches. Without it, a graph with no edges would make every upper bound on matching size empty of content.

Girls arrive one at a time. On arrival, the algorithm learns the edges incident to that girl, chooses an as-yet-unmatched adjacent boy, or declines to match her. A choice cannot later be changed. The paper labels columns so that they arrive in the order n,n−1,…,1n,n-1,\ldots,1n,n−1,…,1. A deterministic algorithm can base its action on the neighborhoods revealed so far. A randomized algorithm can additionally use internal random choices. Its performance p(A)p(A)p(A) is the minimum, over a graph with a perfect matching and a preselected arrival order, of its expected number of matches; the expectation is over its own randomness (Karp, Vazirani, and Vazirani, 1990, p. 352).

The paper's hard family begins with the complete upper-triangular graph TnT_nTn​, where row iii is adjacent to column jjj exactly when i≤ji\le ji≤j. Its columns arrive from largest to smallest. Relabeling the rows by a permutation π\piπ gives an instance TπT_\piTπ​ that still contains a perfect matching. The algorithm RANDOM takes an eligible boy uniformly whenever at least one is available. Write VTn(n,∅)V_{T_n}(n,\varnothing)VTn​​(n,∅) for RANDOM's expected number of matches on TnT_nTn​, beginning with no matched rows (Karp, Vazirani, and Vazirani, 1990, p. 357).

Formalization targets

Main theorem

For every ε>0\varepsilon>0ε>0, there is one threshold NNN such that, for every n≥Nn\ge Nn≥N and every randomized online algorithm AAA on nnn rows and columns, a graph GGG with a perfect matching satisfies

EA[∣M(G)∣]≤(1−e−1+ε)n.\mathbb E_A[|M(G)|]\le (1-e^{-1}+\varepsilon)n.EA​[∣M(G)∣]≤(1−e−1+ε)n.

The threshold is uniform over algorithms. The graph may depend on AAA. This is the quantified upper-bound reading of Theorem 2's p(A)≤n(1−1/e)+o(n)p(A)\le n(1-1/e)+o(n)p(A)≤n(1−1/e)+o(n) (Karp, Vazirani, and Vazirani, 1990, Theorem 2).

Milestone statements

Lemma 13 identifies the average size produced by any deterministic greedy algorithm on a uniformly permuted TnT_nTn​ with RANDOM's expected size on TnT_nTn​. Lemma 14 bounds the worst-case performance of every randomized algorithm, including non-greedy algorithms, by that value. Lemma 16 identifies the value itself by the two-sided limit

lim⁡n→∞VTn(n,∅)n=1−e−1.\lim_{n\to\infty}\frac{V_{T_n}(n,\varnothing)}{n}=1-e^{-1}.n→∞lim​nVTn​​(n,∅)​=1−e−1.

Together these statements give the finite-instance comparison and the asymptotic value named in the goal (Karp, Vazirani, and Vazirani, 1990, Lemmas 13, 14, 16).

Significance

Theorem 2 limits what any randomized online matching rule can guarantee under the paper's oblivious-adversary performance measure. Since the hard graph always admits a perfect matching, the gap from nnn is caused by the order of information and the required irrevocable decisions. The statement applies to the entire algorithm class, rather than comparing two selected procedures. It therefore provides the ceiling against which the paper's lower guarantee for RANKING is measured.

Formalizing the result requires a reusable description of finite online algorithms, their visible histories, and expected matching size. It also requires a definition of uniform random choice from currently eligible rows with an explicit empty-choice case. The 1990 result is proved in the paper; the statements in this mission are targets for machine-checked proofs, not claims of an existing formal proof. The model can support later statements about other matching rules and hard-input distributions without changing the meaning of an online decision.

Difficulty

For a fixed graph, an algorithm can be designed around the graph's particular perfect matching, so one hard graph cannot simply be announced in advance for all deterministic algorithms. The theorem instead has to bound each randomized algorithm against a graph selected for that algorithm. Even then, checking the triangular graph against one rule does not establish a universal bound: different rules can respond differently to the same revealed neighborhoods. The crucial mathematical obstacle is a comparison across all such rules while preserving the restriction that future columns remain unseen. A further asymptotic step is needed to determine RANDOM's value on the triangular family, including both sides of the o(n)o(n)o(n) claim in Lemma 16.

Formalization scope

Both sides are Fin n. An edge set is a subset of row-column pairs. The published perfect-matching predicate is reused: it asks for a bijection whose every selected pair is an edge. The paper's one-based column order n,…,1n,\ldots,1n,…,1 is represented by zero-based order n−1,…,0n-1,\ldots,0n−1,…,0; arrival number zero is the largest column. A deterministic rule receives only the neighborhoods of columns that have arrived, including the current column. A randomized rule is a probability mass function on the finite set of deterministic rules, allowing arbitrary correlations among its choices. The expected size is a finite sum over that mass function.

RANDOM is represented by the conditional-expectation recursion for uniform choice among eligible rows. If none is eligible, the column remains unmatched and there is no division by zero. For n=0n=0n=0, the matching and RANDOM value are zero; the main theorem is eventual in nnn and Lemma 16's value at zero does not affect the limit. Lemma 16 uses a two-sided limit. The paper's o(n)o(n)o(n) in Theorem 2 is stated as ∀ε>0,∃N,∀n≥N\forall\varepsilon>0,\exists N,\forall n\ge N∀ε>0,∃N,∀n≥N with NNN before the algorithm, so the error is uniform across algorithms. A hard graph must contain a perfect matching; dropping this condition would let the empty graph satisfy the inequality trivially.

Useful contributions include finite probabilistic averaging, the greedy comparison, analysis of the recursive RANDOM value, and a proof joining Lemmas 14 and 16 into the uniform goal. The history and matching-run definitions are reusable for other finite online bipartite matching statements. The two numbered claims inside the paper's proof of Lemma 13 and its later remarks are outside the mission's curated milestone list.

Selected references

  • Richard M. Karp, Umesh V. Vazirani, and Vijay V. Vazirani, An Optimal Algorithm for On-line Bipartite Matching, Proceedings of the 22nd Annual ACM Symposium on Theory of Computing, pp. 352–358, 1990. DOI: 10.1145/100216.100262.
9 thms2 active usersReviewed
Algorithmic Game TheoryComplexity TheoryLinear Optimization+1·Captain: mikedeng1

The Polynomial Hierarchy and a Simple Model for Competitive Analysis: Every Optimum of the (p+1)-Level Linear Game J'(F) Is Binary, with x(F) = 1 iff the Σ_p Sentence (3.3) HoldsResearch Paper

Why multi-level programs are hard

Multi-level programs model a hierarchy of decision makers: a leader commits to a decision, a follower optimises given it, a follower of the follower optimises given both, and so on. Bilevel programs are the standard model of Stackelberg competition, toll setting, network interdiction and many other leader–follower problems in operations research (Candler and Townsley 1982; Bard and Falk 1982). When every level has a linear criterion and the constraints are linear, each player's problem looks like a linear program, and it is natural to hope that the whole hierarchy is solvable in polynomial time.

R. G. Jeroslow's 1985 paper (Math. Programming 32, 146–164) shows that this hope fails at every level of the polynomial hierarchy: a (p+1)(p+1)(p+1)-level linear program with fixed criteria can encode the truth of a Σp\Sigma_pΣp​ quantified Boolean sentence. The result places multi-level linear programming in the polynomial hierarchy and is widely cited for the Σp\Sigma_pΣp​-hardness of such programs; NP-hardness of bilevel linear programs (Corollary 4.6) is its special case p=1p=1p=1.

Setting

A multi-level program has real variables x=(x1,…,xp)x=(x^1,\dots,x^p)x=(x1,…,xp), a feasible set S0S_0S0​ (a polyhedron {x:∑iAixi≥b}\{x: \sum_i A^ix^i\ge b\}{x:∑i​Aixi≥b} in the linear case), and players p,p−1,…,1p,p-1,\dots,1p,p−1,…,1 who move in that order; player iii controls xix^ixi and minimises a fixed linear criterion cixc^ixcix. The solution sets are defined from the last mover upwards: S1S_1S1​ is the set of x∈S0x\in S_0x∈S0​ at which player 1's criterion is minimal given the choices of all earlier movers, and in general SjS_{j}Sj​ keeps the points of Sj−1S_{j-1}Sj−1​ minimising cjxc^jxcjx among the points of Sj−1S_{j-1}Sj−1​ that agree with xxx on xj+1,…,xpx^{j+1},\dots,x^pxj+1,…,xp. The value is cpxc^pxcpx on SpS_pSp​, when Sp≠∅S_p\neq\emptysetSp​=∅. The sets SjS_jSj​ can be empty even when S0S_0S0​ is a nonempty polytope: the paper's four-level Example has S4=∅S_4=\emptysetS4​=∅ because S3S_3S3​ is not closed.

A propositional formula FFF over blocks of atoms X1,…,XpX_1,\dots,X_pX1​,…,Xp​ (block XkX_kXk​ has nkn_knk​ atoms) is encoded by the linear system LFL_FLF​: one variable x(G)∈[0,1]x(G)\in[0,1]x(G)∈[0,1] per non-atomic subformula, with the inequalities (3.1a)–(3.1c) for ∨\vee∨, ∧\wedge∧, ¬\neg¬. The quantifier of block XkX_kXk​ is Qk=∃Q_k=\existsQk​=∃ when p−kp-kp−k is even, so

(∃Xp)(∀Xp−1)⋯(Q1X1) [F(X1,…,Xp)=1](3.3)(\exists X_p)(\forall X_{p-1})\cdots(Q_1X_1)\,[F(X_1,\dots,X_p)=1] \qquad (3.3)(∃Xp​)(∀Xp−1​)⋯(Q1​X1​)[F(X1​,…,Xp​)=1](3.3)

is a Σp\Sigma_pΣp​ sentence.

The game J′(F)J'(F)J′(F) adds a bookkeeper, player 000, who moves last. Player k≥1k\ge1k≥1 controls the atoms of XkX_kXk​, and player 111 also controls auxiliary variables yyy; the bookkeeper controls the x(G)x(G)x(G), a variable uuu fixed to 111, and auxiliary variables zzz. Two gadgets, (4.1) and (4.6), let the bookkeeper and player 1 turn the linear criteria into the piecewise-linear functions Zk=1−x(F)+2∑jP(xkj)Z_k=1-x(F)+2\sum_jP(x_{kj})Zk​=1−x(F)+2∑j​P(xkj​) (or x(F)+…x(F)+\dotsx(F)+… for universal QkQ_kQk​) and Z1=(1−x(F))+fr(1−x(F))+10L∑jfr(x1j)+…Z_1=(1-x(F))+fr(1-x(F))+10L\sum_j fr(x_{1j})+\dotsZ1​=(1−x(F))+fr(1−x(F))+10L∑j​fr(x1j​)+…, where LLL is the length of FFF, fr(x)=min⁡{x,1−x}fr(x)=\min\{x,1-x\}fr(x)=min{x,1−x}, and P(x)=1P(x)=1P(x)=1 at x∈{0,1}x\in\{0,1\}x∈{0,1}, 222 otherwise.

Formalization targets

Goal: Theorem 4.5 (p≥2p\ge2p≥2)

Let SSS be the set of optimal solutions Sp+1S_{p+1}Sp+1​ of J′(F)J'(F)J′(F). Then S≠∅S\neq\emptysetS=∅; at every optimum all atom variables and all x(G)x(G)x(G) are binary; and

x(F)=1  ⟺  (3.3) holds,value(J′(F))=2np+1−x(F).x(F)=1 \iff (3.3)\ \text{holds},\qquad \text{value}(J'(F)) = 2n_p+1-x(F).x(F)=1⟺(3.3) holds,value(J′(F))=2np​+1−x(F).

Moreover, when (3.3) holds, v∈Rnpv\in\mathbb R^{n_p}v∈Rnp​ is player ppp's block in some optimum iff vvv is binary and the Πp−1\Pi_{p-1}Πp−1​ sentence (4.17) holds at the truth valuation of vvv.

Milestones

In attack order: Lemma 3.1 (correctness of LFL_FLF​ on binary inputs); Lemma 4.1 (robustness of LFL_FLF​ near binary inputs); Lemmas 4.2 and 4.3 (the bottom two levels of a bounded linear multi-level program are solvable, via LP duality); the bookkeeper identities z=∣2y−x∣z=|2y-x|z=∣2y−x∣ and z=fr(x)z=fr(x)z=fr(x) in S1S_1S1​; (4.2) and (4.3) (player 1's and player kkk's responses on the gadgets); Lemma 4.4 (the induction on kkk with the higher blocks fixed). Companion theorems: Corollary 4.6 (the bilevel case: value 000 iff (∃X1)F(\exists X_1)F(∃X1​)F), the §2 Example, and Proposition 3.2 (the pure binary game J(F)J(F)J(F)).

Significance

The theorem shows that deciding the value of a (p+1)(p+1)(p+1)-level linear program with fixed criteria is at least as hard as deciding Σp\Sigma_pΣp​ sentences, so known exact algorithms for multi-level linear programs cannot be expected to run in polynomial time once p≥2p\ge2p≥2, and even recognising an optimal move is Πp−1\Pi_{p-1}Πp−1​-hard. The bilevel case is an early NP-hardness proof for bilevel linear programming, and the construction (a bookkeeper player and absolute-value gadgets that force binary choices) is a template for hardness reductions to leader–follower problems.

The result is proved in the paper but, to our knowledge, has no machine-checked formalization. A formal development makes precise the solution concept (conditional rather than lexicographic minimisation), which the literature states in several inequivalent ways, and checks a proof whose printed version leaves cases to the reader (the ∧\wedge∧, ¬\neg¬ cases of Lemma 4.1, the universal cases of Lemma 4.4) and applies Lemma 4.3 to a feasible set that is unbounded (the (4.1) variable zzz has no upper bound).

Difficulty

The obvious argument, "each existential player picks a satisfying assignment and each universal player a counterexample", works for the pure binary game J(F)J(F)J(F) (Proposition 3.2) but not for continuous variables: a player may choose fractional values, and the solution sets of a multi-level program need not exist (the §2 Example). The work is in showing that every player is forced to binary choices. Player 1's fractional choices are ruled out only through the robustness estimate of Lemma 4.1 with the weight 10L10L10L, and the existence of optimal solutions at every level has to be established along the induction, since it fails for general three-level programs.

Formalization scope

  • Players are indexed from 000; player iii optimises at level i+1i+1i+1. In J′(F)J'(F)J′(F) the players are 0,…,p0,\dots,p0,…,p as in the paper; in the §2 Example and in J(F)J(F)J(F) the paper's player iii is index i−1i-1i−1. Blocks are 0-based: block k : Fin p is the paper's Xk+1X_{k+1}Xk+1​, owned by player k+1k+1k+1 in J′(F)J'(F)J′(F).
  • solSet encodes the conditional minimisation of (3.7), p. 152. HasValue N w requires SN≠∅S_N\neq\emptysetSN​=∅; the value +∞+\infty+∞ is not modelled.
  • Formulas use ¬,∧,∨\neg,\wedge,\vee¬,∧,∨ (the paper rewrites →\to→ as ¬G1∨G2\neg G_1\vee G_2¬G1​∨G2​); the length counts atoms and connectives. Data are real; the paper's rationality assumption plays no role in the statements.
  • The bookkeeper controls the x(G)x(G)x(G) and has criterion "+z+z+z on (4.1) gadgets, −z-z−z on (4.6) gadgets, and +x(G)+x(G)+x(G) for each non-atomic subformula," as stated on p. 155. The (4.1) variable zzz has no upper bound. The constant 111 of (4.4)/(4.7) is the variable uuu with u=1u=1u=1.
  • "Binary value, zero iff F∈BpF\in B_pF∈Bp​" is stated exactly: the value is 2np+1−x(F)2n_p+1-x(F)2np​+1−x(F), and x(F)=1x(F)=1x(F)=1 iff (3.3). "All optimal solutions are binary" covers the atom variables and the x(G)x(G)x(G), not the gadget variable zzz, which equals 222 at y=1y=1y=1, ξ=0\xi=0ξ=0.
  • The goal is not trivialisable: its first conjunct asserts that the optimal set is nonempty, which fails for general multi-level programs (the §2 Example), so the remaining conjuncts are not vacuous.

A complete development needs a linear-programming duality argument for Lemma 4.2 (Mathlib has IsExtreme and Set.extremePoints; LP duality is on the platform as a single-level theorem) and an induction over the levels of the game. The definitions MultilevelProgram, solSet and the formula encoding LSys are reusable for other complexity results on hierarchical optimisation. Proofs of any milestone, including the generic Lemmas 4.2–4.3 and the formula Lemmas 3.1 and 4.1, are welcome independently.

Selected references

  • R. G. Jeroslow, The polynomial hierarchy and a simple model for competitive analysis, Mathematical Programming 32 (1985) 146–164. https://doi.org/10.1007/BF01586088
  • W. Candler and R. Townsley, A linear two-level programming problem, Computers & Operations Research 9 (1982) 59–76. https://doi.org/10.1016/0305-0548(82)90006-5
  • J. F. Bard and J. E. Falk, An explicit solution to the multi-level programming problem, Computers & Operations Research 9 (1982) 77–100. https://doi.org/10.1016/0305-0548(82)90007-7
  • L. J. Stockmeyer, The polynomial-time hierarchy, Theoretical Computer Science 3 (1976) 1–22. https://doi.org/10.1016/0304-3975(76)90061-X
13 thms1 active userReviewed

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