Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Collections

Discrete Convex Analysis

Murota's Discrete Convex Analysis, chapter by chapter: L-convex and M-convex functions, conjugacy, duality, and discrete separation.

24 completed missions

Missions

1–20 of 24
OpenCompletedAll
🏆Completed
CombinatoricsDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis I: Valuated MatroidsTextbook

Motivation

Matroids abstract the combinatorial content of linear independence: which sets of columns of a matrix are independent, which are maximal (bases), and how bases relate to each other. This abstraction, isolated independently by Whitney (1935) and van der Waerden's school, turned out to be exactly the right level of generality for a large family of greedy and augmenting-path algorithms — a base of a matroid can always be reached from another by a sequence of single-element swaps, and this exchange property is what makes local search on bases correct and efficient.

A natural question, raised in the 1980s once matroid-based combinatorial optimization was mature, is what happens when bases are not merely present or absent but carry real-valued weights that must interact well with the exchange structure. Dress and Wenzel answered this with the notion of a valuated matroid: a real-valued function on the bases of a matroid satisfying a weighted strengthening of the exchange axiom. Their motivation was explicitly algorithmic — valuated matroids are exactly the structures for which a greedy algorithm computes an optimal basis under linear objectives, and more generally under the family of "tilted" objectives obtained by adding an arbitrary linear functional. Independently, valuated matroids arise from the classical Grassmann–Plücker relation applied to matrices over a field with a valuation (hence the name), connecting them to tropical geometry.

This mission formalizes the two theorems of Murota's Discrete Convex Analysis (2003, §2.4) that make this story precise: the classical correspondence between a matroid's base family and its rank function (Theorem 2.29), and the characterization of valuations by a perturbation-robustness property (Theorem 2.32). Theorem 2.32 is also historically the entry point of the book's central theme — it is the special case, for the two-valued lattice {0,1}V\{0,1\}^V{0,1}V, of the general local-exchange criterion for M-convex functions that occupies chapters 6 and 7.

Setting

Let VVV be a finite set (the ground set). A matroid on VVV is a pair (V,B)(V, \mathcal B)(V,B) where B\mathcal BB, the base family, is a nonempty family of subsets of VVV satisfying the simultaneous exchange axiom (B): for every J,J′∈BJ, J' \in \mathcal BJ,J′∈B and every i∈J∖J′i \in J \setminus J'i∈J∖J′, there exists j∈J′∖Jj \in J' \setminus Jj∈J′∖J such that both

J−i+j:=(J∖{i})∪{j}∈BandJ′+i−j:=(J′∖{j})∪{i}∈B.J - i + j := (J \setminus \{i\}) \cup \{j\} \in \mathcal B \quad\text{and}\quad J' + i - j := (J' \setminus \{j\}) \cup \{i\} \in \mathcal B.J−i+j:=(J∖{i})∪{j}∈BandJ′+i−j:=(J′∖{j})∪{i}∈B.

Equivalently (Theorem 2.29 below), a matroid can be described by its rank function ρ:2V→Z\rho : 2^V \to \mathbb Zρ:2V→Z, a set function satisfying:

  • (R1) 0≤ρ(X)≤∣X∣0 \le \rho(X) \le |X|0≤ρ(X)≤∣X∣ for every X⊆VX \subseteq VX⊆V;
  • (R2) monotonicity: X⊆Y  ⟹  ρ(X)≤ρ(Y)X \subseteq Y \implies \rho(X) \le \rho(Y)X⊆Y⟹ρ(X)≤ρ(Y);
  • (R3) submodularity: ρ(X)+ρ(Y)≥ρ(X∪Y)+ρ(X∩Y)\rho(X) + \rho(Y) \ge \rho(X \cup Y) + \rho(X \cap Y)ρ(X)+ρ(Y)≥ρ(X∪Y)+ρ(X∩Y).

A valuation of a base family B\mathcal BB is a function ω:B→R\omega : \mathcal B \to \mathbb Rω:B→R satisfying the axiom (VM): for every J,J′∈BJ, J' \in \mathcal BJ,J′∈B and i∈J∖J′i \in J \setminus J'i∈J∖J′, there is j∈J′∖Jj \in J' \setminus Jj∈J′∖J with J−i+j,J′+i−j∈BJ - i + j, J' + i - j \in \mathcal BJ−i+j,J′+i−j∈B and

ω(J)+ω(J′)≤ω(J−i+j)+ω(J′+i−j).\omega(J) + \omega(J') \le \omega(J - i + j) + \omega(J' + i - j).ω(J)+ω(J′)≤ω(J−i+j)+ω(J′+i−j).

The pair (V,ω)(V, \omega)(V,ω) is then a valuated matroid. For p:V→Rp : V \to \mathbb Rp:V→R, the perturbation of ω\omegaω by ppp is

ω[−p](J)=ω(J)−∑j∈Jp(j).\omega[-p](J) = \omega(J) - \sum_{j \in J} p(j).ω[−p](J)=ω(J)−j∈J∑​p(j).

Formalization targets

Goal: Theorem 2.32 (the valuated matroid characterization)

ω is a valuation of B  ⟺  ∀ p:V→R, {J∈B:ω[−p](J′)≤ω[−p](J) ∀J′∈B} is a nonempty family satisfying (B).\omega \text{ is a valuation of } \mathcal B \iff \forall\, p : V \to \mathbb R,\ \{J \in \mathcal B : \omega[-p](J') \le \omega[-p](J)\ \forall J' \in \mathcal B\} \text{ is a nonempty family satisfying (B)}.ω is a valuation of B⟺∀p:V→R, {J∈B:ω[−p](J′)≤ω[−p](J) ∀J′∈B} is a nonempty family satisfying (B).

The right-hand side says: for every linear perturbation ppp, the set of ω[−p]\omega[-p]ω[−p]-maximal bases is again the base family of a matroid. The universal quantifier over ppp is not optional — a version of this statement quantified over a single fixed ppp is either vacuous or false, and does not capture what makes valuated matroids useful.

Milestone: Theorem 2.29 (the base-family / rank-function correspondence)

The maps

ρ(X)=max⁡{∣X∩J∣:J∈B},B={J⊆V:ρ(J)=∣J∣=ρ(V)}\rho(X) = \max\{|X \cap J| : J \in \mathcal B\}, \qquad \mathcal B = \{J \subseteq V : \rho(J) = |J| = \rho(V)\}ρ(X)=max{∣X∩J∣:J∈B},B={J⊆V:ρ(J)=∣J∣=ρ(V)}

are mutually inverse bijections between nonempty families satisfying (B) and set functions satisfying (R1)-(R3). This is weaker groundwork than the goal, stated first because it fixes the exact axiomatic vocabulary — (B) and (R) — that Theorem 2.32 is built on.

Significance

The result itself. Theorem 2.32 is the reason valuated matroids are the right object for weighted combinatorial optimization on matroids: it says a function on bases behaves correctly under every linear re-weighting of the ground set exactly when it satisfies the local exchange inequality (VM). This is what guarantees, for instance, that a greedy algorithm which is correct for the unweighted matroid extends correctly to families of tilted objectives, and it is the germ of the general local-optimality criterion for M-convex functions (chapters 6–7), which underlies most of the algorithmic content of the rest of the book. Theorem 2.29 is the classical result — due jointly to the development of matroid theory from the 1930s onward — that the base-exchange and rank-submodularity axiomatizations of a matroid carry the same information; it is the finite, unweighted precursor of Theorem 2.32.

Formalizing it. Neither theorem has a machine-checked proof on the platform prior to this mission (see Formalization scope for the prior-art check). Theorem 2.29's own proof is elementary but has two independent halves (each map preserves its target axiom class, and the two maps compose to the identity in both directions) that must all be established; Theorem 2.32's proof, as given in the source, defers entirely to a later, more general chapter-6 theorem, so a solver working only from this mission must either reconstruct a direct combinatorial argument for this special case or await chunk 06 (DiscreteConvex.MConvexFunctions, a separate mission) and specialize its main theorem.

Difficulty

The obvious approach to Theorem 2.32 — fix an optimal basis JJJ for ω[−p]\omega[-p]ω[−p] and try to show the exchange condition on maximizers directly from (VM) — proves one direction (VM implies the maximizer property) in a few lines, since perturbing does not change which exchange moves are available. The converse is the substantial direction: from "the maximizer set is always a matroid, for every ppp," one must recover the single global inequality (VM) that must hold for all pairs J,J′∈BJ, J' \in \mathcal BJ,J′∈B, not just optimal ones. The standard argument constructs, for a given non-optimal pair, a perturbation ppp under which that specific pair becomes simultaneously optimal, and this construction is exactly the step the book skips by citing chapter 6's general theorem. A formalization attempting to bypass this by only checking the maximizer property for a finite or generic sample of perturbations would trivialize the statement to something false or vacuous — a pitfall the goal's explicit ∀ p is designed to prevent.

Formalization scope

The ground set VVV is a Fintype with DecidableEq; 2V2^V2V is represented as Finset (Finset V), and V→RV \to \mathbb RV→R as a plain function type. The rank function is Z\mathbb ZZ-valued (matching the book's own convention for matroid rank, as opposed to the R\mathbb RR-valued conventions used from chapter 6 onward for general M-convex functions); RankOfFamily is implemented with Finset.sup over N\mathbb NN rather than a partial max', so that it is a total function — its junk value at the empty family is never invoked, since every hypothesis in this mission supplies nonemptiness explicitly, matching the book's own phrasing.

A trivializing formalization of the goal is one that quantifies over a single fixed ppp, or allows B\mathcal BB to be empty; both are explicitly excluded by keeping B.Nonempty\mathcal B.\text{Nonempty}B.Nonempty a hypothesis and ppp universally quantified inside the theorem statement itself.

Checked against Mathlib (commit 0df444a360eaa60ab8c11dca51a86af692955474): Mathlib's Matroid structure is axiomatized via the single-element (asymmetric) exchange property, classically but not definitionally equivalent to Murota's simultaneous axiom (B) used throughout this book, and Mathlib provides no constructor recovering a base family or a Matroid from a bare rank function satisfying (R1)-(R3). Theorem 2.29 is therefore genuine, reusable infrastructure, not a restatement of existing Mathlib API. No reference item was found on the platform for either theorem (GET /theorems?q=matroid, q=valuated matroid return only unrelated tropical-geometry and k-server results). Contributions to a shared DiscreteConvex.Combinatorial definitions layer (the exchange and rank axioms) are welcome from later chunks of this series that build on matroid or base-polyhedron structure.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • H. Whitney, "On the abstract properties of linear dependence," American Journal of Mathematics, 57(3), 1935, pp. 509–533.
  • A. W. M. Dress, W. Wenzel, "Valuated matroids," Advances in Mathematics, 93(2), 1992, pp. 214–250.
  • R. A. Brualdi, "Comments on bases in dependence structures," Bulletin of the Australian Mathematical Society, 1(2), 1969, pp. 161–167.
13 thms4 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis XV: Conjugacy of Quadratic Forms and Symmetric M-MatricesTextbook

Motivation

Quadratic minimization problems with a combinatorial sign pattern in their Hessian arise throughout applied mathematics: discretizations of elliptic boundary-value problems such as the Poisson equation, resistor-network energy functionals, and the Dirichlet forms of Markov-process potential theory all produce a symmetric matrix whose off-diagonal entries are nonpositive and whose rows are diagonally dominant (Fukushima, Oshima, and Takeda, Dirichlet Forms and Symmetric Markov Processes, De Gruyter, 1994). Such matrices are exactly the diagonally dominant symmetric M-matrices of classical numerical linear algebra (Berman and Plemmons, Nonnegative Matrices in the Mathematical Sciences, SIAM, 1994). Murota's Discrete Convex Analysis (SIAM, 2003) identifies the combinatorial content of this sign pattern with a discrete convexity property — submodularity, and its strengthening translation submodularity — of the associated quadratic form, and shows that passing to the Legendre-Fenchel conjugate of such a quadratic form (i.e., inverting the matrix) transports this property to a dual combinatorial property, an exchange axiom, on the conjugate side. This mission formalizes that correspondence for the special, matrix-algebraic case of quadratic forms — the case in which Murota's book gives a self-contained proof using only the classical Farkas lemma, before generalizing the same conjugacy to a much broader class of functions in Chapter 8.

Setting

Let VVV be a finite ground set (identified with {1,…,n}\{1,\dots,n\}{1,…,n} in the book) and let L=(ℓij)i,j∈VL = (\ell_{ij})_{i,j\in V}L=(ℓij​)i,j∈V​ be a symmetric real matrix. LLL has off-diagonal nonpositivity if ℓij≤0\ell_{ij}\le 0ℓij​≤0 for all i≠ji\ne ji=j, and diagonal dominance if ∑jℓij≥0\sum_{j} \ell_{ij}\ge 0∑j​ℓij​≥0 for every row iii. The associated quadratic form is g(p)=12p⊤Lpg(p) = \tfrac12 p^\top L pg(p)=21​p⊤Lp for p∈RVp \in \mathbb R^Vp∈RV. For p,q∈RVp,q\in\mathbb R^Vp,q∈RV write p∨qp\vee qp∨q, p∧qp\wedge qp∧q for the componentwise maximum and minimum. A function g:RV→Rg:\mathbb R^V\to\mathbb Rg:RV→R is submodular if g(p)+g(q)≥g(p∨q)+g(p∧q)g(p)+g(q)\ge g(p\vee q)+g(p\wedge q)g(p)+g(q)≥g(p∨q)+g(p∧q) for all p,qp,qp,q, and has translation submodularity if the stronger inequality g(p)+g(q)≥g((p−α1)∨q)+g(p∧(q+α1))g(p)+g(q)\ge g((p-\alpha\mathbf 1)\vee q)+g(p\wedge(q+\alpha\mathbf 1))g(p)+g(q)≥g((p−α1)∨q)+g(p∧(q+α1)) holds for every α≥0\alpha \ge 0α≥0, where 1\mathbf 11 is the all-ones vector (ordinary submodularity is the case α=0\alpha=0α=0).

On the conjugate side, for x∈RVx\in\mathbb R^Vx∈RV write supp⁡+(x)={i:xi>0}\operatorname{supp}^+(x)=\{i : x_i>0\}supp+(x)={i:xi​>0}, supp⁡−(x)={i:xi<0}\operatorname{supp}^-(x)=\{i:x_i<0\}supp−(x)={i:xi​<0}, and let χi\chi_iχi​ denote the iii-th unit vector (χ0\chi_0χ0​ denotes the zero vector). A function f:RV→Rf:\mathbb R^V\to\mathbb Rf:RV→R has the exchange property if for all x,y∈RVx,y\in\mathbb R^Vx,y∈RV and i∈supp⁡+(x−y)i\in\operatorname{supp}^+(x-y)i∈supp+(x−y) there exist j∈supp⁡−(x−y)∪{0}j \in \operatorname{supp}^-(x-y)\cup\{0\}j∈supp−(x−y)∪{0} and α0>0\alpha_0>0α0​>0 such that f(x)+f(y)≥f(x−α(χi−χj))+f(y+α(χi−χj))f(x)+f(y)\ge f(x-\alpha(\chi_i-\chi_j))+f(y+\alpha(\chi_i-\chi_j))f(x)+f(y)≥f(x−α(χi​−χj​))+f(y+α(χi​−χj​)) for every α∈[0,α0]\alpha\in[0,\alpha_0]α∈[0,α0​]. The Legendre-Fenchel conjugate of fff is f∙(p)=sup⁡x{⟨p,x⟩−f(x)}f^\bullet(p) = \sup_x\{\langle p,x\rangle - f(x)\}f∙(p)=supx​{⟨p,x⟩−f(x)}; two functions g,fg,fg,f are conjugate to each other when g=f∙g=f^\bulletg=f∙ and f=g∙f=g^\bulletf=g∙. For positive-definite symmetric M,LM,LM,L, the quadratic forms f(x)=12x⊤Mxf(x)=\tfrac12x^\top Mxf(x)=21​x⊤Mx and g(p)=12p⊤Lpg(p)=\tfrac12p^\top Lpg(p)=21​p⊤Lp are conjugate to each other exactly when MMM and LLL are matrix inverses of one another.

Formalization targets

Goal (Theorem 2.11). For conjugate strictly convex quadratic forms ggg and fff as above,

g has translation submodularity  ⟺  f has the exchange property.g \text{ has translation submodularity} \iff f \text{ has the exchange property.}g has translation submodularity⟺f has the exchange property.

This is the mission's capstone: the statement leaves the correspondence at the level of the two named combinatorial properties, without hard-coding which of the two properties is verified in a given application, so it survives exactly as strongly as the underlying conjugacy fact does.

Supporting milestones, in the order the book develops them: Proposition 2.4 (off-diagonal nonpositivity plus diagonal dominance implies positive semidefiniteness); Proposition 2.6 (off-diagonal nonpositivity is equivalent to plain submodularity of ggg); Theorem 2.7 (the full sign pattern is equivalent to translation submodularity of ggg); Proposition 2.9 (conjugate quadratic forms correspond exactly to inverse matrix pairs); Theorem 2.12 (a nine-way equivalence, for a nonsingular symmetric MMM, among membership in the matrix class L−1\mathcal L^{-1}L−1, two sign-consistency inequalities on the columns of MMM together with their strict forms, two directional-derivative reformulations of the exchange property together with their strict forms, and the exchange property itself together with its strict form); Proposition 2.13 (the Farkas lemma in equality form, together with the strict variant valid for a nonsingular coefficient matrix); and Proposition 2.14 (the class L−1\mathcal L^{-1}L−1 is closed under taking principal submatrices).

Significance

The M-natural exchange property is the function-level analogue of the base-exchange axiom for matroids, and translation submodularity is the analogue, on the "primal" side, of ordinary submodularity for set functions; Chapter 2's quadratic-form case is the historical and pedagogical entry point for the general conjugacy Chapter 8 proves for the full M-convex/ L-convex function classes. Establishing it here, in the self-contained matrix-algebraic setting, isolates exactly which properties of a quadratic form are combinatorial (tied to the coordinate axes) rather than purely convex-analytic (rotation-invariant): submodularity and the exchange property are not preserved by an orthogonal change of variables, in contrast to ordinary convexity, which Proposition 2.4 shows the same sign pattern also implies.

Formalizing this mission produces the first Lean statement, in this project's namespace, of a genuine conjugacy theorem between a primal-side and a dual-side combinatorial convexity property for a concrete function class; nothing of this kind is yet proved (or, so far as the platform's own search shows, formalized at all) elsewhere on the platform. The nine-way equivalence of Theorem 2.12 is a substantial independent contribution beyond the goal itself, since it is what makes the goal's proof possible via elementary linear algebra rather than the general convex-analytic machinery Chapter 8 needs.

Difficulty

The naive approach to Theorem 2.11 tries to derive the exchange property for fff directly from the defining supremum in the conjugate relation f=g∙f = g^\bulletf=g∙, differentiating under the sup; this fails because the exchange property compares fff along a specific combinatorial direction χi−χj\chi_i - \chi_jχi​−χj​ tied to two coordinates, not along an arbitrary direction, and no naive first-order argument isolates the right pair (i,j)(i,j)(i,j) without already knowing the sign pattern of M=L−1M = L^{-1}M=L−1. The book's actual route is Theorem 2.12: it reduces the exchange property to a column-wise sign-consistency statement on MMM itself (conditions (b)/(c)) via the identity f′(x;d)=x⊤Mdf'(x;d) = x^\top Mdf′(x;d)=x⊤Md, and closes the loop back to membership in L−1\mathcal L^{-1}L−1 using the Farkas lemma applied to the linear system ML=IML = IML=I — a genuinely matrix-algebraic argument that does not generalize verbatim to non-quadratic M-/L-convex functions, which is exactly why Chapter 8 needs a different (convex-analytic) proof for the general case.

Formalization scope

Vectors and matrices are indexed by a general finite type V ([Fintype V] [DecidableEq V]) rather than a fixed Fin n, matching this project's convention elsewhere and letting Proposition 2.14's principal-submatrix statement reuse the class predicate at the restricted index type directly. Quadratic forms are real-valued ((V → ℝ) → ℝ, using Matrix.mulVec and dotProduct) since this chapter's functions are always finite everywhere; the Legendre-Fenchel conjugate is EReal-valued via sSup, since a supremum over an infinite domain need not be finite in general even though it is finite here. Every min(0, \dots)-based condition in Theorem 2.12 and the exchange axioms is unfolded as the logically equivalent disjunction over the finitely many terms achieving the minimum, rather than reified via Finset.inf/WithTop machinery — a faithful, checked-equivalent simplification, not a narrowing (see MODERATION_NOTES.md). "Nonsingular" is Matrix.det ≠ 0. No numeric constant needs instantiation anywhere in this mission. The formalization does not trivialize: the goal's exchange property is stated for the specific combinatorial direction χi−χj\chi_i - \chi_jχi​−χj​ with i∈supp⁡+(x−y)i\in \operatorname{supp}^+(x-y)i∈supp+(x−y), j∈supp⁡−(x−y)∪{0}j \in \operatorname{supp}^-(x-y)\cup\{0\}j∈supp−(x−y)∪{0} — not an arbitrary direction, which would reduce the exchange property to a restatement of ordinary convexity and discard the entire combinatorial content the mission is about.

Infrastructure needed: Matrix.PosDef/Matrix.PosSemidef/Matrix.IsSymm (present in Mathlib); everything else (submodularity, translation submodularity, the exchange axioms, the sign-consistency conditions) is defined fresh in DiscreteConvex.CombinatorialB. A solution to the goal will likely want Proposition 2.9, Theorem 2.12, and the Farkas lemma (Proposition 2.13) as lemmas; contributions completing any of the seven milestones independently, or supplying the Schur-complement induction behind Proposition 2.4, are welcome.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003, DOI 10.1137/1.9780898718508, Chapter 2.
  • A. Berman, R. J. Plemmons, Nonnegative Matrices in the Mathematical Sciences, SIAM, 1994.
  • M. Fukushima, Y. Oshima, M. Takeda, Dirichlet Forms and Symmetric Markov Processes, De Gruyter, 1994.
  • J. Farkas, Theorie der einfachen Ungleichungen, J. Reine Angew. Math. 124 (1902), 1–27.
28 thms2 active usersReviewed
🏆Completed
Convex OptimizationLinear OptimizationOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis II: Local Optimality for Integrally Convex FunctionsTextbook

Motivation

For a convex function on Rn\mathbb R^nRn, a point is a global minimizer as soon as it is a local minimizer — this is one of the earliest and most consequential facts of convex analysis, and it underlies why local-search and gradient methods can certify global optimality in convex programs. The discrete analogue is not automatic: a function on the integer lattice Zn\mathbb Z^nZn can be "locally optimal" with respect to any fixed finite neighborhood system and still fail to be a global minimizer, unless the function's discrete structure is compatible with that neighborhood in the right way. Identifying exactly which classes of lattice functions admit a local-to-global optimality principle, and with respect to which neighborhood, is one of the organizing questions of discrete convex analysis.

Integrally convex functions, introduced by Favati and Tardella (1990) and developed systematically by Murota, are the most general class of Zn\mathbb Z^nZn-valued functions for which such a principle holds. They are defined purely in terms of the classical convex closure of a real relaxation, which lets one import theorems from ordinary convex analysis, but the resulting notion of local optimality — checking only the 3n−13^n - 13n−1 neighbors obtained by independently nudging each coordinate by −1-1−1, 000, or +1+1+1 (excluding the trivial no-change case) — is a genuinely discrete, dimension-independent statement about functions whose domain can be arbitrarily large. Almost every discrete convex function class studied later in the book, including M-convex and L-convex functions, is a special case of integral convexity, and this mission's goal theorem is the direct ancestor of the optimality criteria (Theorems 6.26 and 7.14) that drive the algorithms in the rest of the book.

Setting

Let f:Zn→R∪{+∞}f : \mathbb Z^n \to \mathbb R \cup \{+\infty\}f:Zn→R∪{+∞} be a function with nonempty effective domain dom⁡Zf={x∈Zn:f(x)≠+∞}\operatorname{dom}_{\mathbb Z} f = \{x \in \mathbb Z^n : f(x) \ne +\infty\}domZ​f={x∈Zn:f(x)=+∞}. The convex closure of fff is

fˉ(x)=sup⁡p∈Rn, α∈R{⟨p,x⟩+α:⟨p,y⟩+α≤f(y) ∀y∈Zn}(x∈Rn),\bar f(x) = \sup_{p \in \mathbb R^n,\, \alpha \in \mathbb R} \{\langle p,x\rangle + \alpha : \langle p,y\rangle + \alpha \le f(y)\ \forall y \in \mathbb Z^n\} \qquad (x \in \mathbb R^n),fˉ​(x)=p∈Rn,α∈Rsup​{⟨p,x⟩+α:⟨p,y⟩+α≤f(y) ∀y∈Zn}(x∈Rn),

the pointwise supremum of every affine function minorizing fff on all of Zn\mathbb Z^nZn. If fˉ\bar ffˉ​ agrees with fff on integer points, fff is convex extensible. The integral neighborhood of x∈Rnx \in \mathbb R^nx∈Rn is

N(x)={y∈Zn:⌊xi⌋≤yi≤⌈xi⌉, 1≤i≤n},N(x) = \{y \in \mathbb Z^n : \lfloor x_i \rfloor \le y_i \le \lceil x_i \rceil,\ 1 \le i \le n\},N(x)={y∈Zn:⌊xi​⌋≤yi​≤⌈xi​⌉, 1≤i≤n},

and the local convex extension f~\tilde ff~​ relaxes fˉ\bar ffˉ​'s definition by requiring the affine minorant condition only on N(x)N(x)N(x) rather than on all of Zn\mathbb Z^nZn. Always f~≥fˉ\tilde f \ge \bar ff~​≥fˉ​ pointwise, and the two agree on Zn\mathbb Z^nZn. A function fff is integrally convex if f~=fˉ\tilde f = \bar ff~​=fˉ​ everywhere on Rn\mathbb R^nRn — equivalently, if f~\tilde ff~​ is a convex function on all of Rn\mathbb R^nRn (it is automatically convex on every unit cube [z,z+1]n[z, z+1]^n[z,z+1]n with z∈Znz \in \mathbb Z^nz∈Zn, but need not be convex globally without this extra condition).

A discrete set S⊆ZnS \subseteq \mathbb Z^nS⊆Zn is hole free if SSS equals the set of integer points in its own real convex hull, and arg⁡min⁡f[−p]\arg\min f[-p]argminf[−p] denotes the minimizer set, over Zn\mathbb Z^nZn, of the linearly perturbed function f[−p](x)=f(x)−⟨p,x⟩f[-p](x) = f(x) - \langle p,x\ranglef[−p](x)=f(x)−⟨p,x⟩.

Formalization targets

Goal: Theorem 3.21 (local optimality characterizes global optimality)

For integrally convex fff and x∈dom⁡Zfx \in \operatorname{dom}_{\mathbb Z} fx∈domZ​f:

f(x)≤f(y) (∀y∈Zn)  ⟺  f(x)≤f(x+χY−χZ) (∀ Y,Z⊆{1,…,n}),f(x) \le f(y)\ (\forall y \in \mathbb Z^n) \iff f(x) \le f(x + \chi_Y - \chi_Z)\ (\forall\, Y, Z \subseteq \{1,\dots,n\}),f(x)≤f(y) (∀y∈Zn)⟺f(x)≤f(x+χY​−χZ​) (∀Y,Z⊆{1,…,n}),

where χY∈{0,1}n\chi_Y \in \{0,1\}^nχY​∈{0,1}n is the indicator vector of YYY. The right-hand side is a check over at most 3n−13^n - 13n−1 points (each coordinate independently unchanged, incremented, or decremented), regardless of how large dom⁡Zf\operatorname{dom}_{\mathbb Z} fdomZ​f is; this uniform, dimension-only bound is the entire content of the theorem, and is the weakest correct formulation — restricting to a single (Y,Z)(Y,Z)(Y,Z) or letting the right-hand side range over all of Zn\mathbb Z^nZn would trivialize or falsify the equivalence.

Milestones: Propositions 3.18 and 3.19

Proposition 3.18: fff convex extensible   ⟹  \implies⟹ arg⁡min⁡f[−p]\arg\min f[-p]argminf[−p] hole free for every ppp (and conversely, when dom⁡Zf\operatorname{dom}_{\mathbb Z} fdomZ​f is bounded). Proposition 3.19: fff is integrally convex if and only if every restriction f[a,b]f_{[a,b]}f[a,b]​ to a finite integer interval is integrally convex — integral convexity is detectable by looking at bounded pieces of fff one at a time.

Significance

The result itself. Theorem 3.21 is what makes integrally convex functions tractable: without it, verifying global optimality on an infinite or exponentially large integer domain would require checking every point. The theorem reduces this to a check whose size depends only on the dimension nnn, not on the size of the domain, and it does so for the widest class of lattice functions for which such a reduction is possible — the class is defined precisely so that this property holds and no wider natural class enjoys it. Every specialized local-optimality theorem later in the book (for M-convex, M♮^\natural♮-convex, L-convex, and L♮^\natural♮-convex functions) restricts this same neighborhood-checking principle to a class where the local check can be made even smaller (a single-element exchange rather than a full sign pattern) precisely because those classes are integrally convex plus more.

Formalizing it. No matching item exists on the platform: a direct search for "integrally convex" returns no results, and the theorem's own proof leans on results (Theorem 1.1's local-to-global principle for ordinary convex functions on Rn\mathbb R^nRn, and an LP-duality-based alternate formula for f~\tilde ff~​) that are either classical convex analysis or belong to a different chapter of this same book. The remaining work is therefore to give a complete, correct account of the definitional chain — convex closure, local convex extension, integral convexity — in a form a solver can build a proof from directly, and to state the finite local-check equivalence itself exactly at the strength the book proves it, not a plausible-looking weakening of it.

Difficulty

The natural first attempt is to try to prove the "⇐\Leftarrow⇐" direction of Theorem 3.21 by a direct induction on the ℓ1\ell^1ℓ1-distance to a global minimizer, moving one coordinate at a time. This fails in general lattice functions (a function that is only "coordinatewise convex" can have strict local minima that are not global), and the theorem's actual proof instead routes through the real relaxation: it shows the neighborhood-check hypothesis forces xxx to be a local minimizer of the local convex extension f~\tilde ff~​ restricted to the unit ball around xxx, then invokes ordinary convex analysis (local minimality implies global minimality for a convex function on Rn\mathbb R^nRn) to conclude xxx globally minimizes fˉ\bar ffˉ​, and finally uses integral convexity (f~=fˉ\tilde f = \bar ff~​=fˉ​) to transfer this back to fff on Zn\mathbb Z^nZn. The identification of fff's local behavior with f~\tilde ff~​'s convexity on a single unit cube — rather than any coordinatewise or separable argument — is the step that makes the class of integrally convex functions exactly the right one for this theorem, and is where a naive combinatorial argument breaks down.

Formalization scope

The ground set is Zn\mathbb Z^nZn, represented as Fin n → ℤ; fff's codomain is WithTop ℝ (exactly R∪{+∞}\mathbb R \cup \{+\infty\}R∪{+∞}), while the convex closure fˉ\bar ffˉ​ and local convex extension f~\tilde ff~​ take values in EReal (exactly R∪{±∞}\mathbb R \cup \{\pm\infty\}R∪{±∞}, a complete lattice, so their defining suprema are total functions with no side conditions). A trivializing formalization of the goal would quantify the right-hand side over a single fixed (Y,Z)(Y,Z)(Y,Z) pair, or over all of Zn\mathbb Z^nZn instead of the sign-pattern neighbors; both are excluded by keeping Y,ZY, ZY,Z universally quantified Finset (Fin n) ranging over the full 3n3^n3n sign-pattern space (minus the trivial case, which the equivalence still holds through vacuously).

Checked against the platform (GET /theorems?q=integrally convex, 0 hits) and against Mathlib's Analysis/Convex/ for the classical facts this chapter's proof would eventually need (ordinary convex-function local-to-global optimality, LP duality): these are broadly available in Mathlib's convex-analysis library in some form, but none of them is imported here, since none appears in the statement of any item this mission drafts — they belong to a proof this pass does not attempt. Contributions to a shared DiscreteConvex.IntegralConvexity definitions layer are welcome from chunks 06–09, which specialize integral convexity to M-convex and L-convex functions and will need the same convex-closure/local-extension vocabulary.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • P. Favati, F. Tardella, "Convexity in nonlinear integer programming," Ricerca Operativa, 53, 1990, pp. 3–44.
14 thms2 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis XVIII: Integral Convexity of Minimizer SetsTextbook

Motivation

A classical convex function's global minimality is equivalent to its local minimality — the single fact that makes convex optimization tractable, since checking a small neighborhood suffices to certify a global guarantee. The discrete analogue is not automatic: a function on the integer lattice can fail to have any well-behaved "local" notion at all, and even when a discrete convexity-like property is imposed, the naive candidate (a function's values agreeing with its own convex-hull interpolation) does not by itself guarantee that local optimality implies global optimality. Murota's Discrete Convex Analysis (SIAM, 2003) isolates exactly the extra condition — integral convexity — that restores this guarantee, and shows it is general enough to contain every other discrete convexity notion the book studies (M-convex, L-convex, and their variants), making it the common ancestor of the book's entire hierarchy of classes. This mission completes the chapter's account of integral convexity: how it behaves under sums, restrictions, and linear perturbations, how it transfers between a function and its domain or minimizer sets, and a companion fact about hole-freeness under intersection and Minkowski sums that motivates why integral convexity, not mere hole-freeness, is the right notion to use.

Setting

For f:Zn→R∪{+∞}f : \mathbb Z^n \to \mathbb R \cup \{+\infty\}f:Zn→R∪{+∞} with nonempty effective domain dom⁡Zf\operatorname{dom}_{\mathbb Z} fdomZ​f, the convex closure is fˉ(x)=sup⁡p,α{⟨p,x⟩+α:⟨p,y⟩+α≤f(y) ∀y∈Zn}\bar f(x) = \sup_{p,\alpha} \{\langle p,x\rangle + \alpha : \langle p,y\rangle+\alpha \le f(y)\ \forall y \in \mathbb Z^n\}fˉ​(x)=supp,α​{⟨p,x⟩+α:⟨p,y⟩+α≤f(y) ∀y∈Zn}. The integral neighborhood of x∈Rnx \in \mathbb R^nx∈Rn is N(x)={y∈Zn:⌊xi⌋≤yi≤⌈xi⌉}N(x) = \{y \in \mathbb Z^n : \lfloor x_i \rfloor \le y_i \le \lceil x_i \rceil\}N(x)={y∈Zn:⌊xi​⌋≤yi​≤⌈xi​⌉}, and the local convex extension f~\tilde ff~​ replaces "for all y∈Zny \in \mathbb Z^ny∈Zn" in fˉ\bar ffˉ​'s definition with "for all y∈N(x)y \in N(x)y∈N(x)". A function is integrally convex if f~=fˉ\tilde f = \bar ff~​=fˉ​ everywhere on Rn\mathbb R^nRn. A set S⊆ZnS \subseteq \mathbb Z^nS⊆Zn is integrally convex if its indicator function is; it is hole free if S=Sˉ∩ZnS = \bar S \cap \mathbb Z^nS=Sˉ∩Zn, where Sˉ\bar SSˉ is the real convex hull of SSS. The discrete Minkowski sum is S1+S2={x1+x2:x1∈S1,x2∈S2}S_1+S_2 = \{x_1+x_2 : x_1 \in S_1, x_2 \in S_2\}S1​+S2​={x1​+x2​:x1​∈S1​,x2​∈S2​}. A function is separable convex if f(x)=∑ifi(x(i))f(x) = \sum_i f_i(x(i))f(x)=∑i​fi​(x(i)) for univariate functions fif_ifi​ satisfying the discrete convexity inequality fi(t−1)+fi(t+1)≥2fi(t)f_i(t-1)+f_i(t+1) \ge 2f_i(t)fi​(t−1)+fi​(t+1)≥2fi​(t). For p∈Rnp \in \mathbb R^np∈Rn, f[−p](x)=f(x)−⟨p,x⟩f[-p] (x) = f(x) - \langle p,x \ranglef[−p](x)=f(x)−⟨p,x⟩ and arg⁡min⁡f[−p]\arg\min f[-p]argminf[−p] is its minimizer set.

Formalization targets

Goal (Theorem 3.29). For fff with nonempty bounded effective domain,

f is integrally convex  ⟺  arg⁡min⁡f[−p] is an integrally convex set for every p∈Rn.f \text{ is integrally convex} \iff \arg\min f[-p] \text{ is an integrally convex set for every } p \in \mathbb R^n.f is integrally convex⟺argminf[−p] is an integrally convex set for every p∈Rn.

This leaves the characterization at the level of the two named properties (integral convexity of the function versus of every minimizer set), the strongest statement of this kind that holds without extra hypotheses beyond boundedness of the domain.

Supporting milestones. Proposition 3.17 (four basic containment/equality relations between hole-free sets' intersections, Minkowski sums, and their real closures); Proposition 3.22 (for a periodic integrally convex function, global optimality reduces to a one-sided local check); Proposition 3.24 (an integrally convex function plus a separable convex function is integrally convex); Proposition 3.25 (separable convex functions are integrally convex, and integral convexity survives linear perturbation); Proposition 3.26 (an integrally convex set is hole free); Proposition 3.28 (the effective domain and every minimizer set of an integrally convex function are integrally convex sets — the forward direction of the goal); and Proposition 3.30 (for integer-valued integrally convex functions, a finite infimum is always attained).

Significance

Theorem 3.29 turns a statement about a function on all of Rn\mathbb R^nRn (integral convexity, a condition on f~\tilde ff~​ and fˉ\bar ffˉ​ that is a priori about uncountably many points) into a statement about a countable family of discrete sets (the minimizer sets arg⁡min⁡f[−p]\arg\min f[-p]argminf[−p]), giving a genuinely different and often more tractable way to certify or refute integral convexity. Propositions 3.24–3.25 are the closure properties that make integral convexity useful in practice: without them, verifying integral convexity of a function built from simpler pieces (a sum with a separable cost, a linearly reweighted objective) would require re-deriving the property from scratch each time. Proposition 3.17, by contrast, is a cautionary result: Note 3.27 and Example 3.15 (the two hole-free sets whose Minkowski sum has a hole) show that hole-freeness alone does not inherit good behavior under set operations, which is exactly the gap integral convexity's stronger, locally-checkable condition is built to close — this mission's Proposition 3.17 documents the "obvious"/general-purpose relations that hold regardless, so that the reader can see precisely which inclusion is automatic and which requires more.

Difficulty

The naive approach to Theorem 3.29's converse direction (integral convexity of every minimizer set implies integral convexity of fff) tries to check f~(x)=fˉ(x)\tilde f(x) = \bar f(x)f~​(x)=fˉ​(x) directly at an arbitrary x∈dom⁡fx \in \operatorname{dom} fx∈domf; this is circular, since f~\tilde ff~​ and fˉ\bar ffˉ​ are themselves defined via suprema over affine minorants, not via minimizer sets. The book's actual proof instead sets up a primal-dual pair of linear programs whose optimal solutions witness fˉ(x)\bar f(x)fˉ​(x) and f~(x)\tilde f(x)f~​(x) respectively, uses LP duality's complementary slackness to show the dual optimal solution can be chosen supported inside N(x)N(x)N(x), and only then concludes f~(x)=fˉ(x)\tilde f(x) = \bar f(x)f~​(x)=fˉ​(x) — routing the entire argument through the integral convexity of the specific minimizer set S=arg⁡min⁡f[−p∗]S = \arg\min f[-p^*]S=argminf[−p∗] at the optimal dual price p∗p^*p∗. This is why Theorem 3.29's proof needs LP duality (Theorem 3.10, formalized in the previous mission in this series) as an ingredient, not just the closure-property machinery of Propositions 3.24–3.28.

Formalization scope

All apparatus (ConvexClosure, IntegralNeighborhood, LocalConvexExtension, IntegrallyConvex, ArgMinPerturbed, HoleFree, IntegrallyConvexSet, SeparableConvex, MinkowskiSumZ) is redeclared fresh in DiscreteConvex.IntegralConvexityC, mirroring chunk 03-integral-convexity's already-established constructions (Fin n-indexed, WithTop ℝ-valued functions, EReal-valued convex closures via sSup), since a draft mission cannot import another draft's definitions. IntegrallyConvexSet is defined via the book's own primary definition (indicator function integrally convex) rather than either of its two stated equivalent reformulations, since no result in this mission needs those forms as a named predicate. An integer-valued function (Proposition 3.30) is represented as Zⁿ → WithTop ℤ and cast to WithTop ℝ via a small casting map wherever the real-valued apparatus is needed — a faithful embedding. Boundedness of a discrete set is containment in a finite integer interval. The formalization does not trivialize: Theorem 3.29's hypothesis is exactly "nonempty bounded effective domain", not further restricted to, say, a fixed small dimension or a finite ground set with a fixed cardinality bound, and every milestone is stated at the same generality as Propositions 3.24–3.28 and 3.30 give it (arbitrary nnn, arbitrary integrally convex function). Infrastructure needed beyond Mathlib: all definitions are fresh; a contribution completing any milestone, or the LP-duality-based proof of Theorem 3.29's converse direction, would be a natural entry point, alongside chunk 03's Theorem 3.21 as background.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003, DOI 10.1137/1.9780898718508, Chapter 3.
  • K. Murota, A. Shioura, "M-convex function on generalized polymatroid," Mathematics of Operations Research 24 (1999), 95–105 (Lemma 6.13, cited for Proposition 3.30).
28 thms3 active usersReviewed
🏆Completed
CombinatoricsDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis III: Edmonds's Intersection TheoremTextbook

Motivation

Matroid intersection is one of the founding results of combinatorial optimization: given two matroids on a common ground set, the largest common independent set can be found in polynomial time, and its size equals the minimum of a natural upper bound ranging over all subsets — a min-max theorem in the spirit of König's theorem and Menger's theorem, but for a strictly richer combinatorial structure. Jack Edmonds proved this in 1970, and Jack Edmonds and Rick Giles's subsequent generalization to submodular flows, together with André Frank's discrete separation theorem for submodular and supermodular set functions (1982), placed matroid intersection inside a single unifying framework: submodular function duality. This framework explains, in one stroke, matroid intersection, the base-exchange structure of matroids, and a family of other combinatorial min-max theorems that had previously seemed unrelated.

Murota's Discrete Convex Analysis develops this framework as the theory of M-convex sets: sets of integer vectors satisfying a lattice-exchange axiom that turns out to be exactly equivalent to being the integer points of a base polyhedron of an integer-valued submodular set function. This mission formalizes the chapter's central results: the equivalence of four variant forms of the exchange axiom (Theorem 4.3), the M-convex set / submodular function correspondence (Theorem 4.15), Frank's discrete separation theorem (Theorem 4.17), and Edmonds's intersection theorem itself (Theorem 4.18) — the deepest duality result in the theory of submodular functions and the historical origin of the M-convexity concept that the rest of the book generalizes to real-valued functions.

Setting

Let VVV be a finite ground set. A set function ρ:2V→R∪{+∞}\rho : 2^V \to \mathbb R \cup \{+\infty\}ρ:2V→R∪{+∞} with ρ(∅)=0\rho(\emptyset) = 0ρ(∅)=0 and ρ(V)<+∞\rho(V) < +\inftyρ(V)<+∞ is submodular (the class S[R]S[\mathbb R]S[R]) if

ρ(X)+ρ(Y)≥ρ(X∪Y)+ρ(X∩Y)(X,Y⊆V).\rho(X) + \rho(Y) \ge \rho(X \cup Y) + \rho(X \cap Y) \qquad (X, Y \subseteq V).ρ(X)+ρ(Y)≥ρ(X∪Y)+ρ(X∩Y)(X,Y⊆V).

Its base polyhedron and submodular polyhedron are

B(ρ)={x∈RV:x(X)≤ρ(X) (∀X⊆V), x(V)=ρ(V)},P(ρ)={x∈RV:x(X)≤ρ(X) (∀X⊆V)},B(\rho) = \{x \in \mathbb R^V : x(X) \le \rho(X)\ (\forall X \subseteq V),\ x(V) = \rho(V)\}, \qquad P(\rho) = \{x \in \mathbb R^V : x(X) \le \rho(X)\ (\forall X \subseteq V)\},B(ρ)={x∈RV:x(X)≤ρ(X) (∀X⊆V), x(V)=ρ(V)},P(ρ)={x∈RV:x(X)≤ρ(X) (∀X⊆V)},

where x(X)=∑v∈Xx(v)x(X) = \sum_{v \in X} x(v)x(X)=∑v∈X​x(v); a supermodular function μ\muμ is one with −μ-\mu−μ submodular. A nonempty set B⊆ZVB \subseteq \mathbb Z^VB⊆ZV is an M-convex set if it satisfies the exchange axiom (B-EXC[Z]): for x,y∈Bx, y \in Bx,y∈B and uuu in the positive support of x−yx-yx−y, there is vvv in the negative support of x−yx-yx−y with both x−χu+χv∈Bx - \chi_u + \chi_v \in Bx−χu​+χv​∈B and y+χu−χv∈By + \chi_u - \chi_v \in By+χu​−χv​∈B, where χu\chi_uχu​ is the characteristic vector of uuu. A polyhedron P⊆RVP \subseteq \mathbb R^VP⊆RV is integral if P=conv⁡(P∩ZV)P = \operatorname{conv}(P \cap \mathbb Z^V)P=conv(P∩ZV).

Formalization targets

Goal: Theorem 4.18 (Edmonds's intersection theorem)

For submodular set functions ρ1,ρ2∈S[R]\rho_1, \rho_2 \in S[\mathbb R]ρ1​,ρ2​∈S[R],

max⁡{x(V):x∈P(ρ1)∩P(ρ2)}=min⁡{ρ1(X)+ρ2(V∖X):X⊆V},\max\{x(V) : x \in P(\rho_1) \cap P(\rho_2)\} = \min\{\rho_1(X) + \rho_2(V \setminus X) : X \subseteq V\},max{x(V):x∈P(ρ1​)∩P(ρ2​)}=min{ρ1​(X)+ρ2​(V∖X):X⊆V},

with both sides attained. If ρ1,ρ2\rho_1, \rho_2ρ1​,ρ2​ are integer valued, P(ρ1)∩P(ρ2)P(\rho_1) \cap P(\rho_2)P(ρ1​)∩P(ρ2​) is an integral polyhedron and the maximum is attained at an integer point. Dropping the integrality clause and stating only the real max-min equality would leave ordinary LP duality with no discrete content at all; this mission keeps it in the goal at every strength the book proves it.

Milestones: Theorems 4.3, 4.15, 4.17

Theorem 4.3: the exchange axiom (B-EXC[Z]) is equivalent to three variants that impose the exchange condition asymmetrically or only for distinct vectors — groundwork establishing that M-convexity does not depend on which variant is taken as primitive. Theorem 4.15: BBB is M-convex if and only if B=B(ρ)∩ZVB = B(\rho) \cap \mathbb Z^VB=B(ρ)∩ZV for some integer-valued submodular ρ\rhoρ — M-convex sets and integer-valued submodular set functions are two descriptions of the same combinatorial object. Theorem 4.17 (Frank): if a submodular ρ\rhoρ dominates a supermodular μ\muμ pointwise, a single vector x∗x^*x∗ separates them (ρ≥x∗≥μ\rho \ge x^* \ge \muρ≥x∗≥μ pointwise on every subset), integrally when ρ,μ\rho, \muρ,μ are integer valued — derived, in the book, as a direct corollary of the goal theorem.

Significance

The result itself. Edmonds's intersection theorem is the min-max theorem underlying polynomial-time matroid intersection (a matroid's rank function is submodular, so the classical matroid intersection theorem is the special case ρ1,ρ2\rho_1, \rho_2ρ1​,ρ2​ both matroid rank functions), and its generality — arbitrary submodular set functions, not just matroid ranks — is what lets Frank's discrete separation theorem, and through it a wide range of combinatorial duality results in network flows, scheduling, and matroid theory, be derived as corollaries rather than proved from scratch each time. The integrality clause specifically is the fact that makes these duality theorems combinatorial: it guarantees that optimal fractional solutions to the underlying linear program can always be taken integral, without which the connection to discrete optimization would be lost.

Formalizing it. No matching item exists on the platform (searches for "submodular set function", "base polyhedron", "matroid intersection" return no relevant hits; Mathlib's Combinatorics/Matroid/ develops matroid rank functions, a special case, but not general submodular set functions or their polyhedra). This mission gives the first formal statement of the theorem at its natural generality, together with the M-convex-set viewpoint that motivates the rest of the book, and Frank's separation theorem as an explicit worked corollary.

Difficulty

The real-valued half of Theorem 4.18 is ordinary LP duality applied to a cleverly chosen primal program (maximize ⟨p,x⟩\langle p, x\rangle⟨p,x⟩ over P(ρ1)∩P(ρ2)P(\rho_1) \cap P(\rho_2)P(ρ1​)∩P(ρ2​)) and its dual — routine once the right LP is written down. The integrality half is where the combinatorics enters: an optimal dual solution can always be chosen supported on a chain in each ρi\rho_iρi​'s effective domain (an extremal argument maximizing a strictly convex potential over the optimal dual face), and the incidence matrix of a chain of subsets is totally unimodular — this is the fact, external to ordinary LP theory, that forces an integral optimal solution to exist whenever the data (ρ1,ρ2\rho_1, \rho_2ρ1​,ρ2​) are integral. A proof that stops at real-valued LP duality, however carefully done, misses this step entirely and cannot produce the integrality clause; total unimodularity of a chain's incidence matrix is the one piece of combinatorics doing all the discrete work in an otherwise classical convex-duality argument.

Formalization scope

The ground set VVV is a Fintype with DecidableEq; subsets are Finset V, vectors are V → ℝ/V → ℤ. Submodular functions take values in WithTop ℝ (exactly R∪{+∞}\mathbb R \cup \{+\infty\}R∪{+∞}); supermodular functions in WithBot ℝ; comparisons across the two use an explicit embedding into EReal. The max/min in the goal are stated via IsGreatest/IsLeast sharing a common EReal witness, so that "both sides attained, at the same value" — not merely "sup equals inf" — is what the Lean statement asserts, which is essential since the integrality clause's whole content is about which point attains the maximum.

A trivializing formalization of the goal would drop the integrality clause (leaving unqualified LP duality) or replace IsGreatest/IsLeast with a bare supremum/infimum equality (losing the "is attained" content the second half of the theorem needs); both are avoided. Theorem 4.15 is stated as the existential "iff" (some integer submodular ρ\rhoρ realizes BBB) rather than reifying the book's own named bijection Φ,Ψ\Phi, \PsiΦ,Ψ explicitly — a deliberate, documented scope reduction of that one milestone (see MODERATION_NOTES.md), not of the goal. Contributions building the explicit Φ\PhiΦ map, the Lovász extension (needed for Theorem 4.16, not drafted here), or M-convex-set infrastructure reusable by chunks 06–07 (M-convex functions, which build on this chapter's vocabulary) are welcome.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • J. Edmonds, "Submodular functions, matroids, and certain polyhedra," in Combinatorial Structures and Their Applications, Gordon and Breach, 1970, pp. 69–87.
  • A. Frank, "An algorithm for submodular functions on graphs," Annals of Discrete Mathematics, 16, 1982, pp. 97–120.
20 thms2 active usersReviewed
🏆Completed
CombinatoricsDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis IV: Discrete Separation for L-Convex SetsTextbook

Motivation

The classical separating hyperplane theorem says that any two disjoint convex sets in Rn\mathbb R^nRn can be separated by a hyperplane with an arbitrary real normal vector. When the sets in question are not arbitrary convex sets but the integer points of specially structured discrete sets, one can sometimes ask for much more: not merely that a separator exists, but that it can be chosen from a small, structured, dimension-independent family regardless of the size or shape of the sets being separated. Results of this kind — "discrete separation theorems" — are a recurring and often surprising theme in combinatorial optimization, playing the role that the ordinary separation theorem plays in continuous convex analysis, but with genuinely combinatorial content beyond it.

L-convex sets, introduced by Murota as part of the discrete convex analysis framework, are one of the two dual families of well-behaved discrete convex sets studied in the book (the other being M-convex sets, chunk 04 of this series). They are defined by a lattice-closure axiom together with translation invariance, and they correspond one-to-one to integer-valued distance functions satisfying the triangle inequality — objects long familiar from network flow theory and shortest-path duality, even though the L-convexity terminology is not traditionally used there. This mission formalizes the chapter's central results, culminating in Theorem 5.9: two disjoint L-convex sets can always be separated by a vector with entries in {−1,0,1}\{-1, 0, 1\}{−1,0,1}, no matter how large or complicated the sets are.

Setting

Let VVV be a finite ground set. A nonempty set D⊆ZVD \subseteq \mathbb Z^VD⊆ZV is an L-convex set if it satisfies the sublattice axiom (SBS[Z]) — p,q∈D  ⟹  p∨q, p∧q∈Dp, q \in D \implies p \vee q,\ p \wedge q \in Dp,q∈D⟹p∨q, p∧q∈D, where ∨,∧\vee, \wedge∨,∧ are componentwise maximum and minimum — and the translation axiom (TRS[Z]) — p∈D  ⟹  p±1∈Dp \in D \implies p \pm \mathbf 1 \in Dp∈D⟹p±1∈D, where 1\mathbf 11 is the all-ones vector. A distance function γ:V×V→R∪{+∞}\gamma : V \times V \to \mathbb R \cup \{+\infty\}γ:V×V→R∪{+∞} satisfies γ(v,v)=0\gamma(v,v) = 0γ(v,v)=0 for every vvv; it satisfies the triangle inequality if γ(v1,v2)+γ(v2,v3)≥γ(v1,v3)\gamma(v_1,v_2) + \gamma(v_2,v_3) \ge \gamma(v_1,v_3)γ(v1​,v2​)+γ(v2​,v3​)≥γ(v1​,v3​) for all v1,v2,v3v_1, v_2, v_3v1​,v2​,v3​. The admissible-potential polyhedron of γ\gammaγ is

D(γ)={p∈RV:p(v)−p(u)≤γ(u,v) (∀u≠v)}.D(\gamma) = \{p \in \mathbb R^V : p(v) - p(u) \le \gamma(u,v)\ (\forall u \ne v)\}.D(γ)={p∈RV:p(v)−p(u)≤γ(u,v) (∀u=v)}.

The convex hull of a discrete set D⊆ZVD \subseteq \mathbb Z^VD⊆ZV is written Dˉ⊆RV\bar D \subseteq \mathbb R^VDˉ⊆RV.

Formalization targets

Goal: Theorem 5.9 (discrete separation for L-convex sets)

If D1,D2⊆ZVD_1, D_2 \subseteq \mathbb Z^VD1​,D2​⊆ZV are disjoint L-convex sets, there exists x∗∈{−1,0,1}Vx^* \in \{-1,0,1\}^Vx∗∈{−1,0,1}V such that

inf⁡{⟨p,x∗⟩:p∈D1}−sup⁡{⟨p,x∗⟩:p∈D2}≥1.\inf\{\langle p, x^*\rangle : p \in D_1\} - \sup\{\langle p, x^*\rangle : p \in D_2\} \ge 1.inf{⟨p,x∗⟩:p∈D1​}−sup{⟨p,x∗⟩:p∈D2​}≥1.

Dropping the {−1,0,1}V\{-1,0,1\}^V{−1,0,1}V restriction and allowing an arbitrary real separator would recover the classical separation theorem for convex sets, which holds regardless of L-convexity and carries no discrete-convexity content; the three-valued restriction is the weakest correct strengthening and is kept in full.

Milestones: Theorems 5.2, 5.5, 5.7

Theorem 5.2: an L-convex set is hole free (D=Dˉ∩ZVD = \bar D \cap \mathbb Z^VD=Dˉ∩ZV) — its integer points are exactly the integer points of its own convex hull. Theorem 5.5: DDD is L-convex if and only if D=D(γ)∩ZVD = D(\gamma) \cap \mathbb Z^VD=D(γ)∩ZV for some integer-valued distance function γ\gammaγ satisfying the triangle inequality — L-convex sets and such distance functions are two descriptions of the same object, the discrete analogue of chunk 04's M-convex-set / submodular- function correspondence. Theorem 5.7 (parts (1), (4)): L-convex sets are closed under intersection in the strongest sense — the convex hulls intersect exactly where the sets do, and a nonempty intersection of L-convex sets is again L-convex.

Significance

The result itself. Theorem 5.9 packs two claims into one, as the book itself points out: the separator is forced into {−1,0,1}V\{-1,0,1\}^V{−1,0,1}V (explicit in the statement), and disjoint L-convex sets satisfy "convexity in intersection" — their convex hulls are already disjoint whenever the sets themselves are (implicit, and necessary for the stated inequality to be possible at all). The {−1,0,1}\{-1,0,1\}{−1,0,1} structure connects directly to combinatorial duality in network flows: L-convex polyhedra are, without the name, a familiar object there, and a {−1,0,1}\{-1,0,1\}{−1,0,1}-separator corresponds to a signed cut or a negative-cost cycle in an associated graph. Theorem 5.5's correspondence is the L-convex mirror of chunk 04's M-convex/submodular correspondence, and the book explicitly flags that the two will be unified into a single conjugacy relationship in a later chapter (Note 5.6) — this mission's formalization of the L-side is a prerequisite for that later unification.

Formalizing it. No matching item exists on the platform (searches for "L-convex", "distance function", and "negative cycle" return only unrelated results — number-theoretic distance estimates, polytope graph metrics, shortest-path graph structures — none matching the combinatorial L-convexity/discrete-separation content here). This mission gives the first formal statement of L-convex sets and their central separation theorem. Notably, Theorem 5.9's own statement — unlike the analogous M-convex Theorem 4.18 — needs none of the distance-function machinery that its proof uses; only the L-convexity axiom itself appears in the goal, making its formal statement comparatively lean even though the underlying mathematics is just as deep.

Difficulty

The natural first attempt at Theorem 5.9 is to try to construct x∗x^*x∗ directly from the structure of D1,D2D_1, D_2D1​,D2​ — for instance, from a normal vector to a real separating hyperplane, rounded coordinatewise. This does not work: rounding an arbitrary real separator gives no control over its entries, and there is no reason a rounded vector should still separate. The book's actual proof instead represents D1,D2D_1, D_2D1​,D2​ via distance functions γ1,γ2\gamma_1, \gamma_2γ1​,γ2​ (Theorem 5.5), combines them into γ12=min⁡(γ1,γ2)\gamma_{12} = \min(\gamma_1, \gamma_2)γ12​=min(γ1​,γ2​), and extracts the separator from a shortest negative cycle in the associated graph: the vertices of the cycle alternate between the two sets' "tight" arcs, and the alternating ±1\pm 1±1 pattern around the cycle is exactly the {−1,0,1}\{-1,0,1\}{−1,0,1} vector x∗x^*x∗ — with the cycle's negativity translating directly into the required gap of at least 111. Locating the right combinatorial object (a shortest negative cycle, not an arbitrary one) is what pins the separator down to a vector supported on a single alternating cycle rather than an arbitrary {−1,0,1}\{-1,0,1\}{−1,0,1} pattern, and is the step a naive rounding or linear-algebra argument has no analogue of.

Formalization scope

The ground set VVV is a Fintype with DecidableEq; D⊆ZVD \subseteq \mathbb Z^VD⊆ZV is a Set (V → ℤ). Distance functions take values in WithTop ℝ; the goal's infimum and supremum are taken in EReal (a complete lattice), since L-convex sets are always infinite (translation invariance along the all-ones direction), so an ℝ-valued supremum/infimum would silently return a junk value on an unbounded set. The conclusion is stated as sup⁡D2⟨p,x∗⟩+1≤inf⁡D1⟨p,x∗⟩\sup_{D_2}\langle p,x^*\rangle + 1 \le \inf_{D_1}\langle p,x^*\ranglesupD2​​⟨p,x∗⟩+1≤infD1​​⟨p,x∗⟩, an addition-based reformulation of the book's subtraction inequality that avoids EReal's ⊤ - ⊤ ambiguity while remaining equivalent whenever both sides are finite.

A trivializing formalization of the goal would drop the {−1,0,1}V\{-1,0,1\}^V{−1,0,1}V constraint on x∗x^*x∗ (recovering the classical, L-convexity-independent separation theorem) or fix a single coordinate pattern rather than asserting existence over the full three-valued family; neither is done here. Theorem 5.5 is stated existentially rather than via the book's named bijection Φ,Ψ\Phi, \PsiΦ,Ψ (a documented scope reduction, parallel to chunk 04's treatment of Theorem 4.15), and Theorem 5.7 is drafted with only its two representation-independent clauses (parts (1) and (4); see MODERATION_NOTES.md). Contributions building the distance-function/admissible- potential apparatus needed for Theorem 5.7's remaining clauses, or the L-convex/integrally-convex bridge (Theorem 5.10, needing chunk 03's vocabulary), are welcome.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
12 thms2 active usersReviewed
🏆Completed
CombinatoricsConvex OptimizationDiscrete Geometry+2·Captain: Shuze Chen

Discrete Convex Analysis XX: Integral Convexity of L-Convex SetsTextbook

Motivation

Shortest-path distances and network potentials are among the oldest objects in combinatorial optimization: a directed graph with arc lengths, its shortest-path distances, and the "feasible potentials" (vertex labels consistent with those lengths) underlie duality in min-cost flow, scheduling, and difference-constraint systems. Murota's Discrete Convex Analysis (SIAM, 2003) isolates the abstract structure behind these objects — distance functions satisfying the triangle inequality, and their associated sets of admissible potentials — and shows it is governed by exactly the same discrete-convexity machinery as submodular set functions: a one-to-one correspondence with a second family of well-behaved integer point sets, the L-convex sets. Where an M-convex set (chapter 4) is defined by an exchange axiom generalizing matroid base exchange, an L-convex set is defined by closure under coordinatewise lattice operations (∨, ∧) and translation by the all-ones vector — a genuinely different axiom system that nonetheless produces a parallel structural theory: hole-freeness, a polyhedral description via an induced distance function, and integral convexity.

Companion mission 05-lconvex-sets (Discrete Convex Analysis IV) covers this chapter's other half: the hole-free property (Theorem 5.2), the one-to-one correspondence between L-convex sets and integer-valued triangle-inequality distance functions (Theorem 5.5), the intersection properties (Theorem 5.7), and the chapter's discrete separation theorem (Theorem 5.9, its goal). This mission builds the vocabulary those results also need (redeclared here, since sibling drafts cannot yet import one another) and proves the results the chapter leaves for its second half: the fundamental facts connecting a distance function to its admissible potentials (Proposition 5.1), the two-way polyhedral correspondence's supporting propositions (5.3-5.4), Minkowski-sum convexity (Theorem 5.8), and — this mission's goal — the explicit description of an L-convex set's convex hull that establishes its integral convexity (Theorem 5.10).

Setting

Fix a finite ground set VVV. A distance function is a map γ:V×V→R∪{+∞}\gamma : V \times V \to \mathbb R \cup \{+\infty\}γ:V×V→R∪{+∞} with γ(v,v)=0\gamma(v,v) = 0γ(v,v)=0; it may take negative finite values and need not be symmetric. It defines a directed graph Gγ=(V,Aγ)G_\gamma = (V, A_\gamma)Gγ​=(V,Aγ​) with Aγ={(u,v):γ(u,v)<+∞}A_\gamma = \{(u,v) : \gamma(u,v) < +\infty\}Aγ​={(u,v):γ(u,v)<+∞}, arc (u,v)(u,v)(u,v) having length γ(u,v)\gamma(u,v)γ(u,v). Write γˉ(u,v)\bar\gamma(u,v)γˉ​(u,v) for the shortest-path length from uuu to vvv in GγG_\gammaGγ​ (+∞+\infty+∞ if none exists); γ\gammaγ is well defined (γˉ\bar\gammaγˉ​ finite-valued wherever a path exists) exactly when GγG_\gammaGγ​ has no negative cycle. The triangle inequality γ(v1,v2)+γ(v2,v3)≥γ(v1,v3)\gamma(v_1,v_2) + \gamma(v_2,v_3) \ge \gamma(v_1,v_3)γ(v1​,v2​)+γ(v2​,v3​)≥γ(v1​,v3​) defines the class T[R]T[\mathbb R]T[R] (or T[Z]T[\mathbb Z]T[Z] when integer-valued). A vector p∈RVp \in \mathbb R^Vp∈RV is an admissible potential of γ\gammaγ if p(v)−p(u)≤γ(u,v)p(v) - p(u) \le \gamma(u,v)p(v)−p(u)≤γ(u,v) for all u≠vu \ne vu=v; write D(γ)D(\gamma)D(γ) for the set of all such potentials.

A nonempty set D⊆ZVD \subseteq \mathbb Z^VD⊆ZV is L-convex if it satisfies (SBS[Z]): p,q∈D  ⟹  p∨q, p∧q∈Dp, q \in D \implies p \vee q,\ p \wedge q \in Dp,q∈D⟹p∨q, p∧q∈D (coordinatewise max/min), and (TRS[Z]): p∈D  ⟹  p±1∈Dp \in D \implies p \pm \mathbf 1 \in Dp∈D⟹p±1∈D. A set S⊆ZVS \subseteq \mathbb Z^VS⊆ZV is integrally convex if every point of its convex hull S‾\overline SS lies in the convex hull of SSS restricted to that point's integral neighborhood N(p)={y∈ZV:⌊p⌋≤y≤⌈p⌉ coordinatewise}N(p) = \{y \in \mathbb Z^V : \lfloor p \rfloor \le y \le \lceil p \rceil\text{ coordinatewise}\}N(p)={y∈ZV:⌊p⌋≤y≤⌈p⌉ coordinatewise} — a strong, local form of "no holes" saying every real point of the hull is explained by nearby integer points alone.

Formalization targets

Goal: integral convexity of L-convex sets

For an L-convex set D⊆ZVD \subseteq \mathbb Z^VD⊆ZV, writing a=p−⌊p⌋a = p - \lfloor p \rfloora=p−⌊p⌋ for the fractional part of p∈RVp \in \mathbb R^Vp∈RV, α1>⋯>αm\alpha_1 > \cdots > \alpha_mα1​>⋯>αm​ for the distinct nonzero values of aaa, and Ui(p)={v:a(v)≥αi}U_i(p) = \{v : a(v) \ge \alpha_i\}Ui​(p)={v:a(v)≥αi​} (with U0=∅U_0 = \emptysetU0​=∅):

D‾={p∈RV:⌊p⌋+χUi(p)∈D  (i=0,1,…,m)},hence D is integrally convex.\overline D = \{p \in \mathbb R^V : \lfloor p \rfloor + \chi_{U_i(p)} \in D\ \ (i = 0, 1, \ldots, m)\}, \qquad \text{hence } D \text{ is integrally convex}.D={p∈RV:⌊p⌋+χUi​(p)​∈D  (i=0,1,…,m)},hence D is integrally convex.

This is the weakest stable form available: it exhibits an explicit, finite set of at most ∣V∣+1|V|+1∣V∣+1 integer witnesses for every point of the hull, which is what "integrally convex" asserts abstractly, rather than a numerical bound that a sharper construction could later shrink.

Supporting structural targets

Four further results build the correspondence this goal uses: the basic duality between a distance function's admissible potentials, its shortest-path closure, and negative-cycle freedom (Prop. 5.1); the induced-distance-function construction recovering a triangle-inequality distance function from any integer point set, and the convex hull of an L-convex set as its associated polyhedron (Prop. 5.3); the converse construction recovering an L-convex set from an integer-valued distance function (Prop. 5.4); and convexity in Minkowski sum (Thm. 5.8).

Significance

Theorem 5.10 is what makes "L-convex" a genuinely convex-analytic notion rather than a combinatorial curiosity: it shows the convex hull of an L-convex set is not merely a polyhedron (already known from the chapter's polyhedral-description results) but one with the strongest local integrality property discrete convex analysis considers, integral convexity — every real point's hull membership is certified by a small, explicitly constructed set of nearby lattice points, uniformly across the whole set. This is the L-convex counterpart of the corresponding M-convex fact (chapter 4's Theorem 4.24) and is used later in the book wherever L-convex functions (chapter 7) need their epigraphs' local structure. Proposition 5.1 is the combinatorial engine underneath: it is exactly the LP-duality statement between shortest paths and feasible potentials that appears, in various guises, throughout network flow theory, made precise here as the base case the L-convex correspondence rests on.

None of these results are open — Murota presents them as, in his own words, "fundamental facts well known in network flow theory" (Proposition 5.1) systematized into the discrete convex analysis framework. What this mission contributes is a faithful, machine-checked formal statement of each, in the shared Lean vocabulary (LConvexSet, AdmissiblePotentials, ShortestDist) the rest of the Discrete Convex Analysis series can build on; no comparable formalization exists on the platform (see Formalization scope).

Difficulty

The shortest-path closure γˉ\bar\gammaγˉ​ is not a bookkeeping convenience but genuinely graph-theoretic content: proving Proposition 5.1 requires constructing an admissible potential from a shortest-path labeling and, conversely, deriving the negative-cycle-freeness of GγG_\gammaGγ​ from the mere existence of one admissible potential — a min-cost-flow-style LP duality argument, not a direct combinatorial check. Theorem 5.10's difficulty sits in a different place: the naive approach to "DDD is integrally convex" would attempt an inductive argument peeling off one coordinate at a time, but the actual proof constructs a single, uniform family of m+1m+1m+1 witness points from the sorted fractional values of ppp — a Carathéodory-style representation (Eq. (5.11)) that must simultaneously stay inside the integral neighborhood N(p)N(p)N(p) and land in DDD itself via the triangle inequality of DDD's induced distance function, a construction with no one-coordinate-at-a-time shortcut.

Formalization scope

Ground-set elements are a Fintype V with DecidableEq; L-convex sets are Set (V → ℤ); distance functions are V → V → WithTop ℝ; admissible-potential sets are Set (V → ℝ). The shortest-path closure is formalized directly from finite walks (Fin (k+1) → V) rather than via a graph-library shortest-path predicate, matching the book's own construction. The Eq. (5.11) witnesses are built exactly as the book describes them — sorted distinct nonzero fractional values and their level sets — mirroring the Lovász-extension construction of the companion mission 20-ch04b-mconvexsets. No numeric constants are hard-coded anywhere in this mission (rule 7 is vacuous). The goal's explicit witness set (at most ∣V∣+1|V|+1∣V∣+1 points) is not a trivializing special case: it holds for every L-convex set and every point of its hull, with no extra hypothesis narrowing the class. This mission's definitions (LConvexSet, AdmissiblePotentials, DistanceFunction, IsIntegrallyConvex) are redeclared from chunk 05-lconvex-sets (and, for IsIntegrallyConvex/IntegralNeighborhood, from chapter 3's own definitions) rather than imported, since sibling drafts in this series cannot yet reference one another; a later, published version of this book's namespace should consolidate them. Contributions completing any of the five sorrys are welcome; Proposition 5.1's LP-duality argument and the goal's Carathéodory-style construction are the two with the most independent proof content.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • A. J. Hoffman, "On abstract dual linear programs," Naval Research Logistics Quarterly, 10 (1963), pp. 369-373 (feasible-potential duality in network flow theory).
23 thms3 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis V: The M-Optimality Criterion and the Proximity TheoremTextbook

Motivation

Scaling algorithms are one of the standard techniques for solving discrete optimization problems efficiently: instead of searching a huge integer domain directly, an algorithm first solves a coarsened version of the problem — checking optimality only against neighbors reached by a large step size α\alphaα — and then refines the resulting approximate solution down to the true optimum. This strategy is only as good as the guarantee that a coarse-scale local optimum is provably close to a true, fine-scale global optimum; without such a guarantee, refinement could require an unbounded number of steps. Results providing this guarantee are called proximity theorems, and they are a standard tool across combinatorial optimization, from network flow scaling algorithms to submodular function minimization.

M-convex functions, the subject of this chapter, are exactly the class of discrete convex functions for which the classical local-optimality test of chapter 3 (checking a full neighborhood of up to 3n−13^n-13n−1 sign patterns) sharpens to a much smaller, purely pairwise test: checking f(x)≤f(x−χu+χv)f(x) \le f(x - \chi_u + \chi_v)f(x)≤f(x−χu​+χv​) for every pair of coordinates u,vu, vu,v. This mission formalizes the chapter's central definitional equivalence (Theorem 6.2), this pairwise optimality criterion (Theorem 6.26), a structural minimizer-cut lemma (Theorem 6.28), and the chapter's capstone, the M-proximity theorem (Theorem 6.37) — the result that makes M-convex scaling algorithms provably correct, with an explicit, dimension-and-scale-only distance bound between a coarse-scale local optimum and a true global minimizer.

Setting

Let VVV be a finite ground set. A function f:ZV→R∪{+∞}f : \mathbb Z^V \to \mathbb R \cup \{+\infty\}f:ZV→R∪{+∞} with nonempty effective domain dom⁡f\operatorname{dom} fdomf is an M-convex function if it satisfies the exchange axiom (M-EXC[Z]): for x,y∈dom⁡fx, y \in \operatorname{dom} fx,y∈domf and uuu in the positive support of x−yx - yx−y, there is vvv in the negative support of x−yx-yx−y with

f(x)+f(y)≥f(x−χu+χv)+f(y+χu−χv).f(x) + f(y) \ge f(x - \chi_u + \chi_v) + f(y + \chi_u - \chi_v).f(x)+f(y)≥f(x−χu​+χv​)+f(y+χu​−χv​).

An M♮^\natural♮-convex function is one whose lift f~\tilde ff~​ to the extended ground set V~={0}∪V\tilde V = \{0\} \cup VV~={0}∪V — defined by f~(x0,x)=f(x)\tilde f(x_0, x) = f(x)f~​(x0​,x)=f(x) when x0=−x(V)x_0 = -x(V)x0​=−x(V), and +∞+\infty+∞ otherwise — is M-convex; equivalently (Theorem 6.2, below) fff satisfies the axiom (M♮^\natural♮-EXC[Z]), a variant of (M-EXC[Z]) that additionally allows a single-coordinate move (uuu alone, with no compensating vvv). Every M-convex function is M♮^\natural♮-convex, but not conversely. For α\alphaα a positive integer, a point satisfies the scaled local optimality condition at scale α\alphaα if f(xα)≤f(xα+α(χv−χu))f(x_\alpha) \le f(x_\alpha + \alpha(\chi_v - \chi_u))f(xα​)≤f(xα​+α(χv​−χu​)) for all relevant u,vu, vu,v — a check against neighbors α\alphaα steps away rather than adjacent ones.

Formalization targets

Goal: Theorem 6.37 (the M-proximity theorem)

Assume α\alphaα is a positive integer and n=∣V∣n = |V|n=∣V∣.

(1) f M-convex, f(xα)≤f(xα+α(χv−χu)) ∀u,v  ⟹  ∃x∗∈arg⁡min⁡f, ∥xα−x∗∥∞≤(n−1)(α−1),\text{(1) } f \text{ M-convex, } f(x_\alpha) \le f(x_\alpha + \alpha(\chi_v-\chi_u))\ \forall u,v \implies \exists x^* \in \arg\min f,\ \|x_\alpha - x^*\|_\infty \le (n-1)(\alpha-1),(1) f M-convex, f(xα​)≤f(xα​+α(χv​−χu​)) ∀u,v⟹∃x∗∈argminf, ∥xα​−x∗∥∞​≤(n−1)(α−1), (2) f M♮-convex, same hypothesis over u,v∈V∪{0}  ⟹  ∃x∗∈arg⁡min⁡f, ∥xα−x∗∥∞≤n(α−1).\text{(2) } f \text{ M}^\natural\text{-convex, same hypothesis over } u,v \in V \cup \{0\} \implies \exists x^* \in \arg\min f,\ \|x_\alpha - x^*\|_\infty \le n(\alpha-1).(2) f M♮-convex, same hypothesis over u,v∈V∪{0}⟹∃x∗∈argminf, ∥xα​−x∗∥∞​≤n(α−1).

Both bounds are exact and specific to their hypothesis class; replacing either with an unspecified function of nnn and α\alphaα would discard exactly the content chapter 10's algorithms rely on.

Milestones: Theorems 6.2, 6.26, 6.28

Theorem 6.2: M♮^\natural♮-convexity (defined via the lift) is equivalent to the direct exchange axiom (M♮^\natural♮-EXC[Z]) — the chapter's central definitional theorem, needed to work with M♮^\natural♮-convex functions without repeatedly invoking the lift construction. Theorem 6.26 (the M-optimality criterion): global optimality of fff at xxx is equivalent to a purely pairwise local check, f(x)≤f(x−χu+χv)f(x) \le f(x-\chi_u+\chi_v)f(x)≤f(x−χu​+χv​) for all u,vu,vu,v (plus, in the M♮^\natural♮ case, f(x)≤f(x±χv)f(x) \le f(x\pm\chi_v)f(x)≤f(x±χv​)). Theorem 6.28 (the M-minimizer cut): from any point and any coordinate pair minimizing a one-step exchange, one can certify a coordinate-wise bound that some global minimizer must satisfy — the structural fact underlying both the domain-reduction algorithm and, via the same proof technique, the proximity theorem itself.

Significance

The result itself. Theorem 6.26 already sharpens chapter 3's local-to-global criterion (checking a full 3n−13^n-13n−1-point neighborhood) to an O(n2)O(n^2)O(n2)-size pairwise check — the minimum spanning tree optimality criterion is a direct special case. The proximity theorem builds on this to control what happens when the local check is only performed at a coarse scale α\alphaα: it guarantees that scaling-based algorithms, which alternate between coarse-scale local search and scale reduction, terminate with a guaranteed-close approximation at every stage, with an explicit linear-in-nnn, linear-in-α\alphaα error bound rather than a qualitative "eventually converges" guarantee.

Formalizing it. No matching item exists on the platform for M-convex functions, the exchange axiom, or a discrete proximity theorem of this kind. This mission gives the first formal statement of the M-optimality criterion and the M-proximity theorem, together with the exchange-axiom / lift-based-definition equivalence (Theorem 6.2) that the rest of the M-convex function theory (chunks 07, and indirectly 10–14) is built on.

Difficulty

The natural first attempt at Theorem 6.37 is to try a direct coordinatewise argument: since the scaled hypothesis holds for every pair u,vu, vu,v, one might hope to bound ∣xα(v)−x∗(v)∣|x_\alpha(v) - x^*(v)|∣xα​(v)−x∗(v)∣ coordinate by coordinate independently. This does not work, because a single application of the exchange axiom only ever improves fff by trading one coordinate down and one other coordinate up simultaneously — there is no way to move a single coordinate toward a minimizer in isolation without accounting for where the compensating mass goes. The actual proof instead fixes a target coordinate vvv, constructs a chain of strictly decreasing function values y0=xα,y1,…,yky_0 = x_\alpha, y_1, \ldots, y_ky0​=xα​,y1​,…,yk​ by repeatedly applying (M-EXC[Z]) against a fixed near-optimal point x∗x^*x∗ (exactly the technique of Theorem 6.28's proof), and then bounds how far each other coordinate can move along this chain using the scaled hypothesis itself, before summing those bounds via the M-convex domain's hyperplane constraint x(V)=x(V) = x(V)= constant to recover the bound on vvv. The chain construction, not a per-coordinate estimate, is what makes the linear-in-nnn bound provable at all.

Formalization scope

The ground set VVV is a Fintype with DecidableEq; f:ZV→R∪{+∞}f : \mathbb Z^V \to \mathbb R \cup \{+\infty\}f:ZV→R∪{+∞} is (V → ℤ) → WithTop ℝ. M♮^\natural♮-convexity is represented via an explicit lift to Option V (none standing for the extended ground set's new element 000), matching the book's own primary definition; the direct exchange-axiom form is a separate predicate related to it by Theorem 6.2, not conflated with it. `‖x_\alpha - x^*|_\infty \le c$ is stated pointwise.

A trivializing formalization of the goal would replace either exact bound, (n−1)(α−1)(n-1)(\alpha-1)(n−1)(α−1) or n(α−1)n(\alpha-1)n(α−1), with an unspecified asymptotic bound, or merge the two hypothesis classes into a single weaker statement; neither is done here. Propositions establishing dom f as an M-convex set, the M/M♮^\natural♮ relationship (Theorem 6.3), and several structural closure properties are cut from this mission's scope (not needed by the chosen items' statements — see MODERATION_NOTES.md) but are natural targets for a follow-on mission or for chunk 07, which builds directly on this chunk's exchange-axiom vocabulary. Contributions building the arg min f M-convexity corollary (Proposition 6.29) or the scaled minimizer cut (Theorem 6.39, the direct generalization of Theorem 6.28 drafted here) are welcome.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • D. S. Hochbaum, "Lower and upper bounds for the allocation problem and other nonlinear optimization problems," Mathematics of Operations Research, 19(2), 1994, pp. 390–409.
14 thms2 active usersReviewed
🏆Completed
CombinatoricsConvex OptimizationDiscrete Geometry+2·Captain: Shuze Chen

Discrete Convex Analysis XXIII: Directional Derivatives and Subdifferentials of M-Convex FunctionsTextbook

Motivation

An M-convex function is defined on the integer lattice, but chapter 6's earlier results (companion missions 06-mconvex-functions-i, 22-ch06b-mconvexfunctions, 23-ch06c-mconvexfunctions) show it always extends to a genuine convex function on real space. Once that extension exists, every tool of classical convex analysis — directional derivatives, subdifferentials, positive homogeneity — becomes available, and the natural question is whether these classical objects remain combinatorially special when applied to an M-convex function's extension. This mission answers that question at its sharpest: the directional derivative of an M-convex function at any point is again a positively homogeneous M-convex function, its subdifferential is exactly the admissible-potential set of a distance function satisfying the triangle inequality, and this correspondence between positively homogeneous M-convex functions and triangle-inequality distance functions is itself a clean one-to-one correspondence. This closes the loop between chapters 4-5 (M-convex and L-convex sets, distance functions) and the continuous convex-analytic machinery chapter 8 needs for its duality theory.

Companion missions 06-mconvex-functions-i, 22-ch06b-mconvexfunctions, and 23-ch06c-mconvexfunctions cover this chapter's optimality theory, algebraic toolkit, and convex-extensibility characterization. This mission builds the vocabulary those results also need (redeclared here, since sibling drafts cannot yet import one another) and proves the chapter's real-variable capstones: the transfer of M-convexity's basic operations, optimality criterion, and supermodularity to the polyhedral (real-variable) setting, the identification of positively homogeneous M-convex functions with distance functions satisfying the triangle inequality, and — this mission's goal — the full directional-derivative/subdifferential correspondence.

Setting

Fix a finite ground set VVV. A polyhedral convex function g:RV→R∪{+∞}g : \mathbb R^V \to \mathbb R \cup \{+\infty\}g:RV→R∪{+∞} is (polyhedral) M-convex if it satisfies the real-variable exchange axiom (M-EXC[R]): for x,y∈dom⁡Rgx,y \in \operatorname{dom}_{\mathbb R} gx,y∈domR​g and u∈supp⁡+(x−y)u \in \operatorname{supp}^+(x-y)u∈supp+(x−y), some v∈supp⁡−(x−y)v \in \operatorname{supp}^-(x-y)v∈supp−(x−y) and α0>0\alpha_0 > 0α0​>0 make the exchange inequality hold on α∈[0,α0]\alpha \in [0,\alpha_0]α∈[0,α0​]; M♮-convex if its lift to one extra coordinate is M-convex. The directional derivative of ggg at x∈dom⁡Rgx \in \operatorname{dom}_{\mathbb R} gx∈domR​g in direction ddd is g′(x;d)=inf⁡t>0(g(x+td)−g(x))/tg'(x;d) = \inf_{t>0} (g(x+td) - g(x))/tg′(x;d)=inft>0​(g(x+td)−g(x))/t. A function is positively homogeneous if g(tx)=t⋅g(x)g(tx) = t \cdot g(x)g(tx)=t⋅g(x) for all t>0t > 0t>0; write 0M[R→R]0M[\mathbb R \to \mathbb R]0M[R→R] for the positively homogeneous polyhedral M-convex functions. A distance function γ\gammaγ satisfying the triangle inequality and its set of admissible potentials D(γ)D(\gamma)D(γ) were introduced in chapter 5; the subdifferential ∂Rf(x)={p:f(y)−f(x)≥⟨p,y−x⟩ ∀y}\partial_{\mathbb R} f(x) = \{p : f(y) - f(x) \ge \langle p, y-x \rangle\ \forall y\}∂R​f(x)={p:f(y)−f(x)≥⟨p,y−x⟩ ∀y} generalizes this to any function fff at a point xxx in its domain.

Formalization targets

Goal: the directional-derivative/subdifferential correspondence

For f∈M[R→R]f \in M[\mathbb R \to \mathbb R]f∈M[R→R] and x∈dom⁡Rfx \in \operatorname{dom}_{\mathbb R} fx∈domR​f, setting γf,x(u,v)=f′(x;−χu+χv)\gamma_{f,x}(u,v) = f'(x;-\chi_u+\chi_v)γf,x​(u,v)=f′(x;−χu​+χv​):

γf,x satisfies the triangle inequality,∂Rf(x)=D(γf,x)≠∅,f′(x;⋅)=γf,x^(⋅),\gamma_{f,x} \text{ satisfies the triangle inequality}, \quad \partial_{\mathbb R} f(x) = D(\gamma_{f,x}) \ne \emptyset, \quad f'(x;\cdot) = \widehat{\gamma_{f,x}}(\cdot),γf,x​ satisfies the triangle inequality,∂R​f(x)=D(γf,x​)=∅,f′(x;⋅)=γf,x​​(⋅),

with the analogous statement for f∈M[Z→R]f \in M[\mathbb Z \to \mathbb R]f∈M[Z→R] at an integer point xxx, using γf,x(u,v)=f(x−χu+χv)−f(x)\gamma_{f,x}(u,v) = f(x-\chi_u+\chi_v)-f(x)γf,x​(u,v)=f(x−χu​+χv​)−f(x) (Theorem 6.61). This is the weakest stable form: it identifies the subdifferential exactly, as a set, rather than bounding its size or complexity, and holds at every point of the domain uniformly.

Supporting structural targets

Ten further results build the real-variable toolkit and the positive-homogeneity correspondence this goal completes: the transfer of M♮-convexity, the basic operations, the optimality criterion, supermodularity, and weighted-minimizer polyhedrality to the real-variable setting (Theorems 6.48-6.52, Proposition 6.53), the identification of the classes 0M[Z∣R→R]0M[\mathbb Z|\mathbb R \to \mathbb R]0M[Z∣R→R] and 0M[R→R]0M[\mathbb R \to \mathbb R]0M[R→R] and the compatibility of convex extension with positive homogeneity (Proposition 6.56), the two directions of the correspondence between positively homogeneous M-convex functions and triangle-inequality distance functions (Propositions 6.57-6.58, Theorem 6.59), and the fact that a directional derivative of an M-convex function is itself positively homogeneous and M-convex (Proposition 6.60).

Significance

Theorem 6.61 is the technical bridge that lets discrete convex analysis borrow the entire apparatus of classical convex duality: because the subdifferential of an M-convex function is always the admissible-potential set of a chapter-5 distance function, every fact already proved about D(γ)D(\gamma)D(γ) (its polyhedral structure, its own L-convexity, its relationship to shortest paths) transfers immediately to subdifferentials of M-convex functions. This is exactly the mechanism the book calls out as essential for Chapter 8's separation theorem for M♮-convex functions. The 0M↔T0M \leftrightarrow T0M↔T correspondence (Theorem 6.59) is independently significant: it says the positively homogeneous special case of M-convex function theory — which is what directional derivatives of any M-convex function reduce to, by Proposition 6.60 — is exactly as rich as ordinary shortest-path distance function theory, no more and no less, so nothing new needs to be built to understand local behavior at a point.

None of these results are open — they are Murota's account of how the discrete exchange axiom interacts with directional differentiation and subgradients, a bridge chapter between the purely combinatorial theory of chapters 4-6 and the duality theory of chapter 8. What this mission contributes is a faithful, machine-checked formal statement of each, extending the shared Lean vocabulary (MExchangeAxiomR, DirDeriv, GammaHat) the Discrete Convex Analysis series builds on; no comparable formalization exists on the platform (see Formalization scope).

Difficulty

The naive approach to Theorem 6.61 would try to compute ∂Rf(x)\partial_{\mathbb R} f(x)∂R​f(x) directly from the definition of subgradient and separately verify it happens to equal some D(γ)D(\gamma)D(γ); the book's actual proof instead derives the equality of sets from the M-optimality criterion (Theorem 6.52) applied pointwise: p∈∂Rf(x)p \in \partial_{\mathbb R} f(x)p∈∂R​f(x) is shown, via a chain of logical equivalences, to be exactly the condition defining D(γf,x)D(\gamma_{f,x})D(γf,x​), so no separate verification of polyhedrality or nonemptiness is needed beyond what Theorem 6.52 and Proposition 6.60 already supply. The genuine difficulty is upstream, in Proposition 6.60 itself: showing a directional derivative is M-convex requires exploiting the local validity of the identity f(x+d)−f(x)=f′(x;d)f(x+d)-f(x) = f'(x;d)f(x+d)−f(x)=f′(x;d) for small ∥d∥1\|d\|_1∥d∥1​ (Eq. (6.85)) and then extending the exchange property from that neighborhood to all of RV\mathbb R^VRV using positive homogeneity — a two-step argument with no single-step shortcut, since the exchange axiom's defining inequality is not obviously homogeneous-invariant on its own.

Formalization scope

Ground-set elements are a Fintype V with DecidableEq; real-domain functions are (V→ℝ)→WithTop ℝ. The directional derivative is built directly as an infimum of difference quotients over t>0t>0t>0, matching the book's own local characterization (Eq. (6.85)) without a separate limit construction. Positive homogeneity and the classes 0M[R→R]/0M[Z→R] are stated exactly as the book defines them (the latter via positive homogeneity of the convex extension, not of f itself, since f is undefined off Zⱽ). Theorems 6.49-6.50 restate 4 of their 8 operations (matching the identical scope decision for chunk 22-ch06b-mconvexfunctions's Theorem 6.13); Theorem 6.61 omits the dual-integral refinement clauses for the M[R→R|Z]/ M[Z→Z] sub-classes. Both reductions are documented, not trivializing omissions — see Difficulty above and HARD.md/MODERATION_NOTES.md. No numeric constants are hard-coded anywhere in this mission. This mission's definitions are redeclared from chunks 06-mconvex-functions-i, 21-ch05b-lconvexsets (for the distance-function/admissible-potential vocabulary), 22-ch06b-mconvexfunctions, and 23-ch06c-mconvexfunctions rather than imported, since sibling drafts in this series cannot yet reference one another. Contributions completing any of the twelve sorrys are welcome; the goal and Proposition 6.60 carry the most independent proof content.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • K. Murota and A. Shioura, "M-convex function on generalized polymatroid," Mathematics of Operations Research, 24 (1999), pp. 95-105 (the polyhedral M-convex function theory this mission's real-variable results are drawn from).
56 thms2 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis XXIV: Quasi M-Convex FunctionsTextbook

Motivation

Every characterization of M-convexity so far in this series — the exchange axiom, the optimality criterion, convex extensibility, the directional-derivative/subdifferential correspondence — has been an equivalence with M-convexity itself: a function either is M-convex or it is not. Section 6.14 asks a different question: what happens when the exchange axiom's defining inequality is relaxed to only the sign patterns it actually forces? The answer is a hierarchy of "quasi M-convex" conditions — weaker than M-convexity, strong enough to keep the optimality criterion and the proximity/minimizer-cut theorems intact — and, at the top of that hierarchy, a genuinely new characterization of M-convexity itself: a function is M-convex if and only if every one of its linear perturbations is quasi M-convex in the weakest sense. This mission formalizes that entire hierarchy and its capstone, plus two further characterizations of polyhedral M-convexity (via directional derivatives, subdifferentials, and weighted-minimizer polyhedra) that complete the real-variable theory chunk 24-ch06d-mconvexfunctions began.

Setting

Fix a finite ground set VVV and f:ZV→R∪{+∞}f : \mathbb Z^V \to \mathbb R \cup \{+\infty\}f:ZV→R∪{+∞}. Write Δf(z;v,u)=f(z+χv−χu)−f(z)\Delta f(z;v,u) = f(z + \chi_v - \chi_u) - f(z)Δf(z;v,u)=f(z+χv​−χu​)−f(z). Relaxing the exchange axiom's inequality Δf(x;v,u)+Δf(y;u,v)≤0\Delta f(x;v,u) + \Delta f(y;u,v) \le 0Δf(x;v,u)+Δf(y;u,v)≤0 to the sign patterns it forces gives quasi M-convexity (QM) and semistrict quasi M-convexity (SSQM), and requiring only some pair (u,v)(u,v)(u,v) rather than every uuu gives their weaker variants (QMw), (SSQMw); the minimization-only variants (SSQM≠\ne=), (SSQM≠w\ne_w=w​) replace "x,y∈dom⁡fx,y \in \operatorname{dom} fx,y∈domf" with "f(x)≠f(y)f(x) \ne f(y)f(x)=f(y)". The set-level analogue (Q-EXC)/(Q-EXCw) relaxes the M-convex-set exchange axiom the same way. For α∈R\alpha \in \mathbb Rα∈R, the level set L(f,α)={x∈ZV:f(x)≤α}L(f,\alpha) = \{x \in \mathbb Z^V : f(x) \le \alpha\}L(f,α)={x∈ZV:f(x)≤α}. A polyhedral convex function f:RV→R∪{+∞}f : \mathbb R^V \to \mathbb R \cup \{+\infty\}f:RV→R∪{+∞} is (real-variable) M-convex, f∈M[R→R]f \in M[\mathbb R \to \mathbb R]f∈M[R→R], if it satisfies the real exchange axiom (M-EXC[R]) from chunk 24-ch06d-mconvexfunctions; L0[R]L_0[\mathbb R]L0​[R] denotes polyhedra realized as D(γ)D(\gamma)D(γ) for a triangle-inequality distance function γ\gammaγ, and M0[R]M_0[\mathbb R]M0​[R] denotes real M-convex polyhedral cones.

Formalization targets

Goal: the quasi M-convexity hierarchy (Theorem 6.68)

For f:ZV→R∪{+∞}f : \mathbb Z^V \to \mathbb R \cup \{+\infty\}f:ZV→R∪{+∞}: (1) the implication diagram (M-EXC[Z]) ⇒\Rightarrow⇒ (SSQM) ⇒\Rightarrow⇒ (QM), (M-EXCw[Z]) ⇒\Rightarrow⇒ (SSQMw) ⇒\Rightarrow⇒ (QMw), (M-EXC[Z]) ⇔\Leftrightarrow⇔ (M-EXCw[Z]), (SSQM) ⇒\Rightarrow⇒ (SSQMw), (QM) ⇒\Rightarrow⇒ (QMw); (2) fff satisfies (M-EXC[Z]) if and only if f[p]f[p]f[p] satisfies (QMw) for every p∈RVp \in \mathbb R^Vp∈RV. This theorem is not named in the chunk's own extraction table — the automated extractor, which requires a result's label to start a text line, misses it because it opens mid-paragraph directly after an ASCII-rendered implication diagram — but it is the capstone of the section: its own proof is "combining Theorems 6.72 and 6.74," both formalized here as milestones, and part (2) answers exactly the question a reader of this stretch would ask: what does the entire apparatus of quasi M-convexity ultimately say about M-convexity itself?

Supporting structural targets

Fourteen further results build the hierarchy and complete the real-variable theory. Proposition 6.62 checks the two halves of chunk 24-ch06d-mconvexfunctions's Theorem 6.61 agree at integer points. Theorems 6.63–6.64 add two more characterizations of polyhedral M-convexity — via positively homogeneous directional derivatives and L0[R]L_0[\mathbb R]L0​[R]-valued subdifferentials, and via M0[R]M_0[\mathbb R]M0​[R]/M0[Z∣R]M_0[\mathbb Z|\mathbb R]M0​[Z∣R]-valued weighted-minimizer polyhedra — to the convex-extensibility characterization chunk 23-ch06c-mconvexfunctions proved. Theorem 6.67 gives two equivalent reformulations of (QMw) as pointwise inequalities; Propositions 6.69–6.70 and Theorems 6.72–6.74 build the level-set/perturbation machinery the goal needs. Theorem 6.75 is the (SSQM≠w\ne_w=w​) analogue of Theorem 6.67. Theorem 6.76 is the quasi M-optimality criterion (optimality still characterized by local non-improvement, under only the weak quasi-convexity hypotheses). Theorems 6.77–6.79 show the M-minimizer-cut and M-proximity theorems (from missions 06-mconvex-functions-i and 23-ch06c-mconvexfunctions) hold verbatim under the strictly weaker (SSQM≠\ne=) hypothesis.

Significance

The hierarchy's practical payoff is immediate: Theorems 6.77–6.79 mean the algorithms of chapter 10 that rely on minimizer cuts and proximity bounds do not actually need the full exchange axiom to run correctly on nonlinearly rescaled M-convex functions (Example 6.66 shows any nondecreasing scaling ϕ∘f\phi \circ fϕ∘f of an M-convex fff is quasi M-convex, yet nonlinear scalings are common in practice and destroy M-convexity itself). The goal, Theorem 6.68, is significant independently: it says the exchange axiom — a condition that looks irreducibly combinatorial, quantifying over pairs of points and directions — is equivalent to a purely ordinal, perturbation-based condition (every linear tilt of fff has no strict local improvement that a level set can't witness), giving a genuinely different lens on why M-convexity is the right discrete analogue of convexity. Theorems 6.63–6.64 close out chunk 24-ch06d-mconvexfunctions's program of characterizing polyhedral M-convexity in every classical convex-analytic vocabulary at once (directional derivatives, subdifferentials, weighted minimizers), completing the bridge to Chapter 8's duality theory that chunk builds toward.

None of these results are open — they are Murota's account of how far the exchange axiom's defining inequality can be relaxed while keeping optimization theory intact. What this mission contributes is a faithful, machine-checked formal statement of each, including four theorems (6.68, 6.76, 6.77, 6.78) the platform's own automated extractor missed entirely, extending the shared Lean vocabulary (DeltaF, QMw, LevelSet) the Discrete Convex Analysis series builds on; no comparable formalization exists on the platform (see Formalization scope).

Difficulty

The naive approach to the goal would try to prove the implication diagram's six arrows and the perturbation equivalence as six independent facts; the book's own proof of part (2) instead derives it in one step from Theorems 6.72 and 6.74 — themselves nontrivial (Theorem 6.74's proof strengthens the local exchange axiom equivalence (Theorem 6.4) to hold whenever the domain merely satisfies (Q-EXCw), then runs a bipartite-matching argument on a 4-point neighborhood to verify the resulting local condition). The difficulty is genuinely upstream of the goal's own statement: everything the goal needs is already proved by the time Theorem 6.68 is reached, so the formalization work is in stating the sixteen distinct axioms and their level-set reformulations precisely enough that "combining 6.72 and 6.74" is literally how a Lean proof would proceed.

Formalization scope

Ground-set elements are a Fintype V with DecidableEq; integer-domain functions are (V→ℤ)→WithTop ℝ, real-domain ones (V→ℝ)→WithTop ℝ. All fifteen numbered results — including the four (6.68, 6.76, 6.77, 6.78) the automated extractor missed because their labels open mid-paragraph — are placed as milestone or goal, and every clause of every one is stated in full; no partial-coverage scope reduction was needed in this chunk (contrast chunks 22-ch06b-mconvexfunctions/24-ch06d-mconvexfunctions, which restated 4-of-8-part operations theorems). Two formalization choices are recorded in HARD.md/MODERATION_NOTES.md: "inf⁡f[−p]>−∞\inf f[-p] > -\inftyinff[−p]>−∞" is replaced by the equivalent (ArgMinOn ...).Nonempty hypothesis (WithTop ℝ has no −∞-\infty−∞ element), matching mission 23-ch06c-mconvexfunctions's identical substitution; and M0[R]M_0[\mathbb R]M0​[R]/M0[Z∣R]M_0[\mathbb Z|\mathbb R]M0​[Z∣R] are realized via the book's own indicator-function device rather than a freestanding cone axiom. This mission's definitions are redeclared from chunks 06-mconvex-functions-i, 22–24-ch06*-mconvexfunctions rather than imported, since sibling drafts in this series cannot yet reference one another. Contributions completing any of the fifteen sorrys are welcome; the goal and Theorem 6.74 carry the most independent proof content.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • M. Avriel, W. E. Diewert, S. Schaible, and I. Zang, Generalized Concavity, Plenum Press, 1988 (the continuous quasi-convexity theory this chapter's discrete analogue generalizes).
56 thms3 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis VII: The L-Optimality Criterion and the Proximity TheoremTextbook

Motivation

Submodularity — the diminishing-returns property g(p)+g(q)≥g(p∨q)+g(p∧q)g(p) + g(q) \ge g(p \vee q) + g(p \wedge q)g(p)+g(q)≥g(p∨q)+g(p∧q) on a lattice — is one of the most useful structural hypotheses in combinatorial optimization, underlying efficient algorithms for network flows, matroid theory, and set-function minimization. Chapter 7 studies L-convex functions: functions on the integer lattice ZV\mathbb Z^VZV that are submodular and linear along the all-ones direction. This is the "dual" notion, under the conjugacy developed later in the book, to chunk 06's M-convex functions, and it inherits the same strong minimization theory — a purely local optimality criterion and a proximity theorem with an explicit distance bound — while additionally supporting a genuinely new characterization with no M-convex counterpart: discrete midpoint convexity, the direct lattice analogue of the classical real-valued midpoint convexity condition. This mission formalizes the chapter's definitional theorem, its midpoint-convexity characterization, the L-optimality criterion, and the L-proximity theorem itself.

Setting

Let VVV be a finite ground set. A function g:ZV→R∪{+∞}g : \mathbb Z^V \to \mathbb R \cup \{+\infty\}g:ZV→R∪{+∞} with nonempty effective domain is an L-convex function if it satisfies (SBF[Z]): g(p)+g(q)≥g(p∨q)+g(p∧q)g(p) + g(q) \ge g(p \vee q) + g(p \wedge q)g(p)+g(q)≥g(p∨q)+g(p∧q) for all p,qp, qp,q (∨,∧\vee, \wedge∨,∧ componentwise max/min), and (TRF[Z]): there is r∈Rr \in \mathbb Rr∈R with g(p+1)=g(p)+rg(p + \mathbf 1) = g(p) + rg(p+1)=g(p)+r for all ppp, where 1\mathbf 11 is the all-ones vector. An L♮^\natural♮-convex function is one whose lift to the extended ground set {0}∪V\{0\} \cup V{0}∪V is L-convex; equivalently (Theorem 7.1), ggg satisfies the translation-submodularity axiom (SBF♮^\natural♮[Z]): g(p)+g(q)≥g((p−α1)∨q)+g(p∧(q+α1))g(p) + g(q) \ge g((p - \alpha\mathbf 1) \vee q) + g(p \wedge (q + \alpha\mathbf 1))g(p)+g(q)≥g((p−α1)∨q)+g(p∧(q+α1)) for all p,qp, qp,q and all nonnegative integers α\alphaα. Discrete midpoint convexity asks g(p)+g(q)≥g(⌈(p+q)/2⌉)+g(⌊(p+q)/2⌋)g(p) + g(q) \ge g(\lceil (p+q)/2 \rceil) + g(\lfloor (p+q)/2 \rfloor)g(p)+g(q)≥g(⌈(p+q)/2⌉)+g(⌊(p+q)/2⌋) componentwise. For α\alphaα a positive integer, a point satisfies scaled local optimality if g(pα)≤g(pα±αχY)g(p_\alpha) \le g(p_\alpha \pm \alpha \chi_Y)g(pα​)≤g(pα​±αχY​) for every Y⊆VY \subseteq VY⊆V.

Formalization targets

Goal: Theorem 7.18 (the L-proximity theorem)

Assume α\alphaα is a positive integer and n=∣V∣n = |V|n=∣V∣. (1) If ggg is L-convex with g(p)=g(p+1)g(p) = g(p+\mathbf 1)g(p)=g(p+1) for all ppp, and pα∈dom⁡gp_\alpha \in \operatorname{dom} gpα​∈domg satisfies g(pα)≤g(pα+αχY)g(p_\alpha) \le g(p_\alpha + \alpha\chi_Y)g(pα​)≤g(pα​+αχY​) for all Y⊆VY \subseteq VY⊆V, then arg⁡min⁡g≠∅\arg\min g \ne \emptysetargming=∅ and there is p∗∈arg⁡min⁡gp^* \in \arg\min gp∗∈argming with the componentwise bound

pα≤p∗≤pα+(n−1)(α−1)1.p_\alpha \le p^* \le p_\alpha + (n-1)(\alpha-1)\mathbf 1.pα​≤p∗≤pα​+(n−1)(α−1)1.

(2) If ggg is L♮^\natural♮-convex and pαp_\alphapα​ satisfies the two-sided version, then there is p∗p^*p∗ with pα−n(α−1)1≤p∗≤pα+n(α−1)1p_\alpha - n(\alpha-1)\mathbf 1 \le p^* \le p_\alpha + n(\alpha-1)\mathbf 1pα​−n(α−1)1≤p∗≤pα​+n(α−1)1. The bound is a genuine vector (lattice-order) inequality, not an ℓ∞\ell^\inftyℓ∞-norm bound — the form later chapters' applications need.

Milestones: Theorems 7.1, 7.7, 7.14

Theorem 7.1: L♮^\natural♮-convexity (defined via the lift) is equivalent to the direct translation-submodularity axiom. Theorem 7.7: this same class is also characterized by discrete midpoint convexity — a three-way equivalence with the approach property (L♮^\natural♮-APR[Z]) as a bridge — giving L-convexity a genuinely different, more geometric face than anything available on the M-convex side. Theorem 7.14 (the L-optimality criterion): global optimality reduces to a purely local check against the sign-pattern neighbors p±χYp \pm \chi_Yp±χY​, mirroring chunk 06's Theorem 6.26 but with the plain L-convex case additionally requiring the periodicity condition g(p)=g(p+1)g(p) = g(p+\mathbf 1)g(p)=g(p+1).

Significance

The result itself. Discrete midpoint convexity (Theorem 7.7) is philosophically important: it shows the lattice-submodularity definition of L-convexity is not an arbitrary discretization choice but coincides exactly with the most direct discrete analogue of ordinary midpoint convexity, the classical characterization of convex functions via f((p+q)/2)≤(f(p)+f(q))/2f((p+q)/2) \le (f(p)+f(q))/2f((p+q)/2)≤(f(p)+f(q))/2. The L-optimality criterion and L-proximity theorem give L-convex minimization the same algorithmic footing as M-convex minimization (chunk 06): scaling algorithms for L-convex objectives — which arise naturally from network flow and submodular-function duality — inherit a provable, dimension-and-scale-explicit distance guarantee between a coarse-scale local optimum and the true minimizer.

Formalizing it. No matching item exists on the platform for L-convex functions, discrete midpoint convexity, or the L-optimality/proximity theorems. This mission gives the first formal statement of these results, completing (alongside chunk 06's M-convex-function results) both halves of the exchange-axiom-based theory that chapter 8's conjugacy duality later unifies.

Difficulty

A natural shortcut, given the structural parallel to chunk 06, is to assume the L-proximity theorem's proof is a mechanical relabeling of the M-proximity theorem's proof. It is not: the M-convex proof (chunk 06) crucially uses the exchange axiom's additive four-term inequality to build a chain of strictly improving points, whereas the L-convex proof instead exploits (TRF[Z])'s periodicity directly — it reduces to the case pα=0p_\alpha = 0pα​=0 using translation invariance, then constructs a minimal (with respect to the lattice order) point among all sufficiently good solutions and shows this minimality, combined with submodularity (SBF[Z]), forces the componentwise bound. The vector (rather than norm) form of the conclusion is not cosmetic: it is exactly what this lattice-order argument naturally produces, and is the form needed by later chapters' applications.

Formalization scope

The ground set VVV is a Fintype with DecidableEq; g:ZV→R∪{+∞}g : \mathbb Z^V \to \mathbb R \cup \{+\infty\}g:ZV→R∪{+∞} is (V → ℤ) → WithTop ℝ. Unlike chunk 06's M-convex axiom, (SBF[Z]), (TRF[Z]), and (SBF♮^\natural♮[Z]) are stated for all of ZV\mathbb Z^VZV, not restricted to dom⁡g\operatorname{dom} gdomg, so no explicit import of chunk 05's L-convex-set vocabulary was needed for dom g's structure (unlike the corresponding note in chunk 06's BRIEF.md, which flagged the same concern for dom f). L♮^\natural♮-convexity is represented via an explicit lift to Option V, matching the book's own primary definition, with the direct axiom (SBF♮^\natural♮[Z]) kept as a separate object related to it by Theorem 7.1.

A trivializing formalization of the goal would convert its componentwise vector bound into an ℓ∞\ell^\inftyℓ∞-norm bound (losing the direction-of-approach information the vector form carries) or drop Part (1)'s periodicity hypothesis g(p)=g(p+1)g(p) = g(p+\mathbf 1)g(p)=g(p+1); neither is done here. Propositions establishing dom g as an L-convex set, the L/L♮^\natural♮ relationship (Theorem 7.3), the submodular-set-function embedding (Proposition 7.4), and several structural closure properties are cut from this mission's scope (see MODERATION_NOTES.md) but are natural targets for a follow-on mission or for chunk 09, which builds directly on this chunk's exchange-axiom vocabulary, mirroring chunks 06→07.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
18 thms2 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis VIII: Quasi L-Convex Functions and the Quasi-Proximity TheoremTextbook

Motivation

Milgrom and Shannon's theory of quasi-supermodularity, developed for monotone comparative statics in economics, showed that many of the consequences of lattice submodularity survive under a much weaker, purely ordinal relaxation of the defining inequality. Chapter 7's final section imports this idea into discrete convex analysis: does L-convexity's optimality and proximity theory survive when the additive submodularity inequality is relaxed to an ordinal condition on the sign pattern of the two relevant differences, rather than their sum? This mission formalizes the chapter's answer for the strongest of the relevant relaxations, (SSQSB) (semistrict quasi submodularity): yes, and the class is large enough to include every strictly increasing rescaling of an L-convex function — exactly mirroring chunk 07's result for the M-convex side, and completing the "quasi" theory on both halves of the exchange-axiom framework before chapter 8 unifies them under conjugacy.

Setting

Let VVV be a finite ground set and g:ZV→R∪{+∞}g : \mathbb Z^V \to \mathbb R \cup \{+\infty\}g:ZV→R∪{+∞}. Building on chunk 08's submodularity axiom (SBF[Z]), this section introduces four ordinal relaxations. ggg is quasi submodular, satisfying (QSB), if for every p,q∈ZVp, q \in \mathbb Z^Vp,q∈ZV, g(p∧q)≤g(p)g(p \wedge q) \le g(p)g(p∧q)≤g(p) or g(p∨q)≤g(q)g(p \vee q) \le g(q)g(p∨q)≤g(q). ggg is semistrictly quasi submodular, satisfying (SSQSB), if additionally g(p∨q)≥g(q)  ⟹  g(p∧q)≤g(p)g(p \vee q) \ge g(q) \implies g(p \wedge q) \le g(p)g(p∨q)≥g(q)⟹g(p∧q)≤g(p) and symmetrically. The weak variants (QSBw) and (SSQSBw) restrict attention to points of the effective domain and compare max⁡(g(p),g(q))\max(g(p), g(q))max(g(p),g(q)) against min⁡(g(p∧q),g(p∨q))\min(g(p \wedge q), g(p \vee q))min(g(p∧q),g(p∨q)) directly, with (SSQSBw) additionally allowing the four-way tie g(p)=g(q)=g(p∧q)=g(p∨q)g(p) = g(q) = g(p\wedge q) = g(p \vee q)g(p)=g(q)=g(p∧q)=g(p∨q). The linear perturbation of ggg by x:V→Rx : V \to \mathbb Rx:V→R is g[x](p)=g(p)+⟨p,x⟩g[x](p) = g(p) + \langle p, x \rangleg[x](p)=g(p)+⟨p,x⟩.

Formalization targets

Goal: Theorem 7.54 (the quasi L-proximity theorem)

Let ggg satisfy (SSQSB) and g(p)=g(p+1)g(p) = g(p + \mathbf 1)g(p)=g(p+1) for all ppp, n=∣V∣n = |V|n=∣V∣, α\alphaα a positive integer. If pα∈dom⁡gp_\alpha \in \operatorname{dom} gpα​∈domg satisfies g(pα)≤g(pα+αχY)g(p_\alpha) \le g(p_\alpha + \alpha \chi_Y)g(pα​)≤g(pα​+αχY​) for all Y⊆VY \subseteq VY⊆V, then arg⁡min⁡g≠∅\arg\min g \ne \emptysetargming=∅ and there is p∗∈arg⁡min⁡gp^* \in \arg\min gp∗∈argming with the componentwise bound pα≤p∗≤pα+(n−1)(α−1)1p_\alpha \le p^* \le p_\alpha + (n-1)(\alpha-1) \mathbf 1pα​≤p∗≤pα​+(n−1)(α−1)1 — verbatim the same conclusion, and the same exact bound, as chunk 08's Theorem 7.18(1), now established for the strictly larger class satisfying (SSQSB) rather than (SBF[Z]).

Milestones: Theorems 7.49, 7.53

Theorem 7.49: the full nesting chain (SBF[Z]) ⇒\Rightarrow⇒ (SSQSB) ⇒\Rightarrow⇒ (QSB), (SSQSB) ⇒\Rightarrow⇒ (SSQSBw) ⇒\Rightarrow⇒ (QSBw), together with the collapse theorem that (SBF[Z]) holds if and only if every linear perturbation of ggg satisfies (QSBw) — precisely quantifying how weak (QSBw) is pointwise and how the classes reunite under universal perturbation. Theorem 7.53 (the quasi L-optimality criterion): the direct analogue of chunk 08's Theorem 7.14, showing that global (or, for the weaker (QSBw) case, unique-up-to-translation) optimality still reduces to a purely local check against the 2n−22^n - 22n−2 nontrivial sign-pattern neighbors p+χXp + \chi_Xp+χX​.

Significance

The result itself. As with the M-convex case (chunk 07), the proximity theorem is what algorithms actually need: an L-convex-flavored objective transformed by any strictly increasing scalar rescaling (a common device — expressing a network-flow cost in a different currency, or applying a monotone risk adjustment) retains a scaling algorithm's correctness guarantee with exactly the same distance bound, even though the rescaled function is generally no longer L-convex itself.

Formalizing it. No matching item exists on the platform for quasi submodularity or quasi L-convexity in any form. Together with chunk 07 (the M-side quasi-convexity theory), this mission completes the "quasi" relaxation on both halves of the exchange-axiom framework the book develops, immediately before chapter 8 unifies M-convexity and L-convexity under a single conjugacy relationship.

Difficulty

As with chunk 07's quasi M-proximity theorem, the temptation is to imitate chunk 08's L-proximity proof line by line. The overall architecture does survive — translate so pα=0p_\alpha = 0pα​=0, find a lattice-minimal sufficiently-good point, and bound the gap using submodularity — but chunk 08's proof uses (SBF[Z])'s additive inequality directly to compare four function values at once, while this proof must instead route every such comparison through (SSQSB)'s two one-directional implications (Proposition 7.50's quasi-version of the same two-sided inequality), which only ever license moving in one direction at a time depending on which side of a comparison is tight. The book's proof handles this by working with the specific implications (7.43)–(7.44) in place of the L♮-approach property used in chunk 08's proof — an ordinal substitute for the same additive step, at the cost of a case analysis chunk 08's proof did not need.

Formalization scope

This mission builds directly on chunk 08's published items (SBF, DomZ, ArgMin, IndicatorVec), per the platform's textbook convention that a later chapter section of the same book imports an earlier one's definitions; its own namespace DiscreteConvex.LConvexFunctions.Quasi nests under chunk 08's DiscreteConvex.LConvexFunctions accordingly. Note the sign convention of the linear perturbation here, g[x](p)=g(p)+⟨p,x⟩g[x](p) = g(p) + \langle p,x\rangleg[x](p)=g(p)+⟨p,x⟩, is the opposite of the M-side's f[p](x)=f(x)−⟨p,x⟩f[p](x) = f(x) - \langle p,x\ranglef[p](x)=f(x)−⟨p,x⟩ (chunks 06–07) — verified against the book's own formula rather than assumed by analogy.

A trivializing formalization of the goal would silently strengthen (SSQSB) back to plain (SBF[Z]) (making this mission redundant with chunk 08's Theorem 7.18) or loosen the exact bound (n−1)(α−1)(n-1)(\alpha-1)(n−1)(α−1); neither is done. Only (QSB), (SSQSB), (QSBw), (SSQSBw) are drafted, matching exactly what the chosen three items need; the polyhedral L-convex-function bridge (§7.8–7.9, Theorems 7.40–7.46) and the level-set characterizations (Theorems 7.51–7.52) are left for a follow-on mission. Contributions building the 0L ↔ S correspondence (Theorem 7.40, a bridge back to chunk 04's submodular-set-function vocabulary) or the scaled quasi L-minimizer-cut analogue are welcome.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • P. Milgrom, C. Shannon, "Monotone comparative statics," Econometrica, 62(1), 1994, pp. 157–180.
8 thms3 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis XXV: L-Convex Functions via Minimizer PolyhedraTextbook

Motivation

Chapter 5 characterized L-convex sets — sublattice-closed, translation-periodic subsets of ZV\mathbb Z^VZV — and showed they interact cleanly with integral convexity. Chapter 7 asks the functional analogue: which functions on the integer lattice deserve to be called convex in the "L" sense, and how do they relate back to L-convex sets? This mission (continuing mission 08-lconvex-functions-i, which built the axioms (SBF[Z])/(TRF[Z])/(SBF♮^\natural♮[Z]) and proved the L-optimality and L-proximity theorems) answers the second question at its sharpest: an L-convex function is exactly a function whose every weighted-minimizer set is an L-convex polyhedron — the discrete analogue of the fact that a convex function is determined by the convex geometry of its sublevel sets. Along the way it settles the chapter's basic toolkit: operations that preserve L-convexity, the local nature of submodularity, the correspondence with ordinary submodular set functions, and the first two structural facts about the convex extension every L-convex function admits.

Setting

Fix a finite ground set VVV. A function g:ZV→R∪{+∞}g : \mathbb Z^V \to \mathbb R \cup \{+\infty\}g:ZV→R∪{+∞} with nonempty effective domain is L-convex, g∈L[Z→R]g \in L[\mathbb Z \to \mathbb R]g∈L[Z→R], if it satisfies submodularity (SBF[Z]): g(p)+g(q)≥g(p∨q)+g(p∧q)g(p)+g(q) \ge g(p\vee q) + g(p\wedge q)g(p)+g(q)≥g(p∨q)+g(p∧q), and translation invariance (TRF[Z]): ∃r∈R\exists r \in \mathbb R∃r∈R, g(p+1)=g(p)+rg(p+\mathbf 1) = g(p) + rg(p+1)=g(p)+r for all ppp. It is L♮^\natural♮-convex if its lift to one extra coordinate (Eq. (7.2)) is L-convex. A set function ρ:2V→R∪{+∞}\rho : 2^V \to \mathbb R \cup \{+\infty\}ρ:2V→R∪{+∞} is submodular if ρ(X)+ρ(Y)≥ρ(X∪Y)+ρ(X∩Y)\rho(X)+\rho(Y) \ge \rho(X\cup Y) + \rho(X\cap Y)ρ(X)+ρ(Y)≥ρ(X∪Y)+ρ(X∩Y); it corresponds to an L♮^\natural♮-convex function supported on {0,1}V\{0,1\}^V{0,1}V via g(χX)=ρ(X)g(\chi_X) = \rho(X)g(χX​)=ρ(X) (Eq. (7.5)). The convex closure gˉ\bar ggˉ​ of ggg is its extension to RV\mathbb R^VRV by finite convex combinations.

Formalization targets

Goal: L-convexity via minimizer polyhedra (Theorem 7.17)

For g:ZV→R∪{+∞}g : \mathbb Z^V \to \mathbb R \cup \{+\infty\}g:ZV→R∪{+∞} with bounded nonempty effective domain: ggg is L-convex if and only if arg⁡min⁡g[−x]\arg\min g[-x]argming[−x] is an L-convex set for every x∈RVx \in \mathbb R^Vx∈RV; the L♮^\natural♮ analogue holds with L♮^\natural♮-convex sets. This is the direct mirror of mission 23-ch06c-mconvexfunctions's own goal (Theorem 6.43, characterizing M-convex functions via M-convex weighted-minimizer polyhedra) — the book's text calls it exactly "how the concept of L-convex functions can be defined from that of L-convex sets."

Supporting structural targets

Eleven further results build the chapter's basic vocabulary. Theorem 7.2 strengthens translation submodularity to allow negative shifts; Theorem 7.3 places L-convexity inside L♮^\natural♮- convexity; Proposition 7.4 identifies submodular set functions with a subclass of L♮^\natural♮-convex functions via the indicator embedding, and Theorem 7.15 derives the classical submodular-minimizer local-optimality criterion as its corollary; Proposition 7.5 shows submodularity is a local property, needing only unit-distance pairs; Proposition 7.8 transfers L-(natural-)convexity from functions to their effective domains; Proposition 7.9 and Theorem 7.10–7.11 give the chapter's basic examples (univariate and pairwise-difference functions) and its six-operation closure toolkit (scaling, affine reparametrization, linear perturbation, projection, infimal convolution with a separable function, and sums), both for L-convex and L♮^\natural♮- convex functions, the latter also admitting interval and coordinate restrictions; Proposition 7.16 shows minimizer sets of L-convex functions are themselves L-convex, the special case (x=0x=0x=0) the goal generalizes to every linear perturbation; Theorem 7.19 begins the convex-extension program this chapter's next chunk completes, establishing that the convex closure agrees with ggg on ZV\mathbb Z^VZV and inherits its translation constant.

Significance

The goal is significant for the same structural reason as its M-side counterpart: it says L-convexity is not merely a combinatorial condition on lattice differences but is equivalent to a purely polyhedral-geometric one, closing the loop between chapters 5 and 7 the way Theorem 6.43 closes the loop between chapters 4 and 6. Theorem 7.15's corollary status is itself instructive: the well-known fact that a submodular set function's global minimizer needs only local verification against comparable sets — the theoretical basis of every submodular-minimization algorithm in chapter 10 — falls out of the L-optimality criterion (mission 08-lconvex-functions-i's Theorem 7.14) applied to the indicator embedding, rather than needing an independent proof. Theorem 7.10–7.11's six operations are the toolkit every later construction in this chapter and chapter 9's network transformations builds new L-convex functions from old.

None of these results are open — they are Murota's account of the basic function-level theory of L-convexity, mirroring chapter 6's M-convex function theory chunk-by-chunk. What this mission contributes is a faithful, machine-checked formal statement of each, extending the shared Lean vocabulary (SBF, TRF, LNaturalConvex, LConvexSet) that missions 08-lconvex-functions-i and 21-ch05b-lconvexsets began; no comparable formalization exists on the platform (see Formalization scope).

Difficulty

The naive approach to the goal would try to verify L-convexity's submodularity inequality directly against the definition of an L-convex set applied to each minimizer family; the book's actual proof instead routes through Theorem 7.10 (3)'s closure of L-convexity under linear perturbation and Proposition 7.16's minimizer-is-L-convex-set fact for the forward direction, and defers the converse entirely to a later note (Note 7.47, outside this chunk and mission 08's combined range) proved via the integral-convexity machinery of section 7.7 onward. The genuine combinatorial difficulty in this block is upstream, in Theorem 7.10 (5)'s infimal-convolution operation: proving L-convexity of the perturbed function requires a four-term submodularity inequality assembled from the separable function's own convexity and ggg's submodularity applied at the optimal q1,q2q_1,q_2q1​,q2​ simultaneously — a genuine two-hypothesis combination with no single-inequality shortcut.

Formalization scope

Ground-set elements are a Fintype V with DecidableEq; L-(natural-)convex functions are (V→ℤ)→WithTop ℝ. All twelve numbered results found in this chunk's page range are placed, with one documented scope reduction: Theorem 7.19 states only parts (3)-(4) (that the convex closure agrees with ggg on ZV\mathbb Z^VZV and inherits its translation constant), not the explicit Lovász-extension-formula construction of parts (1)-(2) and (5), which needs a sorted-distinct- component apparatus no other result in this chunk requires — see HARD.md. "gU>−∞g_U > -\inftygU​>−∞" and its variants are replaced by the equivalent (DomZ ...).Nonempty hypothesis throughout, matching mission 24-ch06d-mconvexfunctions's identical substitution (WithTop ℝ has no −∞-\infty−∞ element). This mission's base vocabulary (SBF, TRF, LNaturalConvex, etc.) is redeclared verbatim from mission 08-lconvex-functions-i rather than imported, since sibling drafts in this series cannot yet reference one another; LConvexSet is likewise redeclared from mission 21-ch05b-lconvexsets. Contributions completing any of the twelve sorrys are welcome; the goal and Theorem 7.10 carry the most independent proof content.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • K. Murota, "Discrete convex analysis," Mathematical Programming, 83 (1998), pp. 313–371 (the original account of L-convex functions this chapter's basic theory is drawn from).
34 thms3 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis XXVI: Polyhedral L-Convex FunctionsTextbook

Motivation

L-convex functions were defined purely combinatorially, on the integer lattice. Chapter 6's M-convex theory showed that combinatorial definition always extends to a genuine convex function on real space (missions 23-ch06c-mconvexfunctions/24-ch06d-mconvexfunctions); this mission carries out the identical program for the L side. It first shows L♮^\natural♮-convexity is exactly integral convexity plus ordinary submodularity — a clean synonym that also explains why submodular set functions are a natural special case — then builds the entire polyhedral (real-variable) theory of L-convex functions: the axioms (SBF[R])/(TRF[R]), two practical local criteria for verifying submodularity without checking every pair of points, the fact that the classical Lovász extension of a submodular set function is itself a polyhedral L-convex function, the six-operation closure toolkit, and — this mission's goal — the L-optimality criterion in its full polyhedral generality, characterizing global optimality by finitely many directional derivatives.

Setting

Fix a finite ground set VVV. A polyhedral convex function g:RV→R∪{+∞}g : \mathbb R^V \to \mathbb R \cup \{+\infty\}g:RV→R∪{+∞} with nonempty effective domain is polyhedral L-convex, g∈L[R→R]g \in L[\mathbb R \to \mathbb R]g∈L[R→R], if it satisfies (SBF[R]): g(p)+g(q)≥g(p∨q)+g(p∧q)g(p)+g(q) \ge g(p\vee q)+g(p\wedge q)g(p)+g(q)≥g(p∨q)+g(p∧q), and (TRF[R]): ∃r∈R\exists r \in \mathbb R∃r∈R, g(p+α1)=g(p)+αrg(p+\alpha\mathbf 1) = g(p)+\alpha rg(p+α1)=g(p)+αr for all p∈RVp \in \mathbb R^Vp∈RV, α∈R\alpha \in \mathbb Rα∈R; it is polyhedral L♮^\natural♮-convex if its lift to one extra real coordinate is polyhedral L-convex. A function g:ZV→R∪{+∞}g : \mathbb Z^V \to \mathbb R \cup \{+\infty\}g:ZV→R∪{+∞} is integrally convex if its convex closure agrees, at every real point, with the closure taken using only that point's integral neighborhood. The Lovász extension ρ^\hat\rhoρ^​ of a submodular set function ρ\rhoρ is the piecewise-linear interpolation built from the sorted distinct components of p∈RVp \in \mathbb R^Vp∈RV. The directional derivative g′(p;d)g'(p;d)g′(p;d) is inf⁡t>0(g(p+td)−g(p))/t\inf_{t>0}(g(p+td)-g(p))/tinft>0​(g(p+td)−g(p))/t.

Formalization targets

Goal: the polyhedral L-optimality criterion (Theorem 7.33)

For a polyhedral L-convex function ggg and p∈dom⁡Rgp \in \operatorname{dom}_{\mathbb R} gp∈domR​g: g(p)≤g(q)g(p) \le g(q)g(p)≤g(q) for all qqq if and only if g′(p;χY)≥0g'(p;\chi_Y) \ge 0g′(p;χY​)≥0 for every Y⊆VY \subseteq VY⊆V and g′(p;1)=0g'(p;\mathbf 1)=0g′(p;1)=0; for polyhedral L♮^\natural♮-convex ggg, the criterion simplifies to g′(p;±χY)≥0g'(p;\pm\chi_Y) \ge 0g′(p;±χY​)≥0 for every YYY. This is the direct L-side mirror of mission 24-ch06d-mconvexfunctions's M-optimality criterion (Theorem 6.52) and the polyhedral generalization of mission 08-lconvex-functions-i's integer-domain L-optimality criterion (Theorem 7.14): checking global optimality against exponentially many points reduces to ∣V∣+1|V|+1∣V∣+1 (or 2∣V∣2|V|2∣V∣) directional-derivative inequalities.

Supporting structural targets

Twelve further results build the polyhedral theory from the ground up. Theorems 7.20-7.21 identify L♮^\natural♮-convexity with the conjunction of ordinary submodularity and integral convexity — a genuinely different, function-analytic characterization from the exchange-axiom- style definitions used so far. Propositions 7.23-7.24 give two practical sufficient conditions for verifying (SBF[R]) locally, at a single scale, rather than globally. Proposition 7.25 shows the Lovász extension of any submodular set function is automatically polyhedral L-convex — not in this chunk's own extraction table (its label is preceded by an unlabeled restatement of the same fact, which evidently confused the extractor), found and placed by direct reading. Theorem 7.26 shows an L-convex function's convex extension, when polyhedral, inherits polyhedral L-convexity, continuing mission 26-ch07b-lconvexfunctions's Theorem 7.19. Theorems 7.28-7.32 restate the discrete theory's core equivalences (translation submodularity, the L/L♮^\natural♮ correspondence, the six basic operations, restrictions) in the polyhedral setting, and Proposition 7.34 shows minimizer sets of linearly-perturbed polyhedral L-convex functions are themselves L-convex polyhedra — flagged by the book itself as a partial result whose full converse characterization (Theorem 7.45) lies beyond this chunk's range.

Significance

Theorems 7.20-7.21's synonym is structurally important: it means every algorithm and theorem already known for submodular-function minimization over {0,1}V\{0,1\}^V{0,1}V-type domains applies, after a midpoint-convexity check, to the vastly larger class of integer-lattice L♮^\natural♮-convex functions, with no new proof technique required. Proposition 7.25 is the bridge that lets the combinatorial Lovász extension — the workhorse of submodular optimization for forty years — be recognized as a special case of the polyhedral L-convex function theory this mission builds, explaining why algorithms for one transfer so readily to the other. The goal, Theorem 7.33, is the precise tool chapter 10's continuous-relaxation algorithms for L-convex-function minimization actually verify against: a scaling algorithm's claimed optimum is confirmed correct exactly by checking the criterion's finitely many directional-derivative inequalities.

None of these results are open — they are Murota's account of how the integer-lattice theory of L-convexity survives, result by result, the passage to polyhedral convex functions on RV\mathbb R^VRV, mirroring chapter 6's identical program for M-convexity. What this mission contributes is a faithful, machine-checked formal statement of each, including one result (Proposition 7.25) the platform's own automated extractor missed, extending the shared Lean vocabulary (SBFR, TRFR, LovaszExtension, DirDeriv) this series builds on; no comparable formalization exists on the platform (see Formalization scope).

Difficulty

The naive approach to the goal would try to verify g(p)≤g(q)g(p) \le g(q)g(p)≤g(q) directly against every q∈RVq \in \mathbb R^Vq∈RV; the book's actual proof instead reduces this to the finite family of directional derivatives via Theorem 7.20's integral-convexity fact (an L♮^\natural♮-convex function's local behavior determines its global behavior) applied to the polyhedral setting through Theorem 3.21's general optimality criterion for integrally convex functions — a two-layer reduction (polyhedral →\to→ integral-convexity →\to→ finite local check) with no direct one-step argument. The genuine combinatorial content in this block is in Proposition 7.25's proof: showing the Lovász extension is submodular requires the finite-valued case (a direct calculation split on whether the two perturbed coordinates land in the same or different threshold sets) and then a limiting argument over a sequence of finite-valued truncations ρk→ρ\rho_k \to \rhoρk​→ρ for the general, possibly-infinite case — a genuine two-step argument, not a single inequality chase.

Formalization scope

Ground-set elements are a Fintype V with DecidableEq; polyhedral L-(natural-)convex functions are (V→ℝ)→WithTop ℝ. All thirteen numbered results found in this chunk's page range are placed, with one documented scope reduction: Proposition 7.24 states only part (1) (the unconditional-on-magnitude sufficient condition), not part (2)'s sharper, sorted-index-restricted version, which needs the same SortedValues apparatus a second time for no other result's benefit — see HARD.md. "gU>−∞g_U > -\inftygU​>−∞" and its variants are replaced by the equivalent (DomR ...).Nonempty hypothesis throughout, matching mission 24-ch06d-mconvexfunctions's identical substitution. "Domain is closed"/"domain is an interval" (Propositions 7.23-7.24) are stated via Mathlib's IsClosed and Set.OrdConnected respectively, the latter being the precise order-theoretic notion of "interval" in a pointwise-ordered space. This mission's base vocabulary is redeclared from missions 08-lconvex-functions-i, 20-ch04b-mconvexsets (for the Lovász extension machinery), 21-ch05b-lconvexsets, 23-ch06c-mconvexfunctions/ 24-ch06d-mconvexfunctions, and 26-ch07b-lconvexfunctions rather than imported, since sibling drafts in this series cannot yet reference one another. Contributions completing any of the thirteen sorrys are welcome; the goal and Proposition 7.25 carry the most independent proof content.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • K. Murota and A. Shioura, "Extreme points of a generalized polymatroid," Discrete Applied Mathematics, 152 (2005), pp. 268-278 [152] (the polyhedral L-convex function theory this mission's real-variable results are drawn from).
50 thms3 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis XXVII: Directional Derivatives and Quasi L-Convex FunctionsTextbook

Motivation

This mission completes chapter 7's L-convex function theory and closes the book on it. It first finishes the theory of positively homogeneous L-convex functions — showing they coincide exactly with the classical Lovász extensions of submodular set functions, a one-to-one correspondence that recognizes forty years of submodular-optimization machinery as a special case of L-convex function theory. It then proves the L-side capstone this series has been building toward since mission 26-ch07b-lconvexfunctions: polyhedral L-convexity is characterized simultaneously by directional derivatives, subdifferentials, and weighted-minimizer polyhedra — the exact mirror of what mission 24-ch06d-mconvexfunctions proved for M-convex functions. Finally it develops quasi L-convex functions, the L-side analogue of mission 25-ch06e-mconvexfunctions's quasi M-convex functions, ending exactly where chapter 7 itself ends.

Setting

Fix a finite ground set VVV. The class 0L[R→R]0L[\mathbb R \to \mathbb R]0L[R→R] consists of polyhedral L-convex functions that are positively homogeneous; 0L[Z→Z]0L[\mathbb Z \to \mathbb Z]0L[Z→Z], its integer-valued integer-domain analogue. A positively homogeneous L-convex function ggg induces a submodular set function ρg(X)=g(χX)\rho_g(X) = g(\chi_X)ρg​(X)=g(χX​); conversely the Lovász extension ρ^\hat\rhoρ^​ of a submodular set function is positively homogeneous L-convex. The base polyhedron B(ρ)={x∈RV:x(X)≤ρ(X) ∀X, x(V)=ρ(V)}B(\rho) = \{x \in \mathbb R^V : x(X) \le \rho(X)\ \forall X,\ x(V) = \rho(V)\}B(ρ)={x∈RV:x(X)≤ρ(X) ∀X, x(V)=ρ(V)} of a submodular ρ\rhoρ is always an M-convex polyhedron. A function g:ZV→R∪{+∞}g : \mathbb Z^V \to \mathbb R \cup \{+\infty\}g:ZV→R∪{+∞} is quasi submodular (QSB) if g(p∧q)≤g(p)g(p\wedge q) \le g(p)g(p∧q)≤g(p) or g(p∨q)≤g(q)g(p\vee q) \le g(q)g(p∨q)≤g(q) for all p,qp,qp,q; the weaker (QSBw) requires only max⁡{g(p),g(q)}≥min⁡{g(p∧q),g(p∨q)}\max\{g(p),g(q)\} \ge \min\{g(p\wedge q), g(p\vee q)\}max{g(p),g(q)}≥min{g(p∧q),g(p∨q)}.

Formalization targets

Goal: the four characterizations of polyhedral L-convexity (Theorem 7.45)

For a polyhedral convex function ggg with dom⁡Rg≠∅\operatorname{dom}_{\mathbb R} g \ne \emptysetdomR​g=∅: ggg is L-convex if and only if every directional derivative g′(p;⋅)g'(p;\cdot)g′(p;⋅) is 0L[R→R]0L[\mathbb R \to \mathbb R]0L[R→R], if and only if every subdifferential ∂Rg(p)\partial_{\mathbb R} g(p)∂R​g(p) is an M-convex polyhedron, if and only if every weighted minimizer set is an L-convex polyhedron. Not present in this chunk's own extraction table (its label opens mid-paragraph, missed by the same extractor failure already documented for mission 25-ch06e-mconvexfunctions's Theorem 6.68), found by direct reading and chosen as goal because it is the exact L-side mirror of mission 24-ch06d-mconvexfunctions's own milestone Theorem 6.63, and its proof is assembled entirely from results already in this series (Theorem 7.43, Proposition 7.34, Theorem 7.40).

Supporting structural targets

Fifteen further results build the two remaining pieces of chapter 7's theory. Propositions 7.37-7.39 and Theorem 7.40 establish the one-to-one correspondence between positively homogeneous L-convex functions and submodular set functions via the Lovász extension; Proposition 7.41 gives a minimizer-polyhedron characterization of this class, and Proposition 7.42 shows directional derivatives of L-convex functions automatically land in it. Theorem 7.43 — the L-side mirror of mission 24-ch06d-mconvexfunctions's own goal, Theorem 6.61 — proves the directional- derivative/subdifferential correspondence via the induced submodular set function's base polyhedron; Proposition 7.44 checks consistency at integer points, and Theorem 7.46 refines Theorem 7.45 to the integral case. Theorem 7.49 (also missed by the extractor) gives the quasi-submodularity implication hierarchy and its perturbation-equivalence capstone, mirroring mission 25-ch06e-mconvexfunctions's Theorem 6.68 exactly. Proposition 7.50 and Theorems 7.51-7.52 build the level-set/perturbation machinery quasi submodularity needs; Theorems 7.53-7.54 (both missed by the extractor, the latter's full statement requiring one page beyond this chunk's nominal range, at the very end of chapter 7) give the quasi L-optimality and quasi L-proximity theorems.

Significance

The 0L/submodular correspondence (Theorem 7.40) is the precise sense in which L-convex function theory generalizes submodular set function theory rather than merely resembling it: every submodular set function is literally the restriction to {0,1}V\{0,1\}^V{0,1}V of a positively homogeneous L-convex function, and every algorithm for one transfers to the other through this exact dictionary. The goal, Theorem 7.45, completes the parallel structure this series has built since chapter 6: M-convexity and L-convexity are now each characterized in the same four convex-analytic vocabularies, setting up chapter 8's conjugacy theorem, which will show these two characterizations are not merely analogous but literally dual to each other under the Legendre-Fenchel transform. The quasi-submodularity results matter for the same reason as their M-side counterparts: chapter 10's algorithms for L-convex-function minimization remain correct under nonlinear rescalings that destroy L-convexity itself but preserve quasi submodularity.

None of these results are open — they are Murota's account of how far L-convex function theory extends beyond the polyhedral case (to positive homogeneity and its submodular-function incarnation) and how far its exchange-style inequality can be relaxed while preserving optimization theory (to quasi submodularity), mirroring chapter 6's identical two-part program for M-convex functions. What this mission contributes is a faithful, machine-checked formal statement of each, including four theorems (7.45, 7.49, 7.53, 7.54) the platform's own automated extractor missed entirely — one of them requiring a page beyond this chunk's own nominal range to complete, since chapter 7 ends there and this is the last mission covering it — extending the shared Lean vocabulary (ZeroLR, BasePolyhedron, QSBw) this series builds on; no comparable formalization exists on the platform (see Formalization scope).

Difficulty

The naive approach to the goal would try to prove all six pairwise implications among its four conditions independently; the book's own proof instead chains through results already established: (a)⇒(b) is Proposition 7.42, (a)⇒(c) is Theorem 7.43, (a)⇒(d) is Proposition 7.34, (b)⇔(c) uses the 0L/M0[R] correspondence, and (d)⇒(b) is the genuinely hard direction, requiring Proposition 7.41 applied to the directional derivative itself (showing arg⁡min⁡(g′(p;⋅)[−x])\arg\min(g'(p;\cdot)[-x])argmin(g′(p;⋅)[−x]) is an L-convex cone by an explicit description via the admissible-potential set of the distance function underlying arg⁡min⁡g[−x]\arg\min g[-x]argming[−x]). The remaining combinatorial difficulty in this block is in Theorem 7.43's proof: identifying ∂Rg(p)\partial_{\mathbb R} g(p)∂R​g(p) with the base polyhedron B(ρg,p)B(\rho_{g,p})B(ρg,p​) requires the L-optimality criterion (Theorem 7.33, mission 27-ch07c-lconvexfunctions) applied pointwise, a chain of logical equivalences with no single-step shortcut, exactly mirroring how mission 24-ch06d-mconvexfunctions's Theorem 6.61 needed the M-optimality criterion.

Formalization scope

Ground-set elements are a Fintype V with DecidableEq; L-convex functions are (V→ℝ)→WithTop ℝ (polyhedral) or (V→ℤ)→WithTop ℝ (integer-domain). All sixteen numbered results found in this chunk's page range — the twelve in BRIEF.md's own table plus four the extractor missed (Theorems 7.45, 7.49, 7.53, 7.54) — are placed, with two documented, content-preserving scope decisions: Theorem 7.43 omits the dual-integral refinement clauses for L[R→R|Z]/L[Z→Z] (the same decision mission 24-ch06d-mconvexfunctions's Theorem 6.61 made), and Proposition 7.50 states the general inequalities without restating their "In particular" specializations, which add no independent content — see HARD.md. "inf⁡g[−x]>−∞\inf g[-x]>-\inftyinfg[−x]>−∞" is replaced by the equivalent (ArgMinR ...).Nonempty hypothesis throughout, matching mission 25-ch06e-mconvexfunctions's identical substitution. This mission's base vocabulary is redeclared from missions 20-ch04b-mconvexsets, 21-ch05b-lconvexsets, 23-24-ch06*-mconvexfunctions, and 26-27-ch07*-lconvexfunctions rather than imported, since sibling drafts in this series cannot yet reference one another. Contributions completing any of the sixteen sorrys are welcome; the goal and Theorem 7.43 carry the most independent proof content.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • K. Murota and A. Shioura, "M-convex function on generalized polymatroid," Mathematics of Operations Research, 24 (1999), pp. 95-105 [152] (the polyhedral theory Theorems 7.26-7.46 are drawn from).
  • P. Milgrom and C. Shannon, "Monotone comparative statics," Econometrica, 62 (1994), pp. 157-180 [129] (the origin of the quasi-submodularity condition (SSQSB)).
68 thms3 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis X: The Lagrangian Saddle-Point TheoremTextbook

Motivation

Chunk 10 formalized the discrete conjugacy theorem — the Legendre-Fenchel transform's bijection between the classes of M-convex and L-convex functions — and, along the way, a function-level generalization of Edmonds's intersection theorem. Section 8.4 turns that machinery toward a different question: not "how are two convexity classes related," but "when does a discrete optimization problem have a dual that meets it with equality." The classical route to such strong-duality results in continuous convex programming — Lagrangian relaxation, an embedding of the problem in a family of perturbed problems, and a saddle-point characterization of when primal and dual values coincide — has a discrete analogue that needs no continuity, no differentiability, and no convexity in the classical sense at all: only the elementary fact that the Legendre-Fenchel transform, once discretized, is still an involution on the right class of functions. This mission formalizes that discrete Lagrangian duality framework and its central saddle-point theorem, in full generality — before the book specializes it, in the section that follows, to the specific M-convex perturbation that gives the chapter's headline strong-duality result for M-convex programs.

Setting

Let VVV and UUU be finite ground sets. A perturbation of an optimization problem min⁡{f(x):x∈ZV}\min\{f(x) : x \in \mathbb Z^V\}min{f(x):x∈ZV} is a function F:ZV×ZU→Z∪{+∞}F : \mathbb Z^V \times \mathbb Z^U \to \mathbb Z \cup \{+\infty\}F:ZV×ZU→Z∪{+∞} such that F(x,0)=f(x)F(x,0) = f(x)F(x,0)=f(x) for all xxx (Eq. (8.54)) and, for each fixed xxx, F(x,⋅)F(x,\cdot)F(x,⋅) is self-biconjugate: F(x,⋅)∙∙=F(x,⋅)F(x,\cdot)^{\bullet\bullet} = F(x,\cdot)F(x,⋅)∙∙=F(x,⋅) under the discrete Legendre-Fenchel transform of chunk 10 (Eq. (8.55)). The Lagrangian function is K(x,y)=inf⁡{F(x,u)+⟨u,y⟩:u∈ZU}K(x,y) = \inf\{F(x,u) + \langle u,y\rangle : u \in \mathbb Z^U\}K(x,y)=inf{F(x,u)+⟨u,y⟩:u∈ZU} (Eq. (8.58)), valued in Z∪{±∞}\mathbb Z \cup \{\pm\infty\}Z∪{±∞} (formalized in EReal, since both the infimum and the supremum below can be genuinely unbounded). The dual objective is g(y)=inf⁡{K(x,y):x∈ZV}g(y) = \inf\{K(x,y) : x \in \mathbb Z^V\}g(y)=inf{K(x,y):x∈ZV} (Eq. (8.60)). Writing inf⁡(P)=inf⁡xf(x)\inf(P) = \inf_x f(x)inf(P)=infx​f(x), sup⁡(D)=sup⁡yg(y)\sup(D) = \sup_y g(y)sup(D)=supy​g(y), opt⁡(P)={x:f(x)=inf⁡(P)}\operatorname{opt}(P) = \{x : f(x) = \inf(P)\}opt(P)={x:f(x)=inf(P)}, opt⁡(D)={y:g(y)=sup⁡(D)}\operatorname{opt}(D) = \{y : g(y) = \sup(D)\}opt(D)={y:g(y)=sup(D)}, the primal problem PPP is to minimize fff over ZV\mathbb Z^VZV and the dual problem DDD is to maximize ggg over ZU\mathbb Z^UZU.

Formalization targets

Goal: Theorem 8.54 (the saddle-point theorem)

Assuming FFF is self-biconjugate (Eq. (8.55)): both inf⁡(P)\inf(P)inf(P) and sup⁡(D)\sup(D)sup(D) are finite and min⁡(P)=max⁡(D)\min(P) = \max(D)min(P)=max(D) if and only if there exist xˉ∈ZV\bar x \in \mathbb Z^Vxˉ∈ZV, yˉ∈ZU\bar y \in \mathbb Z^Uyˉ​∈ZU with K(xˉ,yˉ)K(\bar x,\bar y)K(xˉ,yˉ​) finite and K(x,yˉ)≤K(xˉ,yˉ)≤K(xˉ,y)K(x,\bar y) \le K(\bar x,\bar y) \le K(\bar x,y)K(x,yˉ​)≤K(xˉ,yˉ​)≤K(xˉ,y) for all x,yx,yx,y — a saddle point of the Lagrangian kernel. When this holds, xˉ∈opt⁡(P)\bar x \in \operatorname{opt}(P)xˉ∈opt(P) and yˉ∈opt⁡(D)\bar y \in \operatorname{opt}(D)yˉ​∈opt(D).

Milestones: Theorem 8.52, Proposition 8.51(1)-(2)

Theorem 8.52 (weak duality): inf⁡(P)≥sup⁡(D)\inf(P) \ge \sup(D)inf(P)≥sup(D) always, with no biconjugacy hypothesis on FFF at all — the baseline the saddle-point theorem sharpens to equality. Proposition 8.51(1)-(2): under self-biconjugacy, the perturbation FFF (and hence the primal objective fff) is itself recoverable from the Lagrangian kernel KKK by a supremum, F(x,u)=sup⁡y{K(x,y)−⟨u,y⟩}F(x,u) = \sup_y\{K(x,y) - \langle u,y\rangle\}F(x,u)=supy​{K(x,y)−⟨u,y⟩} and f(x)=sup⁡yK(x,y)f(x) = \sup_y K(x,y)f(x)=supy​K(x,y) — the algebraic identity the saddle-point theorem's proof turns on directly.

Significance

The result itself. The saddle-point theorem is the general-purpose engine behind every strong-duality result the book proves for specific classes of discrete optimization problems: the book's own next section specializes it (via a particular choice of FFF built from an M-convex regularizer rrr) to obtain strong duality for M-convex programs, but the theorem itself needs no M-convexity, no submodularity, and no exchange axiom — only the elementary self-biconjugacy of a perturbation under the discrete Legendre-Fenchel transform. It is, in that sense, the most general and most reusable strong-duality statement in the book: any future mission proving strong duality for a specific class of discrete programs (M-convex, M2-convex, network flow, or otherwise) by exhibiting a self-biconjugate perturbation can cite this theorem directly rather than reproving the saddle-point argument from scratch.

Formalizing it. No matching item exists on the platform for a discrete Lagrangian saddle- point theorem, discrete weak duality, or this perturbation-based duality framework. (A prior-art search turned up an unrelated continuous Lagrangian saddle-point theorem for convex cones, Luenberger's Chapter 8 §8.4, formalized as VectorSpaceOpt.lagrangian_saddle_sufficient_pointed — a genuinely different setting: no discreteness, no biconjugacy hypothesis, and a one-directional sufficiency statement rather than this mission's iff. Not reused.) This mission gives the first formal statement of discrete Lagrangian duality, and directly reuses chunk 10's ConvexConjugate apparatus (self-biconjugacy is stated using chunk 10's own conjugate-of-conjugate composition), demonstrating exactly the kind of shared-substrate payoff the discrete conjugacy theorem was built to provide.

Difficulty

The saddle-point theorem's "only if" direction is not a routine unwinding of definitions: given min⁡(P)=max⁡(D)\min(P) = \max(D)min(P)=max(D) at finite common value, one must construct the saddle point (xˉ,yˉ)(\bar x,\bar y)(xˉ,yˉ​) — the book's proof takes xˉ∈opt⁡(P)\bar x \in \operatorname{opt}(P)xˉ∈opt(P), yˉ∈opt⁡(D)\bar y \in \operatorname{opt}(D)yˉ​∈opt(D) (which exist because the infimum/supremum are attained at a finite optimum) and verifies the sandwiching inequality using Proposition 8.51(2)'s identity f(x)=sup⁡yK(x,y)f(x) = \sup_y K(x,y)f(x)=supy​K(x,y) together with weak duality, rather than by any direct algebraic manipulation of KKK alone. Skipping straight to a "trivial" biconditional that never invokes Proposition 8.51 would misrepresent the actual proof structure the book relies on for exactly this direction.

Formalization scope

V,UV, UV,U are Fintype ground types; LagrangianKernel, DualObjective, InfP, SupD are EReal-valued to keep both the defining infima/suprema total (a complete lattice) without an artificial finiteness side-condition; OptP, OptD compare PrimalValue/DualObjective against InfP/SupD after casting through chunk 10's ToEReal, mirroring that chunk's own round-trip convention. "Finite" throughout is formalized as ≠ ⊤ ∧ ≠ ⊥ in EReal. The perturbation FFF itself is left fully abstract (an arbitrary function satisfying the self-biconjugacy hypothesis where needed) — this mission does not draft the specific M-convex perturbation FrF_rFr​ (Eq. (8.61)) that the book's next subsection (§8.4.3) uses to specialize this framework to M-convex programs, nor Theorem 8.59 (the resulting M-convex strong-duality theorem) itself, which needs that specific perturbation plus its own regularity conditions (REG)/(OBJ) and a chain of M-convex-specific propositions (8.55–8.58) beyond what the general framework built here provides. A trivializing formalization would state the saddle-point theorem's sandwiching inequality with a weaker order (e.g., only one of the two directions) or would omit the "xˉ∈opt⁡(P),yˉ∈opt⁡(D)\bar x \in \operatorname{opt}(P), \bar y \in \operatorname{opt}(D)xˉ∈opt(P),yˉ​∈opt(D)" consequence clause; neither is done — both inequalities and the full consequence clause are included exactly as the book states them.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
10 thms3 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis XXXI: The Potential Criterion for Network FlowsTextbook

Motivation

Chapter 9 is where discrete convex analysis meets classical network flow theory: the minimum cost flow problem's three hallmark properties — an optimality criterion by potentials, an optimality criterion by negative cycles, and integrality of optimal solutions — are shown to survive, in a precise and increasingly general form, first for arbitrary polyhedral convex costs (MCFP3), then for the M-convex submodular flow problem (MSFP2/MSFP3), the chapter's own combinatorial generalization of the classical problem. This mission places the potential criterion (Theorem 9.4) and its cascade of six corollaries and generalizations, the block of results this book's own text uses to carry every other result in the chapter.

Setting

A digraph G = (V,A) with tail/head maps ∂⁺,∂⁻ : A → V. A flow ξ : A → R has boundary ∂ξ(v) = Σ{ξ(a) : ∂⁺a=v} − Σ{ξ(a) : ∂⁻a=v}. A potential p : V → R has coboundary δp(a) = p(∂⁺a) − p(∂⁻a). The minimum cost flow problem MCFP3 minimizes Γ₃(ξ) = Σₐ fₐ(ξ(a)) + f(∂ξ) over flows, for polyhedral convex arc costs fₐ : R → R∪{+∞} and boundary cost f : Rⱽ → R∪{+∞}; MCFP0 is its linear-cost, fixed-supply special case. The M-convex submodular flow problem MSFP3 is MCFP3 with f additionally M-convex; MSFP2 is its linear-arc-cost special case.

Formalization targets

Goal: The potential criterion for MCFP3 (Theorem 9.4)

For a feasible flow ξ, ξ is optimal for MCFP3 iff there is a potential p with ξ(a) a minimizer of the reduced arc cost fₐ[δp(a)] for every arc and ∂ξ a minimizer of the reduced boundary cost f[−p]; and any such optimal potential characterizes optimality of every feasible flow. This is the hub result of the whole chunk: the book states Theorem 9.14 is "immediate" from it, and every other placed result either specializes it directly or builds on that specialization.

Supporting structural targets

Theorem 9.5 reformulates MCFP0's optimality as the absence of a negative cycle in an auxiliary network; Theorem 9.6 gives MCFP0's primal and dual integrality, the latter identifying the optimal-potential set as an L-convex polyhedron. Theorem 9.14 specializes the goal to MSFP3; Theorem 9.15 upgrades this to a full polyhedral and integrality structure theorem for MSFP3's optimal-flow-boundary and optimal-potential sets (M2-convex and L-convex polyhedra respectively); Theorem 9.16 is the integer-flow analogue, with the boundary set now literally M2-convex and the integer-optimal-potential set literally L-convex. Theorems 9.18 and 9.20 give the negative-cycle reformulation for MSFP2, real and integer flows respectively, generalizing Theorem 9.5 by admitting a third class of auxiliary arcs governed by the M-convex boundary cost's directional derivative (or its discrete difference, in the integer case).

Significance

This is the chapter's demonstration that M-convexity is not merely an abstract combinatorial axiom but the exact structural hypothesis under which classical network-flow duality survives intact: every one of the four "nice properties" the book opens the chapter with (potentials, negative cycles, integrality, efficient algorithms) is preserved verbatim in the M-convex generalization, and this mission's eight results are the proof of that claim for the first three. The chunk's own internal dependency structure — one foundational theorem (9.4) from which every other placed result descends by specialization or direct generalization — is itself characteristic of how this book organizes its combinatorial machinery around a single convex- analytic core.

None of these results are open — they are Murota's own account of network flow duality under M-convexity (sections 9.1, 9.4, and 9.5). What this mission contributes is a faithful, machine-checked formal statement of each, extending the platform's coverage of chapter 9 begun in mission 12-network-flows (which covered §9.1.1-9.1.2 and §9.3, the feasibility and max-flow min-cut results, deliberately leaving this block for later apparatus); no comparable formalization exists on the platform (see Formalization scope).

Difficulty

The eight results span real- and integer-flow versions of two nested problem hierarchies (MCFP0 ⊂ MCFP3, MSFP2 ⊂ MSFP3) and two distinct optimality certificates (potentials, negative cycles), which this mission handles by building one shared apparatus — FeasibleFlowMCFP3, Gamma3, OptimalFlowMCFP3, IsOptimalPotential — that MCFP0 and MSFP3 both instantiate (MCFP0 literally as the linear-cost/singleton-boundary special case of Eq. (9.11)), and one shared generic cycle/negative-cycle apparatus (IsCycle, CycleLength, HasNegativeCycle) instantiated three times with different auxiliary-arc types (A⊕A for MCFP0, A⊕A⊕(V×V) for MSFP2's extra Cξ arcs governed by the boundary cost's directional derivative). "Primal integral" and "dual integral" polyhedral convex functions (the book's own C[Z|R→R]/C[R→R|Z] notation, used in Theorem 9.15) needed a modeling decision, since the book's own definition of these classes lies outside this chunk's page range; see Formalization scope.

Formalization scope

Ground-set vertices V and arcs A are Fintype with DecidableEq. All base M-/L-convexity vocabulary is redeclared from prior missions in this series. "Primal integral" (C[Z|R→R], M[Z|R→R]) is formalized as integer effective domain (IsDomainIntegerArc/IsDomainIntegerR); "dual integral" (C[R→R|Z], M[R→R|Z]) is formalized as the existence of an integer subgradient at every domain point (IsDualIntegralArc/IsDualIntegralR) — a standard equivalent characterization for polyhedral convex functions, and a deliberate modeling choice recorded in MODERATION_NOTES.md rather than a literal transcription of the book's own (out-of-range) definition of these two notation classes. All eight numbered results found in this chunk's page range are placed in full, with no partial-coverage scope reduction. Contributions completing any of the eight sorrys are welcome; the goal and Theorem 9.15 carry the most independent proof content.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • R. T. Rockafellar, Network Flows and Monotropic Optimization, Wiley, 1984 [178] (the classical potential/Fenchel-duality framework this mission's Theorem 9.4 adapts).
  • K. Murota, "Discrete convex analysis," Mathematical Programming, 83 (1998), pp. 313-371 [140] (the Lagrange duality and negative-cycle theory of section 9.5 this mission's Theorems 9.18 and 9.20 draw from).
88 thms2 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis XXXII: Network DualityTextbook

Motivation

Mission 32-ch09b-networkflows established the potential and negative-cycle optimality criteria for M-convex submodular flow problems. This mission finishes chapter 9 with the two topics that close it out: the constructive engine behind the negative-cycle criterion — cycle cancellation, which actually improves a nonoptimal flow rather than merely detecting suboptimality, resting on a delicate "unique-min condition" for bipartite matchings — and network duality, the chapter's capstone structural theorem showing that M-convexity and L-convexity are preserved (and their conjugacy is preserved) under transformation by an arbitrary network.

Setting

For a feasible integer flow ξ in the M-convex submodular flow problem MSFP2, a negative cycle in the auxiliary network (Gξ,ℓξ) witnesses suboptimality (mission 32's Theorem 9.20); cycle cancellation modifies ξ along a smallest such cycle to produce a strictly better flow ξ̄ (Eq. (9.75)). The unique-min condition for a pair (x,y) of integer vectors with ‖x-y‖∞=1 asks whether the bipartite graph G(x,y) — vertices the positive/negative supports of x-y, weights the M-convex exchange values Δf(x;v,u) — has a unique minimum-weight perfect matching; when it does, the M-convex exchange inequality of Proposition 6.25 becomes an equality. Separately, a network G=(V,A;S,T) with entrance set S and exit set T transforms a pair of functions f,g on Zˢ into induced functions f̃,g̃ on Zᵀ (Eqs. (9.81)-(9.82)), the minimum cost to meet a boundary specification at the exit given a production cost at the entrance and a transportation cost along arcs.

Formalization targets

Goal: Network duality for Z→Z functions (Theorem 9.26)

M-(resp. M♮^\natural♮-)convexity and integer-valuedness of f transfer to the induced f̃; L-(resp. L♮^\natural♮-)convexity and integer-valuedness of g transfer to g̃; and if f is M♮^\natural♮-convex, g is its L♮^\natural♮-conjugate, and each arc cost ga is the conjugate of fa, then g̃ is the conjugate of f̃. Chosen as goal: the book calls this "the harmonious relationship between network flow and M-/L-convexity", its own proof runs roughly six pages (the longest argument in this chunk), and it is the general fact from which Theorems 9.27-9.28 (analogues for other type combinations) and Notes 9.29-9.30 (the aggregation and infimal- convolution closure properties of M-convex functions, already placed in mission 22-ch06b-mconvexfunctions's own Theorem 6.13) all descend.

Supporting structural targets

Theorem 9.22 shows cycle cancellation strictly improves the objective; Propositions 9.23-9.25 are "the key ingredient" behind it: Proposition 9.23 shows the unique-min condition upgrades the M-convex exchange inequality to an equality, Proposition 9.24 gives a checkable characterization of when a bipartite weighted graph has a unique minimum-weight perfect matching, and Proposition 9.25 is the fact that makes the machine run — the specific pair (∂ξ,∂ξ̄) arising from cycle cancellation always satisfies the unique-min condition. Theorems 9.27 and 9.28 are the network duality theorem's own analogues for Z→R and R→R functions, the second restoring the conjugacy assertion (missing for Z→R) via the ordinary real Legendre-Fenchel transform.

Significance

Cycle cancellation is this book's constructive answer to the negative-cycle criterion: not just a certificate of suboptimality, but an actual improvement step, the combinatorial core of the cycle-canceling algorithm explained in section 10.4.3 (mission 35-ch10c-algorithms). Its correctness proof is one of the most intricate combinatorial arguments in the entire book — a proof by contradiction using a multiset-union identity (Eq. (9.80)) to derive a smaller negative cycle from an assumed non-uniqueness, itself resting on Proposition 9.24's Monge-like characterization of unique bipartite matchings. Network duality, meanwhile, is the theorem that explains why discrete convex analysis and network flow theory are so tightly intertwined throughout this book: it is the general mechanism (matroid induction, min-max relations, the M-convex aggregation and infimal-convolution closure properties) underlying nearly every construction chapter 2 introduced informally and chapter 6 proved piecemeal.

None of these results are open — they are Murota's own account of cycle cancellation (section 9.5.2) and network duality (section 9.6). What this mission contributes is a faithful, machine-checked formal statement of each, completing the platform's coverage of chapter 9 begun in missions 12-network-flows and 32-ch09b-networkflows; no comparable formalization exists on the platform (see Formalization scope).

Difficulty

Theorem 9.26's own proof needs the full weight of everything chapter 9 has built (Theorem 9.16's potential criterion for integer flows, the conjugacy theorem of chapter 8), which this mission does not re-prove (proofs are sorry throughout, per this pass's scope) but whose statement still needs the induced-function machinery built faithfully: since the general framework's optimal-value-type quantities can genuinely be -∞ (the book's own blanket hypothesis f̃ > -∞ acknowledges this), InducedFTilde/InducedGTilde are EReal-valued, following the same soundness discipline established in mission 31-ch08d-conjugacyduality for Lagrangian duality's derived quantities.

Formalization scope

Ground-set vertices V and arcs A are Fintype with DecidableEq. All base M-/L-convexity vocabulary is redeclared from prior missions. Functions "on Zˢ" for S a proper subset of the ground set are represented as ordinary V→Z functions required to vanish outside S (SupportedOn), rather than as functions on a dependent subset type — a padding-with-zero encoding consistent with this whole series' preference for a single ambient ground-set domain. C[R→R] (univariate real polyhedral convex functions, needed only for Theorem 9.28's arc costs) is formalized as ordinary midpoint-style convexity (IsConvexUnivariateR) rather than the book's own polyhedral characterization, since polyhedrality plays no role in Theorem 9.28's conclusion beyond ensuring the induced functions are well-behaved. All seven numbered results found in this chunk's page range are placed in full, with no partial-coverage scope reduction. Contributions completing any of the seven sorrys are welcome; the goal and Proposition 9.25 carry the most independent proof content.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • K. Murota, "Valuated matroid intersection," SIAM Journal on Discrete Mathematics, 9 (1996), pp. 545-561 [135] (the unique-max lemma Proposition 9.23 reformulates, and the proof technique behind Proposition 9.25).
  • K. Murota, "Discrete convex analysis," Mathematical Programming, 83 (1998), pp. 313-371 [140] (network duality and cycle cancellation for the M-convex submodular flow problem).
74 thms2 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchOptimization·Captain: Shuze Chen

Discrete Convex Analysis XII: Steepest Descent for M-Convex Function MinimizationTextbook

Motivation

Chapters 6 through 9 characterized minimality for M-convex functions structurally (the M-optimality criterion, Theorem 6.26: a point is a global minimizer iff no local swap improves it) without saying how to find one. Chapter 10 turns that structural fact into an algorithm: the local characterization is the termination test of the simplest possible minimization procedure, steepest descent by coordinate swaps. This mission formalizes that algorithm and its two complexity bounds, plus a structural min-max identity for submodular base polyhedra that the chapter's heavier submodular-minimization algorithms build on. It is the first mission in this series whose goal names a method, not just a property of a class of functions — formalizing it faithfully means giving the algorithm itself a Lean representation that the complexity theorem then quantifies over, not just describing its output.

Setting

Let f:ZV→R∪{+∞}f : \mathbb Z^V \to \mathbb R \cup \{+\infty\}f:ZV→R∪{+∞} be an M-convex function (MExchangeAxiom, chunk 06), n=∣V∣n = |V|n=∣V∣. The steepest descent algorithm repeatedly replaces the current point xxx by x−χu+χvx - \chi_u + \chi_vx−χu​+χv​ for a pair u≠vu \ne vu=v minimizing f(x−χu+χv)f(x - \chi_u + \chi_v)f(x−χu​+χv​), stopping when no such swap improves on f(x)f(x)f(x) — at which point, by the M-optimality criterion, xxx is a global minimizer. This mission represents a run of the algorithm as a sequence x:N→ZVx : \mathbb N \to \mathbb Z^Vx:N→ZV satisfying these step and termination relations directly, so that "the number of iterations" is a genuine property of any such run, not an informal gloss. Separately, for a submodular set function ρ:2V→R\rho : 2^V \to \mathbb Rρ:2V→R (chunk 04's SubmodularSetFunction), the base polyhedron B(ρ)B(\rho)B(ρ) (chunk 04's BasePolyhedron) is the polytope {x∈RV:x(X)≤ρ(X) ∀X, x(V)=ρ(V)}\{x \in \mathbb R^V : x(X) \le \rho(X)\ \forall X,\ x(V) = \rho(V)\}{x∈RV:x(X)≤ρ(X) ∀X, x(V)=ρ(V)}.

Formalization targets

Goal: Proposition 10.2 (iteration bound with tie-breaking)

For an M-convex function fff with finite ℓ1\ell^1ℓ1-diameter K1=max⁡{∥x−y∥1:x,y∈dom⁡f}K_1 = \max\{\|x-y\|_1 : x,y \in \operatorname{dom} f\}K1​=max{∥x−y∥1​:x,y∈domf} (Eq. (10.1)), the number of iterations in the steepest descent algorithm using the tie-breaking rule (10.2) — a fixed lexicographic rule for choosing among tied steepest pairs, based on an arbitrary but fixed ordering φ\varphiφ of VVV — is bounded by K1/2K_1/2K1​/2.

Milestones: Proposition 10.1, Proposition 10.8

Proposition 10.1: an unconditional warm-up — if fff has a unique minimizer x∗x^*x∗, any run of the plain (untied) algorithm from x0x^0x0 terminates within ∥x0−x∗∥1/2\|x^0 - x^*\|_1/2∥x0−x∗∥1​/2 iterations, with no tie-breaking rule needed. Proposition 10.8: a structural min-max identity, max⁡{x−(V):x∈B(ρ)}=min⁡{ρ(X):X⊆V}\max\{x^-(V) : x \in B(\rho)\} = \min\{\rho(X) : X \subseteq V\}max{x−(V):x∈B(ρ)}=min{ρ(X):X⊆V}, for a submodular set function ρ\rhoρ — chosen deliberately as a milestone that names no algorithm at all, in contrast to this mission's other two items.

Significance

The result itself. Proposition 10.2 is the complexity backbone of §10.1: it is what turns "steepest descent terminates" (an easy monotonicity observation) into a genuine polynomial bound, and it is the base case the chapter's more elaborate scaling and domain-reduction algorithms (not drafted here) improve on. Proposition 10.8, though algorithm-free, is the structural fact ("verifying membership in B(ρ)B(\rho)B(ρ) seems to need a submodular minimization procedure — but demonstrating optimality of a cut XXX only needs a base with x−(V)=ρ(X)x^-(V) = \rho(X)x−(V)=ρ(X)") that makes Schrijver's and the IFF algorithms' correctness proofs possible, and specializes chunk 04's Edmonds's intersection theorem to the two-function case ρ1=ρ,ρ2=0\rho_1 = \rho, \rho_2 = 0ρ1​=ρ,ρ2​=0.

Formalizing it. A prior-art search (q=steepest descent, q=submodular minimization, q=base polyhedron) found no existing platform items for any of this chapter's results. This mission gives the first formal statement of an M-convex minimization algorithm's complexity, and directly reuses chunk 06's MExchangeAxiom/ArgMin/CharVec/DomZ and chunk 04's SubmodularSetFunction/BasePolyhedron — genuine cross-chunk substrate reuse spanning two different chapters' worth of prior missions.

Difficulty

The chapter's own framing (quoted in BRIEF.md) is that every other chapter's theorems state a property of a class of functions or sets, while this chapter's theorems state properties of a named algorithm run on such an object. Formalizing "the number of iterations in the steepest descent algorithm is bounded by ..." faithfully means giving the algorithm's steps (S0-S3) and termination test a Lean representation that the bound then quantifies over — stating the bound about "the minimizer" alone, with the algorithm silently dropped, would misrepresent the theorem as a fact about minimizers rather than about a procedure that finds one. This mission represents a run of the algorithm as an abstract sequence satisfying the book's own step-transition relations (documented in full in MODERATION_NOTES.md), letting the complexity theorems quantify over any valid run rather than committing to one executable implementation.

A second difficulty is scope: the chapter's recommended primary goal, Proposition 10.18 (Schrijver's algorithm's complexity), needs a scaling procedure with several auxiliary data structures — a materially larger definitional undertaking than steepest descent's simple greedy-swap loop. Per BRIEF.md's own explicit fallback authorization, this mission takes Proposition 10.2 as its goal instead; see HARD.md and MODERATION_NOTES.md for the full reasoning.

Formalization scope

Runs of the algorithm are represented as x : ℕ → V → ℤ satisfying IsSteepestDescentRun (plain) or IsSteepestDescentRunTieBreak (with the tie-breaking rule) — a step relation plus a termination test, with an explicit iteration count N the theorems bound. The tie-breaking key Φ(u,v)\Phi(u,v)Φ(u,v) (Eq. (10.2)) and its lexicographic order are formalized directly (a manual three-way comparison, not Mathlib's default componentwise Prod order). Both iteration bounds are stated as 2 * N ≤ k rather than N ≤ k / 2, avoiding natural-number division. Not drafted: the derived "hence ... in O(F⋅n2K1)O(F \cdot n^2 K_1)O(F⋅n2K1​) time" corollary of Proposition 10.2 (needs a cost-model primitive for "time" and "FFF" this mission does not otherwise use — the iteration-count bound itself, this proposition's genuine combinatorial content, is drafted in full); Schrijver's algorithm and everything in §10.2.2 onward; the steepest descent scaling algorithm, the domain reduction algorithm and its scaling variant (§10.1.2-10.1.3, structurally different algorithms); the quasi-M-convex extension mentioned immediately after Proposition 10.2. A trivializing formalization would state the iteration bound as an unconditional fact about "a" minimizer-finding procedure, or would drop the tie-breaking rule from Proposition 10.2 and thereby understate what the bound actually requires; neither is done.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
16 thms2 active usersReviewed
🏆Completed
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis XXXIII: Near-Optimality for Submodular MinimizationTextbook

Motivation

Chapter 10 turns from structure theory to algorithms: efficient methods for minimizing M-convex functions (via domain reduction) and submodular set functions (via Schrijver's and the Iwata-Fleischer-Fujishige scaling algorithms). Most of chapter 10's numbered results are asymptotic running-time bounds for specific procedural algorithms — a genuinely different kind of claim from the rest of this book (see Formalization scope). This mission places the results of this block that ARE ordinary mathematical propositions: correctness certificates, min-max theorems, and structural facts the algorithms rely on and produce.

Setting

For an M-convex set B ⊆ Z^V, the central part B° (the vectors of B lying away from its boundary, defined via per-coordinate bounds ℓ°_B, u°_B) is what the domain reduction algorithm searches from. For a submodular set function ρ : 2^V → R, the base polyhedron B(ρ) and its extreme bases (one per linear ordering of V, via Eq. (10.12)) let any base be written as a convex combination of finitely many extreme bases (Eq. (10.13)); a candidate minimizer W is certified via the linear orderings representing an optimal base. The Iwata-Fleischer-Fujishige (IFF) scaling algorithm relaxes this problem with a flow-augmentation parameter δ, maintaining a δ-feasible flow φ and vector z = x + ∂φ; near the end of a scaling phase, no augmenting path and no "active triple" together certify near-optimality.

Formalization targets

Goal: Near-optimality from the absence of augmenting paths (Proposition 10.20)

If S ⊆ W ⊆ V∖T, no arc of the auxiliary network leaves W, and no active triple exists, then z⁻(V) ≥ ρ(W)-nδ and x⁻(V) ≥ ρ(W)-n²δ; moreover W exactly minimizes ρ once δ is small enough relative to the smallest positive gap between two values of ρ. Chosen as goal: the book calls this "a key property of the scaling algorithm" and "a relaxation version of the min-max relation in Proposition 10.8", its own proof is the most substantial argument among this chunk's placed results, and Proposition 10.23 is a direct corollary of it.

Supporting structural targets

Proposition 10.8 is the min-max relation underlying the whole of section 10.2 (an Edmonds- intersection-theorem consequence, found by direct reading — the extractor's table missed it). Proposition 10.9 gives the three-part sufficient condition for optimality, in terms of the linear orderings representing an optimal base, that both Schrijver's algorithm and the IFF algorithm use as their termination criterion (also found by direct reading). Propositions 10.5-10.6 establish that the domain reduction algorithm's central part B° is always nonempty, via an explicit vector-extension step (Proposition 10.5 likewise missing from the extractor's table). Proposition 10.23 fixes individual coordinates once a scaling phase ends, and Proposition 10.24 gives the termination certificate for the IFF fixing algorithm's own separate graph-contraction procedure.

Significance

Chapter 10 is where this book cashes out its structure theory as algorithms with provable running times, and this mission places every result of that chapter's first two sections that is a mathematical proposition rather than a runtime bound: two min-max/optimality-certificate theorems (10.8-10.9) that are the combinatorial core making the following two strongly polynomial algorithms (Schrijver's, and Iwata-Fleischer-Fujishige's) correct, one central-part nonemptiness fact (10.5-10.6) underlying the domain reduction algorithm, and the two fixing/termination certificates (10.23-10.24) that let the scaling algorithms actually output a minimizer with a proof of optimality attached, not just a numerical answer.

None of these results are open — they are Murota's own account of submodular-function- minimization algorithms (sections 10.1-10.2). What this mission contributes is a faithful, machine-checked formal statement of each, including two results (Propositions 10.5 and 10.8) the platform's own automated extractor missed; no comparable formalization exists on the platform (see Formalization scope).

Difficulty

Six numbered results in this block (Propositions 10.4, 10.7, 10.17, 10.18, 10.21, 10.22) are excluded as hard: each states a Big-O asymptotic bound on the running time, function-evaluation count, or internal-procedure-call count of a specific iterative algorithm (the domain reduction algorithm, its scaling variant, Schrijver's algorithm, the IFF scaling algorithm). Faithfully stating "this algorithm runs in O(g(n)) time" requires a cost-tracked operational semantics for that specific algorithm — a well-founded recursive procedure with an oracle for evaluating the input function, threading a step/evaluation counter, instantiated over an unbounded family of ground-set sizes n and numeric parameters (K∞, M) — which is a fundamentally different kind of formalization task (computational complexity theory) from every one of the roughly 280 other numbered results in this book, none of which require modeling the cost of computing them. See HARD.md.

Formalization scope

Ground-set elements are a Fintype V with DecidableEq. Base M-convex-set vocabulary is redeclared from prior missions. Linear orderings of V are represented as bijections V ≃ Fin (Fintype.card V) rather than as lists, matching this series' established preference for order-indexed families over sequential data structures. Proposition 10.5's witness vector is stated as an existence claim (the mathematical content of the proposition), rather than by reconstructing the specific recursive modification procedure the book uses to produce it — a choice consistent with how this series has always formalized "the algorithm produces X" claims where X is a mathematical property, by asserting X's existence rather than executing the algorithm (see, e.g., mission 33-ch09c-networkflows's cycle-cancellation theorem). Six numbered results (Propositions 10.4, 10.7, 10.17, 10.18, 10.21, 10.22) are hard; see HARD.md. Contributions completing any of the seven sorrys are welcome; the goal and Proposition 10.9 carry the most independent proof content.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • S. Iwata, L. Fleischer, and S. Fujishige, "A combinatorial strongly polynomial algorithm for minimizing submodular functions," Journal of the ACM, 48 (2001), pp. 761-777 [102] (the IFF scaling algorithm this mission's Proposition 10.20 certifies).
  • A. Schrijver, "A combinatorial algorithm minimizing submodular functions in strongly polynomial time," Journal of Combinatorial Theory, Series B, 80 (2000), pp. 346-355 [182] (Schrijver's algorithm, whose termination criterion is Proposition 10.9).
41 thms2 active usersReviewed
Next

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