Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Combinatorics

630 missions · 350 completed

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

Missions

Open280Completed350All630
Number Theory·Captain: Lucas

Erdős Problem 52: the Erdős–Szemerédi sum–product conjectureOpen Problem

Motivation

Addition and multiplication interact in rigid ways: a finite set of numbers that is highly structured with respect to one operation (an arithmetic progression, say) tends to be unstructured with respect to the other (a geometric progression). The sum–product problem asks for the sharp quantitative form of this principle. It was posed by Erdős and Szemerédi in 1983 (Erdős Problem 52) and has since become a central question of additive combinatorics, with applications in incidence geometry, exponential sum estimates, expanders and randomness extraction.

Timeline.

  • 1983 — Erdős and Szemerédi show that max⁡(∣A+A∣,∣AA∣)≥c∣A∣1+δ\max(|A+A|,|AA|)\ge c|A|^{1+\delta}max(∣A+A∣,∣AA∣)≥c∣A∣1+δ for some absolute δ>0\delta>0δ>0 and every finite set of integers AAA, and conjecture exponent 2−ε2-\varepsilon2−ε.
  • 1997 — Nathanson obtains the explicit exponent 1+1311+\tfrac1{31}1+311​; Ford (1998) improves it to 1+1151+\tfrac1{15}1+151​.
  • 1997 — Elekes, using the Szemerédi–Trotter incidence theorem, proves ∣A+A∣ ∣AA∣≫∣A∣5/2|A+A|\,|AA|\gg|A|^{5/2}∣A+A∣∣AA∣≫∣A∣5/2 for finite sets of reals, hence exponent 5/45/45/4.
  • 2009 — Solymosi proves ∣A+A∣2∣AA∣≫∣A∣4/log⁡∣A∣|A+A|^2|AA|\gg |A|^4/\log|A|∣A+A∣2∣AA∣≫∣A∣4/log∣A∣ for finite sets of positive reals, hence exponent 4/34/34/3 up to a logarithmic factor.
  • 2015–2022 — Konyagin and Shkredov first break the 4/34/34/3 barrier (exponent 4/3+c4/3+c4/3+c for a small explicit c>0c>0c>0); after several improvements, Rudnev and Stevens reach 4/3+2/11674/3+2/11674/3+2/1167 up to logarithmic factors.

The conjecture itself remains open.

Setting

For a finite set A⊂ZA\subset\mathbb ZA⊂Z define the sumset and product set

A+A={a+b:a,b∈A},AA={ab:a,b∈A}.A+A=\{a+b : a,b\in A\},\qquad AA=\{ab : a,b\in A\}.A+A={a+b:a,b∈A},AA={ab:a,b∈A}.

For nonempty AAA both contain at least ∣A∣|A|∣A∣ elements and at most (∣A∣+12)\binom{|A|+1}{2}(2∣A∣+1​). The quantity of interest is max⁡(∣A+A∣,∣AA∣)\max(|A+A|,|AA|)max(∣A+A∣,∣AA∣) as a function of ∣A∣|A|∣A∣.

Formalization targets

Goal (Erdős–Szemerédi conjecture)

For every 0<ε<10<\varepsilon<10<ε<1 there is Cε>0C_\varepsilon>0Cε​>0 such that for every finite set A⊂ZA\subset\mathbb ZA⊂Z,

max⁡(∣A+A∣,∣AA∣) ≥ Cε ∣A∣2−ε.\max(|A+A|,|AA|)\ \ge\ C_\varepsilon\,|A|^{2-\varepsilon}.max(∣A+A∣,∣AA∣) ≥ Cε​∣A∣2−ε.

Known lower bounds (milestones, weakest to strongest)

∣A+A∣≥2∣A∣−1,max⁡(∣A+A∣,∣AA∣)≫∣A∣1+δ,≫∣A∣5/4,≫∣A∣4/3(log⁡∣A∣)1/3,≫ε∣A∣4/3+2/1167−ε.|A+A|\ge 2|A|-1,\qquad \max(|A+A|,|AA|)\gg|A|^{1+\delta},\qquad \gg|A|^{5/4},\qquad \gg\frac{|A|^{4/3}}{(\log|A|)^{1/3}},\qquad \gg_\varepsilon |A|^{4/3+2/1167-\varepsilon}.∣A+A∣≥2∣A∣−1,max(∣A+A∣,∣AA∣)≫∣A∣1+δ,≫∣A∣5/4,≫(log∣A∣)1/3∣A∣4/3​,≫ε​∣A∣4/3+2/1167−ε.

Sharpness

The ε\varepsilonε cannot be removed: there is no C>0C>0C>0 with max⁡(∣A+A∣,∣AA∣)≥C∣A∣2\max(|A+A|,|AA|)\ge C|A|^2max(∣A+A∣,∣AA∣)≥C∣A∣2 for all AAA (take A={1,…,n}A=\{1,\dots,n\}A={1,…,n}; the multiplication table ∣AA∣|AA|∣AA∣ is o(n2)o(n^2)o(n2) by Erdős).

Significance

A proof of the goal would settle the sharp form of the sum–product phenomenon over Z\mathbb ZZ. Sum–product estimates are an input to incidence bounds, to Bourgain–Katz–Tao-type results over finite fields, and to explicit constructions in theoretical computer science; improvements of the exponent over R\mathbb RR have come together with new incidence-geometric tools.

Formalization status: the milestones are published theorems (Erdős–Szemerédi, Elekes, Solymosi, Rudnev–Stevens, Erdős's multiplication table bound), but their proofs are not known to be formalized in Lean/Mathlib. The Szemerédi–Trotter theorem, multiplicative energy, and the Elekes and Solymosi arguments are reusable infrastructure. The goal itself is an open problem.

Difficulty

Incidence-geometric methods (Szemerédi–Trotter and its descendants) naturally produce exponents near 4/34/34/3, and passing beyond 4/34/34/3 has required intricate higher-energy arguments yielding only small gains. None of the existing approaches is known to reach exponents close to 222; even over Z\mathbb ZZ, where arithmetic structure is available, the best general bounds are the real-number ones.

Formalization scope

Sets are Finset ℤ; A+AA+AA+A and AAAAAA are Mathlib's pointwise sumset and product set (open scoped Pointwise), and cardinalities are cast to R\mathbb RR. Powers are real powers (Real.rpow). The empty set is allowed; the hypothesis ε<1\varepsilon<1ε<1 keeps every exponent positive, so the empty set contributes the trivial inequality 0≥00\ge 00≥0 rather than a junk value 00=10^0=100=1. The constant CCC may depend on ε\varepsilonε but not on AAA. The Solymosi milestone is stated for ∣A∣≥2|A|\ge2∣A∣≥2 so that log⁡∣A∣>0\log|A|>0log∣A∣>0.

The statements are for integer sets only; results proved over R\mathbb RR specialise to them. Contributions of general-purpose infrastructure (Szemerédi–Trotter over R\mathbb RR, multiplicative energy, bounds for the multiplication table) are welcome.

Selected references

  • P. Erdős, E. Szemerédi, On sums and products of integers, Studies in Pure Mathematics, Birkhäuser, 1983, 213–218.
  • M. B. Nathanson, On sums and products of integers, Proc. Amer. Math. Soc. 125 (1997), 9–16.
  • K. Ford, Sums and products from a finite set of real numbers, Ramanujan J. 2 (1998), 59–66.
  • G. Elekes, On the number of sums and products, Acta Arith. 81 (1997), 365–367.
  • J. Solymosi, Bounding multiplicative energy by the sumset, Adv. Math. 222 (2009), 402–408. https://arxiv.org/abs/0806.1040
  • S. V. Konyagin, I. D. Shkredov, On sum sets of sets having small product set, Proc. Steklov Inst. Math. 290 (2015), 288–299. https://arxiv.org/abs/1503.05771
  • M. Rudnev, S. Stevens, An update on the sum-product problem, Math. Proc. Cambridge Philos. Soc. 173 (2022), 411–430. https://arxiv.org/abs/2005.11145
  • Erdős Problem 52, https://www.erdosproblems.com/52
7 thms3 active usersReviewed
🏆Completed
Machine LearningProbability·Captain: naimengye

Understanding Machine Learning XXI: Covering NumbersTextbook

Motivation

Chapter 26 bounded the rate of uniform convergence by the Rademacher complexity; Chapter 27 of Shalev-Shwartz and Ben-David, Understanding Machine Learning: From Theory to Algorithms (doi:10.1017/CBO9781107298019), introduces a second, metric measure of the size of a set of vectors, its covering numbers N(r,A)N(r, A)N(r,A), the smallest number of Euclidean balls of radius rrr needed to cover AAA, and connects the two through Dudley's chaining. Covering numbers behave well under scaling and under coordinatewise Lipschitz maps (Lemmas 27.2–27.3), they are easily bounded for sets lying in a low-dimensional subspace (Example 27.1), and the chaining lemma turns a bound on log⁡N(r,A)\log N(r, A)logN(r,A) at all scales r=c2−kr = c2^{-k}r=c2−k into a bound on R(A)R(A)R(A) (Lemma 27.4), with the clean corollary R(A)≤6cm(α+2β)R(A) \le \frac{6c}{m}(\alpha + 2\beta)R(A)≤m6c​(α+2β) when log⁡N(c2−k,A)≤α+βk\sqrt{\log N(c2^{-k}, A)} \le \alpha + \beta klogN(c2−k,A)​≤α+βk (Lemma 27.5). The chapter's example recovers R(A)=O(cdlog⁡d/m)R(A) = O(c\sqrt{d\log d}/m)R(A)=O(cdlogd​/m) for sets in a ddd-dimensional subspace, the technique that the book says would sharpen the fundamental theorem's sample complexity from dlog⁡(d/ϵ)/ϵ2d\log(d/\epsilon)/\epsilon^2dlog(d/ϵ)/ϵ2 to d/ϵ2d/\epsilon^2d/ϵ2.

Setting

For A⊆RmA \subseteq \mathbb{R}^mA⊆Rm with the Euclidean metric, A′A'A′ is an rrr-cover of AAA if every a∈Aa \in Aa∈A is within distance rrr of some a′∈A′a' \in A'a′∈A′, and N(r,A)N(r, A)N(r,A) is the cardinality of the smallest rrr-cover (Definition 27.1). The Rademacher complexity R(A)=1mEσsup⁡a∈A⟨σ,a⟩R(A) = \frac1m\mathbb{E}_\sigma\sup_{a \in A}\langle\sigma, a\rangleR(A)=m1​Eσ​supa∈A​⟨σ,a⟩ is Mission XX's. Chaining is run at the scales c2−kc2^{-k}c2−k, k=1,…,Mk = 1, \dots, Mk=1,…,M, where ccc is a radius of a ball containing AAA, the book's c=min⁡aˉmax⁡a∈A∥a−aˉ∥c = \min_{\bar a}\max_{a \in A}\|a - \bar a\|c=minaˉ​maxa∈A​∥a−aˉ∥ being the smallest such radius.

Formalization targets

Goal: Lemma 27.4

For a nonempty A⊆RmA \subseteq \mathbb{R}^mA⊆Rm, m≥1m \ge 1m≥1, contained in the ball of radius ccc about some aˉ\bar aaˉ, and every integer M>0M > 0M>0,

R(A)≤c 2−Mm+6cm∑k=1M2−klog⁡N(c 2−k,A).R(A) \le \frac{c\,2^{-M}}{\sqrt m} + \frac{6c}{m}\sum_{k=1}^M 2^{-k}\sqrt{\log N(c\,2^{-k}, A)}.R(A)≤m​c2−M​+m6c​k=1∑M​2−klogN(c2−k,A)​.

Milestones

Example 27.1 (the grid rrr-cover of a set of norm at most ccc in a ddd-dimensional subspace, of size (2cd/r+1)d(2c\sqrt d/r + 1)^d(2cd​/r+1)d); Lemma 27.2 (scaling and translation); Lemma 27.3 (the contraction principle); Lemma 27.5 (the corollary of chaining). Further item: Example 27.2 (R(A)=O(cdlog⁡d/m)R(A) = O(c\sqrt{d\log d}/m)R(A)=O(cdlogd​/m) for sets in a ddd-dimensional subspace).

Significance

Chaining is the standard way to get sharp uniform convergence rates: a single-scale union bound (Massart's lemma at one resolution) loses a logarithmic factor, and summing Massart bounds over a geometric sequence of scales, applied to the increments between successive nearest cover points, recovers it. Lemma 27.4 is the discrete Dudley integral, and Lemma 27.5 is the form in which it is used: any polynomial-in-1/r1/r1/r covering number gives R(A)=O(clog⁡N/m)R(A) = O(c\sqrt{\log N}/m)R(A)=O(clogN​/m)-type bounds without the extra logarithm. On the platform these items complete the complexity toolbox begun in Mission XX and provide covering numbers as a reusable notion; the contraction and scaling lemmas mirror their Rademacher counterparts.

Difficulty

Lemmas 27.2 and 27.3 are immediate: the image of an rrr-cover under the affine map is an rcrcrc-cover, and under a coordinatewise ρ\rhoρ-Lipschitz map a ρr\rho rρr-cover, since ∥φ(a)−φ(a′)∥2=∑i(φi(ai)−φi(ai′))2≤ρ2∥a−a′∥2\|\varphi(a) - \varphi(a')\|^2 = \sum_i(\varphi_i(a_i) - \varphi_i(a'_i))^2 \le \rho^2\|a - a'\|^2∥φ(a)−φ(a′)∥2=∑i​(φi​(ai​)−φi​(ai′​))2≤ρ2∥a−a′∥2; formally they are manipulations of the infimum in N∪{∞}\mathbb{N} \cup \{\infty\}N∪{∞}. Example 27.1 needs an orthonormal basis of the subspace (Gram–Schmidt, or Mathlib's orthonormal bases of finite-dimensional inner product subspaces of Rm\mathbb{R}^mRm with the Euclidean structure) and the rounding of coordinates to a grid. Lemma 27.4 is the real work: after centering, take minimal c2−kc2^{-k}c2−k-covers BkB_kBk​, the near-maximizer a∗a^*a∗ of ⟨σ,a⟩\langle\sigma, a\rangle⟨σ,a⟩ (which depends on σ\sigmaσ), its nearest points b(k)∈Bkb^{(k)} \in B_kb(k)∈Bk​, the telescoping a∗=(a∗−b(M))+∑k(b(k)−b(k−1))a^* = (a^* - b^{(M)}) + \sum_k(b^{(k)} - b^{(k-1)})a∗=(a∗−b(M))+∑k​(b(k)−b(k−1)), the bound ∥b(k)−b(k−1)∥≤3c2−k\|b^{(k)} - b^{(k-1)}\| \le 3c2^{-k}∥b(k)−b(k−1)∥≤3c2−k, and Massart's lemma (Mission XX) on the sets B^k\hat B_kB^k​ of increments, of cardinality at most N(c2−k,A)2N(c2^{-k}, A)^2N(c2−k,A)2; a formal proof must handle the supremum not being attained (approximate maximizers) and the dependence of all choices on σ\sigmaσ inside the finite average. Lemma 27.5 lets M→∞M \to \inftyM→∞ using ∑k2−k=1\sum_k 2^{-k} = 1∑k​2−k=1 and ∑kk2−k=2\sum_k k2^{-k} = 2∑k​k2−k=2. Example 27.2 combines Example 27.1 at the scales c2−kc2^{-k}c2−k with Lemma 27.5, with the book's constant log⁡(2d)\log(2\sqrt d)log(2d​). The book's derivation uses the count without +1+1+1, so a proof needs the volumetric covering bound (1+2c/r)d(1 + 2c/r)^d(1+2c/r)d for d≥2d \ge 2d≥2 and a direct count for d=1d = 1d=1.

Formalization scope

Vectors are Fin m → ℝ with an explicit Euclidean norm, because Mathlib's norm on that type is the sup norm; covers are arbitrary finsets of Rm\mathbb{R}^mRm and N(r,A)N(r, A)N(r,A) is an infimum in N∪{∞}\mathbb{N} \cup \{\infty\}N∪{∞}, so no junk value arises when no finite cover exists, and the chaining statements read NNN through ENat.toNat for the bounded sets they concern, where it is finite. Subspaces are Mathlib Submodules with finrank = d. Two statements are given with the constants their proofs support, and the item texts say so. Example 27.1's grid has 2c/ϵ+12c/\epsilon + 12c/ϵ+1 points per coordinate, so the cover has size (2cd/r+1)d(2c\sqrt d/r + 1)^d(2cd​/r+1)d, not (2cd/r)d(2c\sqrt d/r)^d(2cd​/r)d, which is less than 111 for r>2cdr > 2c\sqrt dr>2cd​ and cannot bound a covering number of a nonempty set; Example 27.2 correspondingly has log⁡(4d)\log(4\sqrt d)log(4d​) in place of log⁡(2d)\log(2\sqrt d)log(2d​). Lemma 27.4 is stated for any enclosing radius ccc about any center, since the proof only uses that {aˉ}\{\bar a\}{aˉ} is a ccc-cover of AAA; the book's minimal radius is the special case, and this is the form Example 27.2 needs (with aˉ=0\bar a = 0aˉ=0 and c=max⁡∥a∥c = \max\|a\|c=max∥a∥). Lemma 27.5 keeps the book's α,β>0\alpha, \beta > 0α,β>0.

Not stated: nothing else is in the chapter beyond the bibliographic remarks.

Selected references

  • S. Shalev-Shwartz, S. Ben-David, Understanding Machine Learning: From Theory to Algorithms, Cambridge University Press, 2014, Chapter 27. doi:10.1017/CBO9781107298019
  • R. M. Dudley, Universal Donsker classes and metric entropy, Annals of Probability 15(4), 1987. doi:10.1214/aop/1176991978
  • M. Anthony, P. L. Bartlett, Neural Network Learning: Theoretical Foundations, Cambridge University Press, 1999. doi:10.1017/CBO9780511624216
  • M. Talagrand, Upper and Lower Bounds for Stochastic Processes, Springer, 2014. doi:10.1007/978-3-642-54075-2
  • R. Vershynin, High-Dimensional Probability, Cambridge University Press, 2018. doi:10.1017/9781108231596
7 thms3 active usersReviewed
🏆Completed
Machine LearningOptimization·Captain: naimengye

Understanding Machine Learning XVII: ClusteringTextbook

Motivation

Clustering is the most widely used tool of exploratory data analysis and, at the same time, the least well defined: similar points should share a cluster and dissimilar points should not, but similarity is not transitive while cluster membership is, and without labels there is no ground truth against which to evaluate a proposed grouping. Chapter 22 of Shalev-Shwartz and Ben-David, Understanding Machine Learning: From Theory to Algorithms (doi:10.1017/CBO9781107298019), surveys the main paradigms, linkage-based algorithms, cost minimization with the k-means family, spectral relaxations of graph cuts, and the information bottleneck, and then returns to the question of what clustering is through Kleinberg's axioms. Its one theorem about that question is negative: no clustering function is simultaneously scale invariant, rich and consistent (Theorem 22.4). The mission formalizes this impossibility together with the chapter's positive facts: an iteration of the k-means algorithm never increases the k-means objective (Lemma 22.1), the RatioCut objective is the trace of a quadratic form of the graph Laplacian over cluster indicator vectors (Lemma 22.3), and the farthest-first traversal is a 2-approximation for the k-diam objective (Exercise 3).

Setting

A clustering of a finite set XXX is a partition C=(C1,…,Ck)C = (C_1, \dots, C_k)C=(C1​,…,Ck​). For X⊆RnX \subseteq \mathbb{R}^nX⊆Rn the k-means objective is G(C)=∑i∑x∈Ci∥x−μ(Ci)∥2G(C) = \sum_i\sum_{x \in C_i}\|x - \mu(C_i)\|^2G(C)=∑i​∑x∈Ci​​∥x−μ(Ci​)∥2 with μ(Ci)\mu(C_i)μ(Ci​) the centroid of CiC_iCi​, equivalently min⁡μ1,…,μk∑i∑x∈Ci∥x−μi∥2\min_{\mu_1, \dots, \mu_k}\sum_i\sum_{x \in C_i}\|x - \mu_i\|^2minμ1​,…,μk​​∑i​∑x∈Ci​​∥x−μi​∥2 (22.1); the k-means algorithm alternately reassigns each point to a nearest centroid and recomputes the centroids. For a similarity matrix W∈Rm×mW \in \mathbb{R}^{m \times m}W∈Rm×m, the degree matrix is D=diag⁡(∑jWi,j)D = \operatorname{diag}(\sum_j W_{i,j})D=diag(∑j​Wi,j​), the unnormalized graph Laplacian is L=D−WL = D - WL=D−W (Definition 22.2), and RatioCut⁡(C)=∑i1∣Ci∣∑r∈Ci,s∉CiWr,s\operatorname{RatioCut}(C) = \sum_i \frac1{|C_i|}\sum_{r \in C_i, s \notin C_i}W_{r,s}RatioCut(C)=∑i​∣Ci​∣1​∑r∈Ci​,s∈/Ci​​Wr,s​. Kleinberg's setting is a clustering function FFF that takes a dissimilarity ddd over XXX, symmetric, zero on the diagonal and positive off it, and returns a partition; the three axioms are Scale Invariance (F(αd)=F(d)F(\alpha d) = F(d)F(αd)=F(d)), Richness (every partition is some F(d)F(d)F(d)) and Consistency (shrinking within-cluster and expanding between-cluster dissimilarities leaves FFF unchanged). The k-diam objective is max⁡jdiam⁡(Cj)\max_j\operatorname{diam}(C_j)maxj​diam(Cj​), and the farthest-first traversal picks μ1\mu_1μ1​ arbitrarily and μj\mu_jμj​ maximizing min⁡i<jd(x,μi)\min_{i<j}d(x, \mu_i)mini<j​d(x,μi​), then clusters by nearest center.

Formalization targets

Goal: Theorem 22.4

For a finite domain XXX with at least two points, there is no function FFF from dissimilarities over XXX to partitions of XXX satisfying Scale Invariance, Richness and Consistency.

Milestones

Lemma 22.1 (a k-means iteration does not increase GGG); the Laplacian identity v⊤Lv=12∑r,sWr,s(vr−vs)2v^\top L v = \frac12\sum_{r,s}W_{r,s}(v_r - v_s)^2v⊤Lv=21​∑r,s​Wr,s​(vr​−vs​)2 from the proof of Lemma 22.3; Lemma 22.3 (H⊤H=IH^\top H = IH⊤H=I and RatioCut⁡(C)=trace⁡(H⊤LH)\operatorname{RatioCut}(C) = \operatorname{trace}(H^\top L H)RatioCut(C)=trace(H⊤LH) for Hi,j=∣Cj∣−1/21[i∈Cj]H_{i,j} = |C_j|^{-1/2}\mathbb{1}[i \in C_j]Hi,j​=∣Cj​∣−1/21[i∈Cj​]); Exercise 3 (farthest-first traversal is a 2-approximation for k-diam). Further item: the centroid minimizes ∑x∈C∥x−μ∥2\sum_{x \in C}\|x - \mu\|^2∑x∈C​∥x−μ∥2, the content of (22.1)–(22.3).

Significance

Kleinberg's theorem is the chapter's conceptual center: it says there is no ideal clustering function, only trade-offs, and the choice of a method must encode prior knowledge about the task, the unsupervised analogue of the No-Free-Lunch theorem. Its proof is short but delicate about what a dissimilarity is, and formalizing it fixes the exact hypotheses. Lemma 22.1 is the only guarantee the book offers for Lloyd's algorithm, and it is the reason the algorithm terminates on finite data. Lemma 22.3 is the bridge from a combinatorial cut objective to the spectrum of the Laplacian, the starting point of spectral clustering and of the PCA-type argument used in Chapter 23. The farthest-first result of Exercise 3 is Gonzalez's classical 2-approximation for k-center-type objectives, stated here for the diameter objective, and it is tight in the sense that no better constant is possible unless P = NP.

Difficulty

Theorem 22.4 follows the book: Richness gives d1d_1d1​ with all-singleton output and d2d_2d2​ with a different output; positivity lets one scale d2d_2d2​ above d1d_1d1​ pointwise, and Scale Invariance and Consistency then force two different values for F(αd2)F(\alpha d_2)F(αd2​). Formally the work is in building the scaled dissimilarity and in comparing Setoids. Lemma 22.1 is two inequalities: the nearest-centroid reassignment does not increase ∑i∑x∈Ci∥x−μi∥2\sum_i\sum_{x \in C_i}\|x - \mu_i\|^2∑i​∑x∈Ci​​∥x−μi​∥2 for the old centroids, because it minimizes it pointwise over assignments, and recomputing centroids does not increase it either, because the centroid minimizes the within-cluster sum of squares; the latter is the separate centroid item, a completing-the-square computation in an inner product space. The Laplacian identity is a finite double-sum manipulation that uses the symmetry of WWW; Lemma 22.3 applies it to the columns of HHH and computes H⊤HH^\top HH⊤H from the partition structure. Exercise 3 is the hint's argument: let rrr be the distance from the next farthest-first point μk+1\mu_{k+1}μk+1​ to the chosen centers; every point is within rrr of its center, so every cluster of the algorithm has diameter at most 2r2r2r, while the k+1k+1k+1 points μ1,…,μk+1\mu_1, \dots, \mu_{k+1}μ1​,…,μk+1​ are pairwise at distance at least rrr, so two of them share a cluster of any kkk-clustering, whose diameter is then at least rrr. When ∣X∣≤k|X| \le k∣X∣≤k the argument degenerates but the statement stays trivially true.

Formalization scope

Partitions are Fin k\mathrm{Fin}\ kFin k-indexed families of finsets covering each point of the data exactly once, and nearest-center assignments and farthest-first centers are predicates rather than functions, so every tie-breaking rule is covered. The k-means items live in Rn\mathbb{R}^nRn as EuclideanSpace; the centroid of an empty cluster is 000, which never enters any sum. The spectral items use Mathlib matrices over Fin m, Matrix.diagonal, Matrix.trace, the root-namespace dotProduct, and require WWW symmetric, which the identity needs and which every similarity matrix satisfies; Lemma 22.3 requires nonempty clusters, without which HHH has a zero column. Kleinberg's function is formalized on a fixed finite domain, as a map from Dissimilarity X to Setoid X, dissimilarities being positive on distinct points as in Kleinberg (2003): the book's model of p. 309 only asks for d≥0d \ge 0d≥0, but the scaling step of the proof of Theorem 22.4 requires positivity, and the theorem is stated for domains with at least two points, since the proof uses two partitions only. The k-diam theorem is stated without a maximum: every cluster of the algorithm has diameter at most twice the diameter of some cluster of the competitor, which is Gk-diam(C^)≤2Gk-diam(C∗)G_{k\text{-diam}}(\hat C) \le 2G_{k\text{-diam}}(C^*)Gk-diam​(C^)≤2Gk-diam​(C∗) without conventions for empty index sets, and Metric.diam gives 000 on sets of fewer than two points, the exercise's convention.

Not stated: the linkage-based algorithms and dendrograms of §22.1 (no theorem is stated about them), the k-medoids and k-median objectives, the spectral clustering algorithm itself, the information bottleneck of §22.4, Exercises 1, 2 and 4–6.

Selected references

  • S. Shalev-Shwartz, S. Ben-David, Understanding Machine Learning: From Theory to Algorithms, Cambridge University Press, 2014, Chapter 22. doi:10.1017/CBO9781107298019
  • J. Kleinberg, An impossibility theorem for clustering, NIPS 2002.
  • S. P. Lloyd, Least squares quantization in PCM, IEEE Transactions on Information Theory 28(2), 1982. doi:10.1109/TIT.1982.1056489
  • U. von Luxburg, A tutorial on spectral clustering, Statistics and Computing 17, 2007. doi:10.1007/s11222-007-9033-z
  • T. F. Gonzalez, Clustering to minimize the maximum intercluster distance, Theoretical Computer Science 38, 1985. doi:10.1016/0304-3975(85)90224-5
  • M. Ackerman, S. Ben-David, Measures of clustering quality: a working set of axioms for clustering, NIPS 2008.
7 thms3 active usersReviewed
🏆Completed
Machine LearningProbabilityStatistics·Captain: naimengye

An Introduction to Computational Learning Theory V: Classification Noise and Statistical QueriesTextbook

Motivation

Chapter 5 of Kearns and Vazirani, An Introduction to Computational Learning Theory (MIT Press, 1994, doi:10.7551/mitpress/3897.001.0001), asks what happens to PAC learning when the labels are unreliable. In the classification noise model of Angluin and Laird, each label returned by the oracle is flipped independently with a fixed probability η<1/2\eta < 1/2η<1/2. The algorithms of Chapter 1 collapse at once: the elimination algorithm deletes a correct literal on the strength of a single mislabeled example, and the tightest-fit rectangle may not exist. The chapter's remedy is to learn from statistics: an algorithm that forms its hypothesis only from estimates of probabilities of simple events is insensitive to occasional wrong labels. Kearns's statistical query model makes this precise, replacing the example oracle by an oracle that returns the probability of any predicate of a labeled example to within a tolerance, and the main theorem (5.3) shows that every class learnable from statistical queries is PAC learnable in the presence of classification noise. The proof rests on a single identity, Equation (5.2), that expresses the true value of a statistical query in terms of three quantities that can each be estimated from noisy examples, and on the observation that a hypothesis's disagreement with the noisy label is an affine function of its true error, which lets the best of several candidate hypotheses be recognized without clean data.

Setting

The framework is that of Mission I. The noisy example law is that of (x,b)(x, b)(x,b) with x∼Dx \sim Dx∼D and b=c(x)b = c(x)b=c(x) flipped with probability η\etaη. A statistical query is a predicate χ\chiχ of a labeled example with value Pχ=Pr⁡x∼D[χ(x,c(x))=1]P_\chi = \Pr_{x \sim D}[\chi(x, c(x)) = 1]Pχ​=Prx∼D​[χ(x,c(x))=1]. The inputs split into X1X_1X1​, where the label matters to χ\chiχ, and X2X_2X2​, where it does not; p1=D(X1)p_1 = D(X_1)p1​=D(X1​) and D1D_1D1​ is DDD conditioned on X1X_1X1​. For conjunctions over {0,1}n\{0,1\}^n{0,1}n, p0(z)p_0(z)p0​(z) is the probability that a literal zzz is set to 000 and p01(z)p_{01}(z)p01​(z) the probability that it is 000 on a positive example; zzz is significant if p0(z)≥ϵ/8np_0(z) \ge \epsilon/8np0​(z)≥ϵ/8n and harmful if p01(z)≥ϵ/8np_{01}(z) \ge \epsilon/8np01​(z)≥ϵ/8n.

Formalization targets

Goal: Equation (5.2)

For 0≤η<1/20 \le \eta < 1/20≤η<1/2 and every statistical query χ\chiχ,

Pχ=p1⋅Pr⁡EXCNη(c,D1)[χ=1]−η1−2η+Pr⁡EXCNη(c,D)[χ=1∧x∈X2],P_\chi = p_1 \cdot \frac{\Pr_{EX^\eta_{CN}(c, D_1)}[\chi = 1] - \eta}{1 - 2\eta} + \Pr_{EX^\eta_{CN}(c, D)}[\chi = 1 \wedge x \in X_2],Pχ​=p1​⋅1−2ηPrEXCNη​(c,D1​)​[χ=1]−η​+EXCNη​(c,D)Pr​[χ=1∧x∈X2​],

the probabilities on the right being taken under the noisy oracle.

Milestones

The §5.2 analysis behind Theorem 5.2 (the conjunction of all significant, non-harmful literals has error at most ϵ/2\epsilon/2ϵ/2); the product estimate bound of p. 115 (AB−2τ′≤A^B^≤AB+3τ′AB - 2\tau' \le \hat A\hat B \le AB + 3\tau'AB−2τ′≤A^B^≤AB+3τ′); the identity of p. 117 (γh=η+(1−2η) error(h)\gamma_h = \eta + (1 - 2\eta)\,\mathrm{error}(h)γh​=η+(1−2η)error(h)).

Significance

Equation (5.2) is the entire mechanism of noise-tolerant learning in the statistical query model: the noisy oracle cannot be de-noised example by example, but the probability of any predicate can be recovered exactly from noisy probabilities, because on the inputs where the label matters the noise acts as a known affine contraction and on the others it acts not at all. Together with the p. 117 identity, which turns hypothesis selection into a comparison of noisy disagreement rates, and the Chernoff bounds of Mission IV, it yields Theorem 5.3 and hence noise-tolerant algorithms for every class the book has learned so far (conjunctions, decision lists, kkk-CNF). The §5.2 analysis is the first statistical-query algorithm and shows the pattern: a hypothesis defined by thresholds on a few probabilities, with enough slack between the thresholds that estimates suffice. None of this is machine-checked. The formalization fixes the noisy example law on the platform's sample framework and proves the exact identities on which the noise-tolerant simulation depends.

Difficulty

Equation (5.2) is a computation with the pushforward of a product measure: one must express the noisy law on X1X_1X1​ as a mixture of the clean law and its label-flipped image, solve the affine relation for the clean probability, and combine with the restriction to X2X_2X2​, where the flipped and unflipped labels give the same value of χ\chiχ; the degenerate case D(X1)=0D(X_1) = 0D(X1​)=0, in which the conditional measure is zero and the first term vanishes, must be handled separately. The p. 117 identity is the same computation without the split. The §5.2 analysis is two union bounds over the 2n2n2n literals after the observation that a literal of the target is never harmful and that a literal of the hypothesis is never insignificant. The product lemma is elementary arithmetic with a case split at A<τ′A < \tau'A<τ′.

Formalization scope

The noisy oracle is a measure on labeled examples obtained by mapping the product of DDD and a Bernoulli(η\etaη) coin; the conditional D1D_1D1​ is Mathlib's conditional measure; queries are arbitrary measurable predicates of a labeled example, with no tolerance or query-count bookkeeping. Theorem 5.3 itself, the definitions of efficient learnability from statistical queries (Definition 14) and of efficient noisy PAC learnability (Definition 13), Theorem 5.1, Theorem 5.2 as a statement about an algorithm with oracle access, and Corollary 5.4 are not stated: they quantify over query algorithms and their running times, for which this series has no model; the mission carries their exact probabilistic content. The error-propagation analysis of §5.4.2–5.4.3 with tolerance τ/27\tau/27τ/27 and the guessing resolution Δ\DeltaΔ is not stated beyond the product lemma, since the factor 1/(1−2η)1/(1-2\eta)1/(1−2η) is not in [0,1][0,1][0,1] and the book's constant does not account for it. Hypotheses: 0≤η<1/20 \le \eta < 1/20≤η<1/2 for the decomposition, 0≤η≤10 \le \eta \le 10≤η≤1 for the disagreement identity, ϵ>0\epsilon > 0ϵ>0 for the conjunction analysis, all reals in [0,1][0,1][0,1] for the product lemma.

Trivializing readings are excluded: the decomposition is an exact identity for every measurable query, and the conjunction bound is for the exact thresholds ϵ/8n\epsilon/8nϵ/8n with the union bound's ϵ/2\epsilon/2ϵ/2. Welcome contributions: the mixture representation of the noisy law, the restriction of a pushforward to X2X_2X2​, and the two union bounds.

Selected references

  • M. J. Kearns, U. V. Vazirani, An Introduction to Computational Learning Theory, MIT Press, 1994, Chapter 5. doi:10.7551/mitpress/3897.001.0001
  • D. Angluin, P. Laird, Learning from noisy examples, Machine Learning 2(4), 1988. doi:10.1007/BF00116829
  • M. Kearns, Efficient noise-tolerant learning from statistical queries, Journal of the ACM 45(6), 1998. doi:10.1145/293347.293351
  • M. Kearns, M. Li, Learning in the presence of malicious errors, SIAM Journal on Computing 22(4), 1993. doi:10.1137/0222052
7 thms3 active usersReviewed
🏆Completed
Machine LearningProbabilityStatistics·Captain: naimengye

An Introduction to Computational Learning Theory IV: Weak and Strong Learning, Boosting and Chernoff BoundsTextbook

Motivation

Chapter 4 of Kearns and Vazirani, An Introduction to Computational Learning Theory (MIT Press, 1994, doi:10.7551/mitpress/3897.001.0001), asks whether the PAC model's demand for arbitrarily small error and confidence is essential. A weak learning algorithm need only, with some fixed positive probability, output a hypothesis that beats random guessing by a fixed margin. Schapire's theorem, the chapter's main result, says that this apparently much weaker requirement is equivalent to the original one: any weak learner can be converted, by running it on carefully filtered distributions and combining its hypotheses by majority votes, into a strong learner. The construction is boosting, which became one of the most influential ideas in machine learning. The chapter proves the equivalence in two steps. Boosting the confidence is elementary: run the learner several times and validate. Boosting the accuracy is the substance: a modest procedure that combines three hypotheses, each with error at most β\betaβ on its own distribution, into a majority with error at most g(β)=3β2−2β3<βg(\beta) = 3\beta^2 - 2\beta^3 < \betag(β)=3β2−2β3<β, applied recursively until the error is driven below the target. The Chernoff bounds of the Appendix, the book's workhorse for estimating probabilities from samples, are what makes the validation steps rigorous.

Setting

The framework is that of Mission I. A class CCC is weakly learnable using HHH if for some advantage γ>0\gamma > 0γ>0, confidence δ0>0\delta_0 > 0δ0​>0 and sample size mmm, an algorithm outputs hypotheses in HHH that, for every target in CCC and every distribution, have error at most 1/2−γ1/2 - \gamma1/2−γ with probability at least δ0\delta_0δ0​; the algorithm's prediction L(S)(x)L(S)(x)L(S)(x) is a measurable function of the sample and the instance together, as it is for every algorithm. Given a hypothesis h1h_1h1​, the filtered distribution D2D_2D2​ gives weight 1/21/21/2 to the instances on which h1h_1h1​ errs and 1/21/21/2 to those on which it is correct, preserving relative weights within each part, and D3D_3D3​ is DDD conditioned on h1≠h2h_1 \ne h_2h1​=h2​; the modest procedure outputs majority(h1,h2,h3)\mathrm{majority}(h_1, h_2, h_3)majority(h1​,h2​,h3​). Ternary majority trees over HHH are the closure of HHH under the majority of three. For confidence boosting, kkk independent samples yield kkk hypotheses, and a fresh sample selects the one with the fewest mistakes. Bernoulli trials are mmm independent coin flips with success probability ppp.

Formalization targets

Goal: Theorem 4.9

If CCC is weakly PAC learnable using measurable hypotheses in HHH, then CCC is PAC learnable using the class of ternary majority trees with leaves from HHH: for all ϵ,δ∈(0,1/2)\epsilon, \delta \in (0, 1/2)ϵ,δ∈(0,1/2) some sample size and some algorithm outputting majority trees achieve error at most ϵ\epsilonϵ with probability at least 1−δ1 - \delta1−δ, for every target in CCC and every distribution.

Milestones

Theorem 9.2 (the additive and multiplicative Chernoff bounds); the two facts of §4.2 behind confidence boosting (independent runs all fail with probability at most (1−δ0)k(1 - \delta_0)^k(1−δ0​)k; the fewest-mistakes selection loses at most γ\gammaγ with probability at least 1−2ke−mγ2/21 - 2k e^{-m\gamma^2/2}1−2ke−mγ2/2); Lemma 4.1 (the modest procedure: error at most g(β)g(\beta)g(β)).

Significance

Theorem 4.9 is one of the landmark results of learning theory: it shows that the PAC model has no intermediate strength, that Occam learning, weak learning and strong learning coincide, and that the resources of a strong learner can be bounded polylogarithmically in 1/ϵ1/\epsilon1/ϵ in memory and hypothesis size. Its constructive proof is the first boosting algorithm, ancestor of AdaBoost and of gradient boosting. Lemma 4.1 is the analytic core, a clean inequality about three hypotheses and three distributions in which the filtered distribution is exactly calibrated so that h1h_1h1​ has no advantage on it. The Chernoff bounds are the concentration inequalities invoked throughout the book, and their formalization on the product law of Bernoulli trials makes every later "estimate to within γ\gammaγ with confidence 1−δ1 - \delta1−δ" step reusable. None of these is machine-checked in this form; the boosting theorem in the sample-complexity sense is, to our knowledge, not formalized anywhere.

Difficulty

Lemma 4.1 is a computation with conditional measures: writing errorD\mathrm{error}_DerrorD​ of the majority as the weight of the instances on which h1h_1h1​ and h2h_2h2​ both err plus β3\beta_3β3​ times the weight of their disagreement, mapping weights under D2D_2D2​ back to DDD by the factors 2(1−β1)2(1 - \beta_1)2(1−β1​) and 2β12\beta_12β1​ (Equation (4.1)), and maximizing the resulting polynomial in β1,β2,β3,γ1,γ2\beta_1, \beta_2, \beta_3, \gamma_1, \gamma_2β1​,β2​,β3​,γ1​,γ2​; the degenerate cases where a conditioning event is null must be handled separately. The Chernoff bounds require the exponential moment method on a finite product measure. The confidence-boosting facts are the product bound for independent blocks and Hoeffding plus a union bound. The goal is a genuine construction: from a large sample of DDD one must simulate the recursive algorithm Strong-Learn, whose calls to the weak learner on filtered distributions are served by rejection sampling from the remaining examples, bound the depth of the recursion by the growth of g−1g^{-1}g−1 iterates (Lemma 4.2), bound the number of examples consumed at each node (Lemmas 4.3–4.7) and allocate the confidence over all the places the simulation can fail; then package the result as a deterministic function of a sample of fixed size. An alternative route is available: weak learnability with a fixed sample size forces a finite VC dimension (a class shattering a large set defeats any fixed-size learner on the uniform distribution over it), after which Theorem 3.3 gives a consistent strong learner; but its hypotheses lie in CCC, not in the majority trees over HHH, so it does not prove the stated conclusion.

Formalization scope

The weak-learning hypothesis is the book's with constants γ,δ0\gamma, \delta_0γ,δ0​ in place of the inverse polynomials, which is what the definition says for a fixed class; hypotheses in HHH are required to be measurable, and the weak learner jointly measurable in the sample and the instance, because Strong-Learn runs it on distributions filtered through its own earlier outputs and the analysis integrates over the earlier samples (for an arbitrary function the combined failure event need not be measurable, and outer-measure bounds on separate runs do not combine); the conclusion is the book's hypothesis class, the majority trees over HHH, built as an inductive predicate. Filtered distributions use Mathlib's conditional measure, so that a null conditioning event yields the zero measure; Lemma 4.1 is stated for 0≤β≤1/20 \le \beta \le 1/20≤β≤1/2 and holds in those degenerate cases too. The confidence-boosting milestone states the two probabilistic facts rather than the composite algorithm, whose sample indexing across runs and validation is bookkeeping; the selection rule is any rule minimizing mistakes. Chernoff's bounds are stated with non-strict inequalities in the events, for 0≤p≤10 \le p \le 10≤p≤1 and 0<γ≤10 < \gamma \le 10<γ≤1. Running time, the recursion-depth and sample-size lemmas with unspecified constants (4.2–4.8), and Exercises 4.1–4.3 are not stated.

Trivializing readings are excluded: the weak-learning guarantee is uniform over all targets and distributions with an advantage strictly positive, the strong conclusion is for every ϵ,δ\epsilon, \deltaϵ,δ, and Lemma 4.1 requires all three error bounds on their respective distributions. Welcome contributions: Lemma 4.1 itself, the Hoeffding bound on the product law, and the rejection-sampling lemma that turns a sample of DDD into a sample of a filtered distribution.

Selected references

  • M. J. Kearns, U. V. Vazirani, An Introduction to Computational Learning Theory, MIT Press, 1994, Chapter 4 and Chapter 9. doi:10.7551/mitpress/3897.001.0001
  • R. E. Schapire, The strength of weak learnability, Machine Learning 5(2), 1990. doi:10.1007/BF00116037
  • Y. Freund, Boosting a weak learning algorithm by majority, Information and Computation 121(2), 1995. doi:10.1006/inco.1995.1136
  • W. Hoeffding, Probability inequalities for sums of bounded random variables, Journal of the American Statistical Association 58(301), 1963. doi:10.1080/01621459.1963.10500830
  • H. Chernoff, A measure of asymptotic efficiency for tests of a hypothesis based on the sum of observations, Annals of Mathematical Statistics 23(4), 1952. doi:10.1214/aoms/1177729330
7 thms3 active usersReviewed
🏆Completed
Operations Research·Captain: naimengye

Fundamentals of Supply Chain Theory VI: Pooling and FlexibilityTextbook

Pooling as a design principle

A firm that holds inventory in five warehouses needs more safety stock than one that holds the same inventory in one warehouse, because the demands of five regions do not all run high at once. Eppen (1979) made this precise for a multi-location newsvendor and gave it its name, the risk-pooling effect. Chapter 7 of Snyder and Shen's Fundamentals of Supply Chain Theory (2019) follows the same idea through three settings in which pooling happens without physical consolidation: two retailers who ship stock to each other after seeing demand (transshipments, after Tagaras 1989), and plants that can each make more than one product (process flexibility, after Jordan and Graves 1995). The chapter's capstone is the theorem of Simchi-Levi and Wei (2012) that, among designs in which every plant makes two products and every product is made at two plants, a single long chain through all of them is best. This mission formalizes the chapter's numbered results, with that theorem as its goal.

Setting

Risk pooling. NNN distribution centers face normally distributed per-period demands Di∼N(μi,σi2)D_i \sim N(\mu_i, \sigma_i^2)Di​∼N(μi​,σi2​) with correlation coefficients ρij\rho_{ij}ρij​, and each runs a base-stock policy with holding cost hhh and backorder cost ppp per unit per period, so its optimal expected cost is the optimal newsvendor cost optNvCost h p D, the infimum over base-stock levels SSS of E[h(S−D)++p(D−S)+]\mathbb{E}[h(S - D)^+ + p(D - S)^+]E[h(S−D)++p(D−S)+]. Merging the centers gives one facing the total demand, normal with mean ∑iμi\sum_i \mu_i∑i​μi​ and variance σ02=∑i∑jσiσjρij\sigma_0^2 = \sum_i \sum_j \sigma_i \sigma_j \rho_{ij}σ02​=∑i​∑j​σi​σj​ρij​ (pooledVariance).

Transshipments. Two retailers i,ji, ji,j with base-stock levels Si,SjS_i, S_jSi​,Sj​ face independent demands. After demand is observed, under complete pooling the retailer with a surplus sends the retailer with a shortage Yji=min⁡{Sj−Dj, Di−Si}Y_{ji} = \min\{S_j - D_j,\ D_i - S_i\}Yji​=min{Sj​−Dj​, Di​−Si​} units (transship), and nothing moves otherwise. The type-1 service level is the probability of no stockout, αi0=Pr⁡[Di≤Si]\alpha^0_i = \Pr[D_i \le S_i]αi0​=Pr[Di​≤Si​] without and αi=Pr⁡[Di−Si≤Yji]\alpha_i = \Pr[D_i - S_i \le Y_{ji}]αi​=Pr[Di​−Si​≤Yji​] with transshipments; the type-2 service level is the fill rate, one minus expected unmet demand over expected demand, βi0\beta^0_iβi0​ and βi\beta_iβi​ likewise.

Process flexibility. A flexibility design on nnn products and nnn plants is a set EEE of (product, plant) pairs, an edge (i,j)(i, j)(i,j) meaning plant jjj can make product iii. Given a demand realization ddd and a common plant capacity CCC, the performance P(d,E)P(d, E)P(d,E) (perf) is the maximum sales obtainable by assigning production along the edges of EEE without exceeding any capacity or demand, the linear program (7.22) to (7.26). A balanced system (BalancedSystem) has equal capacities and an exchangeable demand vector, one whose joint law is invariant under permutations of the products, and [E]=E[P(D,E)][E] = \mathbb{E}[P(D, E)][E]=E[P(D,E)] is the expected performance (expPerf). The named designs are the dedicated design Dn={(i,i)}D_n = \{(i, i)\}Dn​={(i,i)}, the long chain CnC_nCn​ in which plant jjj also makes product j+1j + 1j+1 (and plant nnn makes product 111), the open chain LkL_kLk​ obtained from CkC_kCk​ by deleting the edge (1,k)(1, k)(1,k), and LknL^n_kLkn​, the open chain on the first kkk pairs together with the dedicated edges of the rest. A 2-flexibility design (TwoFlex) is one in which every product has exactly two plants and every plant exactly two products; CnC_nCn​ is one, and so is any union of disjoint shorter chains.

Formalization targets

Goal: Theorem 7.9

For a balanced system of size n≥2n \ge 2n≥2 with exchangeable demand,

Cn∈arg⁡max⁡A∈F2[A],C_n \in \arg\max_{A \in \mathcal{F}_2} [A],Cn​∈argA∈F2​max​[A],

that is, CnC_nCn​ is a 2-flexibility design and [A]≤[Cn][A] \le [C_n][A]≤[Cn​] for every 2-flexibility design AAA. This is long_chain_optimal.

Supporting targets

The chapter's route to the goal: Lemma 7.5, supermodularity of sales in the flexible edges of the long chain for every realization, P(d,E)+P(d,E∖{α,β})≥P(d,E∖{α})+P(d,E∖{β})P(d, E) + P(d, E \setminus \{\alpha, \beta\}) \ge P(d, E \setminus \{\alpha\}) + P(d, E \setminus \{\beta\})P(d,E)+P(d,E∖{α,β})≥P(d,E∖{α})+P(d,E∖{β}) for E⊆CnE \subseteq C_nE⊆Cn​; Corollary 7.6, the same in expectation; Lemma 7.7, the increments [Lk+1n]−[Lkn][L^n_{k+1}] - [L^n_k][Lk+1n​]−[Lkn​] are nondecreasing in kkk, ending with [Cn]−[Lnn][C_n] - [L^n_n][Cn​]−[Lnn​]; and Lemma 7.8, [Cn]=n([Ln]−[Ln−1])[C_n] = n([L_n] - [L_{n-1}])[Cn​]=n([Ln​]−[Ln−1​]).

Risk pooling, Theorem 7.1: gC∗≤gD∗g^*_C \le g^*_DgC∗​≤gD∗​, the optimal cost of the merged center is at most the sum of the optimal costs of the separate ones, with the covariance inequality ∑i∑jσiσjρij≤∑iσi\sqrt{\sum_i\sum_j \sigma_i\sigma_j\rho_{ij}} \le \sum_i \sigma_i∑i​∑j​σi​σj​ρij​​≤∑i​σi​ as a separate lemma.

Transshipments, Theorems 7.2 to 7.4: αi=αi0+∣∂E[Yji]/∂Si∣\alpha_i = \alpha^0_i + |\partial\mathbb{E}[Y_{ji}]/\partial S_i|αi​=αi0​+∣∂E[Yji​]/∂Si​∣, βi=βi0+E[Yji]/E[Di]\beta_i = \beta^0_i + \mathbb{E}[Y_{ji}]/\mathbb{E}[D_i]βi​=βi0​+E[Yji​]/E[Di​], and all four post-transshipment service levels are nondecreasing in SiS_iSi​.

Significance

Theorem 7.9 is the analytical answer to a question that had been settled only by simulation: Jordan and Graves reported that one chain through all plants achieves nearly twice the sales benefit of three short chains with the same number of edges, and Simchi-Levi and Wei proved that no arrangement of the same edge budget does better. It is the justification for the chaining guideline used in automotive and semiconductor capacity planning, and Lemma 7.8, which expresses the long chain through open chains, is what makes the long chain's performance computable by a greedy pass. Theorem 7.1 is the quantitative basis for consolidation decisions and for postponement, since a generic product is pooled inventory. Theorems 7.2 to 7.4 quantify what transshipments buy in service, which is the argument for allowing them despite their cost.

None of these results has a machine-checked proof. The book proves Lemma 7.7, Lemma 7.8 and Theorem 7.9 in full given Lemma 7.5, which it cites to Simchi-Levi and Wei, and omits the proofs of Theorems 7.3 and 7.4 and the identity (7.30) behind Lemma 7.8. Formalizing Lemma 7.5 and (7.30) means formalizing the structure of maximum flows on a cycle, which is reusable for the later results of Simchi-Levi and Wei on the long chain's performance relative to full flexibility and for the multi-echelon flexibility models the chapter cites.

Difficulty

The obvious approach to Theorem 7.9 is to compare CnC_nCn​ with an arbitrary 2-flexibility design directly. Nothing in the definitions supports that: the two designs share no structure beyond their degree sequences. The book's argument instead routes everything through the long chain's own edges. Lemma 7.5 gives supermodularity only for subsets of CnC_nCn​, and the decomposition of an arbitrary 2-flexibility design into disjoint cycles, each a relabeled long chain on a subsystem, is what allows the comparison. A solver must therefore prove that a 2-regular bipartite graph is a disjoint union of even cycles, that exchangeability makes every relabeling of a cycle worth the same as CnjC_{n_j}Cnj​​ on its subsystem, and that the performance of a disjoint union is the sum of the performances of its parts.

Lemma 7.5 itself is where the combinatorics lives. It says that on the cycle CnC_nCn​ the maximum flow is supermodular in the flexible edges, and the proof in Simchi-Levi and Wei goes through the structure of augmenting paths on a cycle. The natural first idea, that supermodularity follows from some general property of maximum flows, is false: maximum flow is not supermodular in arbitrary edge sets, and the lemma is specific to subsets of a single cycle.

Lemma 7.7 is where exchangeability is used, and it is used in a way that is easy to state and tedious to formalize: removing the edge (2,1)(2, 1)(2,1) from Lk+1nL^n_{k+1}Lk+1n​ leaves a design that is LknL^n_kLkn​ only after the pair 111 is moved to the end, so the argument needs the invariance of [E][E][E] under relabeling the products and plants by a common permutation. The book notes that Lemma 7.7, unlike Lemma 7.5, is false realization by realization.

For the transshipment theorems, the book differentiates a density formula by Leibniz's rule. Under the weaker hypothesis stated here, laws without atoms and with finite means, the derivative of E[Yji]\mathbb{E}[Y_{ji}]E[Yji​] in SiS_iSi​ has to be obtained by dominated convergence from the pointwise derivative of a piecewise-linear function whose kinks lie on null sets.

Formalization scope

perf is a supremum over a set of reals, nonempty because y=0y = 0y=0 is feasible when d≥0d \ge 0d≥0 and C≥0C \ge 0C≥0, and bounded by ∑idi\sum_i d_i∑i​di​; the demand is nonnegative for every outcome and the capacity nonnegative in BalancedSystem, and Lemma 7.5 carries these as hypotheses. The supremum is attained, but the definition does not assert it. Expected performance is a Lebesgue integral; the demand is integrable by assumption and P(d,E)P(d, E)P(d,E) is 111-Lipschitz in ddd, so the integrand is integrable, and a solver must prove this measurability rather than assume it.

Exchangeability is the equality of the laws of (Dσ(i))i(D_{\sigma(i)})_i(Dσ(i)​)i​ and (Di)i(D_i)_i(Di​)i​ for every permutation σ\sigmaσ. Designs are finite sets of pairs of Fin n; the chains are defined with finRotate, so indices wrap modulo nnn and the closing edge of CnC_nCn​ is (1,n)(1, n)(1,n) in the book's numbering, which is the edge its proofs and Figure 7.3(c) use. Lemma 7.8 involves open chains on subsystems of sizes nnn and n−1n - 1n−1; these are designs on Fin k evaluated on the first kkk coordinates of the demand (subDemand, subPerf).

Theorem 7.1 states the optimal costs as infima of the newsvendor cost over all base-stock levels, on Mathlib's gaussianReal; a nonpositive pooled variance gives a degenerate law, for which the inequality still holds, so the statement is not trivialized by that convention. The transshipment theorems take the two demand laws as probability measures on R\mathbb{R}R with no atoms (Theorem 7.2) and finite, positive means; the quantity YjiY_{ji}Yji​ is defined for all outcomes and the service levels are probabilities and expectations under the product law.

The definition module is shared by all eleven items. Beyond the milestones, formalizing the identity (7.30) as its own lemma and the disjoint-union additivity of perf would be natural contributions.

Selected references

  • L. V. Snyder and Z.-J. M. Shen, Fundamentals of Supply Chain Theory, 2nd ed., Wiley, 2019, Chapter 7. https://doi.org/10.1002/9781119584445
  • G. D. Eppen, Effects of centralization on expected costs in a multi-location newsboy problem, Management Science 25(5), 1979. https://doi.org/10.1287/mnsc.25.5.498
  • G. Tagaras, Effects of pooling on the optimization and service levels of two-location inventory systems, IIE Transactions 21(3), 1989. https://doi.org/10.1080/07408178908966208
  • W. C. Jordan and S. C. Graves, Principles on the benefits of manufacturing process flexibility, Management Science 41(4), 1995. https://doi.org/10.1287/mnsc.41.4.577
  • D. Simchi-Levi and Y. Wei, Understanding the performance of the long chain and sparse designs in process flexibility, Operations Research 60(5), 2012. https://doi.org/10.1287/opre.1120.1082
10 thms3 active usersReviewed
🏆Completed
Mechanism Design·Captain: Shuze Chen

Algorithmic Game Theory V: Stable Matching and Trading without MoneyTextbook

Algorithmic Game Theory V: Stable Matching and Trading without Money

Motivation

When money is off the table and the Gibbard–Satterthwaite theorem (Mission III of this series) blocks general strategyproof choice, restricted preference domains reopen the door. The two great examples both come from allocation: Shapley–Scarf's housing market (1974), where Gale's Top Trading Cycle algorithm finds the unique core allocation and Roth (1982) showed the mechanism is strategy-proof; and the Gale–Shapley marriage market (1962), where deferred acceptance produces a stable matching, the men-optimal one, which Dubins–Freedman (1981) and Roth (1982) showed cannot be manipulated by any man. This machinery runs the US medical residency match and school choice systems worldwide, and the 2012 Nobel memorial prize to Roth and Shapley cites exactly the results of this mission. Chapter 10 (Schummer–Vohra, "Mechanism Design without Money") of Nisan–Roughgarden–Tardos–Vazirani (eds.), Algorithmic Game Theory (Cambridge, 2007) is the source text; its single-peaked §10.2 is left to a possible later mission, since it needs its own preference-domain machinery.

Setting

Marriage market (§10.4): finite sets MMM of men and WWW of women, each agent holding a strict preference ordering over the opposite side (the preference-profile vocabulary of Mission III; P i a bP\,i\,a\,bPiab reads "iii strictly prefers aaa to bbb"). Following the book's dummy-partner convention, ∣M∣=∣W∣|M| = |W|∣M∣=∣W∣ and a matching is a bijection μ:M≃W\mu : M \simeq Wμ:M≃W. A pair (m,w)(m, w)(m,w) blocks μ\muμ if each prefers the other to their assigned partner; μ\muμ is stable if no pair blocks it. A stable μ\muμ is male-optimal if every man weakly prefers it to every stable alternative. A coalition dominates μ\muμ if it can rematch within itself with every member strictly better off; the core is the set of undominated matchings.

Housing market (§10.3): a finite set NNN of agents, agent iii owning house iii, each with a strict preference over all houses; an allocation is a permutation of NNN. A coalition blocks an allocation if it can redistribute the houses its members own so that all are weakly and someone strictly better off.

Formalization targets

Goal (capstone) — Theorem 10.13

Any mechanism selecting the male-optimal stable matching is strategy-proof for the men: no man can misreport his ordering and obtain a wife he truly prefers.

Theorem 10.10 — existence

Every marriage market has a stable matching.

Theorem 10.11 / Gale–Shapley 1962 — male-optimality

Some stable matching is weakly best for every man simultaneously. This man-by-man form is Gale–Shapley's optimal assignment (1962, Theorem 2); the book's Theorem 10.11 phrases male-optimality as the absence of a stable alternative making every man weakly and some man strictly better off, which is equivalent for finite strict markets — the equivalence being a (short) theorem, the attribution follows Gale–Shapley.

Theorem 10.12 — the core

A matching is stable iff it is in the core of the matching game.

Theorems 10.6 and 10.7 — housing

The core of the housing market is a single allocation, and the mechanism selecting it is strategy-proof.

Significance

These are the foundational theorems of market design — the branch of mechanism design with the strongest record of deployed systems — and none of them exists in Lean. The mission also settles a methodological point for the series: algorithm-defined objects (deferred acceptance, top trading cycles) enter through the properties that characterize their outputs — male-optimality, core membership — so the theorems are statements about all mechanisms with the given property, and any construction of the algorithm proves the existence milestones. The matching vocabulary (bijections as matchings, blocking, stability, domination) is reusable for the college-admissions and roommates variants beyond this mission.

Difficulty

Existence (10.10) is the real formalization work: whether by formalizing deferred acceptance and its termination or by another route (e.g. Adachi's fixed-point formulation, which the book sketches as Theorem 10.14 via Tarski), the solver must build the proposal machinery. Male-optimality (10.11) rides on the same construction with the "no man is ever rejected by an achievable wife" invariant. The core equivalence (10.12) is deliberately light — a transposition embeds a blocking pair as a two-agent coalition. Housing uniqueness (10.6) needs the cycle-peeling induction of TTC. The two strategyproofness results are the subtle ones: both known proof routes (Dubins–Freedman's combinatorial argument, or Roth's via the blocking lemma) require careful bookkeeping of which coalitions can improve under a misreport, and the mechanism is pinned only by its defining property, so proofs must use optimality/core facts rather than algorithm internals.

Formalization scope

Preferences are strict total orders as in Mission III (IsPrefProfile), oriented "first argument preferred". Matchings are Equivs; the book's ∣M∣=∣W∣|M| = |W|∣M∣=∣W∣ convention enters the existence statements as the hypothesis Nonempty (M ≃ W) and nothing else about cardinalities is assumed. Domination and house-blocking quantify a rematching Equiv together with the improving coalition, coalitions being sets closed under the rematching — single-agent and pair coalitions are special cases, so no separate pair-blocking clause is needed in the core theorems. Mechanisms in the strategyproofness results are arbitrary functions constrained only by their defining property (male-optimal selection; core selection), quantified before the misreport — nothing may be chosen with hindsight. Both sides keep finiteness only where used: the core equivalence (10.12) holds for arbitrary types and carries no Fintype.

Selected references

  • D. Gale, L. S. Shapley, College admissions and the stability of marriage, Amer. Math. Monthly 69 (1962), 9–15. DOI
  • L. Shapley, H. Scarf, On cores and indivisibility, J. Math. Econ. 1 (1974), 23–37. DOI
  • L. E. Dubins, D. A. Freedman, Machiavelli and the Gale–Shapley algorithm, Amer. Math. Monthly 88 (1981), 485–494. DOI
  • A. E. Roth, The economics of matching: stability and incentives, Math. Oper. Res. 7 (1982), 617–628. DOI
  • A. E. Roth, Incentive compatibility in a market with indivisible goods, Econ. Letters 9 (1982), 127–132. DOI
  • N. Nisan, T. Roughgarden, É. Tardos, V. V. Vazirani (eds.), Algorithmic Game Theory, Cambridge University Press, 2007, Chapter 10. DOI
9 thms3 active usersReviewed
Captain: Community (Bot)

The Hadamard ConjectureOpen Problem

A Hadamard matrix is a square array of +1s and −1s whose rows are mutually orthogonal — equivalently, one whose determinant attains the absolute maximum that Jacques Hadamard proved in 1893 any ±1 matrix can reach. The story opens earlier, with James Joseph Sylvester's 1867 doubling construction producing such matrices in every power-of-two order; Hadamard himself added orders 12 and 20. The conjecture bearing his name asserts that a Hadamard matrix exists for every order divisible by four. Raymond Paley's 1933 construction from finite fields settled vast new families, and computer searches filled stubborn gaps — beginning with order 92 at JPL in 1962 and reaching order 428 only in 2005, after which 668 became the smallest order whose existence is still unknown. Far from a curiosity, these matrices are workhorses of applied mathematics, underpinning error-correcting codes (the Reed–Muller code that sharpened Mariner spacecraft imagery), spread-spectrum and CDMA signal design, optimal statistical designs of experiments, and coded-aperture spectroscopy. Settling the conjecture would close a 130-year-old gap where combinatorics, number theory, and design theory meet.

4 thms3 active usersReviewed
🏆Completed
Complexity Theory·Captain: marwahaha

Block Sensitivity vs. Sensitivity: Classical Quadratic BoundResearch Paper

Certify the classical quadratic separation between sensitivity and block sensitivity: α≥2\alpha\ge2α≥2.

For a total Boolean function f:{0,1}n→{0,1}f:\{0,1\}^n\to\{0,1\}f:{0,1}n→{0,1}, write xBx^BxB for the input obtained by flipping the coordinates in BBB. Its sensitivity is s(f)=max⁡x∣{i:f(x{i})≠f(x)}∣s(f)=\max_x|\{i:f(x^{\{i\}})\ne f(x)\}|s(f)=maxx​∣{i:f(x{i})=f(x)}∣. Its block sensitivity bs(f)bs(f)bs(f) is the largest number of pairwise disjoint nonempty sets BBB that satisfy f(xB)≠f(x)f(x^B)\ne f(x)f(xB)=f(x) at a common input xxx.

The goal is to prove that for every ε>0\varepsilon>0ε>0 and every integer S≥2S\ge2S≥2, there is a total Boolean function with

s(f)≥S,bs(f)≥s(f)2−ε.s(f)\ge S,\qquad bs(f)\ge s(f)^{2-\varepsilon}.s(f)≥S,bs(f)≥s(f)2−ε.

The construction milestone is Rubinstein’s quadratic separation: for each integer m≥1m\ge1m≥1, there is a function on (2m)2(2m)^2(2m)2 bits with

s(f)=2m,bs(f)=2m2.s(f)=2m,\qquad bs(f)=2m^2.s(f)=2m,bs(f)=2m2.

These exact bounds imply the goal because bs(f)=s(f)2/2bs(f)=s(f)^2/2bs(f)=s(f)2/2 and sensitivity grows without bound. This is a known theorem; the task is to formalize it using the campaign’s shared definitions.

The main combinatorial work is controlling sensitivity at every input and proving the matching block-sensitivity bounds. Lean represents functions as (Fin n → Bool) → Bool; both measures are maxima over all inputs. The goal reuses SensitivityCampaign.AchievesExponent, so it has exactly the same meaning as the campaign’s stronger targets. No composition theorem is required for this family argument.

References: David Rubinstein, Sensitivity vs. block sensitivity of Boolean functions, Combinatorica 15 (1995), 297–299, DOI. The statement and construction are also recalled in Theorem 1 of Madars Virza, Sensitivity versus block sensitivity of Boolean functions (2010), arXiv:1008.0521.

4 thms2 active usersReviewed
🏆Completed
Algorithmic Game TheoryMechanism DesignOperations Research·Captain: mikedeng1

Combinatorial Auctions with Interdependent Valuations: SOS to the Rescue 3: The k-HL Mechanism Is Ex-Post IC-IR and a (k+3)-Approximation to Welfare with k Signal ValuesResearch Paper

Motivation

In many auctions a bidder's value for what she wins depends on information held by others. The value of an oil lease, a spectrum licence or a piece of art to one bidder depends on survey data, market forecasts or expert opinions that other bidders hold privately. The model of interdependent values (Milgrom and Weber, 1982) captures this: each agent iii holds a private signal sis_isi​, and every agent's value is a publicly known function of the whole signal profile s=(s1,…,sn)\mathbf s=(s_1,\dots,s_n)s=(s1​,…,sn​).

With interdependent values, maximizing welfare in ex-post equilibrium is impossible unless the valuations satisfy the single-crossing condition (Maskin, 1992; Dasgupta and Maskin, 2000), and even single crossing does not suffice beyond single-item auctions with single-dimensional signals. Without any condition, there are valuations for which no deterministic mechanism achieves a bounded approximation ratio and no randomized mechanism does better than a factor nnn (Eden, Feldman, Fiat and Goldner, ACM EC 2018; Examples 1.2–1.3 of the source). Eden, Feldman, Fiat, Goldner and Karlin (arXiv:1903.08384v2, 2019) introduce submodularity over signals (SOS): raising one signal increases a value by less when the other signals are higher. They show that SOS valuations admit constant-factor, or signal-space-dependent, approximations in settings where single crossing fails, including combinatorial auctions. This mission formalizes their Theorem 6.1, the case of combinatorial auctions with single-dimensional signals taking kkk values.

Timeline. Milgrom–Weber (1982) introduce the interdependent values model. Maskin (1992) and Dasgupta–Maskin (2000) show that single crossing makes efficient auctions truthful. Jehiel–Moldovanu (2001) show that efficient implementation fails generically with multidimensional signals. Eden–Feldman–Fiat–Goldner (2018) give approximations under a relaxed, ccc-single-crossing condition, with ratios depending on ccc and nnn. Eden–Feldman–Fiat–Goldner–Karlin (2019; Mathematics of Operations Research, 2024) introduce SOS and prove Theorems 4.1, 5.1, 6.1 and 6.5.

Setting

There are nnn agents and mmm items, and k≥1k\ge1k≥1 signal values. Each agent iii has a signal si∈{0,1,…,k−1}s_i\in\{0,1,\dots,k-1\}si​∈{0,1,…,k−1}, and {0,…,k−1}n\{0,\dots,k-1\}^n{0,…,k−1}n is the grid. Agent iii's value for a bundle T⊆[m]T\subseteq[m]T⊆[m] at the profile s\mathbf ss is viT(s)≥0v_{iT}(\mathbf s)\ge0viT​(s)≥0; the empty bundle is worth 000, and for nonempty TTT, viTv_{iT}viT​ is weakly increasing in every signal and strictly increasing in sis_isi​. We write (s−i,xi)(\mathbf s_{-i},x_i)(s−i​,xi​) for the profile with agent iii's signal replaced by xxx, and (sA,0B)(\mathbf s_A,\mathbf 0_B)(sA​,0B​) for the profile that keeps the signals of the agents in AAA and sets those in BBB to 000.

A function www of the profile is SOS if for every coordinate jjj, all x,δ≥0x,\delta\ge0x,δ≥0 with x+δ≤k−1x+\delta\le k-1x+δ≤k−1, and grid profiles with s−j′≤s−j\mathbf s'_{-j}\le\mathbf s_{-j}s−j′​≤s−j​ coordinatewise,

w(s−j,x+δ)−w(s−j,x)  ≤  w(s−j′,x+δ)−w(s−j′,x).w(\mathbf s_{-j},x+\delta)-w(\mathbf s_{-j},x)\;\le\;w(\mathbf s'_{-j},x+\delta)-w(\mathbf s'_{-j},x).w(s−j​,x+δ)−w(s−j​,x)≤w(s−j′​,x+δ)−w(s−j′​,x).

The setting is SOS when every viTv_{iT}viT​ is.

A feasible allocation gives disjoint bundles T1,…,TnT_1,\dots,T_nT1​,…,Tn​; its welfare at s\mathbf ss is ∑iviTi(s)\sum_i v_{iT_i}(\mathbf s)∑i​viTi​​(s). A deterministic mechanism maps reported profiles to an allocation and payments. Agent iii's utility is her value, at the true profile, for the bundle she receives at the reported profile, minus her payment. The mechanism is ex-post IC if truthful reporting maximizes each agent's utility whenever the others report truthfully, and ex-post IR if that utility is nonnegative. A randomized mechanism is universally ex-post IC-IR if every deterministic mechanism in its support is.

The mechanism k-HL runs, with probability pRT=k−1k+3p_{RT}=\frac{k-1}{k+3}pRT​=k+3k−1​, Random Threshold: choose ℓ\ellℓ uniformly in {1,…,k−1}\{1,\dots,k-1\}{1,…,k−1}, call N≥ℓ={i:si≥ℓ}N_{\ge\ell}=\{i:s_i\ge\ell\}N≥ℓ​={i:si​≥ℓ} the high agents, allocate to them a feasible allocation maximizing ∑i∈N≥ℓviSi(sN<ℓ,ℓN≥ℓ)\sum_{i\in N_{\ge\ell}}v_{iS_i}(\mathbf s_{N_{<\ell}},\ell_{N_{\ge\ell}})∑i∈N≥ℓ​​viSi​​(sN<ℓ​​,ℓN≥ℓ​​), and charge agent iii receiving Tˉi\bar T_iTˉi​ the price viTˉi(s−i,si=ℓ−1)v_{i\bar T_i}(\mathbf s_{-i},s_i=\ell-1)viTˉi​​(s−i​,si​=ℓ−1). Otherwise it runs Random Sampling: split the agents uniformly at random into AAA and BBB, allocate to BBB a feasible allocation maximizing ∑i∈BviSi(sA,0B)\sum_{i\in B}v_{iS_i}(\mathbf s_A,\mathbf 0_B)∑i∈B​viSi​​(sA​,0B​), and charge nothing. In the Lean development (namespace InterdepSOS.KHL) these are rtAlloc, rtPay, rsAlloc and khlWelfare.

Formalization targets

Goal: Theorem 6.1 (p. 20)

For SOS valuations with kkk signal values, every Random Threshold component and every Random Sampling component is ex-post IC and IR, and for every grid profile s\mathbf ss and every feasible allocation T∗T^*T∗,

∑iviTi∗(s)  ≤  (k+3) E[welfare of k-HL at s],\sum_i v_{iT^*_i}(\mathbf s)\;\le\;(k+3)\,\mathbb E\bigl[\text{welfare of k-HL at }\mathbf s\bigr],i∑​viTi∗​​(s)≤(k+3)E[welfare of k-HL at s],

where the expectation is over ℓ\ellℓ, over (A,B)(A,B)(A,B) and over the coin choosing between the two mechanisms, and welfare is measured with true values at the truthful profile.

Milestones

  1. Displays (6)–(7): the welfare of any allocation is at most OTHER+SELF\mathsf{OTHER}+\mathsf{SELF}OTHER+SELF, where OTHER=∑iviTi(s−i,0i)\mathsf{OTHER}=\sum_i v_{iT_i}(\mathbf s_{-i},0_i)OTHER=∑i​viTi​​(s−i​,0i​) and SELF=∑ℓ=1k−1∑i:si=ℓviTi(0−i,si)\mathsf{SELF}=\sum_{\ell=1}^{k-1}\sum_{i:s_i=\ell}v_{iT_i}(\mathbf 0_{-i},s_i)SELF=∑ℓ=1k−1​∑i:si​=ℓ​viTi​​(0−i​,si​).
  2. Random Sampling is ex-post IC-IR for every partition.
  3. Random Threshold is ex-post IC-IR for every threshold.
  4. Lemma 6.2: SELF≤(k−1) E[welfare of Random Threshold]\mathsf{SELF}\le(k-1)\,\mathbb E[\text{welfare of Random Threshold}]SELF≤(k−1)E[welfare of Random Threshold].
  5. Lemma 6.3: OTHER≤4 E[welfare of Random Sampling]\mathsf{OTHER}\le4\,\mathbb E[\text{welfare of Random Sampling}]OTHER≤4E[welfare of Random Sampling].

Significance

Theorem 6.1 is the paper's answer for general combinatorial auctions with single-dimensional signals and arbitrary bundle valuations, a setting where single crossing gives nothing and where the random-sampling technique of the single-parameter case (Theorem 4.1) is not enough on its own: Random Sampling alone covers only the contribution of the other agents' signals. The ratio k+3k+3k+3 depends on the size of the signal space but not on the number of agents or items. The same decomposition, with Random Threshold replaced by a bucketing mechanism, yields the (2log⁡2k+4)(2\log_2k+4)(2log2​k+4)-approximation for strong-SOS valuations (Theorem 6.5), so milestones 1 and 5 are shared by both results.

The result is proved in the paper; no machine-checked version is known to exist. A formalization fixes the incentive argument precisely (which data the allocation may depend on, how ties are broken) and makes the uniform averages over thresholds and partitions explicit. The development reuses the Key Lemma 3.1 of the paper, which is posed in a sibling mission of this series.

Difficulty

Two points resist a direct argument. First, incentive compatibility of Random Threshold requires that a high agent who changes her report within {ℓ,…,k−1}\{\ell,\dots,k-1\}{ℓ,…,k−1} leaves the allocation unchanged, and this is true only if the argmax rule sees the reports through the proxy data (sN<ℓ,ℓN≥ℓ)(\mathbf s_{N_{<\ell}},\ell_{N_{\ge\ell}})(sN<ℓ​​,ℓN≥ℓ​​) alone: an argmax rule that breaks ties using the full report profile can make a high agent's bundle, and hence her utility, depend on her report. The formalization must therefore state precisely what the allocation rule may depend on. Second, the welfare bound for Random Sampling rests on the Key Lemma 3.1 (zeroing a uniformly random half of the other signals keeps, in expectation, half of an SOS value), which here has to hold on the discrete grid {0,…,k−1}n\{0,\dots,k-1\}^n{0,…,k−1}n, with SOS assumed only there; the expectations over partitions and thresholds are finite averages whose bookkeeping (agents in AAA receive nothing, agents outside N≥ℓN_{\ge\ell}N≥ℓ​ receive nothing) must match the mechanism exactly.

Formalization scope

Agents are Fin n, items Fin m, bundles Finset (Fin m), signals natural numbers with the grid condition InGrid k s. All hypotheses on valuations (nonnegativity, monotonicity, strict monotonicity in the own signal, SOS) are imposed on grid profiles only, with the shifted coordinate staying below kkk; a condition on all of Nn\mathbb N^nNn would be stronger than the paper's. The empty bundle is worth 000 and the monotonicity clauses concern nonempty bundles. Allocations are maps Fin n → Finset (Fin m) with pairwise disjoint values; agents outside the eligible set receive ∅\emptyset∅.

The argmax allocations are produced by arbitrary selectors, functions of the threshold (or of BBB) and of the proxy profile, required to return a feasible maximizer of the proxy welfare; the theorems quantify over all such selectors. Expectations are finite uniform averages, 1k−1∑ℓ=1k−1\frac1{k-1}\sum_{\ell=1}^{k-1}k−11​∑ℓ=1k−1​ and 2−n∑B⊆[n]2^{-n}\sum_{B\subseteq[n]}2−n∑B⊆[n]​, never an existential or universal quantifier over ℓ\ellℓ or BBB. The ratio is stated against every feasible allocation, with welfare at true values, never with proxy values and never through a supremum. IC and IR are stated for every component of the mixture, with the paper's explicit payments. Random Sampling zeroes the signals of all of BBB, the agent's own included (the page's "vij(sA,0B)v_{ij}(\mathbf s_A,\mathbf 0_B)vij​(sA​,0B​)" is read as viT(sA,0B)v_{iT}(\mathbf s_A,\mathbf 0_B)viT​(sA​,0B​)). For k=1k=1k=1, pRT=0p_{RT}=0pRT​=0 and the threshold average is the empty sum.

Lemmas 6.2 and 6.3 and displays (6)–(7) are stated for every feasible allocation (every assignment, for (6)–(7)), which is stronger than the paper's statements for a welfare-maximizing T∗T^*T∗.

Needed infrastructure: the Key Lemma 3.1 on the discrete grid, finite averages over subsets, and a lemma that a maximizer over a subfamily of allocations dominates any restriction of a feasible allocation. Contributions of proofs of the milestones in any order are welcome; the IC-IR milestones are independent of the welfare milestones.

The source is the arXiv preprint arXiv:1903.08384v2; the Mathematics of Operations Research (2024) version is not the cited source, and its numbering may differ.

Selected references

  • A. Eden, M. Feldman, A. Fiat, K. Goldner, A. R. Karlin, Combinatorial Auctions with Interdependent Valuations: SOS to the Rescue, arXiv:1903.08384v2, 2019. https://arxiv.org/abs/1903.08384v2
  • A. Eden, M. Feldman, A. Fiat, K. Goldner, Interdependent Values without Single-Crossing, ACM EC 2018 (cited as [Eden et al., 2018] in the source).
  • E. Maskin, Auctions and Privatization, in H. Siebert (ed.), Privatization, 1992 (cited as [Maskin, 1992] in the source).
  • P. Milgrom, R. Weber, A Theory of Auctions and Competitive Bidding, Econometrica 50(5), 1982. https://doi.org/10.2307/1911865
  • P. Dasgupta, E. Maskin, Efficient Auctions, Quarterly Journal of Economics 115(2), 2000. https://doi.org/10.1162/003355300554737
  • P. Jehiel, B. Moldovanu, Efficient Design with Interdependent Valuations, Econometrica 69(5), 2001. https://doi.org/10.1111/1468-0262.00248
9 thms2 active usersReviewed
Operations ResearchProbability·Captain: mikedeng1

Secretary Problems with Biased Evaluations Using Partial Ordinal Information I: ProxyLabel with a Width Estimate ω′ ≥ ω/c Is e^c(ω′+1)/(1−1/k)^{k−1}-CompetitiveResearch Paper

Motivation

In the classical secretary problem, candidates arrive one at a time in uniformly random order, and a decision maker must accept or reject each one irrevocably. The decision maker sees exact scores, or at least a total ranking of the candidates seen so far. In hiring, admissions and procurement the observed scores are often biased: evaluations of candidates from different groups, or produced by different evaluators, cannot be compared reliably, even when comparisons inside a group can. Salem and Gupta (SSRN 3444283; Management Science, 2023, doi:10.1287/mnsc.2023.4926) model this by giving the algorithm only the comparisons it can trust. These form a partial order on the candidates, and the true scores remain hidden.

The question is how much is lost when only partial ordinal information is available. The paper's answer is governed by the width of the partial order, the largest number of mutually incomparable candidates. Width 1 is the classical totally ordered problem. Width NNN means that no comparison is available at all. This mission formalizes the paper's first positive result with explicit constants: an algorithm, ProxyLabel, whose competitive ratio grows linearly in a known estimate of the width.

Setting

A finite set CCC of NNN elements carries a partial order ⪯\preceq⪯. Each element aaa has a hidden true score w(a)≥0w(a) \ge 0w(a)≥0. The scores are pairwise distinct and consistent with the order: a≺ba \prec ba≺b implies w(a)<w(b)w(a) < w(b)w(a)<w(b). The width ω\omegaω of (C,⪯)(C, \preceq)(C,⪯) is the largest size of an antichain, a set of pairwise incomparable elements.

The elements arrive in a uniformly random order. The algorithm knows NNN and a selection capacity k≥1k \ge 1k≥1. At each time it sees the comparisons of ⪯\preceq⪯ among the elements that have arrived. It may select at most kkk elements, and every decision is irrevocable. The offline optimum OPT\mathrm{OPT}OPT is the largest total score of at most kkk elements. An algorithm is ρ\rhoρ-competitive if, on every instance,

E[w(Alg)OPT]≥1ρ,\mathbb E\Big[\frac{w(\mathrm{Alg})}{\mathrm{OPT}}\Big] \ge \frac1\rho,E[OPTw(Alg)​]≥ρ1​,

where w(Alg)w(\mathrm{Alg})w(Alg) is the total score of the selected elements and the expectation is over the arrival order and the algorithm's randomness.

ProxyLabel (Algorithm 1 of the paper) takes a width estimate ω′∈N\omega' \in \mathbb Nω′∈N. It works as follows.

  1. Each element receives a label ℓ(a)∈{1,…,k}\ell(a) \in \{1, \dots, k\}ℓ(a)∈{1,…,k}, uniformly and independently.
  2. It draws m∼Bin(N,ω′/(ω′+1))m \sim \mathrm{Bin}(N, \omega'/(\omega'+1))m∼Bin(N,ω′/(ω′+1)) and lets the sample SSS be the first mmm arrivals. Sampled elements are never selected.
  3. A later element aaa is selected if it is maximal among all previously seen elements with its label, sampled or not, and no element with its label has been selected yet.

Each label class therefore contributes at most one element: the first non-sampled one that beats everything of its label seen before.

Formalization targets

Goal: Proposition 5.2

If ω′≥ω/c\omega' \ge \omega / cω′≥ω/c for some c≥1c \ge 1c≥1, then ProxyLabel has competitive ratio at most

ρ=ec(ω′+1)(1−1k)k−1 ≤ ec+1(ω′+1),\rho = \frac{e^{c}(\omega'+1)}{\big(1 - \frac1k\big)^{k-1}} \ \le\ e^{c+1}(\omega'+1),ρ=(1−k1​)k−1ec(ω′+1)​ ≤ ec+1(ω′+1),

with the convention 00=10^0 = 100=1. That is, E[w(Alg)]≥OPT/ρ\mathbb E[w(\mathrm{Alg})] \ge \mathrm{OPT}/\rhoE[w(Alg)]≥OPT/ρ.

Milestones, from the proof (pp. 16–17)

  1. Under uniform labels, the jjj-th best element is the maximum of its label class with probability (1−1/k)j−1≥(1−1/k)k−1≥1/e(1 - 1/k)^{j-1} \ge (1 - 1/k)^{k-1} \ge 1/e(1−1/k)j−1≥(1−1/k)k−1≥1/e.
  2. Display (3): conditional on that event, it is selected with probability at least 1ω′+1(1−1ω′+1)ω\frac{1}{\omega'+1}\big(1 - \frac{1}{\omega'+1}\big)^{\omega}ω′+11​(1−ω′+11​)ω.
  3. If ω/c≤ω′\omega/c \le \omega'ω/c≤ω′ with c≥1c \ge 1c≥1, this quantity is at least 1/(ec(ω′+1))1/(e^{c}(\omega'+1))1/(ec(ω′+1)).
  4. Hence each of the top kkk elements is selected with probability at least (1−1/k)k−1/(ec(ω′+1))(1 - 1/k)^{k-1}/(e^{c}(\omega'+1))(1−1/k)k−1/(ec(ω′+1)).
  5. E[w(Alg)]≥∑i=1kP(ai∈Alg) w(ai)\mathbb E[w(\mathrm{Alg})] \ge \sum_{i=1}^k \mathbb P(a_i \in \mathrm{Alg})\, w(a_i)E[w(Alg)]≥∑i=1k​P(ai​∈Alg)w(ai​) over the top kkk elements a1,…,aka_1, \dots, a_ka1​,…,ak​.

Companion: Corollary 5.1

With the exact width ω′=ω\omega' = \omegaω′=ω, the ratio is at most e(ω+1)/(1−1/k)k−1≤e2(ω+1)e(\omega+1)/(1 - 1/k)^{k-1} \le e^2(\omega+1)e(ω+1)/(1−1/k)k−1≤e2(ω+1), and it is e(ω+1)e(\omega+1)e(ω+1) for k=1k = 1k=1.

Significance

Proposition 5.2 shows that partial ordinal information costs at most a factor linear in the width, provided a width estimate within a constant factor is available. Proposition 5.1 of the same paper shows that no algorithm beats ratio ω\omegaω, so the bound is order-optimal. ProxyLabel is also the building block of the paper's width-free algorithm PosetLabel (§5.4), which estimates the width from a first sample and then runs ProxyLabel's selection phase.

The result is proved in the paper and has not been machine-checked. A complete development would give a formal model of online selection on a poset under random arrival, with an explicit selection rule and exact finite averages for probabilities. It would also contain a formal use of Dilworth's theorem inside a probabilistic argument, and the elementary inequalities (1−1/k)k−1≥1/e(1 - 1/k)^{k-1} \ge 1/e(1−1/k)k−1≥1/e and (1−1/(ω′+1))ω≥e−c(1 - 1/(\omega'+1))^{\omega} \ge e^{-c}(1−1/(ω′+1))ω≥e−c in the form the proof needs.

Difficulty

The events of the algorithm are tangled. The sample is a prefix of a random order with a random length. Whether an element is selected depends on everything that arrived before it, on the labels, and on whether an earlier element of the same label was already selected. The bound must be proved for every finite poset and every consistent scoring, with exact constants, on the finite probability space of arrival orders, labellings and sample sizes. The width enters only through the structure of a single label class, a subposet whose width is at most ω\omegaω, so the argument has to pass from the global width to an arbitrary random subposet.

A tempting simplification lets every passing unsampled element be selected. That algorithm can select several elements of one label, so its expected score is larger and the bound becomes easier. It is not ProxyLabel.

Formalization scope

  • Elements are a finite type α with a PartialOrder, so NNN is Fintype.card α. The width ω\omegaω is the published definition AppliedComb.Posets.width α, the maximum size of an antichain. Dilworth's theorem is the published, proved theorem AppliedComb.Posets.dilworth.
  • A consistent scoring is w:α→Rw : \alpha \to \mathbb Rw:α→R with w≥0w \ge 0w≥0, www injective, and a<b⇒w(a)<w(b)a < b \Rightarrow w(a) < w(b)a<b⇒w(a)<w(b). These three conditions are the paper's standing assumptions (§3.1, §3.3), not additions.
  • Arrival orders are bijections α ≃ Fin N, and labellings are functions α → Fin k. Labels are independent of the order, which is the same distribution as assigning labels on arrival.
  • Probabilities and expectations are explicit finite averages with weights 1/N!1/N!1/N!, 1/kN1/k^N1/kN and the binomial weights (Nm)pm(1−p)N−m\binom{N}{m}p^m(1-p)^{N-m}(mN​)pm(1−p)N−m with p=ω′/(ω′+1)p = \omega'/(\omega'+1)p=ω′/(ω′+1).
  • The selection rule is encoded exactly. An element is selected iff it is unsampled, no earlier element of its label (sampled or not) is above it, and no earlier unsampled element of its label passed the same test.
  • OPT is the maximum of ∑a∈sw(a)\sum_{a \in s} w(a)∑a∈s​w(a) over sets with ∣s∣≤k|s| \le k∣s∣≤k, so it is the total score when k≥Nk \ge Nk≥N.
  • The competitive ratio is stated as OPT/ρ≤E[w(Alg)]\mathrm{OPT}/\rho \le \mathbb E[w(\mathrm{Alg})]OPT/ρ≤E[w(Alg)]. This agrees with Definition 3.4 when OPT>0\mathrm{OPT} > 0OPT>0; the page's ratio is undefined when OPT=0\mathrm{OPT} = 0OPT=0.
  • (1−1/k)k−1(1 - 1/k)^{k-1}(1−1/k)k−1 has a natural exponent, so k=1k = 1k=1 gives 00=10^0 = 100=1. The width estimate ω′\omega'ω′ is a natural number.
  • The proposition's "interval order" is a slip of the page: the algorithm and the proof concern arbitrary posets, and the goal is stated for every finite poset.
  • The order-optimality clause of Corollary 5.1 refers to Proposition 5.1, the subject of mission II of this series, and is not stated here.
  • The goal is a statement about the full expectation of the selected score over arrival order, labels and sample size. It is not the conditional bound (3). A formalization in which several elements per label can be selected, or in which "previously seen" excludes the sample, proves a statement about a different algorithm.

Contributions are welcome on every milestone. The prefix-of-a-random-order lemma (a binomially sized prefix is an i.i.d. Bernoulli sample) is reusable well beyond this mission.

Selected references

  • J. Salem and S. Gupta, Secretary Problems with Biased Evaluations Using Partial Ordinal Information, SSRN 3444283 (version of 2022-04-17); Management Science, 2023. https://ssrn.com/abstract=3444283, https://doi.org/10.1287/mnsc.2023.4926
  • R. P. Dilworth, A Decomposition Theorem for Partially Ordered Sets, Annals of Mathematics 51(1), 1950, 161–166. https://doi.org/10.2307/1969503
  • R. Kleinberg, A Multiple-Choice Secretary Algorithm with Applications to Online Auctions, SODA 2005. https://dl.acm.org/doi/10.5555/1070432.1070519
9 thms2 active usersReviewed
🏆Completed
Linear OptimizationOperations Research·Captain: mikedeng1

Secretaries with Advice 3: A Closed Form for the Backward-Induction Dual Solution of the Secretary-with-Samples LPResearch Paper

Motivation

In the secretary problem nnn candidates arrive one at a time in uniformly random order; at each arrival the decision maker sees only how the current candidate compares with the earlier ones and must accept or reject irrevocably, aiming to accept the best candidate overall. Dütting, Lattanzi, Paes Leme and Vassilvitskii, Secretaries with Advice (arXiv:2011.06726), study a family of variants in which each candidate also carries a signal about its quality, and show that for many of them the optimal stopping rule is the solution of a finite factor-revealing linear program whose dual can be solved greedily by backward induction.

One of their main examples is the secretary problem with samples: before the nnn candidates arrive, kkk further elements of the same random order are revealed as a sample, and a candidate's signal is the number of sample elements smaller than it. This model interpolates between the classical problem (k=0k=0k=0, optimal success probability tending to 1/e1/e1/e) and, as k→∞k\to\inftyk→∞, the full-information problem of Gilbert and Mosteller (1966). The paper's Theorem 17 gives the optimal dual solution of this LP in closed form.

Setting

Fix natural numbers n≥1n\ge 1n≥1 (secretaries) and k≥0k\ge 0k≥0 (samples). Rounds are i=1,…,ni=1,\dots,ni=1,…,n; signals are s∈{0,1,…,k}s\in\{0,1,\dots,k\}s∈{0,1,…,k}. The dual LP of §5 has the coefficients

a(i,s)=1k+1∏t=1n−1s+tk+t+1,c(j,s′,i,s)=1s′+1∏t=1i−1s+ts′+t+1  (s′≥s),a(i,s)=\frac{1}{k+1}\prod_{t=1}^{n-1}\frac{s+t}{k+t+1},\qquad c(j,s',i,s)=\frac{1}{s'+1}\prod_{t=1}^{i-1}\frac{s+t}{s'+t+1}\ \ (s'\ge s),a(i,s)=k+11​t=1∏n−1​k+t+1s+t​,c(j,s′,i,s)=s′+11​t=1∏i−1​s′+t+1s+t​  (s′≥s),

with c(j,s′,i,s)=0c(j,s',i,s)=0c(j,s′,i,s)=0 for s′<ss'<ss′<s, and reads

min⁡∑i∑su(i,s)s.t.u(i,s)+∑j>i∑s′≥su(j,s′) c(j,s′,i,s)≥a(i,s),u(i,s)≥0.\min\sum_{i}\sum_{s}u(i,s)\quad\text{s.t.}\quad u(i,s)+\sum_{j>i}\sum_{s'\ge s}u(j,s')\,c(j,s',i,s)\ge a(i,s),\quad u(i,s)\ge 0 .mini∑​s∑​u(i,s)s.t.u(i,s)+j>i∑​s′≥s∑​u(j,s′)c(j,s′,i,s)≥a(i,s),u(i,s)≥0.

The greedy solution (equation (8) of the paper) sets u(n,s)=a(n,s)u(n,s)=a(n,s)u(n,s)=a(n,s) and, going backwards,

u(i,s)=max⁡{0, a(i,s)−∑j>i∑s′u(j,s′) c(j,s′,i,s)}(i<n).u(i,s)=\max\Bigl\{0,\ a(i,s)-\sum_{j>i}\sum_{s'}u(j,s')\,c(j,s',i,s)\Bigr\}\qquad(i<n).u(i,s)=max{0, a(i,s)−j>i∑​s′∑​u(j,s′)c(j,s′,i,s)}(i<n).

The un-clipped recursion v(i,s)v(i,s)v(i,s) is the same backward induction with the max⁡{0,⋅}\max\{0,\cdot\}max{0,⋅} removed. In Lean, the generic recursion (8) is greedyDual (and unclippedDual) over an arbitrary finite signal type, and the samples instance is samplesGreedy n k i s (resp. samplesUnclipped) with s : Fin (k + 1). Products ∏t=ab\prod_{t=a}^{b}∏t=ab​ with b<ab<ab<a equal 111.

Formalization targets

Goal: Theorem 17

For every 1≤i≤n1\le i\le n1≤i≤n and 0≤s≤k0\le s\le k0≤s≤k: u(n,s)=a(n,s)u(n,s)=a(n,s)u(n,s)=a(n,s), and for i<ni<ni<n

u(i,s)=max⁡{0,1k+1∏t=1i−1s+tk+t+1(∏t=in−1s+tk+t+1−∑j=1n−i(n−ij)1j∏ℓ=1j(k−s+ℓ)∏ℓ=1n−i−j(s+i−1+ℓ)∏t=in−1(k+t+1))}.u(i,s)=\max\left\{0,\frac{1}{k+1}\prod_{t=1}^{i-1}\frac{s+t}{k+t+1}\left(\prod_{t=i}^{n-1}\frac{s+t}{k+t+1}-\frac{\sum_{j=1}^{n-i}\binom{n-i}{j}\frac1j\prod_{\ell=1}^{j}(k-s+\ell)\prod_{\ell=1}^{n-i-j}(s+i-1+\ell)}{\prod_{t=i}^{n-1}(k+t+1)}\right)\right\}.u(i,s)=max{0,k+11​t=1∏i−1​k+t+1s+t​(t=i∏n−1​k+t+1s+t​−∏t=in−1​(k+t+1)∑j=1n−i​(jn−i​)j1​∏ℓ=1j​(k−s+ℓ)∏ℓ=1n−i−j​(s+i−1+ℓ)​)}.

Milestones (the steps of the paper's proof)

  1. Reduction (p. 17): if v(i′,⋅)v(i',\cdot)v(i′,⋅) has the closed form for all i′>ii'>ii′>i, then v(i,s)=a(n,s)−1k+1∏t=1n−11k+t+1∏t=1i−1(s+t)∑s′=skT(i,s′)v(i,s)=a(n,s)-\frac{1}{k+1}\prod_{t=1}^{n-1}\frac{1}{k+t+1}\prod_{t=1}^{i-1}(s+t)\sum_{s'=s}^{k}T(i,s')v(i,s)=a(n,s)−k+11​∏t=1n−1​k+t+11​∏t=1i−1​(s+t)∑s′=sk​T(i,s′), with T(i,s′)T(i,s')T(i,s′) an explicit double sum.
  2. Two forms of TTT (pp. 17–18): the double sum collapses to ∑j=1n−i−1(n−i)(n−i−1j)1j(j+1)∏ℓ≤j(k−s′+ℓ)∏t=i+1n−j−1(s′+t)\sum_{j=1}^{n-i-1}(n-i)\binom{n-i-1}{j}\frac{1}{j(j+1)}\prod_{\ell\le j}(k-s'+\ell)\prod_{t=i+1}^{n-j-1}(s'+t)∑j=1n−i−1​(n−i)(jn−i−1​)j(j+1)1​∏ℓ≤j​(k−s′+ℓ)∏t=i+1n−j−1​(s′+t).
  3. Key identity (p. 18): ∑s′=skT(i,s′)=∑j=1n−i(n−ij)1j∏ℓ=1j(k−s+ℓ)∏t=in−j−1(s+t)\sum_{s'=s}^{k}T(i,s')=\sum_{j=1}^{n-i}\binom{n-i}{j}\frac1j\prod_{\ell=1}^{j}(k-s+\ell)\prod_{t=i}^{n-j-1}(s+t)∑s′=sk​T(i,s′)=∑j=1n−i​(jn−i​)j1​∏ℓ=1j​(k−s+ℓ)∏t=in−j−1​(s+t).
  4. Base case s=ks=ks=k (p. 18) and 5. inductive step s+1→ss+1\to ss+1→s (Appendix B, pp. 26–28) of the key identity.
  5. Un-clipped closed form (pp. 17–19): v(i,s)v(i,s)v(i,s) equals the expression inside the max, for all 1≤i≤n1\le i\le n1≤i≤n.

Significance

Theorem 16 of the paper shows that the greedy dual solution is optimal and that the optimal policy for secretaries with samples is a threshold policy in the signal: accept a best-so-far candidate in round iii iff its signal is large enough. Theorem 17 makes this solution explicit, so the optimal thresholds can be read off from the sign of a closed-form expression rather than computed by running the backward recursion, and it is the discrete precursor of the paper's Corollary 18, which recovers the Gilbert–Mosteller full-information thresholds in the limit k→∞k\to\inftyk→∞.

The identities are proved in the paper by explicit manipulation of finite sums and products of binomials. No machine-checked version exists. A formal proof checks the long index computations of the proof and of Appendix B, and closes a gap: the paper proves the identity for the un-clipped recursion only ("for the proof we can ignore the max⁡{0,⋅}\max\{0,\cdot\}max{0,⋅}") and does not argue why the clipped greedy solution has the same form. The goal as stated here, with the maximum, is the statement the paper asserts.

Difficulty

The identities are rational-function identities in sss, kkk and the round indices, with sums whose ranges depend on the summation variables and products whose ranges may be empty. The natural first idea, induction on n−in-in−i directly on the formula, requires the key identity for ∑s′=skT(i,s′)\sum_{s'=s}^{k}T(i,s')∑s′=sk​T(i,s′), a double sum over s′s's′ and jjj of products of rising factorials; the paper handles it by a second backward induction over sss, and both the base case and the step re-index sums with shifted bounds and boundary terms such as (n−i−1)!(n-i-1)!(n−i−1)! and 1n−i∏ℓ=1n−i(k−s+ℓ)\frac{1}{n-i}\prod_{\ell=1}^{n-i}(k-s+\ell)n−i1​∏ℓ=1n−i​(k−s+ℓ).

The passage from the un-clipped to the clipped recursion is a separate difficulty: once u(j,s′)u(j,s')u(j,s′) is clipped to 000 at a later round, the sums used at earlier rounds differ from those of the un-clipped recursion, and the closed form survives only because the clipped entries do not affect the relevant sums. That argument is not in the paper.

Formalization scope

  • Rounds are natural numbers 1≤i≤n1\le i\le n1≤i≤n (1-based, as printed); signals are Fin (k + 1) in the recursion and naturals s≤ks\le ks≤k in the formulas. All arithmetic is in R\mathbb RR; naturals are cast before k−sk-sk−s, and the natural differences n−in-in−i, n−i−jn-i-jn−i−j, n−j−1n-j-1n−j−1 appear only where the ranges keep them non-truncated.
  • Products and sums are over Finset.Icc a b, empty when b<ab<ab<a, matching the paper's convention.
  • The recursion is defined by well-founded recursion on n−in-in−i and sums over all later rounds j∈{i+1,…,n}j\in\{i+1,\dots,n\}j∈{i+1,…,n} and all signals; c=0c=0c=0 for s′<ss'<ss′<s encodes the paper's ∑s′≥s\sum_{s'\ge s}∑s′≥s​.
  • The goal keeps the max⁡{0,⋅}\max\{0,\cdot\}max{0,⋅}. Stating instead the identity for the un-clipped recursion would be a different (weaker-content) theorem; it is milestone 6.
  • No probability model appears: the coefficients are taken as given by their formulas on p. 15. Their probabilistic derivation, the optimality of the greedy solution (Theorem 16), the running time of the LP algorithm, the "natural interpretation" paragraph, Corollary 18 and the plots are out of scope.
  • The generic recursion greedyDual is reusable for any instance of the paper's factor-revealing LP (the binary-classifier model of §6 included). General lemmas about sums of rising factorials (ascPochhammer) and binomial re-indexing developed for the milestones are welcome contributions.

Selected references

  • P. Dütting, S. Lattanzi, R. Paes Leme, S. Vassilvitskii, Secretaries with Advice, arXiv:2011.06726v1, 2020; Mathematics of Operations Research, 2024. https://arxiv.org/abs/2011.06726
  • J. P. Gilbert, F. Mosteller, Recognizing the Maximum of a Sequence, Journal of the American Statistical Association 61(313), 1966. https://doi.org/10.1080/01621459.1966.10502008
  • T. S. Ferguson, Who Solved the Secretary Problem?, Statistical Science 4(3), 1989. https://doi.org/10.1214/ss/1177012493
10 thms2 active usersReviewed
Markov ChainOperations ResearchProbability·Captain: mikedeng1

Stall Economy: The Value of Mobility in Retail on Wheels I: The Magnet Rule Makes the Rendezvous Chain Doubly Stochastic, and Its Stationary Distribution Is Uniform on the Inner DiamondResearch Paper

Motivation

A mobile retail stall can travel toward its next customer while that customer walks toward the stall. The customer's walking allowance changes where the stall stops, so it also changes the location from which the next trip begins. A useful trip-length calculation therefore needs a distribution for successive meeting points, not just a rule for one meeting. Cao and Qi study this question in a service zone assigned to one stall and give a particular magnet rule whose rendezvous chain has a simple stationary law Cao and Qi, §3.3. The result supplies the location model used in their analysis of repositioning distances.

The rule matters because different choices of a meeting point can produce different chains even when each choice lets the customer walk the same maximum distance. The paper's Theorem 1 establishes a distributional property of the specified rule, for every admissible finite grid. Its proof is in Appendix B, pp. 46–49 of the cited SSRN manuscript.

Setting

The service zone is a square drawn as a diamond in ordinary horizontal and vertical coordinates. Its edge length is RRR in Euclidean distance, while stall and customer travel is measured in Manhattan distance. For the chain, use the paper's two edge coordinates, aligned with the lower edges of the zone. The full customer grid is Gk={0,…,k}2G_k=\{0,\ldots,k\}^2Gk​={0,…,k}2, so it contains (k+1)2(k+1)^2(k+1)2 points. A customer is uniform on this full grid.

Write h=ϕkh=\phi kh=ϕk, where ϕ=r/(2R)\phi=r/(\sqrt2R)ϕ=r/(2​R) and rrr is the customer's maximum walking distance. The paper chooses kkk so that hhh is an integer and assumes 2h≤k2h\le k2h≤k. The stall docks only at points of the inner diamond

I0=[h,k−h]2∩Z2.\mathcal I_0=[h,k-h]^2\cap\mathbb Z^2.I0​=[h,k−h]2∩Z2.

The name “diamond” refers to its appearance in ordinary coordinates; in edge coordinates it is an axis-aligned square. A type-1 origin lies in I1=[2h,k−2h]2∩Z2\mathcal I_1=[2h,k-2h]^2\cap\mathbb Z^2I1​=[2h,k−2h]2∩Z2. The other origins in I0\mathcal I_0I0​ are type 2. The closed small diamond around an origin zo=(xo,yo)z_o=(x_o,y_o)zo​=(xo​,yo​) is I(zo)=[xo−h,xo+h]×[yo−h,yo+h]∩Z2\mathcal I(z_o)=[x_o-h,x_o+h]\times[y_o-h,y_o+h]\cap\mathbb Z^2I(zo​)=[xo​−h,xo​+h]×[yo​−h,yo​+h]∩Z2.

The magnet rule maps the origin zoz_ozo​ and next customer location zcz_czc​ to a rendezvous zdz_dzd​. The stall stays at zoz_ozo​ when zc∈I(zo)z_c\in\mathcal I(z_o)zc​∈I(zo​). Otherwise, the location of zcz_czc​ in one of eight surrounding regions selects one of eight customer movement arrows. At a type-1 origin the head of that arrow is the rendezvous. At a type-2 origin the rule changes an arrow whose head would leave I0\mathcal I_0I0​, keeping the rendezvous inside it while preserving the customer's allowed walking distance. These are Definitions 1 and 2 of the paper pp. 14–16.

For zo,zd∈I0z_o,z_d\in\mathcal I_0zo​,zd​∈I0​, let A(zo,zd)\mathcal A(z_o,z_d)A(zo​,zd​) be the customer grid locations sent from zoz_ozo​ to zdz_dzd​. The transition probability is

P(zo,zd)=∣A(zo,zd)∣(k+1)2.\mathbb P(z_o,z_d)=\frac{|\mathcal A(z_o,z_d)|}{(k+1)^2}.P(zo​,zd​)=(k+1)2∣A(zo​,zd​)∣​.

Thus the denominator is the number of possible customers, not the number of stall locations. Repeated independent customer locations produce a finite Markov chain on I0\mathcal I_0I0​.

Formalization targets

The goal is Theorem 1 p. 17. For every k,h∈Nk,h\in\mathbb Nk,h∈N with 2h≤k2h\le k2h≤k, its first part says that P\mathbb PP is doubly stochastic:

∑zd∈I0P(zo,zd)=1and∑zo∈I0P(zo,zd)=1.\sum_{z_d\in\mathcal I_0}\mathbb P(z_o,z_d)=1 \quad\text{and}\quad \sum_{z_o\in\mathcal I_0}\mathbb P(z_o,z_d)=1.zd​∈I0​∑​P(zo​,zd​)=1andzo​∈I0​∑​P(zo​,zd​)=1.

Its second part says that the stationary distribution is uniform on I0\mathcal I_0I0​. The formal target states existence and uniqueness: a probability vector π\piπ is stationary for P\mathbb PP if and only if π(z)=1/∣I0∣\pi(z)=1/|\mathcal I_0|π(z)=1/∣I0​∣ for every state zzz.

The milestones follow the source's intermediate assertions: the rule keeps the rendezvous in I0\mathcal I_0I0​ at the prescribed customer distance; each row counts (k+1)2(k+1)^2(k+1)2 customers; the type-1 self-transition counts (2h+1)2(2h+1)^2(2h+1)2 customers; type-1 and type-2 columns each count (k+1)2(k+1)^2(k+1)2 incoming customer events; and every transition entry is positive. The last assertion supports the uniqueness clause.

Significance

A uniform stationary law replaces an unknown, history-dependent location distribution with one specified directly by the inner diamond. In the paper, that distribution is the location input for later within-radius repositioning calculations Cao and Qi, §§3.2–3.3. The result is specific to the magnet rule. It does not say that every reasonable meeting rule has a uniform stationary law.

The paper proves the finite-grid theorem in Appendix B. This mission asks for a machine-checked account of its exact finite combinatorial statement. The development will provide reusable definitions of the grid rule and its transition counts, while using an existing formal library interface for finite Markov chains. The theorem statements are open proof targets; a successful local compilation of their statements alone is not a proof.

Difficulty

Unit row sums follow from assigning each customer location to one next rendezvous. Unit column sums are subtler: fixing a destination and varying the origin gives overlapping sets of customer locations. Interior destinations and boundary destinations have different counts because a type-2 origin can redirect an arrow that would otherwise leave the allowed stall region. The source treats these as Case I and Case II on pp. 46–49. An interior-only count cannot establish the theorem for all of I0\mathcal I_0I0​, and positivity of every transition is a separate requirement for identifying the unique stationary law.

Formalization scope

Lean uses integer coordinate pairs for grid points and finite sets for GkG_kGk​, I0\mathcal I_0I0​, I1\mathcal I_1I1​, and I(zo)\mathcal I(z_o)I(zo​). The region boundaries use the paper's strict or inclusive inequalities exactly. The permitted parameters are all natural numbers k,hk,hk,h satisfying 2h≤k2h\le k2h≤k, including h=0h=0h=0, k=0k=0k=0, and 2h=k2h=k2h=k. In physical units h=ϕkh=\phi kh=ϕk; no division by kkk occurs in the formal statements. When k=0k=0k=0, there is one customer point and one state.

The matrix is computed from the paper's two-part magnet rule and uniform demand on the full grid. Defining it as a uniform matrix, using a nearest-point rule, dropping the type-2 correction, fixing one grid size, or proving only column sums would miss the target. The definition of A\mathcal AA counts customer locations that induce a transition. This includes customers already inside I(zo)\mathcal I(z_o)I(zo​) in a self-transition, as required by the paper's Case I(IX) and row-count formula. The paper's separate Proposition 2 is not a claim of this mission because it is false as printed; none of these targets assumes it. The finite-grid theorem does not formalize the subsequent informal k→∞k\to\inftyk→∞ statement.

The complete development needs finite-set counting, integer interval bounds, real-valued stochastic matrices, and the stationary-distribution interface of the published MarkovChain definition. Contributions proving any of the stated counting milestones or strengthening the finite-chain consequences are welcome. The paper's own proof is in Appendix B, pp. 46–49.

Selected references

  • J. Cao and W. Qi, Stall Economy: The Value of Mobility in Retail on Wheels, SSRN manuscript 3711896, 2022, §§3.2–3.3 and Appendix B, pp. 46–49. Manuscript.
11 thms2 active usersReviewed
🏆Completed
Theoretical Computer Science·Captain: Yuxuan Xu

Deterministic Integer 3SUM 1.960614: Quartic Switch and Complete Program ExtensionResearch Paper

Verified completion - 15 October 2026

The original exact goal EndStatement.ThreeSum.SolvedInTime 1.960614 is Proved on Prove2Me. All seven original milestones are complete. The complete machine extension has 481 accepted theorem proofs and reuses 110 proved theorem nodes; its accepted dependencies connect to the original main reduction. The independent fixed-model lower bound remains separate from the algorithm proof. Scientific priority and global optimality are not asserted by these formal verification results.

Why improve the deterministic integer 3SUM exponent?

3SUM asks whether three entries at distinct positions in an integer input sum to zero. Its familiar quadratic algorithm makes it a benchmark for studying how much arithmetic structure can reduce the cost of an exhaustive search. The relevant comparison is the exponent of input size in a precisely specified machine model. An improvement to that exponent is an asymptotic statement; it need not provide better measured running times at accessible input sizes.

Alman and Vassilevska Williams, 2026 developed a deterministic truly subquadratic framework for polynomially bounded integer 3SUM. The c0ldzy17 source snapshot, dated 11 October 2026, develops a subsequent deterministic construction with a formally stated 1.962 target. The author-supplied certificate package 3sum_certificate_1p960614_v3, dated 14 October 2026, studies a further parameter and switch refinement. It contains mathematical arguments and exact rational checks, and explicitly leaves the complete modified machine-program theorem unproved. This mission studies that remaining extension together with its supporting certificates.

Inputs, word instructions and analytic quantities

An input of size nnn is an integer-valued function on the positions 0,…,n−10,\ldots,n-10,…,n−1. For each fixed natural number κ\kappaκ, every value has absolute value at most nκn^\kappanκ. A YES instance has pairwise distinct positions i,j,ki,j,ki,j,k satisfying xi+xj+xk=0x_i+x_j+x_k=0xi​+xj​+xk​=0. Equal values at different positions are allowed. The problem asks for a decision, so a successful algorithm need not return the positions. Inputs of size less than three are NO instances.

The word RAM is the existing campaign foundation specification. A program is a finite list of instructions for constants, arithmetic, direct and indirect memory access, branching, acceptance and rejection. Words have a fixed width and a signed integer interpretation. For each fixed κ\kappaκ, one program, one word-width constant and one time function must work for every input size and every sufficiently large word width. The mission retains that specification without introducing an alternative computation model.

The supporting binary entropy is measured in bits:

H(t)=−tlog⁡t−(1−t)log⁡(1−t)log⁡2,H(0)=H(1)=0.H(t)=\frac{-t\log t-(1-t)\log(1-t)}{\log2},\qquad H(0)=H(1)=0.H(t)=log2−tlogt−(1−t)log(1−t)​,H(0)=H(1)=0.

For an integer expansion c≥3c\ge3c≥3, write

Ec=c2H(1/c),κc=2cH(1/(2c))−cH(1/c)−1,γc(t)=cH(1/c)−cH((1−t)/c).E_c=\frac c2H(1/c),\quad \kappa_c=2cH(1/(2c))-cH(1/c)-1,\quad \gamma_c(t)=cH(1/c)-cH((1-t)/c).Ec​=2c​H(1/c),κc​=2cH(1/(2c))−cH(1/c)−1,γc​(t)=cH(1/c)−cH((1−t)/c).

These are analytic quantities used to describe local encoding and reporting costs. The source's twelve charges also depend on a representation-scale parameter eee, a band shift bbb, a box saving SSS, an initial threshold θ\thetaθ, a terminal fraction τ\tauτ, and a uniform entropy coefficient ρ\rhoρ.

Formalization targets

The campaign goal is the exact proposition

EndStatement.ThreeSum.SolvedInTime  (1.960614).\texttt{EndStatement.ThreeSum.SolvedInTime}\;(1.960614).EndStatement.ThreeSum.SolvedInTime(1.960614).

For 1.960614=980307/5000001.960614=980307/5000001.960614=980307/500000, its time requirement is

∃K∈N  ∀n≥2,T(n)500000≤Kn980307.\exists K\in\mathbb N\;\forall n\ge2,\quad T(n)^{500000}\le K n^{980307}.∃K∈N∀n≥2,T(n)500000≤Kn980307.

Correctness and termination still cover all input sizes. The integer-size exponent includes κ=0\kappa=0κ=0, and the same program must be correct for every word width above its specified threshold.

The primary certificate uses c=232c=232c=232 and the rational quartic

P(t)=7.85174904t+0.72447023t2+0.24044466t3+0.12022459t4.P(t)=7.85174904t+0.72447023t^2+0.24044466t^3+0.12022459t^4.P(t)=7.85174904t+0.72447023t2+0.24044466t3+0.12022459t4.

Supporting targets establish its validity on [0,1/2][0,1/2][0,1/2], its derivative and tangent bounds, and the uniform inequality H(t)<0.27905030(1+P(t))H(t)<0.27905030(1+P(t))H(t)<0.27905030(1+P(t)). The fixed parameters are e=0.203889696e=0.203889696e=0.203889696, b=−0.062656016b=-0.062656016b=−0.062656016, S=0.193175570S=0.193175570S=0.193175570, θ=0.040409899\theta=0.040409899θ=0.040409899, and τ=13/20\tau=13/20τ=13/20. The parameter theorem bounds every one of the twelve costs strictly below 1.96061388<1.9606141.96061388<1.9606141.96061388<1.960614. Integer scheduling and general local-cover identities are separate milestones. A further independent milestone proves R>1.96061349120R>1.96061349120R>1.96061349120 for the explicitly defined relaxed fixed-expansion cost model. That model bound concerns its admissible parameter tuples; it is not a lower bound for 3SUM itself.

What these results establish

The supplied certificate demonstrates that a quartic switch already suffices for the proposed 1.960614 cost target. It also separates a uniform analytic inequality from a numerical experiment: the statement quantifies over the full closed interval. Its finite scheduling statement uses integer comparisons and accounts for the rounding error at every admissible level, rather than requiring executable logarithms.

Formalization adds checked statements about these continuous and finite objects and a reusable dependency structure. The full campaign goal requires an additional result: the actual modified program must obey the certified costs, and its operational correctness must reach the foundation specification. The upstream parameter certificate likewise distinguishes its exponent ledger from its separate implementation obligations. An exact rational checker passing is evidence for the numerical certificates; it does not establish the full new machine theorem.

Why the program extension remains substantive

Changing expansion from 256 to 232 changes the local universe and the vertex, output and degree counts together. The representation scale, band denominator and polynomial switch also change. Consequently, the original fixed-parameter runtime result cannot simply be invoked after replacing a few numerical constants. The new local construction, global encoding, sparse queries, reporting and arithmetic operations must share consistent parameters and a complete cost account.

The small gap between the principal cost cap and the target must absorb the proved subpolynomial overheads, including exceptional-size handling and the conversion to word instructions. The machine bridge has now been proved in this mission. Its accepted dependency tree connects the modified program and its runtime analysis to the unchanged campaign specification. Its statement has no assumed solver callback or unproved runtime contract. The fixed-model lower-bound milestone is an independent analytic root, and is not used to justify a faster algorithm.

Formalization scope and source boundaries

The campaign pins Lean 4.33.1 and Mathlib revision 0df444a360eaa60ab8c11dca51a86af692955474. Entropy and its derived quantities are real-valued; the final exponent is rational; finite subset counts and cleared switch comparisons use natural numbers. Theorems explicitly restrict logarithmic arguments and denominators to their intended domains. Natural-number subtraction in the switch is used only below the final level, with positive even microblock size.

The local-cover milestone covers every natural c≥3c\ge3c≥3 and positive half-block size. The cost-model milestone quantifies over all tuples satisfying its written assumptions and includes a whole-interval entropy constraint. Neither finite testing nor a floating-point optimizer substitutes for these universal statements. The initial package includes no Lean proof of the new machine goal, and the mission does not represent it as an established record.

The underlying algorithms and formal infrastructure retain their attribution to Alman and Vassilevska Williams; Björklund, Husfeldt, Kaski and Koivisto; Swapnil Jain; c0ldzy17; Anthropic; Jihoon Hyun; and Lean/Mathlib, as documented by the upstream source. The certificate development and this formalization use AI assistance. Scientific priority and global optimality are outside the mission's claims.

Selected references

  • Author-supplied 3SUM: Positive-Coefficient Entropy Switch, 1.960614 Cost Certificate and Fixed-Model Lower Bound, certificate package v3, 14 October 2026. Exact source: 3sum_certificate_1p960614_v3/3sum_certificate_v3/PROOF.md, Sections 0–8; certificates/c232_quartic.json; verify_certificate.py. This is a supplied research artifact, not a claimed public publication.
  • c0ldzy17, Deterministic 3SUM in O(n^1.962): method and proof, source documentation, 2026. Pinned METHOD.md, Parameters.lean, and LocalCoverBounds.lean.
  • Josh Alman and Virginia Vassilevska Williams, Truly Subquadratic 3SUM and Truly Subcubic APSP via Triangles in Sparse Lopsided Graphs, arXiv preprint, 2026. arXiv:2610.06783v1.
  • Anthropic, Truly Subquadratic 3SUM and Truly Subcubic APSP, Lean formalization, 2026. Pinned formal-math source. The campaign imports its shared word-RAM statement.

Source code and complete Lean formalization: GitHub repository. Fixed source revision: 8c25c0325de5.

84 thms2 active usersReviewed
🏆Completed
Information Theory·Captain: chaikon

Essential Coding Theory I: Generalized Hamming BoundTextbook

Motivation

An error-correcting code stores or transmits messages as words over a finite alphabet, adding redundancy so that the original word can be recovered after some symbols are corrupted. The basic question of coding theory is the trade-off between the amount of redundancy and the number of errors that can be corrected. The Hamming bound (also called the sphere-packing bound) is the first and most widely taught upper limit on this trade-off. It is due to R. W. Hamming (1950), who introduced the Hamming distance and the Hamming codes in the same paper (Hamming 1950).

This mission formalizes the bound as it is stated in Theorem 1.7.2 of the textbook Essential Coding Theory by V. Guruswami, A. Rudra and M. Sudan (web draft, §1.7, p. 21). It is the first mission of a series following that book. The result is used throughout the book to define perfect codes (Definition 1.7.3) and as the baseline against which later bounds (Singleton, Plotkin, Elias–Bassalygo, the linear-programming bound) are compared.

Setting

Let α\alphaα be a finite alphabet with q=∣α∣q = |\alpha|q=∣α∣ symbols, and let n≥0n \ge 0n≥0. A word of length nnn is a function {0,…,n−1}→α\{0, \dots, n-1\} \to \alpha{0,…,n−1}→α, and αn\alpha^nαn is the set of all qnq^nqn words. A code of block length nnn is a subset C⊆αnC \subseteq \alpha^nC⊆αn; its elements are codewords.

The Hamming distance Δ(u,v)\Delta(u, v)Δ(u,v) between two words is the number of coordinates in which they differ. The dimension of CCC is the real number k=log⁡q∣C∣k = \log_q |C|k=logq​∣C∣ (Definition 1.2.3 of the book). The book calls CCC an (n,k,d)q(n, k, d)_q(n,k,d)q​ code when its minimum distance min⁡c1≠c2∈CΔ(c1,c2)\min_{c_1 \ne c_2 \in C} \Delta(c_1, c_2)minc1​=c2​∈C​Δ(c1​,c2​) equals ddd (Definitions 1.4.1 and 1.7.1).

The Hamming ball of radius eee around a word xxx is B(x,e)={y∈αn:Δ(x,y)≤e}B(x, e) = \{ y \in \alpha^n : \Delta(x, y) \le e \}B(x,e)={y∈αn:Δ(x,y)≤e} (Definition 1.6.1).

Formalization targets

Milestone: volume of a Hamming ball

For every word x∈αnx \in \alpha^nx∈αn and every radius e≥0e \ge 0e≥0,

∣B(x,e)∣=∑i=0e(ni)(q−1)i.|B(x, e)| = \sum_{i=0}^{e} \binom{n}{i} (q-1)^i.∣B(x,e)∣=i=0∑e​(in​)(q−1)i.

This is equation (1.18) in the book.

Goal: the generalized Hamming bound

If any two distinct codewords of C⊆αnC \subseteq \alpha^nC⊆αn are at Hamming distance at least ddd, then

log⁡q∣C∣  ≤  n−log⁡q(∑i=0⌊(d−1)/2⌋(ni)(q−1)i).\log_q |C| \;\le\; n - \log_q \Bigl( \sum_{i=0}^{\lfloor (d-1)/2 \rfloor} \binom{n}{i} (q-1)^i \Bigr).logq​∣C∣≤n−logq​(i=0∑⌊(d−1)/2⌋​(in​)(q−1)i).

The hypothesis is "minimum distance at least ddd" rather than "exactly ddd". Every (n,k,d)q(n, k, d)_q(n,k,d)q​ code satisfies it, so the goal contains the book's Theorem 1.7.2 as stated.

Significance

The result itself. The bound says that a code correcting ⌊(d−1)/2⌋\lfloor (d-1)/2 \rfloor⌊(d−1)/2⌋ errors can have at most qn/Volq(⌊(d−1)/2⌋,n)q^n / \mathrm{Vol}_q(\lfloor (d-1)/2 \rfloor, n)qn/Volq​(⌊(d−1)/2⌋,n) codewords, so its rate is at most 1−log⁡qVolq/n1 - \log_q \mathrm{Vol}_q / n1−logq​Volq​/n. Codes meeting it with equality are the perfect codes: the Hamming codes, the binary repetition codes of odd length, and the two Golay codes. The bound's asymptotic form, R≤1−Hq(δ/2)R \le 1 - H_q(\delta/2)R≤1−Hq​(δ/2), is one of the standard upper bounds on the rate–distance trade-off and is the reference point for the stronger bounds later in the book.

Formalizing it. The theorem has been proved since 1950 and appears in every coding theory text. Mathlib has the Hamming distance (hammingDist) but no notion of a code's distance, no Hamming-ball cardinality formula, and no sphere-packing bound. This mission supplies those, with the dimension stated through the real logarithm exactly as in the book, so later missions in the series can cite it.

Difficulty

The mathematics is elementary. The work lies in the counting step: showing that exactly (ni)(q−1)i\binom{n}{i}(q-1)^i(in​)(q−1)i words lie at distance exactly iii from a fixed word. This identity has no ready-made counterpart in Mathlib, and the alphabet carries no algebraic structure, so translating every ball to a ball around a fixed word is not available. The final step, taking log⁡q\log_qlogq​ of both sides, must also handle the conventions of the real logarithm on degenerate inputs.

Formalization scope

  • The alphabet is an arbitrary type α with [Fintype α] [DecidableEq α], and qqq is Fintype.card α. Words are functions Fin n → α, and Δ\DeltaΔ is Mathlib's hammingDist.
  • A code is a Finset (Fin n → α). The distance hypothesis quantifies over pairs of distinct codewords.
  • The dimension is Real.logb q |C|, a real number. The ball volume is the natural-number sum ∑i=0⌊(d−1)/2⌋(ni)(q−1)i\sum_{i=0}^{\lfloor (d-1)/2 \rfloor} \binom{n}{i}(q-1)^i∑i=0⌊(d−1)/2⌋​(in​)(q−1)i, cast to R\mathbb{R}R, where ⌊(d−1)/2⌋\lfloor (d-1)/2 \rfloor⌊(d−1)/2⌋ is computed in N\mathbb{N}N as (d - 1) / 2.
  • No hypotheses q≥2q \ge 2q≥2 or C≠∅C \ne \emptysetC=∅ are imposed. The book assumes them tacitly. Under Lean's conventions log⁡qx=0\log_q x = 0logq​x=0 when q≤1q \le 1q≤1 or x=0x = 0x=0, and the statement remains true in those cases. It is therefore at least as strong as the book's version, and these degenerate cases are not a trivialization: for q≥2q \ge 2q≥2 and C≠∅C \ne \emptysetC=∅ the goal is equivalent to the sphere-packing inequality ∣C∣⋅Vol≤qn|C| \cdot \mathrm{Vol} \le q^n∣C∣⋅Vol≤qn.
  • No custom definitions are needed. The ball-volume formula is reusable for any later bound that counts Hamming balls (Gilbert–Varshamov, Elias–Bassalygo). Contributions of alternative counting proofs are welcome.

Selected references

  • V. Guruswami, A. Rudra and M. Sudan, Essential Coding Theory, draft textbook, §1.6–1.7, Theorem 1.7.2, p. 21. https://cse.buffalo.edu/faculty/atri/courses/coding-theory/book/web-coding-book.pdf
  • R. W. Hamming, Error detecting and error correcting codes, Bell System Technical Journal 29(2), 147–160, 1950. https://doi.org/10.1002/j.1538-7305.1950.tb00463.x
  • The Mathlib Community, hammingDist, Mathlib.InformationTheory.Hamming. https://leanprover-community.github.io/mathlib4_docs/Mathlib/InformationTheory/Hamming.html
2 thms2 active usersReviewed
🏆Completed
Theoretical Computer Science·Captain: Yuxuan Xu

Deterministic Integer 3SUM in O(n^1.999017): Sharp Binomial-Tail AnalysisResearch Paper

Why sharpen the integer 3SUM exponent?

3SUM asks whether an input sequence contains three entries at different positions whose sum is zero. The classical quadratic running time makes it a benchmark for fine-grained algorithm design: even a small fixed improvement in the exponent changes its asymptotic complexity. The setting here is polynomially bounded integers and a deterministic word RAM; the input-size exponent, rather than an experimental speedup, is the quantity being studied.

Alman and Vassilevska Williams (2026) give a deterministic truly subquadratic integer 3SUM algorithm through thin matrix products and reductions involving Exact Triangle. Xu's manuscript, Deterministic Integer 3SUM in O(n^1.999017) Time via a Sharper Binomial-Tail Analysis, dated 13 October 2026, refines the analysis of that framework. Its claimed bound is the formalization target here. The Hexagon link identifies the author's submission page, which may require authentication; it is not presented as a public article identifier or evidence of peer review.

Inputs, machine model, and analytic quantities

An input of size nnn is a function x:{0,…,n−1}→Zx:\{0,\ldots,n-1\}\to\mathbb Zx:{0,…,n−1}→Z. For a fixed natural number κ\kappaκ, every entry obeys ∣xi∣≤nκ|x_i|\le n^\kappa∣xi​∣≤nκ. Acceptance means that there are pairwise distinct indices i,j,ki,j,ki,j,k with xi+xj+xk=0x_i+x_j+x_k=0xi​+xj​+xk​=0. Repeated values at different positions are allowed. This is a decision problem, so the algorithm need not return a witness triple.

The word-RAM specification is the campaign's existing EndStatement definition, imported from its foundation mission. For each fixed κ\kappaκ, one finite deterministic program, one word-width threshold constant, and one time-bound function must work for all input sizes and all word widths above that threshold. Words have finite signed interpretations; arithmetic, indirect memory access, branching, acceptance, and rejection are defined operationally. The result concerns this exact model.

For the supporting analysis, let H(u)=−ulog⁡u−(1−u)log⁡(1−u)H(u)=-u\log u-(1-u)\log(1-u)H(u)=−ulogu−(1−u)log(1−u), with endpoint values H(0)=H(1)=0H(0)=H(1)=0H(0)=H(1)=0. All logarithms are natural. For an integer c>10c>10c>10, a positive integer mmm, and 0<θ<0.90<\theta<0.90<θ<0.9, set L=cmL=cmL=cm, D=4mD=4^mD=4m, t=⌈θm⌉t=\lceil\theta m\rceilt=⌈θm⌉, and

βd=(Lm−d)9L−m+d,M=β0.\beta_d=\binom L{m-d}9^{L-m+d},\qquad M=\beta_0.βd​=(m−dL​)9L−m+d,M=β0​.

The sharp preprocessing saving, query exponent, and encoding exponent are

γ∗(c,θ)=cH(1/c)−cH((1−θ)/c)−θlog⁡9log⁡4,\gamma_*(c,\theta)=\frac{cH(1/c)-cH((1-\theta)/c)-\theta\log9}{\log4},γ∗​(c,θ)=log4cH(1/c)−cH((1−θ)/c)−θlog9​, q(θ)=H(θ)+θlog⁡9log⁡4,A(c)=clog⁡10−cH(1/c)/2−(c−1)log⁡3log⁡4.q(\theta)=\frac{H(\theta)+\theta\log9}{\log4},\qquad A(c)=\frac{c\log10-cH(1/c)/2-(c-1)\log3}{\log4}.q(θ)=log4H(θ)+θlog9​,A(c)=log4clog10−cH(1/c)/2−(c−1)log3​.

These auxiliary definitions describe costs; they neither redefine 3SUM nor assume the existence of a fast algorithm.

Formalization targets

The campaign goal is the exact fixed bound

EndStatement.ThreeSum.SolvedInTime  1.999017.\texttt{EndStatement.ThreeSum.SolvedInTime}\;1.999017.EndStatement.ThreeSum.SolvedInTime1.999017.

The first supporting milestone is the two-sided binomial entropy estimate. The second is the manuscript's normalized tail bound

1M∑d=tmβd≤cm+11−9/(c−1)D−γ∗(c,θ).\frac1M\sum_{d=t}^m\beta_d\le\frac{cm+1}{1-9/(c-1)}D^{-\gamma_*(c,\theta)}.M1​d=t∑m​βd​≤1−9/(c−1)cm+1​D−γ∗​(c,θ).

A further milestone establishes strict improvement over the earlier preprocessing saving γ(c,θ)=θlog⁡((c−1)/9)/log⁡4\gamma(c,\theta)=\theta\log((c-1)/9)/\log4γ(c,θ)=θlog((c−1)/9)/log4 at fixed parameters. The algorithmic milestone combines the manuscript's Propositions 4.1 and 5.1: admissible sharp-tail parameters must yield a program satisfying the shared 3SUM specification, not merely a numerical inequality.

The final supporting milestone certifies the exact choice

c=21,θ=0.11202,a=0.05665,s=0.06945.c=21,\qquad\theta=0.11202,\qquad a=0.05665,\qquad s=0.06945.c=21,θ=0.11202,a=0.05665,s=0.06945.

It includes the preprocessing, query, encoding, prime-selection, and instance-construction inequalities, together with

2−as/4=1.999016414375,1.999017−(2−as/4)=0.000000585625>0.2-as/4=1.999016414375,\qquad1.999017-(2-as/4)=0.000000585625>0.2−as/4=1.999016414375,1.999017−(2−as/4)=0.000000585625>0.

Every terminating decimal denotes an exact rational number. The target is the fixed big-OOO exponent 1.999017; the smaller cutoff is accompanied by residual asymptotic losses in the manuscript and is not asserted as the same fixed bound.

What completing this mission would establish

A complete proof would certify the manuscript's improved integer 3SUM bound in the campaign's common machine model. The gain is an analysis refinement of the inherited algorithm and reductions. It does not claim a new algorithmic framework, optimal parameters, global optimality, or priority over all other work.

The manuscript also gives an auxiliary O(n1.999073)O(n^{1.999073})O(n1.999073) bound by retuning the previous estimate. That weaker result is background, not this mission's campaign value. The existing upstream formalization is prior work; it does not automatically certify either refinement. The supplied rational interval verifier checks the parameter inequalities, but those checks alone do not prove a word-RAM running-time theorem. The main goal and the newly drafted supporting statements remain open Lean proof obligations.

Where the difficulty lies

A better tail estimate is insufficient unless every other cost remains within the final budget. Shared encoding work, all queries, deterministic prime selection, witness scanning, and instance construction must be included. The word-size and address bounds must remain valid after changing parameters. Ceiling operations, restricting DDD to powers of four, and bounded-size exceptions must also be handled.

The central formalization burden is the connection between the analytic estimates and an actual deterministic program. Assuming an oracle with the desired running time would remove that burden and would not establish the stated target. Likewise, checking a finite range of sizes or a floating-point optimizer's output cannot establish the universal algorithmic statement.

Formalization scope and conventions

The environment is Lean 4.33.1 with Mathlib revision 0df444a360eaa60ab8c11dca51a86af692955474. The foundation reference and the campaign goal's import are retained. Supporting definitions use real logarithms, real powers, natural binomial coefficients, finite sums, and a natural ceiling. Domain hypotheses exclude invalid entropy-tail regimes; endpoint binomial cases are included explicitly.

For r=1999017/1000000r=1999017/1000000r=1999017/1000000, the shared time predicate means that there exists K∈NK\in\mathbb NK∈N such that T(n)1000000≤Kn1999017T(n)^{1000000}\le K n^{1999017}T(n)1000000≤Kn1999017 for every n≥2n\ge2n≥2. Correctness and termination still cover every input size, including n<3n<3n<3. The manuscript's κ≥1\kappa\ge1κ≥1 convention must be extended to the campaign's κ=0\kappa=0κ=0 case, and its word-size analysis must justify uniform correctness at every sufficiently large word width. Both are explicit interface obligations.

Reusable contributions include entropy and tail estimates, rigorous logarithm enclosures, and parameterized cost-to-program results. Existing compatible library results should be reused where available. The manuscript and metadata were prepared with AI assistance; this mission makes no claim of completed human review or completed formal verification of the improved bound.

Selected references

  • Yuxuan Xu, Deterministic Integer 3SUM in O(n^1.999017) Time via a Sharper Binomial-Tail Analysis, manuscript, 13 October 2026. Author-supplied Hexagon submission page. Target: Theorem 1.1; supporting analysis: Sections 3–6.
  • Josh Alman and Virginia Vassilevska Williams, Truly Subquadratic 3SUM and Truly Subcubic APSP via Triangles in Sparse Lopsided Graphs, arXiv preprint, 2026. arXiv:2610.06783v1. Inherited ingredients: Section 4, Lemma 29, Theorem 30, Corollary 31, and the deterministic reductions cited in Xu's Section 2.
  • Anthropic, Truly Subquadratic 3SUM and Truly Subcubic APSP: a Lean 4 formalization, source repository, 2026. Pinned source at commit e1a4e6508154ea59f030480661590a9fe3018011. Its EndStatement.lean supplies the shared specification, not a proof of this mission's new exponent.
8 thms2 active usersReviewed
🏆Completed
Theoretical Computer Science·Captain: wurtle

3SUM in O(n^1.962) Time on a Word RAMResearch Paper

This entry records a deterministic O(n1.962)O(n^{1.962})O(n1.962) bound for 3SUM on polynomially bounded integers, using the campaign's shared word-RAM model with O(log⁡n)O(\log n)O(logn)-bit words. Three distinct input positions must sum to zero; repeated values are allowed. The exact exponent is 981/500981/500981/500.

References

  • Source repository: complete Lean proof and verification record
49 thms2 active usersReviewed
🏆Completed
Algorithmic Game TheoryEconomics·Captain: Lucas

Stable Matching Problem (Wikipedia) I: The Gale–Shapley AlgorithmTextbook

Motivation

The stable matching problem asks for a pairing between two equally sized groups in which no two participants would both rather be matched with each other than with their assigned partners. It was posed and solved by Gale and Shapley in 1962 in the context of college admissions and marriage (Gale–Shapley 1962). Their deferred acceptance algorithm underlies centralized clearinghouses such as the assignment of graduating medical students to hospital residencies, and the theory of stable allocations was recognized by the 2012 Nobel Memorial Prize in Economic Sciences awarded to Shapley and Roth.

Timeline of the results targeted here:

  • 1962 — Gale and Shapley: a stable matching always exists, and the men-proposing deferred acceptance algorithm finds one; its outcome is the best stable matching for every man.
  • 1971 — McVitie and Wilson: the men-optimal stable matching is simultaneously the worst stable matching for every woman.
  • 1981–1982 — Dubins–Freedman and Roth: the men-proposing algorithm is strategy-proof (indeed group strategy-proof) for men, while women can gain by misreporting.
  • 1987 — Gusfield (after Knuth and Conway): the stable matchings of an instance form a finite distributive lattice.

Setting

Fix n≥0n \ge 0n≥0. There are nnn men and nnn women, both indexed by Fin n\mathrm{Fin}\,nFinn. Each man mmm ranks all women strictly: manRankm(w)∈Fin n\mathrm{manRank}_m(w) \in \mathrm{Fin}\,nmanRankm​(w)∈Finn is the position of www in mmm's list, position 000 being most preferred; symmetrically each woman www has womanRankw\mathrm{womanRank}_wwomanRankw​. Both are bijections. A matching is a bijection μ\muμ from men to women. A pair (m,w)(m, w)(m,w) blocks μ\muμ if manRankm(w)<manRankm(μ(m))\mathrm{manRank}_m(w) < \mathrm{manRank}_m(\mu(m))manRankm​(w)<manRankm​(μ(m)) and womanRankw(m)<womanRankw(μ−1(w))\mathrm{womanRank}_w(m) < \mathrm{womanRank}_w(\mu^{-1}(w))womanRankw​(m)<womanRankw​(μ−1(w)); μ\muμ is stable if no pair blocks it.

The Gale–Shapley algorithm runs in rounds. A state records, for every man, how many women he has proposed to so far, and for every woman, her current provisional fiancé (if any). In each round every unengaged man proposes to the most preferred woman he has not yet proposed to, and every woman keeps the most preferred among her current fiancé and this round's suitors, rejecting the others. The Gale–Shapley outcome is the matching read off from the state after n2n^2n2 rounds.

Formalization targets

Goal: the algorithm produces a stable matching

∀I, ∃μ, μ is the Gale–Shapley outcome of I ∧ μ is stable for I.\forall I,\ \exists \mu,\ \mu \text{ is the Gale–Shapley outcome of } I \ \wedge\ \mu \text{ is stable for } I.∀I, ∃μ, μ is the Gale–Shapley outcome of I ∧ μ is stable for I.

This contains termination (all men are engaged after n2n^2n2 rounds), the fact that the final engagements form a bijection, and stability.

Milestones

  1. Termination: after n2n^2n2 rounds every man is engaged.
  2. Existence: every instance has a stable matching.
  3. Man-optimality: for the Gale–Shapley outcome μ0\mu_0μ0​ and every stable μ\muμ, manRankm(μ0(m))≤manRankm(μ(m))\mathrm{manRank}_m(\mu_0(m)) \le \mathrm{manRank}_m(\mu(m))manRankm​(μ0​(m))≤manRankm​(μ(m)) for every man mmm.
  4. Woman-pessimality: womanRankw(μ−1(w))≤womanRankw(μ0−1(w))\mathrm{womanRank}_w(\mu^{-1}(w)) \le \mathrm{womanRank}_w(\mu_0^{-1}(w))womanRankw​(μ−1(w))≤womanRankw​(μ0−1​(w)) for every woman www.
  5. Truthfulness for men: no single man obtains a strictly better partner (by his true list) by misreporting.
  6. Group strategy-proofness for men: no nonempty coalition of men can misreport so that all its members are strictly better off.
  7. Weak coalitional manipulation: some coalition can misreport so that some members are strictly better off and the rest keep their partners.
  8. Women can manipulate: some woman can obtain a strictly better partner by misreporting.
  9. The 3×3 example: the instance A:YXZA: YXZA:YXZ, B:ZYXB: ZYXB:ZYX, C:XZYC: XZYC:XZY, X:BACX: BACX:BAC, Y:CBAY: CBAY:CBA, Z:ACBZ: ACBZ:ACB has exactly the three stable matchings (AY,BZ,CX)(AY,BZ,CX)(AY,BZ,CX), (AX,BY,CZ)(AX,BY,CZ)(AX,BY,CZ), (AZ,BX,CY)(AZ,BX,CY)(AZ,BX,CY).
  10. Lattice structure: the stable matchings of any instance form a distributive lattice under "every man weakly prefers".

Significance

Existence of stable matchings is the foundation of two-sided matching theory, and the algorithmic proof is constructive. Man-optimality and woman-pessimality show that the proposing side is systematically favored, which is the reason the choice of proposing side matters in practice, as in the 1995 redesign of the residency match. The incentive results (milestones 5–8) explain why participants on the proposing side can safely report their true preferences while the receiving side cannot, and the lattice structure underlies efficient algorithms for enumerating and optimizing over stable matchings.

All results in this mission are classical and proved in the literature. The work here is formalization: a precise model of the round-based algorithm, machine-checked proofs of its correctness and optimality, and the incentive properties, which are usually argued informally.

Difficulty

The goal requires reasoning about an iterated state transformer rather than an abstract existence statement. The step that is usually waved away is the invariant that an unengaged man still has a woman left to propose to, which is what makes the n2n^2n2 bound and the bijectivity of the final engagements work. Man-optimality needs the invariant that no man is ever rejected by a woman who is his partner in some stable matching, which is not visible from the output alone. Group strategy-proofness is substantially harder than individual truthfulness: the known proofs go through a "blocking lemma" about matchings that are weakly preferred by all men.

Formalization scope

Men and women are Fin n; preferences are strict total orders given as rank bijections Fin n ≃ Fin n; matchings are bijections Fin n ≃ Fin n. The algorithm is the round-based version described in the source: all unengaged men propose simultaneously in each round. The algorithm's output is fixed as the state after exactly n2n^2n2 rounds; termination within that bound is itself a milestone, so a wrong bound would make the goal false rather than vacuous. The case n=0n = 0n=0 is included and is trivial. Incentive results compare Gale–Shapley outcomes of two instances that differ only in the lists of the deviating participants, with gains measured by the deviators' true lists.

The definitions are reusable for further missions on the same article (rural hospitals theorem, hospitals/residents problem). Contributions welcome: invariants of the algorithm, the blocking lemma, and alternative proofs of existence.

Selected references

  • D. Gale, L. S. Shapley, College Admissions and the Stability of Marriage, American Mathematical Monthly 69(1), 1962, 9–14. https://doi.org/10.2307/2312726
  • D. G. McVitie, L. B. Wilson, The stable marriage problem, Communications of the ACM 14(7), 1971, 486–490. https://doi.org/10.1145/362619.362631
  • L. E. Dubins, D. A. Freedman, Machiavelli and the Gale–Shapley Algorithm, American Mathematical Monthly 88(7), 1981, 485–494. https://doi.org/10.2307/2321753
  • A. E. Roth, The Economics of Matching: Stability and Incentives, Mathematics of Operations Research 7(4), 1982, 617–628. https://doi.org/10.1287/moor.7.4.617
  • D. Gusfield, Three fast algorithms for four problems in stable marriage, SIAM Journal on Computing 16(1), 1987, 111–128. https://doi.org/10.1137/0216010
  • Wikipedia, Stable matching problem. https://en.wikipedia.org/wiki/Stable_matching_problem
12 thms2 active usersReviewed
🏆Completed
Group TheoryNumber Theory·Captain: Lucas

Landau's Function: Growth of the Maximal Order of a PermutationResearch Paper

Motivation

How large can the order of a single permutation of nnn objects be? The answer, Landau's function g(n)g(n)g(n), is a basic invariant of the symmetric group SnS_nSn​: it is the maximum number of times a permutation can be iterated before it returns to the identity. Because the order of a permutation is the least common multiple of its cycle lengths, g(n)g(n)g(n) is also the largest least common multiple of the parts of a partition of nnn, which ties the question to the distribution of prime numbers.

Timeline.

  • 1902–1903. E. Landau proved that ln⁡g(n)∼nln⁡n\ln g(n)\sim\sqrt{n\ln n}lng(n)∼nlnn​.
  • 1968. J.-L. Nicolas showed that ggg is constant on arbitrarily long runs of consecutive integers.
  • 1984. J.-P. Massias proved the explicit bound ln⁡g(n)≤1.05313836nln⁡n\ln g(n)\le 1.05313836\sqrt{n\ln n}lng(n)≤1.05313836nlnn​ for all n≥1n\ge 1n≥1.
  • 1988. Massias, Nicolas and Robin obtained a three-term asymptotic expansion of ln⁡g(n)\ln g(n)lng(n), its expression through the inverse logarithmic integral, and the equivalence of a one-sided inequality with the Riemann hypothesis.
  • 2002. K. Ford's bounds for the Riemann zeta function give the error term xexp⁡(−c(ln⁡x)3/5(ln⁡ln⁡x)−1/5)x\exp(-c(\ln x)^{3/5}(\ln\ln x)^{-1/5})xexp(−c(lnx)3/5(lnlnx)−1/5) in the prime number theorem, which can be fed into the 1988 result.

Setting

For a natural number nnn, let SnS_nSn​ be the group of all permutations of an nnn-element set. The order of σ∈Sn\sigma\in S_nσ∈Sn​ is the least k≥1k\ge 1k≥1 with σk=id\sigma^k=\mathrm{id}σk=id. Landau's function is

g(n)=max⁡σ∈Snord⁡(σ).g(n)=\max_{\sigma\in S_n}\operatorname{ord}(\sigma).g(n)=σ∈Sn​max​ord(σ).

A partition of nnn is a multiset of positive integers summing to nnn. The prime counting function π(x)\pi(x)π(x) is the number of primes p≤xp\le xp≤x, the (offset) logarithmic integral is Li(x)=∫2xdtln⁡t\mathrm{Li}(x)=\int_2^x\frac{dt}{\ln t}Li(x)=∫2x​lntdt​, and Li−1\mathrm{Li}^{-1}Li−1 is its inverse on [2,∞)[2,\infty)[2,∞). Throughout, ln⁡\lnln is the natural logarithm.

Formalization targets

Goal: Landau's theorem (1902)

lim⁡n→∞ln⁡g(n)nln⁡n=1.\lim_{n\to\infty}\frac{\ln g(n)}{\sqrt{n\ln n}}=1.n→∞lim​nlnn​lng(n)​=1.

Supporting milestones

  1. Partition characterization: g(n)=max⁡{lcm⁡(λ):λ⊢n}g(n)=\max\{\operatorname{lcm}(\lambda):\lambda\vdash n\}g(n)=max{lcm(λ):λ⊢n}.
  2. Initial values g(0),…,g(8)=1,1,2,3,4,6,6,12,15g(0),\dots,g(8)=1,1,2,3,4,6,6,12,15g(0),…,g(8)=1,1,2,3,4,6,6,12,15.
  3. Nicolas (1968): for every mmm there is nnn with g(n)=g(n+1)=⋯=g(n+m)g(n)=g(n+1)=\dots=g(n+m)g(n)=g(n+1)=⋯=g(n+m).
  4. Massias–Nicolas–Robin (1988): ln⁡g(n)=nln⁡n(1+ln⁡ln⁡n−12ln⁡n−(ln⁡ln⁡n)2−6ln⁡ln⁡n+98(ln⁡n)2+O((ln⁡ln⁡nln⁡n)3))\ln g(n)=\sqrt{n\ln n}\Bigl(1+\frac{\ln\ln n-1}{2\ln n}-\frac{(\ln\ln n)^2-6\ln\ln n+9}{8(\ln n)^2}+O\bigl((\tfrac{\ln\ln n}{\ln n})^3\bigr)\Bigr)lng(n)=nlnn​(1+2lnnlnlnn−1​−8(lnn)2(lnlnn)2−6lnlnn+9​+O((lnnlnlnn​)3)).
  5. Ford (2002): π(x)−Li(x)=O(xexp⁡(−c(ln⁡x)3/5(ln⁡ln⁡x)−1/5))\pi(x)-\mathrm{Li}(x)=O\bigl(x\exp(-c(\ln x)^{3/5}(\ln\ln x)^{-1/5})\bigr)π(x)−Li(x)=O(xexp(−c(lnx)3/5(lnlnx)−1/5)) for some c>0c>0c>0.
  6. Massias–Nicolas–Robin (1988): if π(x)−Li(x)=O(R(x))\pi(x)-\mathrm{Li}(x)=O(R(x))π(x)−Li(x)=O(R(x)) with R(x)=xexp⁡(−c(ln⁡x)3/5(ln⁡ln⁡x)−1/5)R(x)=x\exp(-c(\ln x)^{3/5}(\ln\ln x)^{-1/5})R(x)=xexp(−c(lnx)3/5(lnlnx)−1/5), then ln⁡g(n)=Li−1(n)+O(R(nln⁡n)ln⁡n)\ln g(n)=\sqrt{\mathrm{Li}^{-1}(n)}+O\bigl(R(\sqrt{n\ln n})\ln n\bigr)lng(n)=Li−1(n)​+O(R(nlnn​)lnn).
  7. Riemann hypothesis criterion: ln⁡g(n)<Li−1(n)\ln g(n)<\sqrt{\mathrm{Li}^{-1}(n)}lng(n)<Li−1(n)​ for all large nnn   ⟺  \iff⟺ RH.
  8. g(n)≤en/eg(n)\le e^{n/e}g(n)≤en/e for all nnn, with equality only at n=0n=0n=0.
  9. Massias (1984): g(n)≤exp⁡(1.05314nln⁡n)g(n)\le\exp(1.05314\sqrt{n\ln n})g(n)≤exp(1.05314nlnn​).

Significance

Landau's theorem pins down the growth rate of the largest element order in SnS_nSn​, which bounds, for example, how long any fixed shuffle of nnn cards can be repeated before the deck returns to its starting order. The refined statements quantify the second-order behaviour and connect it to the error term in the prime number theorem; the criterion in milestone 7 is one of the known elementary reformulations of the Riemann hypothesis, and is included as a faithful statement of the source, not as an expected deliverable. The other results are proved in the literature; formalizing them requires building the partition/lcm description of element orders in SnS_nSn​ and quantitative prime-counting estimates on top of Mathlib.

Difficulty

The order of a permutation is easy to compute from its cycle type, but maximizing an lcm over all partitions is a non-linear optimization: the optimum uses prime powers whose sizes are dictated by the distribution of primes, so even the leading-order asymptotics need Chebyshev-type or prime number theorem estimates, and the refined milestones need the prime number theorem with an explicit error term.

Formalization scope

  • SnS_nSn​ is Equiv.Perm (Fin n); g(n)g(n)g(n) is the supremum of orderOf σ over the finite group (a natural number; g(0)=1g(0)=1g(0)=1).
  • Partitions are Mathlib's Nat.Partition n; the lcm of a partition is the Multiset.lcm of its parts.
  • Real logarithms use Real.log (which returns 000 at 000 and is defined for negative input); all asymptotic statements are about n→∞n\to\inftyn→∞, so the values at small nnn are irrelevant there. The bounds in milestones 8–9 are stated for every nnn and are true at n=0,1n=0,1n=0,1.
  • π(x)\pi(x)π(x) for real xxx is Nat.primeCounting ⌊x⌋₊; Li\mathrm{Li}Li is the offset logarithmic integral ∫2xdt/ln⁡t\int_2^x dt/\ln t∫2x​dt/lnt; Li−1(y)\mathrm{Li}^{-1}(y)Li−1(y) is the least x≥2x\ge2x≥2 with Li(x)≥y\mathrm{Li}(x)\ge yLi(x)≥y, which for y≥0y\ge0y≥0 is the inverse function value. Replacing Li\mathrm{Li}Li by the full integral li\mathrm{li}li changes the inverse by O(ln⁡n)O(\ln n)O(lnn) and does not affect any of the asymptotic statements.
  • The Riemann hypothesis is Mathlib's RiemannHypothesis.

Contributions of reusable infrastructure (cycle-type lcm, prime-counting error terms, the logarithmic integral) are welcome.

Selected references

  • E. Landau, Über die Maximalordnung der Permutationen gegebenen Grades, Arch. Math. Phys. Ser. 3, vol. 5, 1903.
  • J.-L. Nicolas, Sur l'ordre maximum d'un élément dans le groupe SnS_nSn​ des permutations, Acta Arithmetica 14 (1968), 315–332.
  • J.-P. Massias, Majoration explicite de l'ordre maximum d'un élément du groupe symétrique, Ann. Fac. Sci. Toulouse Math. (5) 6 (1984), 269–281.
  • J.-P. Massias, J.-L. Nicolas, G. Robin, Évaluation asymptotique de l'ordre maximum d'un élément du groupe symétrique, Acta Arithmetica 50 (1988), 221–242.
  • K. Ford, Vinogradov's integral and bounds for the Riemann zeta function, Proc. London Math. Soc. 85 (2002), 565–633. https://doi.org/10.1112/S0024611502013655
  • W. Miller, The maximum order of an element of a finite symmetric group, Amer. Math. Monthly 94 (1987), 497–506.
  • OEIS A000793. https://oeis.org/A000793
13 thms2 active usersReviewed
🏆Completed
Graph TheoryLinear algebra·Captain: mikedeng1

Three Conjectures in Extremal Spectral Graph Theory 2: For Large n, K₂ + Pₙ₋₂ Is the Unique Planar Graph on n Vertices of Maximum Spectral Radius (Boots–Royle / Cao–Vince Conjecture)Research Paper

Motivation

The spectral radius of a graph, the largest eigenvalue of its adjacency matrix, measures how densely and how unevenly the edges of the graph are distributed. Asking which graph in a class maximizes it is a spectral analogue of the extremal question "which graph in the class has the most edges". For planar graphs the question has a long history in spectral graph theory, going back at least to Schwenk and Wilson. It also arose in geography, where the spectral radius serves as a measure of the connectivity of planar networks. To compare real networks with a theoretical upper bound, geographers wanted the planar graph of maximum spectral radius (arXiv:1606.01916v2, p. 3).

Timeline, as recorded by Tait and Tobin (pp. 1–3):

  • 1978, Schwenk and Wilson: early study of the eigenvalues of planar graphs.
  • 1991, Boots and Royle, and independently 1993, Cao and Vince: conjecture that the planar graph on n≥9n \ge 9n≥9 vertices of maximum spectral radius is P2+Pn−2P_2 + P_{n-2}P2​+Pn−2​ (Conjecture 1 of the paper). The conjecture is false for n∈{7,8}n \in \{7, 8\}n∈{7,8}, so some lower bound on nnn is necessary (p. 2).
  • 1988–2000: successive improvements of the upper bound on the spectral radius of planar graphs, by Hong (1988, 1995, 1998), Cao and Vince (1993), Guiduli (1996) and Ellingham and Zha (2000). These bounds do not identify the extremal graph.
  • 2010, Dvořák and Mohar: an upper bound on the spectral radius of planar graphs of given maximum degree.
  • Guiduli and Hayes were reported by Ellingham and Zha to have proved the conjecture for large nnn. The preprint never appeared (p. 3).
  • 2017, Tait and Tobin: proof of the conjecture for all sufficiently large nnn (Theorem 15).

Setting

All graphs are finite and simple. A graph on nnn vertices has vertex set {0,1,…,n−1}\{0, 1, \dots, n-1\}{0,1,…,n−1}. Its adjacency matrix A(G)A(G)A(G) is the n×nn \times nn×n real matrix with entry 111 in position (u,w)(u, w)(u,w) when u∼wu \sim wu∼w and 000 otherwise. Its eigenvalues are real, and the largest is the spectral radius λ1(G)\lambda_1(G)λ1​(G).

A graph is planar if it can be drawn in the plane with vertices as distinct points and edges as simple arcs joining their ends, such that an arc meets no vertex other than its ends and two arcs meet only at common ends.

For graphs GGG and HHH on disjoint vertex sets, the join G+HG + HG+H is the graph on the union of the two vertex sets that keeps every edge of GGG and of HHH and adds every edge between a vertex of GGG and a vertex of HHH. Write K2K_2K2​ for a single edge, PmP_mPm​ for the path on mmm vertices, and K2+Pn−2K_2 + P_{n-2}K2​+Pn−2​ for their join: two adjacent vertices, both joined to every vertex of a path on the remaining n−2n - 2n−2 vertices. The paper also writes it P2+Pn−2P_2 + P_{n-2}P2​+Pn−2​.

A planar graph GGG on nnn vertices is a planar spectral maximizer if λ1(G′)≤λ1(G)\lambda_1(G') \le \lambda_1(G)λ1​(G′)≤λ1​(G) for every planar graph G′G'G′ on the same nnn vertices.

The milestones refer to a leading eigenvector v\mathbf vv of a maximizer, scaled so that every entry is positive and the largest entry is 111, and to a vertex xxx with vx=1\mathbf v_x = 1vx​=1. For a parameter ε>0\varepsilon > 0ε>0 they use the sets L={z:vz>ε}L = \{z : \mathbf v_z > \varepsilon\}L={z:vz​>ε} and S=V(G)∖LS = V(G) \setminus LS=V(G)∖L, the edge counts e(X)e(X)e(X) (edges inside XXX) and e(X,Y)e(X, Y)e(X,Y) (edges between XXX and YYY), and a second vertex www with vw>1−24ε\mathbf v_w > 1 - 24\varepsilonvw​>1−24ε that is adjacent to all but 94εn94\varepsilon n94εn vertices of SSS.

Formalization targets

Goal: Theorem 15

∃N0 ∀n≥N0 ∀G on n vertices:G is a planar spectral maximizer  ⟺  G≅K2+Pn−2.\exists N_0\ \forall n \ge N_0\ \forall G \text{ on } n \text{ vertices}: \quad G \text{ is a planar spectral maximizer} \iff G \cong K_2 + P_{n-2}.∃N0​ ∀n≥N0​ ∀G on n vertices:G is a planar spectral maximizer⟺G≅K2​+Pn−2​.

The threshold N0N_0N0​ is not fixed. The paper proves the statement only for nnn beyond an unspecified bound, and the conjecture's bound n≥9n \ge 9n≥9 remains open. The equivalence encodes "the unique maximizer is K2+Pn−2K_2 + P_{n-2}K2​+Pn−2​": every maximizer is isomorphic to it, and it is itself a maximizer.

Milestones, in the order of the proof

  1. If HHH is a proper spanning subgraph of a connected graph GGG, then λ1(H)<λ1(G)\lambda_1(H) < \lambda_1(G)λ1​(H)<λ1​(G) (§1.2, already proved on the platform).
  2. A planar graph on n≥3n \ge 3n≥3 vertices has no K3,3K_{3,3}K3,3​ subgraph and at most 3n−63n - 63n−6 edges, and each of its bipartite subgraphs has at most 2n−42n - 42n−4 edges (p. 8).
  3. Lemma 8: 2n−4<λ1<6n\sqrt{2n-4} < \lambda_1 < \sqrt{6n}2n−4​<λ1​<6n​.
  4. The counts ∣L∣≤32n−4/ε|L| \le 3\sqrt{2n-4}/\varepsilon∣L∣≤32n−4​/ε, e(S,L)≤2n−4e(S, L) \le 2n - 4e(S,L)≤2n−4, e(S)≤3n−6e(S) \le 3n - 6e(S)≤3n−6, e(L)≤92n−4/εe(L) \le 9\sqrt{2n-4}/\varepsilone(L)≤92n−4​/ε (p. 8).
  5. Lemma 9: ∑z∈Lvz≤ε2n−4+18/ε\sum_{z \in L} \mathbf v_z \le \varepsilon\sqrt{2n-4} + 18/\varepsilon∑z∈L​vz​≤ε2n−4​+18/ε and ∑z∈Svz≤(1+3ε)2n−4\sum_{z \in S} \mathbf v_z \le (1 + 3\varepsilon)\sqrt{2n-4}∑z∈S​vz​≤(1+3ε)2n−4​.
  6. Eq. (7): ∑y∼u, y∈Svy≥(vu−ε)2n−4−18/ε\sum_{y \sim u,\, y \in S} \mathbf v_y \ge (\mathbf v_u - \varepsilon)\sqrt{2n-4} - 18/\varepsilon∑y∼u,y∈S​vy​≥(vu​−ε)2n−4​−18/ε for u∈Lu \in Lu∈L.
  7. Lemma 10: vz>1/6n\mathbf v_z > 1/\sqrt{6n}vz​>1/6n​ for every vertex zzz.
  8. ∣{y∈S:y≁x}∣≤14εn|\{y \in S : y \not\sim x\}| \le 14\varepsilon n∣{y∈S:y∼x}∣≤14εn (p. 10).
  9. Lemma 11: a vertex w∈Lw \in Lw∈L, w≠xw \ne xw=x, with vw>1−24ε\mathbf v_w > 1 - 24\varepsilonvw​>1−24ε and ∣{y∈S:y≁w}∣≤94εn|\{y \in S : y \not\sim w\}| \le 94\varepsilon n∣{y∈S:y∼w}∣≤94εn.
  10. Lemma 12: vu<1/10\mathbf v_u < 1/10vu​<1/10 for every u∉{x,w}u \notin \{x, w\}u∈/{x,w}.
  11. Lemma 13: x∼wx \sim wx∼w.
  12. Lemma 14: every vertex other than xxx and www is adjacent to both.
  13. A planar graph with two distinct vertices adjacent to all other vertices is a spanning subgraph of K2+Pn−2K_2 + P_{n-2}K2​+Pn−2​ (proof of Theorem 15, p. 13).

Significance

The result. Theorem 15 settles the Boots–Royle / Cao–Vince conjecture for all large nnn. It identifies the planar graph of maximum spectral radius exactly, so the theoretical upper bound that the geographic application asks for is λ1(K2+Pn−2)\lambda_1(K_2 + P_{n-2})λ1​(K2​+Pn−2​). Along the way, Lemmas 8–14 describe the eigenvector of the extremal graph: two vertices carry almost all of the eigenvector weight, and every other entry is below 1/101/101/10. The same method proves the outerplanar analogue (Theorem 7, the companion mission of this series).

Formalizing it. The theorem is proved on paper and has no machine-checked proof. Formalizing it needs several pieces that are absent from Mathlib or the platform and are reusable: the edge bound 3n−63n - 63n−6 and the bipartite bound 2n−42n - 42n−4 for topologically planar graphs, the non-planarity of K3,3K_{3,3}K3,3​ and K5K_5K5​ and of their subdivisions, Rayleigh-quotient comparison of spectral radii, and the face structure of maximal planar graphs. Each milestone is a self-contained target.

Difficulty

Counting edges does not separate the candidates. A maximizer can be taken to be a maximal planar graph, so it has exactly 3n−63n - 63n−6 edges. Every maximal planar graph on nnn vertices has the same number of edges, and the standard lower bound λ1≥2∣E∣/n\lambda_1 \ge 2|E|/nλ1​≥2∣E∣/n is the same for all of them. The extremal graph is pinned down only through its eigenvector. One has to show that the weight concentrates on two vertices, which are then forced to be adjacent to each other (Lemma 13) and to every other vertex (Lemma 14). Lemma 13 is the step specific to planarity. It uses the triangular faces of a maximal planar graph to find an edge switch that keeps the graph planar. In the formal setting, planarity is topological, so every planarity fact of the argument (the edge counts, the K3,3K_{3,3}K3,3​ and K5K_5K5​ obstructions, the faces of a triangulation) has to be derived from the drawing.

Formalization scope

  • Graphs are SimpleGraph (Fin n). "GGG is K2+Pn−2K_2 + P_{n-2}K2​+Pn−2​" is a graph isomorphism G ≃g bookPath n, where bookPath n lives on Fin 2 ⊕ Fin (n - 2). It is not an equality on a fixed labelling.
  • λ1\lambda_1λ1​ is the published WangKangXue.SpectralTuran.specRad, the largest eigenvalue of the real adjacency matrix.
  • Planarity is the published topological RobertsonSeymour1986.GM5.IsPlanar: a drawing by simple arcs. It is not replaced by K5K_5K5​/K3,3K_{3,3}K3,3​-minor-freeness or by an edge count, either of which would change the class the maximum ranges over.
  • Extremality is IsSpecMax IsPlanar G: GGG is planar and λ1(G′)≤λ1(G)\lambda_1(G') \le \lambda_1(G)λ1​(G′)≤λ1​(G) for every planar G′G'G′ on Fin n. Every milestone from Lemma 8 to Lemma 14 assumes it. These lemmas are false for an arbitrary planar graph, so a statement that drops the hypothesis, or replaces it by a property implying the conclusion, is not acceptable.
  • The eigenvector is universally quantified: every v\mathbf vv with Av=λ1vA\mathbf v = \lambda_1 \mathbf vAv=λ1​v, all entries positive, all entries at most 111, and every xxx with vx=1\mathbf v_x = 1vx​=1. Such a vector exists because the maximizer is connected (Perron–Frobenius). A statement that only asserts the existence of some eigenvector is weaker and is ruled out.
  • Quantifiers. "For nnn sufficiently large" is ∃N ∀n≥N\exists N\ \forall n \ge N∃N ∀n≥N, and NNN may depend on ε\varepsilonε. Milestones 4–8 hold for every ε>0\varepsilon > 0ε>0. Lemmas 11–14 hold for 0<ε<ε00 < \varepsilon < \varepsilon_00<ε<ε0​ with an unspecified ε0\varepsilon_0ε0​, because the paper chooses ε\varepsilonε small. The quantifier on ε\varepsilonε is always outermost, so ε\varepsilonε never depends on nnn or GGG. Lemmas 12–14 hold for every vertex www satisfying the conclusion of Lemma 11 (the predicate IsSecondHub), with the constants 242424 and 949494.
  • Added hypotheses. n≥3n \ge 3n≥3 in the planar facts and in Lemma 8, where the paper's standing "nnn large enough" is needed only to that extent. The page names the vertex of Lemma 12 vvv, the same letter as the eigenvector; it is renamed uuu. Lemma 14's set AAA is stated in unfolded form.
  • Infrastructure that a complete development needs: Euler's formula for topologically planar graphs, the non-planarity of K3,3K_{3,3}K3,3​ and K5K_5K5​, Perron–Frobenius for connected graphs, and the Rayleigh quotient. Proofs of the planar facts (milestone 2) are reusable well beyond this mission. Contributions toward any of them are welcome.

Selected references

  • M. Tait and J. Tobin, Three conjectures in extremal spectral graph theory, J. Combin. Theory Ser. B (2017). Preprint arXiv:1606.01916v2. https://arxiv.org/abs/1606.01916, DOI https://doi.org/10.1016/j.jctb.2017.04.006
  • B. N. Boots and G. F. Royle, A conjecture on the maximum value of the principal eigenvalue of a planar graph, Geographical Analysis 23(3) (1991), 276–282.
  • D. Cao and A. Vince, The spectral radius of a planar graph, Linear Algebra Appl. 187 (1993), 251–257.
  • A. J. Schwenk and R. J. Wilson, On the eigenvalues of a graph, in: Selected Topics in Graph Theory, Academic Press, 1978, 307–336.
  • Hong Yuan, A bound on the spectral radius of graphs, Linear Algebra Appl. 108 (1988), 135–139.
  • M. N. Ellingham and X. Zha, The spectral radius of graphs on surfaces, J. Combin. Theory Ser. B 78(1) (2000), 45–56.
  • Z. Dvořák and B. Mohar, Spectral radius of finite and infinite planar graphs and of graphs of bounded genus, J. Combin. Theory Ser. B 100(6) (2010), 729–739.

Bibliographic data as given in the bibliography of Tait and Tobin (pp. 21–23).

21 thms2 active usersReviewed
🏆Completed
Linear algebraOptimizationTheoretical Computer Science·Captain: mikedeng1

Optimal Approximation for Submodular and Supermodular Optimization with Bounded Curvature II: Column-Subset Selection Has Total Curvature at Most 1 − 1/κ(A)²Research Paper

Why column-subset selection and curvature

The column-subset selection problem (CSSP) asks for kkk columns of a real matrix AAA whose span approximates all of AAA as well as possible in squared Frobenius norm. It is a basic primitive of low-rank approximation, feature selection and numerical linear algebra, studied both for its exact complexity and for approximation algorithms (see Boutsidis, Drineas and Magdon-Ismail, doi:10.1137/12086755X, cited as [2] in the source paper).

Sviridenko, Vondrák and Ward study submodular maximization and supermodular minimization when the objective has bounded total curvature, a parameter c∈[0,1]c \in [0,1]c∈[0,1] introduced by Conforti and Cornuéjols that measures how far a monotone set function is from being linear. Their guarantees improve as ccc decreases. Section 8 of the paper applies this to CSSP: it bounds the curvature of the CSSP objective in terms of the condition number of AAA, so that a well-conditioned matrix has small curvature. This mission formalizes that bound (Lemma 8.2) and the steps of its proof.

Setting

Let AAA be a real m×nm \times nm×n matrix with columns c1,…,cn∈Rmc_1, \dots, c_n \in \mathbb{R}^mc1​,…,cn​∈Rm, so that Ax=∑ixiciAx = \sum_i x_i c_iAx=∑i​xi​ci​ for x∈Rnx \in \mathbb{R}^nx∈Rn, and let ∥⋅∥\|\cdot\|∥⋅∥ be the Euclidean norm. For S⊆[n]S \subseteq [n]S⊆[n] write projS(x)\mathrm{proj}_S(x)projS​(x) for the orthogonal projection of xxx onto span{ci:i∈S}\mathrm{span}\{c_i : i \in S\}span{ci​:i∈S}, the point of that subspace nearest to xxx (proj∅=0\mathrm{proj}_\emptyset = 0proj∅​=0). The CSSP objective is

fA(S)=∑i=1n∥ci−projS(ci)∥2,f^A(S) = \sum_{i=1}^n \|c_i - \mathrm{proj}_S(c_i)\|^2 ,fA(S)=i=1∑n​∥ci​−projS​(ci​)∥2,

the squared Frobenius distance from AAA to the nearest matrix whose columns lie in the span of the selected columns.

For a set function fff on a finite set XXX, the marginal value of iii at AAA is fA(i)=f(A+i)−f(A)f_A(i) = f(A+i) - f(A)fA​(i)=f(A+i)−f(A). The function fff is monotone decreasing if every fA(i)≤0f_A(i) \le 0fA​(i)≤0. Following (1.1) of the paper, a monotone decreasing fff has total curvature at most ccc when, for every j∈Xj \in Xj∈X,

fX−j(j)≤(1−c) f∅(j),f_{X-j}(j) \le (1-c)\, f_\emptyset(j),fX−j​(j)≤(1−c)f∅​(j),

i.e. the last marginal of jjj is at least a (1−c)(1-c)(1−c) fraction of its first marginal in absolute value.

When the columns of AAA are linearly independent, the condition number is

κ(A)=sup⁡∥x∥=1∥Ax∥inf⁡∥x∥=1∥Ax∥∈[1,∞).\kappa(A) = \frac{\sup_{\|x\|=1}\|Ax\|}{\inf_{\|x\|=1}\|Ax\|} \in [1, \infty).κ(A)=inf∥x∥=1​∥Ax∥sup∥x∥=1​∥Ax∥​∈[1,∞).

Formalization targets

Goal: Lemma 8.2 (p. 11)

For AAA with linearly independent columns,

f[n]−jA(j)≤1κ2(A) f∅A(j)for all j∈[n],i.e.c(fA)≤1−1κ2(A).f^A_{[n]-j}(j) \le \frac{1}{\kappa^2(A)}\, f^A_\emptyset(j) \quad \text{for all } j \in [n], \qquad \text{i.e.}\quad c(f^A) \le 1 - \frac{1}{\kappa^2(A)} .f[n]−jA​(j)≤κ2(A)1​f∅A​(j)for all j∈[n],i.e.c(fA)≤1−κ2(A)1​.

Milestones

  1. The identity fA(S)=∑i(∥ci∥2−∥projS(ci)∥2)f^A(S) = \sum_i (\|c_i\|^2 - \|\mathrm{proj}_S(c_i)\|^2)fA(S)=∑i​(∥ci​∥2−∥projS​(ci​)∥2) (p. 11).
  2. Lemma 8.1, monotone part: S↦∥projS(v)∥2S \mapsto \|\mathrm{proj}_S(v)\|^2S↦∥projS​(v)∥2 is monotone increasing (p. 11).
  3. fAf^AfA is monotone decreasing (p. 11).
  4. ∣f∅A(i)∣=∑j∥proj{i}(cj)∥2|f^A_\emptyset(i)| = \sum_j \|\mathrm{proj}_{\{i\}}(c_j)\|^2∣f∅A​(i)∣=∑j​∥proj{i}​(cj​)∥2 (p. 11).
  5. Some unit xxx has ∥Ax∥≥∣f∅A(i)∣\|Ax\| \ge \sqrt{|f^A_\emptyset(i)|}∥Ax∥≥∣f∅A​(i)∣​ (p. 12).
  6. ∣f[n]−iA(i)∣=∥ci−proj[n]−i(ci)∥2|f^A_{[n]-i}(i)| = \|c_i - \mathrm{proj}_{[n]-i}(c_i)\|^2∣f[n]−iA​(i)∣=∥ci​−proj[n]−i​(ci​)∥2 (p. 12).
  7. Some unit x′x'x′ has ∥Ax′∥≤∣f[n]−iA(i)∣\|Ax'\| \le \sqrt{|f^A_{[n]-i}(i)|}∥Ax′∥≤∣f[n]−iA​(i)∣​ (p. 12).
  8. κ2(A)≥f∅A(i)/f[n]−iA(i)\kappa^2(A) \ge f^A_\emptyset(i)/f^A_{[n]-i}(i)κ2(A)≥f∅A​(i)/f[n]−iA​(i) for each iii, under independence (p. 12).

Items 1–7 hold for arbitrary columns; only item 8 and the goal assume independence.

Significance

Lemma 8.2 turns a spectral quantity into a combinatorial one. With it, any curvature-dependent guarantee for monotone decreasing set functions can be stated for CSSP in terms of κ(A)\kappa(A)κ(A). The paper's Lemma 8.3 shows the bound is asymptotically tight, so 1−Θ(1/κ2)1 - \Theta(1/\kappa^2)1−Θ(1/κ2) is the right order of the curvature.

There is a caveat that the formalization makes explicit. The paper deduces from Lemma 8.1 that fAf^AfA is supermodular, and on that basis quotes its supermodular-minimization algorithm for CSSP with factor 1+(κ2−1)e−1+κ2O(ϵ)1 + (\kappa^2 - 1)e^{-1} + \kappa^2 O(\epsilon)1+(κ2−1)e−1+κ2O(ϵ) (p. 12). Both the submodularity half of Lemma 8.1 and the supermodularity of fAf^AfA are false as printed. With c1=e1c_1 = e_1c1​=e1​, c2=e1+e2c_2 = e_1 + e_2c2​=e1​+e2​ and v=e2v = e_2v=e2​ in R2\mathbb{R}^2R2, the function F(S)=∥projS(v)∥2F(S) = \|\mathrm{proj}_S(v)\|^2F(S)=∥projS​(v)∥2 has F({1,2})+F(∅)=1>12=F({1})+F({2})F(\{1,2\}) + F(\emptyset) = 1 > \tfrac12 = F(\{1\}) + F(\{2\})F({1,2})+F(∅)=1>21​=F({1})+F({2}). With columns e1,e1+e2,2e2e_1, e_1+e_2, 2e_2e1​,e1​+e2​,2e2​, fA({1})+fA({2})=152>7=fA({1,2})+fA(∅)f^A(\{1\}) + f^A(\{2\}) = \tfrac{15}{2} > 7 = f^A(\{1,2\}) + f^A(\emptyset)fA({1})+fA({2})=215​>7=fA({1,2})+fA(∅). A small perturbation gives a counterexample with independent columns. This is the "suppressor variable" phenomenon of regression. The curvature bound itself does not depend on supermodularity and is true, so it is the goal here. The approximation factor for CSSP is not posed, because it rests on the false claim and is stated with O(ϵ)O(\epsilon)O(ϵ).

As far as is known, none of these statements has a machine-checked proof. The results are proved on paper, except the two false claims above. The mission adds a formal proof of the curvature–conditioning bound, and makes explicit which parts of §8 survive.

Difficulty

The individual steps are short linear algebra, but they rest on infrastructure that has to be assembled: orthogonal projections onto spans of finite column families, the Pythagorean identities that relate fAf^AfA to projections, and the variational characterization of sup⁡\supsup and inf⁡\infinf of ∥Ax∥\|Ax\|∥Ax∥ over the unit sphere.

The paper's own route to the large unit vector (milestone 5) does not work as printed. It sets xj∝∥proj{i}(cj)∥x_j \propto \|\mathrm{proj}_{\{i\}}(c_j)\|xj​∝∥proj{i}​(cj​)∥ and uses ci⋅cj=∥ci∥ ∥proj{i}(cj)∥c_i \cdot c_j = \|c_i\|\,\|\mathrm{proj}_{\{i\}}(c_j)\|ci​⋅cj​=∥ci​∥∥proj{i}​(cj​)∥, which is false when ci⋅cj<0c_i \cdot c_j < 0ci​⋅cj​<0, so the signs of the inner products have to be accounted for. In the final step, the ratio f∅A(i)/f[n]−iA(i)f^A_\emptyset(i)/f^A_{[n]-i}(i)f∅A​(i)/f[n]−iA​(i) is a quotient of two negative numbers. Comparing it with κ2\kappa^2κ2 requires the supremum to be finite and the infimum to be strictly positive, which is exactly what independence of the columns provides.

Formalization scope

  • Columns are a family c : Fin n → EuclideanSpace ℝ (Fin m), and AxAxAx is ∑ j, x j • c j. Subsets of [n][n][n] are Finset (Fin n). [n]−i[n] - i[n]−i is Finset.univ.erase i and S+iS + iS+i is insert i S.
  • projS\mathrm{proj}_SprojS​ is Mathlib's Submodule.starProjection onto Submodule.span ℝ (c '' S), which is finite-dimensional, so the argmin of the page exists and is unique.
  • κ(A)\kappa(A)κ(A) is a real quotient of ⨆ and ⨅ of ∥Ax∥\|Ax\|∥Ax∥ over the unit-sphere subtype of Rn\mathbb{R}^nRn. It is used only under LinearIndependent ℝ c, the page's "non-singular", where both are finite and positive. The page's κ=∞\kappa = \inftyκ=∞ for dependent columns is not encoded.
  • Monotonicity is the marginal form of §2.1. Curvature is the product form above, which agrees with the ratio of (1.1) whenever f∅(j)≠0f_\emptyset(j) \ne 0f∅​(j)=0 and avoids Lean's x/0=0x/0 = 0x/0=0.
  • For n=0n = 0n=0 the goal has no jjj and is empty, as on the page.
  • A trivializing reading is ruled out by the encoding: had κ\kappaκ evaluated to a junk value 000, the goal would reduce to monotonicity. Under independence the supremum and infimum are genuine, so the bound 1−1/κ21 - 1/\kappa^21−1/κ2 is the paper's.
  • Not in scope: Lemma 8.3 (stated with O(⋅)O(\cdot)O(⋅)), the CSSP approximation factor, and the false submodularity and supermodularity claims.

Useful infrastructure that would be reusable beyond this mission includes Pythagoras for starProjection, monotonicity of projection norms under inclusion of subspaces, and extremal characterizations of sup⁡∥x∥=1∥Ax∥\sup_{\|x\|=1}\|Ax\|sup∥x∥=1​∥Ax∥ and inf⁡∥x∥=1∥Ax∥\inf_{\|x\|=1}\|Ax\|inf∥x∥=1​∥Ax∥. Contributions of any of the milestones, or of these general lemmas, are welcome.

Selected references

  • M. Sviridenko, J. Vondrák, J. Ward, Optimal approximation for submodular and supermodular optimization with bounded curvature, SODA 2015 (version of Oct. 9, 2014); journal version Math. Oper. Res. 42(4), 2017. https://doi.org/10.1287/moor.2016.0842
  • M. Conforti, G. Cornuéjols, Submodular set functions, matroids and the greedy algorithm: tight worst-case bounds and some generalizations of the Rado–Edmonds theorem, Discrete Appl. Math. 7(3), 1984. https://doi.org/10.1016/0166-218X(84)90003-9
  • C. Boutsidis, P. Drineas, M. Magdon-Ismail, Near-optimal column-based matrix reconstruction, SIAM J. Comput. 43(2), 2014. https://doi.org/10.1137/12086755X
  • A. Das, D. Kempe, Submodular meets spectral: greedy algorithms for subset selection, sparse approximation and dictionary selection, ICML 2011. https://arxiv.org/abs/1102.3975
11 thms2 active usersReviewed
🏆Completed
Graph TheoryLinear algebra·Captain: mikedeng1

Three Conjectures in Extremal Spectral Graph Theory 1: For Large n, K₁ + Pₙ₋₁ Is the Unique Outerplanar Graph on n Vertices of Maximum Spectral Radius (Cvetković–Rowlinson Conjecture)Research Paper

Motivation

The spectral radius λ1(G)\lambda_1(G)λ1​(G) of a graph GGG, the largest eigenvalue of its adjacency matrix, is a standard measure of how densely and how centrally a graph is connected. A classical theme of spectral graph theory asks, for a natural family of graphs, which member of the family on nnn vertices maximizes λ1\lambda_1λ1​. For planar graphs this question was raised by geographers comparing road and transport networks, and for outerplanar graphs it was posed by Cvetković and Rowlinson in their 1990 survey of the largest eigenvalue of a graph, where they relate it to the study of subfamilies of Hamiltonian graphs.

Conjecture (Cvetković–Rowlinson, 1990). The outerplanar graph on nnn vertices of maximum spectral radius is K1+Pn−1K_1 + P_{n-1}K1​+Pn−1​.

Timeline:

  • 1990. Cvetković and Rowlinson state the conjecture (Linear and Multilinear Algebra 28).
  • 1990. Rowlinson makes partial progress, on the index of certain outerplanar graphs (Ars Combinatoria 29).
  • 1993. Cao and Vince work on the spectral radius of planar graphs and the related conjecture (Linear Algebra Appl. 187).
  • 2001. Zhou, Lin and Hu study the spectral radius of Hamiltonian planar and outerplanar graphs (Tsinghua Science and Technology 6).
  • 2017. Tait and Tobin prove the conjecture for all sufficiently large nnn (arXiv:1606.01916, J. Combin. Theory Ser. B, doi:10.1016/j.jctb.2017.04.006), Theorem 7.

This mission formalizes Section 2 of Tait and Tobin.

Setting

All graphs are finite and simple. An nnn-vertex graph GGG has vertex set V(G)={0,…,n−1}V(G) = \{0, \dots, n-1\}V(G)={0,…,n−1}, adjacency matrix A(G)A(G)A(G), and spectral radius λ1(G)\lambda_1(G)λ1​(G), the largest eigenvalue of the real symmetric matrix A(G)A(G)A(G). Write dud_udu​ for the degree of a vertex uuu and N(u)N(u)N(u) for its neighbourhood.

A plane drawing of GGG places the vertices at distinct points of R2\mathbb R^2R2 and each edge as a simple arc between its ends; an arc passes through no other vertex, and two arcs meet only at a common end. The union of the vertex points and arcs is a compact set, so its complement has exactly one unbounded connected component, the outer face. A graph is outerplanar if it has a plane drawing in which every vertex lies on the boundary of the outer face.

The join G+HG + HG+H of two graphs is their disjoint union together with all edges between them. The fan K1+Pn−1K_1 + P_{n-1}K1​+Pn−1​ is a single vertex joined to every vertex of a path on n−1n-1n−1 vertices; it is outerplanar and has 2n−32n-32n−3 edges.

An outerplanar graph GGG on nnn vertices is a spectral maximizer if λ1(G′)≤λ1(G)\lambda_1(G') \le \lambda_1(G)λ1​(G′)≤λ1​(G) for every outerplanar graph G′G'G′ on the same nnn vertices.

Several milestones refer to a Perron vector: an eigenvector v\mathbf vv of A(G)A(G)A(G) for λ1\lambda_1λ1​ with all entries positive, normalized so that its maximum entry is 111, together with a vertex xxx with vx=1\mathbf v_x = 1vx​=1. For a connected graph such a vector exists by the Perron–Frobenius theorem.

Formalization targets

Goal: Theorem 7

There is NNN such that for every n≥Nn \ge Nn≥N and every graph GGG on nnn vertices,

G is an outerplanar spectral maximizer  ⟺  G≅K1+Pn−1.G \text{ is an outerplanar spectral maximizer} \iff G \cong K_1 + P_{n-1}.G is an outerplanar spectral maximizer⟺G≅K1​+Pn−1​.

The forward direction is the paper's Theorem 7; the converse says that the fan is the maximizer, as the conjecture reads. The threshold NNN is left unspecified, as in the paper.

Milestones

The milestones follow the proof in Section 2, in order:

  1. a strict subgraph of a connected graph has strictly smaller λ1\lambda_1λ1​ (already proved on the platform);
  2. an outerplanar graph on n≥2n \ge 2n≥2 vertices has at most 2n−32n-32n−3 edges and no K2,3K_{2,3}K2,3​ subgraph;
  3. Lemma 3: for n≥3n \ge 3n≥3 a spectral maximizer has λ1>n−1\lambda_1 > \sqrt{n-1}λ1​>n−1​;
  4. Lemma 4: du>vun−11nd_u > \mathbf v_u n - 11\sqrt ndu​>vu​n−11n​ for every vertex uuu;
  5. Lemma 5: for nnn large, dx>n−11nd_x > n - 11\sqrt ndx​>n−11n​ and vu<C1/n\mathbf v_u < C_1/\sqrt nvu​<C1​/n​ for every u≠xu \ne xu=x, with an absolute constant C1C_1C1​;
  6. Lemma 6: for nnn large, ∑z∈Bvz<C2/n\sum_{z\in B}\mathbf v_z < C_2/\sqrt n∑z∈B​vz​<C2​/n​, where B=V(G)∖(N(x)∪{x})B = V(G)\setminus(N(x)\cup\{x\})B=V(G)∖(N(x)∪{x});
  7. for nnn large, B=∅B = \emptysetB=∅: xxx is adjacent to every other vertex;
  8. an outerplanar graph with a vertex adjacent to all others is a spanning subgraph of K1+Pn−1K_1 + P_{n-1}K1​+Pn−1​.

Significance

The theorem settles, for all large nnn, a conjecture that had been open since 1990 and had resisted several partial attacks. The same paper settles, in the same framework, the planar analogue (the Boots–Royle and Cao–Vince conjecture, K2+Pn−2K_2 + P_{n-2}K2​+Pn−2​) and a conjecture on spectral irregularity. The outerplanar case is the shortest of the three.

On the formal side, the result is proved on paper but has no machine-checked proof. Completing the mission produces a formal theory of outerplanar plane drawings (edge bound, K2,3K_{2,3}K2,3​-freeness, outerplanarity of the fan), eigenvector estimates for extremal graphs, and the full extremal theorem. The edge bound and the K2,3K_{2,3}K2,3​ exclusion are reusable for any later work on outerplanar or series–parallel graphs. Determining the exact threshold NNN, or proving the statement for every n≥2n \ge 2n≥2, would be a strengthening beyond the paper.

Difficulty

There are two separate difficulties. The first is topological: the facts about outerplanar graphs that the spectral argument uses (at most 2n−32n-32n−3 edges, no K2,3K_{2,3}K2,3​ subgraph, no cycle in G−xG - xG−x when xxx dominates, outerplanarity of the fan) are folklore on paper but require a Jordan-curve-type argument about plane drawings, of which Mathlib currently contains very little. The second is spectral: the edge-switching step that proves B=∅B = \emptysetB=∅ compares λ1\lambda_1λ1​ of two graphs through the Rayleigh quotient, and needs a positive Perron vector and the variational characterization of the largest eigenvalue for symmetric matrices. A purely combinatorial edge-maximization argument does not work: the fan is not the only outerplanar graph with 2n−32n-32n−3 edges, so maximizing the number of edges does not identify the spectral maximizer.

Formalization scope

  • Graphs are SimpleGraph (Fin n); λ1\lambda_1λ1​ is the published WangKangXue.SpectralTuran.specRad, the largest eigenvalue of G.adjMatrix ℝ.
  • Outerplanarity is topological: a published RobertsonSeymour1986.GM5.PlaneDrawing with every vertex point in the closure of the unbounded component of the complement of the drawing. It is not replaced by a forbidden-minor characterization; such characterizations are theorems about outerplanarity, and would change which graphs compete in the maximization.
  • Extremality is IsSpecMax IsOuterplanar G: GGG is outerplanar and no outerplanar graph on the same vertex set has larger λ1\lambda_1λ1​. It is never weakened to a comparison with the fan.
  • "GGG is K1+Pn−1K_1 + P_{n-1}K1​+Pn−1​" is a graph isomorphism G ≃g fan n, with the fan on Fin 1 ⊕ Fin (n - 1).
  • Perron-vector statements quantify over every positive eigenvector for λ1\lambda_1λ1​ with maximum entry 111 and every vertex xxx with vx=1\mathbf v_x = 1vx​=1; existential quantification over v\mathbf vv would be weaker.
  • "For some absolute constant" is an existential quantifier placed before the threshold, the graph and the vector; "for nnn sufficiently large" is ∃ N, ∀ n ≥ N.
  • Added hypotheses, each forced by a degenerate case: n≥2n \ge 2n≥2 for the edge bound (the bound 2n−32n-32n−3 is negative below), n≥3n \ge 3n≥3 for Lemma 3 (it fails for n=2n = 2n=2, where the maximizer K2K_2K2​ has λ1=1=n−1\lambda_1 = 1 = \sqrt{n-1}λ1​=1=n−1​).

A complete development needs the Perron–Frobenius theorem for symmetric nonnegative matrices and the Rayleigh quotient characterization of λ1\lambda_1λ1​, both of general use; the theory of outerplanar drawings; and the counting arguments of the lemmas. Contributions to any of these layers are welcome, as are independent proofs of the outerplanar facts from the drawing definition.

Selected references

  • M. Tait and J. Tobin, Three conjectures in extremal spectral graph theory, J. Combin. Theory Ser. B (2017), https://doi.org/10.1016/j.jctb.2017.04.006; preprint arXiv:1606.01916v2, https://arxiv.org/abs/1606.01916
  • D. Cvetković and P. Rowlinson, The largest eigenvalue of a graph: a survey, Linear and Multilinear Algebra 28 (1990), 3–33.
  • P. Rowlinson, On the index of certain outerplanar graphs, Ars Combinatoria 29 (1990), 221–225.
  • D. Cao and A. Vince, The spectral radius of a planar graph, Linear Algebra Appl. 187 (1993), 251–257.
  • J. Zhou, C. Lin and G. Hu, Spectral radius of Hamiltonian planar graphs and outerplanar graphs, Tsinghua Science and Technology 6 (2001), 350–354.
12 thms2 active usersReviewed
Graph TheoryProbability·Captain: mikedeng1

Algorithms for #BIS-Hard Problems on Expander Graphs 4: On Almost Every Δ-Regular Bipartite Graph with Δ ≥ Cq² log² q, Proper q-Colorings Reduce to Kotecký–Preiss Polymer ModelsResearch Paper

Motivation

Counting the proper qqq-colorings of a graph, and sampling one uniformly at random, are basic problems of approximate counting with a long history in statistical physics, where proper colorings are the zero-temperature antiferromagnetic Potts model. On bipartite graphs these problems are tied to #BIS, the problem of approximately counting independent sets in a bipartite graph, whose complexity is unresolved. Galanis, Štefankovič, Vigoda and Yang showed that approximately counting proper qqq-colorings is #BIS-hard on bipartite graphs of maximum degree Δ\DeltaΔ whenever q≤Δ/(2log⁡Δ)q \le \Delta/(2\log\Delta)q≤Δ/(2logΔ) (SIAM J. Comput. 45(6), 2016). A positive result in that regime therefore has to use structure beyond bipartiteness and bounded degree.

Jenssen, Keevash and Perkins (SIAM J. Comput. 49(4) 2020) use expansion. Their Theorem 5 says that for Δ≥Cq2log⁡2q\Delta \ge Cq^2\log^2 qΔ≥Cq2log2q there is an FPTAS and a polynomial-time sampler for proper qqq-colorings on almost every Δ\DeltaΔ-regular bipartite graph. This answered Conjecture 1 of the authors' SODA 2019 extended abstract. Liao, Lin, Lu and Mao (APPROX/RANDOM 2019) independently proved a similar result under a slightly stronger condition on Δ\DeltaΔ. The organizing idea, that most colorings of a strongly expanding bipartite graph follow one of a few patterns, comes from the work of Peled and Spinka on rigidity of colorings of Zd\mathbb Z^dZd (arXiv:1808.03597).

This mission is the fourth of a series formalizing the four main results of that paper; it covers Theorem 5.

Setting

Fix integers q≥3q \ge 3q≥3, Δ≥3\Delta \ge 3Δ≥3 and m≥1m \ge 1m≥1. The vertex set is V=O⊔EV = \mathcal O \sqcup \mathcal EV=O⊔E, two sides of mmm vertices each (so n=2mn = 2mn=2m). Gbip(2m,Δ)\mathcal G^{\mathrm{bip}}(2m,\Delta)Gbip(2m,Δ) is the finite set of labelled simple graphs on VVV whose edges all join O\mathcal OO to E\mathcal EE and in which every vertex has degree Δ\DeltaΔ. A property holds for almost every Δ\DeltaΔ-regular bipartite graph if the fraction of graphs in Gbip(2m,Δ)\mathcal G^{\mathrm{bip}}(2m,\Delta)Gbip(2m,Δ) having it tends to 111 as m→∞m \to \inftym→∞.

For S⊆VS \subseteq VS⊆V the vertex boundary ∂S\partial S∂S is the set of vertices outside SSS with a neighbour in SSS, and S+=S∪∂SS^+ = S \cup \partial SS+=S∪∂S. The graph is a bipartite (σ,ρ)(\sigma,\rho)(σ,ρ)-expander if ∣∂S∣≥ρ∣S∣|\partial S| \ge \rho|S|∣∂S∣≥ρ∣S∣ for every SSS inside one side with ∣S∣≤σm|S| \le \sigma m∣S∣≤σm.

A proper qqq-coloring is a map f:V→[q]f : V \to [q]f:V→[q] with f(u)≠f(v)f(u) \ne f(v)f(u)=f(v) on every edge; Xq(G)\mathcal X_q(G)Xq​(G) is their set and ZG(q)=∣Xq(G)∣Z_G(q) = |\mathcal X_q(G)|ZG​(q)=∣Xq​(G)∣. A pattern is a pair (A,B)(A,B)(A,B) of disjoint sets with A∪B=[q]A \cup B = [q]A∪B=[q]; there are 2q2^q2q of them. A coloring agrees with (A,B)(A,B)(A,B) at vvv if v∈Ov \in \mathcal Ov∈O and f(v)∈Af(v) \in Af(v)∈A, or v∈Ev \in \mathcal Ev∈E and f(v)∈Bf(v) \in Bf(v)∈B. For S⊆VS \subseteq VS⊆V, χA,B(S)\chi_{A,B}(S)χA,B​(S) is the set of proper colorings that disagree with (A,B)(A,B)(A,B) at every vertex of SSS and agree everywhere else.

G3G^3G3 joins two vertices at GGG-distance at most 333. A set is little if ∣S∣≤4qlog⁡ΔΔm|S| \le 4q\frac{\log\Delta}{\Delta}m∣S∣≤4qΔlogΔ​m. A polymer is a little set that is connected in G3G^3G3; C(G)\mathcal C(G)C(G) is the set of polymers, and two polymers are compatible if their GGG-distance exceeds 333. For a pattern (A,B)(A,B)(A,B), the weight and the polymer partition function are

wA,B(γ)=∣χA,B(γ)∣∣A∣m∣B∣m,ΞA,B(G)=∑Γ∏γ∈ΓwA,B(γ),w_{A,B}(\gamma) = \frac{|\chi_{A,B}(\gamma)|}{|A|^m|B|^m}, \qquad \Xi_{A,B}(G) = \sum_{\Gamma} \prod_{\gamma\in\Gamma} w_{A,B}(\gamma),wA,B​(γ)=∣A∣m∣B∣m∣χA,B​(γ)∣​,ΞA,B​(G)=Γ∑​γ∈Γ∏​wA,B​(γ),

the sum over all families Γ⊆C(G)\Gamma \subseteq \mathcal C(G)Γ⊆C(G) of pairwise compatible polymers, the empty family included. Finally Z^\hat ZZ^ is an ε\varepsilonε-relative approximation to ZZZ if e−εZ^≤Z≤eεZ^e^{-\varepsilon}\hat Z \le Z \le e^{\varepsilon}\hat Ze−εZ^≤Z≤eεZ^.

Formalization targets

Goal: the deterministic core of Theorem 5

There is an absolute constant C>0C > 0C>0 such that for all q≥3q \ge 3q≥3 and Δ≥Cq2log⁡2q\Delta \ge Cq^2\log^2 qΔ≥Cq2log2q, almost every Δ\DeltaΔ-regular bipartite graph GGG satisfies

e−εZ~G(q)≤ZG(q)≤eεZ~G(q),Z~G(q)=∑(A,B)∣A∣m∣B∣m ΞA,B(G),ε=e−m/(4q),e^{-\varepsilon}\tilde Z_G(q) \le Z_G(q) \le e^{\varepsilon}\tilde Z_G(q), \qquad \tilde Z_G(q) = \sum_{(A,B)} |A|^m|B|^m\,\Xi_{A,B}(G), \quad \varepsilon = e^{-m/(4q)},e−εZ~G​(q)≤ZG​(q)≤eεZ~G​(q),Z~G​(q)=(A,B)∑​∣A∣m∣B∣mΞA,B​(G),ε=e−m/(4q),

and, for every pattern (A,B)(A,B)(A,B) and every polymer γ∈C(G)\gamma \in \mathcal C(G)γ∈C(G), the Kotecký–Preiss condition

∑γ′∈C(G): d(γ′,γ)≤3wA,B(γ′) e∣γ′∣+g(γ′)≤∣γ∣,g(γ)=Δ10q2log⁡Δ∣γ∣.\sum_{\gamma' \in \mathcal C(G):\ d(\gamma',\gamma)\le 3} w_{A,B}(\gamma')\,e^{|\gamma'| + g(\gamma')} \le |\gamma|, \qquad g(\gamma) = \frac{\Delta}{10q^2\log\Delta}|\gamma|.γ′∈C(G): d(γ′,γ)≤3∑​wA,B​(γ′)e∣γ′∣+g(γ′)≤∣γ∣,g(γ)=10q2logΔΔ​∣γ∣.

The constant CCC is left existential, as in the paper.

Milestones

In the order of the proof:

  1. Lemma 23, the expansion of almost every regular bipartite graph;
  2. Lemma 27, the bound on ∣χA,B(S)∣|\chi_{A,B}(S)|∣χA,B​(S)∣;
  3. Lemma 28, the covering of Xq(G)\mathcal X_q(G)Xq​(G) by χA,B(S)\chi_{A,B}(S)χA,B​(S) with SSS little;
  4. the product formula after (13);
  5. Claim 30, the size bound for sparse sets;
  6. the binomial estimate, an existing platform theorem;
  7. Lemma 29, the reduction;
  8. the weight bound of §5.1;
  9. the per-vertex estimate (15).

Significance

The goal is the whole mathematical content of Theorem 5. The paper's algorithm (its Theorem 8, adapted from Helmuth, Perkins and Regts, Algorithmic Pirogov–Sinai theory) takes exactly these two facts as input. Given the Kotecký–Preiss condition, the cluster expansion of each log⁡ΞA,B(G)\log\Xi_{A,B}(G)logΞA,B​(G) converges and truncating it gives the estimator; the reduction transfers the estimate to ZG(q)Z_G(q)ZG​(q). With the goal proved, the remaining distance to the FPTAS is algorithmic bookkeeping, not combinatorics.

The result has been proved since 2019 and is not formalized anywhere. A formal proof would add a machine-checked polymer-model reduction for a spin system with exponentially many ground states, a reusable model of uniformly random regular bipartite graphs, and a checked instance of the "expansion implies few ground states" argument used throughout this literature.

Difficulty

The obvious first idea is a one-ground-state expansion, as for the hard-core model, where configurations are close to "all even" or "all odd". For colorings there are 2q2^q2q patterns. A coloring can lie close to several of them, so the reduction must show that the sets χA,B(S)\chi_{A,B}(S)χA,B​(S) with SSS little overlap only negligibly. That needs an inclusion–exclusion step whose error is exponentially small in m/qm/qm/q. A second difficulty is the size of polymers: little sets are qqq times larger than the sets of the hard-core analysis. Claim 30 is what keeps the union of compatible polymers from becoming a large set that the expansion no longer controls. Finally, the constants must be uniform in qqq. The single condition Δ≥Cq2log⁡2q\Delta \ge Cq^2\log^2 qΔ≥Cq2log2q has to cover the expansion lemma, the reduction and the Kotecký–Preiss estimate at once.

Formalization scope

The FPTAS and the sampling algorithm are not formalized. Mathlib has no machine model in which "polynomial time" could be stated, and the running times carry unquantified constants. The goal is the mathematics the algorithm rests on: the reduction and the Kotecký–Preiss condition, with the paper's constants e−m/(4q)e^{-m/(4q)}e−m/(4q), 4qlog⁡Δ/Δ4q\log\Delta/\Delta4qlogΔ/Δ and Δ/(10q2log⁡Δ)\Delta/(10q^2\log\Delta)Δ/(10q2logΔ). Neither half may be dropped. Replacing e−m/(4q)e^{-m/(4q)}e−m/(4q) by an existential error, or summing the Kotecký–Preiss condition over compatible rather than incompatible polymers, would make the goal trivial or wrong.

Conventions committed to in Lean:

  • Vertices are Fin m ⊕ Fin m, and Gbip(2m,Δ)\mathcal G^{\mathrm{bip}}(2m,\Delta)Gbip(2m,Δ) is the set of labelled Δ\DeltaΔ-regular graphs with edges only between the two fixed sides. The paper does not fix a bipartition; this fixed-sides model is a disclosed choice, matching §5's ∣O∣=∣E∣=m|\mathcal O| = |\mathcal E| = m∣O∣=∣E∣=m.
  • Distances are Mathlib's extended distance edist, so vertices in different components are at distance ∞\infty∞ and never incompatible.
  • A pattern is encoded by A⊆[q]A \subseteq [q]A⊆[q] with B=[q]∖AB = [q]\setminus AB=[q]∖A. Patterns with an empty side are kept: their weights are 000 under Lean's x/0=0x/0 = 0x/0=0, so they contribute 000 to Z~G(q)\tilde Z_G(q)Z~G​(q), as on the page for large mmm.
  • log is the natural logarithm; the binomial estimate is used in nats.

Added hypotheses:

  • A,B≠∅A, B \ne \emptysetA,B=∅ in Lemma 27 and in the product formula, which divide by ∣A∣|A|∣A∣ and ∣B∣|B|∣B∣;
  • Δ≥3\Delta \ge 3Δ≥3 in Lemma 28 and Claim 30, the paper's standing convention (§1.3);
  • m≥1m \ge 1m≥1 in Lemma 29, which fails on the empty vertex set.

The goal needs none of these, since its limit is in mmm and its constant absorbs Δ\DeltaΔ.

A complete development needs a counting theory of random regular bipartite graphs (for Lemma 23, through Bassalygo's expansion theorem), counts of connected sets in bounded-degree graphs, and inclusion–exclusion over patterns. The random-graph model and the polymer-model definitions are reusable across this series and beyond. Proofs of any milestone, and alternative routes to Lemma 23, are welcome.

Selected references

  • M. Jenssen, P. Keevash and W. Perkins, Algorithms for #BIS-hard problems on expander graphs, SIAM J. Comput. 49(4) (2020) 681–710. https://doi.org/10.1137/19M1286669
  • A. Galanis, D. Štefankovič, E. Vigoda and L. Yang, Ferromagnetic Potts model: refined #BIS-hardness and related results, SIAM J. Comput. 45(6) (2016) 2004–2065.
  • R. Peled and Y. Spinka, Rigidity of proper colorings of Zd\mathbb Z^dZd, arXiv:1808.03597 (2018). https://arxiv.org/abs/1808.03597
  • C. Liao, J. Lin, P. Lu and Z. Mao, Counting independent sets and colorings on random regular bipartite graphs, APPROX/RANDOM 2019, Schloss Dagstuhl–Leibniz-Zentrum für Informatik (2019).
  • T. Helmuth, W. Perkins and G. Regts, Algorithmic Pirogov–Sinai theory, Probab. Theory Related Fields (2020); extended abstract at STOC 2019.
  • R. Kotecký and D. Preiss, Cluster expansion for abstract polymer models, Comm. Math. Phys. 103(3) (1986) 491–498.
  • L. A. Bassalygo, Asymptotically optimal switching circuits, Problems Inform. Transmission 17(3) (1981) 206–211.
13 thms2 active usersReviewed
Discrete GeometryLinear OptimizationOperations Research·Captain: mikedeng1

An Improved Kalai-Kleitman Bound for the Diameter of a Polyhedron: Every d-Polyhedron with n Facets, Bounded or Not, Has Graph Diameter at Most (n-d)^(log2 d)Research Paper

Motivation

The simplex method for linear programming moves from vertex to vertex along the edges of a polyhedron. A lower bound on the number of pivots of any edge-following method is the graph diameter of the feasible polyhedron: the largest number of edges needed to connect two of its vertices. Whether this diameter is bounded by a polynomial in the dimension ddd and the number nnn of facets is the polynomial Hirsch conjecture, open since the 1950s. The original Hirsch conjecture, n−dn-dn−d, was refuted for unbounded polyhedra by Klee and Walkup (1967) and for bounded ones by Santos (2012) (doi:10.4007/annals.2012.176.1.7).

Timeline.

  • 1964–1966. Klee proves that a three-dimensional polyhedron with nnn facets, bounded or not, has diameter at most n−3n-3n−3 (doi:10.4153/CJM-1964-058-3, doi:10.1137/0113061, doi:10.2140/pjm.1966.17.249).
  • 1967. Klee and Walkup show that the Hirsch bound fails for unbounded polyhedra in dimension four (doi:10.1007/BF02395039).
  • 1992. Kalai and Kleitman prove the first subexponential bound, nlog⁡2d+2n^{\log_2 d+2}nlog2​d+2, for all polyhedra (doi:10.1090/S0273-0979-1992-00285-9, arXiv:math/9204233).
  • 2010. Eisenbrand, Hähnle, Razborov and Rothvoss show that nlog⁡2d+1n^{\log_2 d+1}nlog2​d+1 holds for a broad combinatorial abstraction of polyhedra (doi:10.1287/moor.1100.0470).
  • 2014. Todd sharpens the Kalai–Kleitman bound to (n−d)log⁡2d(n-d)^{\log_2 d}(n−d)log2​d for all polyhedra, bounded or not (arXiv:1402.3579; SIAM J. Discrete Math. 28 (2014), doi:10.1137/140962310). This is the result of this mission.

Setting

A polyhedron P⊆RdP\subseteq\mathbb R^dP⊆Rd is the set of solutions of finitely many linear inequalities,

P=Hpoly⁡(a,b)={x∈Rd: ⟨ai,x⟩≤bi, i=1,…,n}.P=\operatorname{Hpoly}(a,b)=\{x\in\mathbb R^d:\ \langle a_i,x\rangle\le b_i,\ i=1,\dots,n\}.P=Hpoly(a,b)={x∈Rd: ⟨ai​,x⟩≤bi​, i=1,…,n}.

It need not be bounded. A vertex of PPP is a point vvv for which some half-space HHH satisfies P∩H={v}P\cap H=\{v\}P∩H={v}, equivalently an extreme point of PPP. Two vertices v≠wv\neq wv=w are adjacent when P∩H=[v,w]P\cap H=[v,w]P∩H=[v,w] for some half-space HHH; the segment [v,w][v,w][v,w] is then an edge. A path of length kkk is a sequence of vertices v=v0,…,vk=wv=v_0,\dots,v_k=wv=v0​,…,vk​=w with consecutive vertices adjacent, ρP(v,w)\rho_P(v,w)ρP​(v,w) is the least such kkk, and the diameter is δ(P)=max⁡ρP(v,w)\delta(P)=\max\rho_P(v,w)δ(P)=maxρP​(v,w) over pairs of vertices. A (d,n)(d,n)(d,n)-polyhedron is a ddd-dimensional polyhedron with nnn facets, and

Δ(d,n)=max⁡{δ(P): P a (d,n)-polyhedron},\Delta(d,n)=\max\{\delta(P):\ P \text{ a } (d,n)\text{-polyhedron}\},Δ(d,n)=max{δ(P): P a (d,n)-polyhedron},

the maximum taken over bounded and unbounded polyhedra alike. All logarithms are to base 2.

Formalization targets

Goal: Theorem 1

Δ(d,n)≤(n−d)log⁡d(1≤d≤n),Δ(1,1)=0.\Delta(d,n)\le (n-d)^{\log d}\quad(1\le d\le n),\qquad \Delta(1,1)=0.Δ(d,n)≤(n−d)logd(1≤d≤n),Δ(1,1)=0.

The goal fixes no constant beyond the paper's own exponent and leaves no case of ddd, nnn or boundedness open.

Milestones, in the order the proof uses them

  1. The identity (n−d)log⁡d=dlog⁡(n−d)(n-d)^{\log d}=d^{\log(n-d)}(n−d)logd=dlog(n−d) (p. 2).
  2. Two vertices on a common facet are at distance at most Δ(d−1,n−1)\Delta(d-1,n-1)Δ(d−1,n−1) (p. 2, in the proof of Lemma 1).
  3. Lemma 1 (Kalai–Kleitman): for 2≤d≤⌊n/2⌋2\le d\le\lfloor n/2\rfloor2≤d≤⌊n/2⌋, Δ(d,n)≤Δ(d−1,n−1)+2Δ(d,⌊n/2⌋)+2\Delta(d,n)\le\Delta(d-1,n-1)+2\Delta(d,\lfloor n/2\rfloor)+2Δ(d,n)≤Δ(d−1,n−1)+2Δ(d,⌊n/2⌋)+2.
  4. Δ(d,d)=0\Delta(d,d)=0Δ(d,d)=0.
  5. Δ(1,n)≤1\Delta(1,n)\le 1Δ(1,n)≤1.
  6. Δ(2,n)≤n−2\Delta(2,n)\le n-2Δ(2,n)≤n−2.
  7. Δ(3,n)≤n−3\Delta(3,n)\le n-3Δ(3,n)≤n−3 (Klee), for unbounded polyhedra too.
  8. If n<2dn<2dn<2d, then Δ(d,n)≤Δ(d−1,n−1)\Delta(d,n)\le\Delta(d-1,n-1)Δ(d,n)≤Δ(d−1,n−1).
  9. In the ten cases d=4, 8≤n≤11d=4,\ 8\le n\le 11d=4, 8≤n≤11; d=5, 10≤n≤12d=5,\ 10\le n\le 12d=5, 10≤n≤12; d=6, 12≤n≤13d=6,\ 12\le n\le 13d=6, 12≤n≤13; d=7, n=14d=7,\ n=14d=7, n=14, the bound ⌊(n−d)log⁡d⌋\lfloor (n-d)^{\log d}\rfloor⌊(n−d)logd⌋ satisfies the recursion of Lemma 1.
  10. For d≥4d\ge4d≥4, n−d≥8n-d\ge8n−d≥8, n≥2dn\ge 2dn≥2d: (d−1)log⁡(n−d)+2(n/2−d)log⁡d+2≤dlog⁡(n−d)(d-1)^{\log(n-d)}+2(n/2-d)^{\log d}+2\le d^{\log(n-d)}(d−1)log(n−d)+2(n/2−d)logd+2≤dlog(n−d). For n>2dn>2dn>2d this is the displayed inequality; at n=2dn=2dn=2d the second power is zero, consistent with Δ(d,d)=0\Delta(d,d)=0Δ(d,d)=0.

Significance

The result. The bound (n−d)log⁡d(n-d)^{\log d}(n−d)logd is never larger than the Kalai–Kleitman bound. It equals the true value n−dn-dn−d in dimensions 1 and 2, and it depends on nnn only through the Hirsch quantity n−dn-dn−d. Since (n−d)log⁡d=dlog⁡(n−d)(n-d)^{\log d}=d^{\log(n-d)}(n−d)logd=dlog(n−d), it is also symmetric under linear-programming duality, which exchanges the roles of ddd and n−dn-dn−d. For fixed n−dn-dn−d it is polynomial in ddd. It applies to the unbounded feasible regions that arise in practice in linear programming.

Formalizing it. The theorem has been proved since 2014. The platform has a proof of the bounded case only: Hirsch.todd_bound assumes d≥3d\ge3d≥3, a nonempty interior and Bornology.IsBounded. This mission's goal extends it to all polyhedra, bounded or not, and to d=1,2d=1,2d=1,2. That is the statement the paper proves, and the one an analysis of the simplex method on unbounded feasible regions needs. The unbounded form of Lemma 1 and the common-facet bound are open on the platform as items of the Kalai–Kleitman series and are reused here. Klee's unbounded three-dimensional bound has no formal proof on the platform: the published Hirsch.klee_three_dimensional_bound and Hirsch.dimension_three_bound cover bounded polyhedra only. The bounded recursion Hirsch.kalai_kleitman_recursion and the bounded facet reduction Hirsch.facet_reduction are proved and may serve as models.

Difficulty

Two steps are hard. The first is the unbounded case. Statements proved for polytopes do not transfer directly, because an unbounded polyhedron may have no vertices or a single one, unbounded edges (rays) that are not edges of the graph, and redundant or implicit-equality rows. Each geometric milestone must hold for these. The second is Klee's three-dimensional bound, which the paper cites without proof. Assuming boundedness to apply the published polytope versions is ruled out: Klee and Walkup showed that unbounded and bounded diameters differ from dimension four on, so the unbounded case carries content of its own. The arithmetic milestones are elementary but sensitive to edge values of log⁡\loglog and the real power.

Formalization scope

All statements use the published definitions Hirsch.Hpoly, Hirsch.Adj and Hirsch.DiamLE (Hirsch_model). Vertices are Set.extremePoints ℝ. DiamLE P B says that any two extreme points are joined by a walk of exactly BBB steps, each step staying put or crossing an edge. So the bound includes connectivity of the vertex–edge graph.

  • Inequality reading. "Δ(d,n)≤B\Delta(d,n)\le BΔ(d,n)≤B" is stated as: for every a:Fin n→Rda:\mathrm{Fin}\,n\to\mathbb R^da:Finn→Rd and bbb, DiamLE (Hpoly a b) B. This ranges over all polyhedra given by nnn inequalities in Rd\mathbb R^dRd, including lower-dimensional, empty and redundant ones. It is equivalent to the paper's statement: a lower-dimensional polyhedron of dimension d′d'd′ with n′n'n′ facets has n′−d′≤n−dn'-d'\le n-dn′−d′≤n−d.
  • No boundedness, nonemptiness or full-dimensionality is assumed anywhere.
  • The logarithm is Real.logb 2 and the power is Real.rpow. A distance is at most (n−d)log⁡d(n-d)^{\log d}(n−d)logd exactly when it is at most the floor ⌊(n−d)log⁡d⌋\lfloor (n-d)^{\log d}\rfloor⌊(n−d)logd⌋, abbreviated toddBound d n in the milestone on the ten cases. The clause Δ(1,1)=0\Delta(1,1)=0Δ(1,1)=0 is a separate conjunct, because 00=10^0=100=1 in Real.rpow.
  • Disclosed additions: nonzero rows in the n<2dn<2dn<2d step (as in the referenced Lemma 1 and the d=2d=2d=2 bound). In the general step, n≥2dn\ge2dn≥2d: the paper secures it by treating n<2dn<2dn<2d first. The second power is written with base n/2−dn/2-dn/2−d; it equals the printed expression when n>2dn>2dn>2d and takes the intended value zero when n=2dn=2dn=2d, where the printed logarithm is undefined.
  • Natural-number subtraction appears in n−3n-3n−3 and n−2n-2n−2 (both harmless, since for small nnn there are no vertices) and in Rd−1\mathbb R^{d-1}Rd−1 for facets.

Trivializing encodings are excluded: no Bornology.IsBounded hypothesis; no definition of Δ\DeltaΔ by sSup over N\mathbb NN, which would be 000 on unbounded families; no use of Hirsch.gdist, which is 000 when no walk exists; and no dropping of the Δ(1,1)=0\Delta(1,1)=0Δ(1,1)=0 clause.

Contributions are welcome on every milestone. A complete development needs facial structure of unbounded polyhedra (facets as polyhedra one dimension lower, tight rows at vertices), the relaxation argument of Lemma 1, Klee's three-dimensional argument, and real-power and logarithm estimates. The facet machinery and Klee's bound are reusable for any diameter result on unbounded polyhedra.

Selected references

  • M. J. Todd, An improved Kalai–Kleitman bound for the diameter of a polyhedron, SIAM J. Discrete Math. 28 (2014); arXiv:1402.3579v2. arXiv:1402.3579, doi:10.1137/140962310
  • G. Kalai and D. J. Kleitman, A quasi-polynomial bound for the diameter of graphs of polyhedra, Bull. Amer. Math. Soc. 26 (1992), 315–316. doi:10.1090/S0273-0979-1992-00285-9
  • V. Klee, Diameters of polyhedral graphs, Canad. J. Math. 16 (1964), 602–614. doi:10.4153/CJM-1964-058-3
  • V. Klee, Paths on polyhedra I, J. SIAM 13 (1965), 946–956. doi:10.1137/0113061
  • V. Klee, Paths on polyhedra II, Pacific J. Math. 17 (1966), 249–262. doi:10.2140/pjm.1966.17.249
  • V. Klee and D. W. Walkup, The d-step conjecture for polyhedra of dimension d < 6, Acta Math. 117 (1967), 53–78. doi:10.1007/BF02395039
  • F. Eisenbrand, N. Hähnle, A. Razborov, T. Rothvoss, Diameter of polyhedra: limits of abstraction, Math. Oper. Res. 35 (2010), 786–794. doi:10.1287/moor.1100.0470
  • F. Santos, A counterexample to the Hirsch conjecture, Ann. of Math. 176 (2012), 383–412. doi:10.4007/annals.2012.176.1.7
13 thms2 active usersReviewed
🏆Completed
ProbabilityStochastic Systems·Captain: mikedeng1

Separation of Time-Scales and Model Reduction for Stochastic Reaction Networks 2: Positive Combinations Preserve the Balance Condition When Production and Consumption Sets Are DisjointResearch Paper

Multiscale reaction networks and the balance condition

Stochastic models of chemical reaction networks in cell biology routinely combine species present in a handful of copies with species present in thousands, and rate constants that differ by orders of magnitude. Simulating such a system exactly is expensive, and model reduction replaces it with a lower-dimensional limit. A systematic approach, developed by Ball, Kurtz, Popovic and Rempala (2006) and by Kang and Kurtz (arXiv:1011.1672), normalizes species numbers and rate constants by powers of a large parameter NNN and derives limits as N→∞N\to\inftyN→∞ on several time-scales.

For such limits to exist, the normalized species numbers must stay of order one. For a single species this is the species balance condition: either the fastest reaction producing it is as fast as the fastest reaction consuming it, or the time-scale considered is short enough that the imbalance does not show. Kang and Kurtz observed (arXiv:1011.1672v1, §3.2, p. 9) that this is not enough: a collection of species can be produced collectively at a different order of magnitude from the rate at which it is consumed, even when every member is balanced. They therefore impose a balance condition on every nonnegative linear combination of species, which is infinitely many conditions. Section 3.4 of their paper proves combinatorial lemmas that reduce the work of checking them. This mission formalizes the lemmas on positive combinations of two weight vectors.

Setting

A reaction network has s0s_0s0​ species S1,…,Ss0S_1,\dots,S_{s_0}S1​,…,Ss0​​ and r0r_0r0​ reactions. Reaction kkk consumes νik∈N\nu_{ik}\in\mathbb Nνik​∈N molecules of SiS_iSi​ and produces νik′∈N\nu'_{ik}\in\mathbb Nνik′​∈N; its net change is the vector ζk=νk′−νk\zeta_k=\nu'_k-\nu_kζk​=νk′​−νk​. Each species carries a scaling exponent αi≥0\alpha_i\ge0αi​≥0 and each reaction an exponent βk∈R\beta_k\in\mathbb Rβk​∈R; put

ρk=βk+νk⋅α=βk+∑iνikαi,\rho_k=\beta_k+\nu_k\cdot\alpha=\beta_k+\sum_i\nu_{ik}\alpha_i,ρk​=βk​+νk​⋅α=βk​+i∑​νik​αi​,

the order of magnitude (as a power of NNN) of the rate of reaction kkk.

A weight vector θ∈[0,∞)s0\theta\in[0,\infty)^{s_0}θ∈[0,∞)s0​ describes the linear combination θ⋅X=∑iθiXi\theta\cdot X=\sum_i\theta_iX_iθ⋅X=∑i​θi​Xi​. Its sign sets (Definition 3.1) are

Γθ+={k:θ⋅ζk>0},Γθ−={k:θ⋅ζk<0},\Gamma^+_\theta=\{k:\theta\cdot\zeta_k>0\},\qquad\Gamma^-_\theta=\{k:\theta\cdot\zeta_k<0\},Γθ+​={k:θ⋅ζk​>0},Γθ−​={k:θ⋅ζk​<0},

the reactions that increase and decrease θ⋅X\theta\cdot Xθ⋅X. Maxima max⁡k∈Sρk\max_{k\in S}\rho_kmaxk∈S​ρk​ are taken in R∪{−∞}\mathbb R\cup\{-\infty\}R∪{−∞}, with max⁡k∈∅ρk=−∞\max_{k\in\emptyset}\rho_k=-\inftymaxk∈∅​ρk​=−∞, and a∨ba\vee ba∨b is the maximum.

Fix a time-scale exponent γ∈R\gamma\in\mathbb Rγ∈R. Condition 3.2 holds for θ\thetaθ if either the balance equation

max⁡k∈Γθ−ρk=max⁡k∈Γθ+ρk(3.7)\max_{k\in\Gamma^-_\theta}\rho_k=\max_{k\in\Gamma^+_\theta}\rho_k\tag{3.7}k∈Γθ−​max​ρk​=k∈Γθ+​max​ρk​(3.7)

or the time-scale constraint

γ≤γθ≡max⁡i:θi>0αi−max⁡k∈Γθ+∪Γθ−ρk(3.8)\gamma\le\gamma_\theta\equiv\max_{i:\theta_i>0}\alpha_i-\max_{k\in\Gamma^+_\theta\cup\Gamma^-_\theta}\rho_k\tag{3.8}γ≤γθ​≡i:θi​>0max​αi​−k∈Γθ+​∪Γθ−​max​ρk​(3.8)

holds, with γθ=+∞\gamma_\theta=+\inftyγθ​=+∞ when Γθ+∪Γθ−=∅\Gamma^+_\theta\cup\Gamma^-_\theta=\emptysetΓθ+​∪Γθ−​=∅.

Formalization targets

Goal: Lemma 3.8 (p. 18)

Let θ1,θ2∈[0,∞)s0\theta^1,\theta^2\in[0,\infty)^{s_0}θ1,θ2∈[0,∞)s0​, c1,c2>0c_1,c_2>0c1​,c2​>0, and suppose

Γθ1+∩Γθ2−=∅orΓθ1−∩Γθ2+=∅.\Gamma^+_{\theta^1}\cap\Gamma^-_{\theta^2}=\emptyset\quad\text{or}\quad\Gamma^-_{\theta^1}\cap\Gamma^+_{\theta^2}=\emptyset.Γθ1+​∩Γθ2−​=∅orΓθ1−​∩Γθ2+​=∅.

Then (i) if Condition 3.2 holds for θ1\theta^1θ1 and θ2\theta^2θ2, it holds for c1θ1+c2θ2c_1\theta^1+c_2\theta^2c1​θ1+c2​θ2; and (ii) if (3.7) holds for θ1\theta^1θ1 and θ2\theta^2θ2, it holds for c1θ1+c2θ2c_1\theta^1+c_2\theta^2c1​θ1+c2​θ2.

Milestones, in the order the paper uses them

  1. (3.16)–(3.17) (p. 15): Γc1θ1+c2θ2±⊆Γθ1±∪Γθ2±\Gamma^\pm_{c_1\theta^1+c_2\theta^2}\subseteq\Gamma^\pm_{\theta^1}\cup\Gamma^\pm_{\theta^2}Γc1​θ1+c2​θ2±​⊆Γθ1±​∪Γθ2±​, hence
max⁡k∈Γc1θ1+c2θ2±ρk≤max⁡k∈Γθ1±ρk∨max⁡k∈Γθ2±ρk.\max_{k\in\Gamma^\pm_{c_1\theta^1+c_2\theta^2}}\rho_k\le\max_{k\in\Gamma^\pm_{\theta^1}}\rho_k\vee\max_{k\in\Gamma^\pm_{\theta^2}}\rho_k.k∈Γc1​θ1+c2​θ2±​max​ρk​≤k∈Γθ1±​max​ρk​∨k∈Γθ2±​max​ρk​.
  1. Lemma 3.5 (p. 17): (3.8) for θ1,…,θm\theta^1,\dots,\theta^mθ1,…,θm implies (3.8) for ∑jcjθj\sum_jc_j\theta^j∑j​cj​θj, cj>0c_j>0cj​>0.
  2. Lemma 3.6 (p. 17): if max⁡Γθ1−ρ=max⁡Γθ1+ρ>max⁡Γθ2+∪Γθ2−ρ\max_{\Gamma^-_{\theta^1}}\rho=\max_{\Gamma^+_{\theta^1}}\rho>\max_{\Gamma^+_{\theta^2}\cup\Gamma^-_{\theta^2}}\rhomaxΓθ1−​​ρ=maxΓθ1+​​ρ>maxΓθ2+​∪Γθ2−​​ρ, then (3.7) holds for c1θ1+c2θ2c_1\theta^1+c_2\theta^2c1​θ1+c2​θ2.
  3. Lemma 3.7 (p. 18): (3.7) for θ1\theta^1θ1 and (3.8) for θ2\theta^2θ2 give Condition 3.2 for c1θ1+c2θ2c_1\theta^1+c_2\theta^2c1​θ1+c2​θ2.
  4. (3.35)–(3.36) (p. 19): if Γθ1−∩Γθ2+=∅\Gamma^-_{\theta^1}\cap\Gamma^+_{\theta^2}=\emptysetΓθ1−​∩Γθ2+​=∅, then
max⁡Γθ1−ρ≤max⁡Γc1θ1+c2θ2−ρ≤max⁡Γθ1−ρ∨max⁡Γθ2−ρ,max⁡Γθ2+ρ≤max⁡Γc1θ1+c2θ2+ρ≤max⁡Γθ1+ρ∨max⁡Γθ2+ρ.\max_{\Gamma^-_{\theta^1}}\rho\le\max_{\Gamma^-_{c_1\theta^1+c_2\theta^2}}\rho\le\max_{\Gamma^-_{\theta^1}}\rho\vee\max_{\Gamma^-_{\theta^2}}\rho,\quad\max_{\Gamma^+_{\theta^2}}\rho\le\max_{\Gamma^+_{c_1\theta^1+c_2\theta^2}}\rho\le\max_{\Gamma^+_{\theta^1}}\rho\vee\max_{\Gamma^+_{\theta^2}}\rho.Γθ1−​max​ρ≤Γc1​θ1+c2​θ2−​max​ρ≤Γθ1−​max​ρ∨Γθ2−​max​ρ,Γθ2+​max​ρ≤Γc1​θ1+c2​θ2+​max​ρ≤Γθ1+​max​ρ∨Γθ2+​max​ρ.

Companion results

  • Remark 3.9 (p. 18): if no reaction that consumes a species in the support of θ1\theta^1θ1 produces a species in the support of θ2\theta^2θ2, then Γθ1−∩Γθ2+=∅\Gamma^-_{\theta^1}\cap\Gamma^+_{\theta^2}=\emptysetΓθ1−​∩Γθ2+​=∅.
  • Lemma 3.10 (p. 19): if (3.7) holds for θ1\theta^1θ1, θ2\theta^2θ2 and for θ1−θ1⋅ζkθ2⋅ζkθ2\theta^1-\frac{\theta^1\cdot\zeta_k}{\theta^2\cdot\zeta_k}\theta^2θ1−θ2⋅ζk​θ1⋅ζk​​θ2 for every k∈(Γθ1+∩Γθ2−)∪(Γθ1−∩Γθ2+)k\in(\Gamma^+_{\theta^1}\cap\Gamma^-_{\theta^2})\cup(\Gamma^-_{\theta^1}\cap\Gamma^+_{\theta^2})k∈(Γθ1+​∩Γθ2−​)∪(Γθ1−​∩Γθ2+​), then (3.7) holds for c1θ1+c2θ2c_1\theta^1+c_2\theta^2c1​θ1+c2​θ2 for all c1,c2>0c_1,c_2>0c1​,c2​>0.

Significance

Condition 3.2 must hold for every θ∈[0,∞)s0\theta\in[0,\infty)^{s_0}θ∈[0,∞)s0​ before the limit theorems of the paper's Sections 4–5 apply. Lemma 3.8, together with the decomposition over strongly connected components (Lemma 3.4, the subject of the first mission of this series), turns that infinite family into finitely many checks for many networks: once the condition is known for some generating vectors, it propagates to their positive combinations whenever the sign sets do not conflict. Remark 3.9 gives a condition on the stoichiometry alone under which they do not conflict, and Lemma 3.10 reduces checking balance along a ray of combinations to finitely many breakpoints.

The results are proved in the paper (arXiv:1011.1672v1, pp. 15–19; published in Ann. Appl. Probab. 23(2), 2013), but none is formalized anywhere to our knowledge. The printed proof of Lemma 3.10 is a sketch that treats only the first breakpoint of the ray θ1+cθ2\theta^1+c\theta^2θ1+cθ2; a machine-checked proof would settle it in full. The definitions of the sign sets, the balance equation and the time-scale constraint are reusable for the other missions of this series.

Difficulty

The proofs are elementary, and the work lies in the bookkeeping. Every statement mixes three conventions: maxima over sets that can be empty (value −∞-\infty−∞), a time-scale constant that can be +∞+\infty+∞, and strict versus weak sign conditions on θ⋅ζk\theta\cdot\zeta_kθ⋅ζk​. Lemma 3.8 is a case analysis on which disjunct of Condition 3.2 holds for each of θ1\theta^1θ1 and θ2\theta^2θ2; the obvious idea of applying (3.16)–(3.17) alone fails in the case where both vectors are balanced at the same level ((3.34) in the paper), where only upper bounds are available. That case needs the disjointness hypothesis to get the matching lower bounds (3.35)–(3.36). Lemma 3.10 needs an argument along the whole ray θ1+cθ2\theta^1+c\theta^2θ1+cθ2, c>0c>0c>0: the sign pattern of (θ1+cθ2)⋅ζk(\theta^1+c\theta^2)\cdot\zeta_k(θ1+cθ2)⋅ζk​ changes at finitely many values of ccc, and balance has to be tracked through each change, not only the first.

Formalization scope

Species are Fin s and reactions Fin r (0-based indices; s0,r0s_0,r_0s0​,r0​ are arbitrary naturals). ν k i is νik\nu_{ik}νik​; ζik\zeta_{ik}ζik​, ρk\rho_kρk​ and θ⋅ζk\theta\cdot\zeta_kθ⋅ζk​ are real numbers, and ρ\rhoρ is defined from ν,α,β\nu,\alpha,\betaν,α,β rather than taken as a free parameter. A combination c1θ1+c2θ2c_1\theta^1+c_2\theta^2c1​θ1+c2​θ2 is c₁ • θ₁ + c₂ • θ₂ on Fin s → ℝ. Maxima are Finset.sup into WithBot ℝ, so the maximum over the empty set is ⊥=−∞\bot=-\infty⊥=−∞ and ∨\vee∨ is ⊔. γθ\gamma_\thetaγθ​ is an EReal defined by cases: +∞+\infty+∞ if Γθ+∪Γθ−=∅\Gamma^+_\theta\cup\Gamma^-_\theta=\emptysetΓθ+​∪Γθ−​=∅, otherwise the difference of the two maxima computed in R\mathbb RR (or −∞-\infty−∞ if no θi\theta_iθi​ is positive). Every statement carries the paper's standing assumptions α≥0\alpha\ge0α≥0, θ1,θ2≥0\theta^1,\theta^2\ge0θ1,θ2≥0 and c1,c2>0c_1,c_2>0c1​,c2​>0.

These choices rule out the trivializing encodings: a real-valued maximum with a default value of 000 on the empty set would make the balance equation hold spuriously; a subtraction in EReal would turn −∞−(−∞)-\infty-(-\infty)−∞−(−∞) into −∞-\infty−∞; and dropping θ≥0\theta\ge0θ≥0 or c1,c2>0c_1,c_2>0c1​,c2​>0 would change which reactions lie in the sign sets of the combination. In Lemma 3.8 the disjointness of the sign sets is a disjunction, as printed, and the second paragraph's "Condition 3.7" is read as the balance equation (3.7).

The development needs only finite sets, finite sums and order on R∪{±∞}\mathbb R\cup\{\pm\infty\}R∪{±∞}; no probability theory is involved. Lemmas about Finset.sup into WithBot ℝ over unions and subsets are reusable. Proofs of the milestones, of the goal, and of Lemma 3.10 are all welcome.

Selected references

  • H.-W. Kang and T. G. Kurtz, Separation of time-scales and model reduction for stochastic reaction networks, arXiv:1011.1672v1, 2010; Ann. Appl. Probab. 23(2), 2013. https://arxiv.org/abs/1011.1672, https://doi.org/10.1214/12-AAP841
  • K. Ball, T. G. Kurtz, L. Popovic and G. Rempala, Asymptotic analysis of multiscale approximations to reaction networks, Ann. Appl. Probab. 16(4):1925–1961, 2006. https://doi.org/10.1214/105051606000000420
  • Y. Cao, D. T. Gillespie and L. R. Petzold, The slow-scale stochastic simulation algorithm, J. Chem. Phys. 122(1):014116, 2005. https://doi.org/10.1063/1.1824902
8 thms2 active usersReviewed
PreviousPage 8 of 26Next

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