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.

35 missions

Missions

1–20 of 35
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
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis XVI: Substitutes and Complements in Network FlowsTextbook

Motivation

In economics, a pair of goods are substitutes if raising the price of one increases demand for the other, and complements if it decreases it; formally, a utility or value function is submodular in the substitutes case and supermodular in the complements case. A natural question is which of these two regimes a given optimization problem's value function falls into, and whether the answer depends on the underlying combinatorial structure of the problem rather than being a coincidence of the particular numbers involved. Murota's Discrete Convex Analysis (SIAM, 2003) answers this question for the maximum-weight circulation problem in a directed network: the value function is submodular in some coordinates and supermodular in others, purely as a consequence of a graph-theoretic distinction — whether the arcs involved are parallel or series — and this chapter shows the distinction is explained precisely by the dual pair of discrete convexity notions (L-natural-convexity and M-natural-convexity) developed elsewhere in the book. This mission also completes the quadratic-forms thread the previous mission in this series (Discrete Convex Analysis XV) began, by formalizing its natural generalization to functions that may take the value +∞+\infty+∞.

Setting

Let G=(V,A)G=(V,A)G=(V,A) be a directed graph with vertex set VVV and arc set AAA; write ∂+a\partial^+a∂+a, ∂−a\partial^-a∂−a for the initial and terminal vertex of arc aaa. For a flow ξ:A→R\xi:A\to\mathbb Rξ:A→R, its boundary is ∂ξ(v)=∑a:∂+a=vξ(a)−∑a:∂−a=vξ(a)\partial\xi(v)=\sum_{a:\partial^+a=v}\xi(a)-\sum_{a:\partial^-a=v}\xi(a)∂ξ(v)=∑a:∂+a=v​ξ(a)−∑a:∂−a=v​ξ(a), the net flow leaving vvv. Given a capacity c:A→R≥0c:A\to\mathbb R_{\ge0}c:A→R≥0​, ξ\xiξ is a feasible circulation for ccc if 0≤ξ(a)≤c(a)0\le\xi(a)\le c(a)0≤ξ(a)≤c(a) for every arc and ∂ξ(v)=0\partial\xi(v)=0∂ξ(v)=0 for every vertex. For a weight w:A→Rw:A\to\mathbb Rw:A→R, F(w,c)=max⁡{⟨w,ξ⟩:ξ feasible for c}F(w,c)=\max\{\langle w,\xi\rangle : \xi\text{ feasible for }c\}F(w,c)=max{⟨w,ξ⟩:ξ feasible for c} is the maximum-weight circulation value, and ξ\xiξ is optimal for www (with capacity ccc) if it is feasible and attains this maximum. A simple cycle is an alternating sequence of pairwise distinct vertices v0,…,vk−1v_0,\dots,v_{k-1}v0​,…,vk−1​ and arcs a1,…,aka_1,\dots,a_ka1​,…,ak​ with {∂+ai,∂−ai}={vi−1,vi}\{\partial^+a_i,\partial^-a_i\}=\{v_{i-1},v_i\}{∂+ai​,∂−ai​}={vi−1​,vi​} (indices mod kkk) and v0=vkv_0=v_kv0​=vk​. Two arcs are parallel if every simple cycle containing both of them orients them oppositely, and series if every such cycle orients them the same way; a set of arcs is parallel (series) if its arcs are pairwise parallel (series). A circuit is a {0,±1}\{0,\pm1\}{0,±1}-valued π:A→R\pi:A\to\mathbb Rπ:A→R with ∂π=0\partial\pi=0∂π=0 whose support forms a simple cycle. For x∈Rnx\in\mathbb R^nx∈Rn, 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}. A function g:Rn→Rg:\mathbb R^n\to\mathbb Rg:Rn→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), supermodular with the reverse inequality, and has translation submodularity (is L-natural-convex) if the stronger inequality g(p)+g(q)≥g((p−α1)∨q)+g(p∧(q+α1))g(p)+g(q)\ge g((p-\alpha\mathbf1)\vee q)+g(p\wedge(q+\alpha\mathbf1))g(p)+g(q)≥g((p−α1)∨q)+g(p∧(q+α1)) holds for every α≥0\alpha\ge0α≥0. A function fff has the M-natural exchange property (is M-natural-convex) if for 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 with 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 α∈[0,α0]\alpha\in[0,\alpha_0]α∈[0,α0​]; a function is M-natural-concave or L-natural-concave if its negation is M-natural- or L-natural-convex.

Formalization targets

Goal (Theorem 2.23). For PPP a parallel arc set and SSS a series arc set,

F is L-natural-convex in wP and M-natural-concave in cP,F\text{ is L-natural-convex in }w_P\text{ and M-natural-concave in }c_P,F is L-natural-convex in wP​ and M-natural-concave in cP​, F is M-natural-convex in wS and L-natural-concave in cS,F\text{ is M-natural-convex in }w_S\text{ and L-natural-concave in }c_S,F is M-natural-convex in wS​ and L-natural-concave in cS​,

where wPw_PwP​, cPc_PcP​ denote FFF's dependence on the coordinates of www, ccc indexed by PPP (resp. SSS) with the remaining coordinates held fixed. This is the mission's capstone: it upgrades the plain submodularity/supermodularity split of Theorem 2.22 to the sharper pair of combinatorial convexity classes that explains it.

Supporting milestones. Proposition 2.21 (the classical fact that FFF is convex in www and concave in ccc, with no combinatorial content — the baseline against which Theorem 2.23's sharper claim is measured); Theorem 2.16 (the general, possibly-+∞+\infty+∞-valued extension of the quadratic-form conjugacy from Discrete Convex Analysis XV's Theorem 2.11, to functions restricted to a linear subspace); Theorem 2.22 (plain submodularity/supermodularity of FFF in wP,cPw_P,c_PwP​,cP​ and wS,cSw_S,c_SwS​,cS​, the result Theorem 2.23 strengthens); and Propositions 2.24–2.28 (the graph-theoretic lemmas — sparse intersection of a circuit's support with a parallel or series arc set, merging two circuits along a series set, and three existence statements for optimality-preserving perturbations — that the book's own proof of Theorem 2.23 is built from).

Significance

Theorem 2.23 gives a structural explanation, rather than a case-by-case verification, for a phenomenon well known in network flow theory: that convexity/concavity and submodularity/supermodularity are independent properties, appearing in all four combinations depending on which side of the problem (weights or capacities) and which graph-theoretic role (parallel or series) is varied. Without it, (2.55)'s four combinations would be four separate facts with no common cause; with it, they are corollaries of two applications of a single pair of dual discrete-convexity notions, the same notions the book uses throughout to unify matroid theory, submodular optimization, and convex analysis. Formalizing this mission produces, so far as a platform search shows, the first Lean statement of a combinatorial-convexity classification result for a network optimization value function, together with the graph-theoretic vocabulary (simple cycles, parallel/series arcs, circuits) needed to state it — infrastructure with no prior formalized counterpart on the platform that a later mission on network flows or matroid union could reuse.

Difficulty

The naive approach to Theorem 2.23 tries to verify translation submodularity or the exchange property directly from the linear-programming definition of FFF as a maximum over a polytope, treating wP↦F(w,c)w_P\mapsto F(w,c)wP​↦F(w,c) as an abstract convex-piecewise-linear function; this loses the graph structure entirely and gives at best the plain submodularity of Theorem 2.22, not the sharper L-natural/M-natural classification, because submodularity alone does not distinguish a combinatorially meaningful discrete convexity from an arbitrary submodular function. The book's actual route instead works with explicit optimal circulations for the two perturbed weight vectors and reconstructs a feasible pair achieving the target inequality by rerouting flow along a circuit — and the existence of a usable circuit (one that touches the perturbed arcs in a way compatible with the parallel or series structure) is exactly what Propositions 2.24–2.28 supply via the conformal decomposition of a difference of two circulations into elementary circuits. This is why those five propositions, although individually narrow existence lemmas, are included as milestones: they are the load-bearing combinatorial content the naive convex-analytic argument cannot reach.

Formalization scope

The graph is {V A : Type*} with src dst : A → V rather than a bundled structure, matching the book's own ∂+,∂−\partial^+,\partial^-∂+,∂− notation directly. F(w,c)F(w,c)F(w,c) is a real sSup over feasible circulations' weights (existence of a maximizer is not asserted, since no proof is attempted this pass); IsOptimalCirc is a separate, directly-stated primitive for "ξ\xiξ is optimal for www", matching the book's own working vocabulary in the propositions that need it. A simple cycle is formalized as an injective cyclically-indexed vertex sequence together with a matching arc sequence, exactly as the book's own footnote defines it; parallel and series arcs are defined by quantifying over every such representation of every simple cycle containing the two arcs, which is checked to be independent of which of a cycle's two traversal directions or starting vertex is chosen. Viewing FFF as a function of wPw_PwP​ alone extends a partial vector by a fixed background vector on the complement of PPP, the same partial-application device the book uses informally. M-natural- and L-natural-concavity are recorded as the corresponding convexity property of the negated function, the standard convention. The formalization does not trivialize: parallel and series arc sets are genuine graph-theoretic hypotheses (not, e.g., specialized to ∣P∣=1|P|=1∣P∣=1 or a graph with no simple cycles, which would make the parallel/series distinction vacuous), and Theorem 2.23's four conclusions are stated with the same combinatorial-convexity predicates (TranslationSubmodular, MNatExchangeR) used for the book's sharpest discrete convexity classes, not weakened to plain submodularity/supermodularity. Theorem 2.16 additionally needs Set (V → ℝ)-valued subspaces K, H (following the book's own set-builder notation for ker M and X⊥ rather than bundling them as Mathlib Submodules) and a WithTop ℝ-valued Legendre- Fenchel conjugate. Infrastructure needed beyond Mathlib: all graph, circulation, and combinatorial-convexity vocabulary is defined fresh in DiscreteConvex.CombinatorialC; a contribution proving any of the five graph-theoretic lemmas (Propositions 2.24–2.28) or the convex/concave halves of Proposition 2.21 independently would be a natural entry point.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003, DOI 10.1137/1.9780898718508, Chapter 2.
  • K. Murota, A. Shioura, "Conjugacy relationship between M-convex and L-convex functions in continuous variables," Mathematical Programming 101 (2004), 415–433.
  • R. T. Rockafellar, Network Flows and Monotropic Optimization, Wiley, 1984.
41 thms3 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
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis XVII: Fenchel Duality and Linear-Programming IntegralityTextbook

Motivation

Duality is the organizing principle of convex optimization: a minimization problem's optimal value equals a maximization problem's optimal value, and this coincidence, rather than being a lucky accident, follows from a separating-hyperplane argument that applies whenever the two problems' feasible regions are shaped compatibly enough. Werner Fenchel formalized this in the 1950s for pairs of convex and concave functions related by the Legendre-Fenchel transform, and the resulting Fenchel duality theorem specializes, for linear objectives over polyhedral feasible regions, to linear programming duality — the fact, central to the entire theory of combinatorial optimization, that a linear program's optimal value can always be certified from above and below by a pair of primal and dual feasible solutions. Murota's Discrete Convex Analysis (SIAM, 2003) collects this classical machinery, together with the integrality theory that lets it produce combinatorial (integer-valued) certificates rather than merely real ones, as the technical foundation the rest of the book builds its discrete theory on top of.

Setting

For f:Rn→R∪{+∞}f : \mathbb R^n \to \mathbb R \cup \{+\infty\}f:Rn→R∪{+∞}, the epigraph is epi⁡f={(x,Y):Y≥f(x)}\operatorname{epi} f = \{(x,Y) : Y \ge f(x)\}epif={(x,Y):Y≥f(x)}, and fff is convex iff epi⁡f\operatorname{epi} fepif is a convex set; fff is proper if additionally its effective domain dom⁡f={x:f(x)<+∞}\operatorname{dom} f = \{x : f(x) < +\infty\}domf={x:f(x)<+∞} is nonempty, and closed if epi⁡f\operatorname{epi} fepif is topologically closed. A function h:Rn→R∪{−∞}h : \mathbb R^n \to \mathbb R \cup \{-\infty\}h:Rn→R∪{−∞} is concave, proper, closed analogously via its hypograph. The convex conjugate 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)}, and the concave conjugate h∘(p)=inf⁡x{⟨p,x⟩−h(x)}h^\circ(p) = \inf_x\{\langle p,x\rangle - h(x)\}h∘(p)=infx​{⟨p,x⟩−h(x)}. The relative interior ri⁡S\operatorname{ri} SriS of a set SSS is the interior of SSS relative to its affine hull. A function is polyhedral if its epigraph (or hypograph) is a finite intersection of half-spaces. Given an m×nm \times nm×n matrix AAA, b∈Rmb \in \mathbb R^mb∈Rm, c∈Rnc \in \mathbb R^nc∈Rn, the primal and dual linear programs are min⁡{c⊤x:Ax=b, x≥0}\min\{c^\top x : Ax=b,\ x\ge0\}min{c⊤x:Ax=b, x≥0} and max⁡{b⊤y:A⊤y≤c}\max\{b^\top y : A^\top y \le c\}max{b⊤y:A⊤y≤c}, with feasible regions PPP, DDD. A matrix is totally unimodular if every square submatrix has determinant 000, 111, or −1-1−1. A discrete set S⊆ZnS \subseteq \mathbb Z^nS⊆Zn is hole free if S=Sˉ∩ZnS = \bar S \cap \mathbb Z^nS=Sˉ∩Zn, where Sˉ\bar SSˉ is the convex hull of SSS's real embedding; 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​}.

Formalization targets

Goal (Theorem 3.6, Fenchel duality). For proper convex fff and proper concave hhh satisfying at least one of four alternative conditions — a relative-interior condition on dom⁡f∩dom⁡h\operatorname{dom} f \cap \operatorname{dom} hdomf∩domh, a polyhedrality condition on the same, or the analogous pair of conditions on dom⁡f∙∩dom⁡h∘\operatorname{dom} f^\bullet \cap \operatorname{dom} h^\circdomf∙∩domh∘ together with closedness of fff, hhh —

inf⁡x{f(x)−h(x)}=sup⁡p{h∘(p)−f∙(p)},\inf_x\{f(x)-h(x)\} = \sup_p\{h^\circ(p)-f^\bullet(p)\},xinf​{f(x)−h(x)}=psup​{h∘(p)−f∙(p)},

with the extremum on the appropriate side attained whenever the common value is finite. This is the mission's capstone: the four alternative hypotheses make it the most broadly applicable statement of the four convex-duality results in this mission, each of the other three being either a special case in substance (Theorem 3.5, separation, which 3.6 is proved from) or a literal specialization to linear data (Theorem 3.10, LP duality).

Supporting milestones. Theorem 3.2 (biconjugation: f∙f^\bulletf∙ is always closed proper convex, and g∙∙=gg^{\bullet\bullet}=gg∙∙=g for closed proper convex ggg); Theorem 3.5 (the separation theorem for convex/concave functions, under two of Theorem 3.6's four hypotheses); Theorem 3.9 (the Farkas lemma, equality form); Theorem 3.10 (LP duality: weak duality, strong duality with attainment, and complementary slackness); Theorem 3.13 (total unimodularity of the constraint matrix guarantees an integral optimal solution whenever an optimal solution exists); Proposition 3.14 (an explicit potential function certifying a minimum-weight bipartite perfect matching, via the totally unimodular incidence-matrix LP); and Proposition 3.16 (for a translation-invariant family of hole-free discrete sets, the property that discrete disjointness implies closure disjointness is equivalent to the discrete Minkowski sum matching the integer points of the closures' Minkowski sum).

Significance

Fenchel duality is the single result from which the separation theorem, LP duality, and (via the totally-unimodular incidence matrix of a bipartite graph) the combinatorial duality underlying weighted bipartite matching all descend, in one unbroken chain of specialization; formalizing this chain in one mission exhibits that structure directly, rather than treating each result as an independent fact. Proposition 3.16 plays a different role: it is the chapter's warning that naive discrete analogues of convexity (hole-freeness) do not automatically inherit convexity's good closure properties under Minkowski sums, which is exactly the gap the book's later M-convexity and L-convexity machinery is built to close — this mission's Proposition 3.16 is therefore the motivating negative result for the rest of the book's positive theory, not a loose end. So far as a platform search shows, no existing formalization matches this chunk's specific combination of extended-valued (possibly ±∞\pm\infty±∞) functions, the four-alternative Fenchel duality hypothesis, or the bipartite-matching-via-total-unimodularity argument; the one related platform result (VectorSpaceOpt.fenchel_duality, from Luenberger) is for real-valued functions on general normed spaces under a single relative-interior-and-solidness hypothesis, a different generality from the extended-valued, four-hypothesis statement here.

Difficulty

The naive approach to Theorem 3.6 tries to prove the duality gap is zero directly from the definitions of the two conjugates, which only gives the easy inequality inf⁡≥sup⁡\inf \ge \supinf≥sup (a one-line computation, shown in the book's own proof in three lines); the substantive content is the reverse inequality, and it genuinely fails without a constraint-qualification hypothesis like (a1)-(b2) — Example 3.8 in the book exhibits a convex/concave pair with inf⁡=0≠−1=sup⁡\inf = 0 \ne -1 = \supinf=0=−1=sup when none of the four conditions hold. The book's actual route reduces Theorem 3.6 to the separation theorem (Theorem 3.5) applied to fff shifted down by the (assumed finite) infimum, which produces the separating affine function directly; this is why Theorem 3.5, although logically a special case in spirit, earns its own milestone rather than being subsumed silently.

Formalization scope

All convex and concave functions are represented uniformly as (V → ℝ) → EReal-valued (Fintype V), rather than mixing WithTop ℝ for convex and WithBot ℝ for concave functions, so that Theorem 3.2's biconjugate — whose properness is a conclusion, not an assumption — has a well-defined codomain without extra casts. Convexity is defined via the epigraph being a convex subset of the ordinary real vector space (V→R)×R(V\to\mathbb R)\times\mathbb R(V→R)×R (Mathlib's Convex ℝ), following the book's own equivalent characterization, rather than unfolding the direct inequality definition, which would require a extended-arithmetic scalar-multiplication convention (0\cdot(+\infty)=0) that Mathlib does not provide for EReal. The relative interior is defined directly from the book's own metric-ball-intersected-with-affine-hull description, since Mathlib has no relative-interior primitive at the pinned revision. Polyhedra are finite intersections of explicit half-spaces. A bipartite perfect matching is represented as a bijection between the two vertex sides restricted to the edge set — a faithful, not narrower, representation since every perfect matching between equal-size parts arises this way. The formalization does not trivialize: Theorem 3.6's four hypotheses are carried in full (not reduced to the easiest single case), and no result is stated only for finite-valued (never ±∞\pm\infty±∞) functions, which would discard the entire point of the extended-value convex-analysis framework this chapter sets up for the rest of the book. Infrastructure needed beyond Mathlib's Convex, Matrix, and EReal API: all epigraph/hypograph, conjugate, relative-interior, and polyhedral apparatus is defined fresh in DiscreteConvex.IntegralConvexityB; a contribution proving any of the seven milestones independently, or supplying Mathlib-quality relative-interior lemmas, would be a natural entry point.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003, DOI 10.1137/1.9780898718508, Chapter 3.
  • R. T. Rockafellar, Convex Analysis, Princeton University Press, 1970.
  • A. Schrijver, Theory of Linear and Integer Programming, Wiley, 1986.
37 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
CombinatoricsConvex OptimizationDiscrete Geometry+2·Captain: Shuze Chen

Discrete Convex Analysis XIX: Discrete Separation for M-Convex SetsTextbook

Motivation

Submodular set functions are the combinatorial stand-in for convexity: a function ρ:2V→R\rho : 2^V \to \mathbb Rρ:2V→R on the subsets of a finite ground set VVV 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), and this single diminishing-returns inequality drives an enormous range of combinatorial optimization — matroid rank functions, graph cut capacities, entropy, coverage functions, and the max-flow min-cut theorem all arise as special or dual cases (Edmonds 1970; Lovász 1983; Fujishige 2005). M-convex sets are the "vector" incarnation of the same idea: subsets BBB of ZV\mathbb Z^VZV satisfying an exchange axiom that generalizes the basis-exchange property of matroids to sets of integer points lying on a common hyperplane. Murota's Discrete Convex Analysis (SIAM, 2003) develops both sides of this correspondence and proves they coincide exactly: M-convex sets are precisely the integer points of the base polyhedra of integer-valued submodular functions. This mission covers the second half of that development — the structural theory (integrality, holes, Minkowski sums) that turns the correspondence into a working calculus, and its capstone, a discrete separation theorem for two disjoint M-convex sets whose separating hyperplane is forced to have {0,1}\{0,1\}{0,1}- or {0,−1}\{0,-1\}{0,−1}-valued coefficients.

Companion mission 04-mconvex-sets (Discrete Convex Analysis III) covers the same chapter's foundational results: the equivalence of the exchange-axiom variants, the one-to-one correspondence between M-convex sets and integer submodular functions (Theorem 4.15), Edmonds's intersection theorem (Theorem 4.18), and Frank's discrete separation theorem for submodular/ supermodular pairs (Theorem 4.17). This mission builds on that vocabulary (redeclared here, since draft missions in the same series cannot yet import one another) and proves the results the chapter leaves for its second half.

Setting

Fix a finite ground set VVV. A vector x∈ZVx \in \mathbb Z^Vx∈ZV assigns an integer x(v)x(v)x(v) to each v∈Vv \in Vv∈V; write x(X)=∑v∈Xx(v)x(X) = \sum_{v \in X} x(v)x(X)=∑v∈X​x(v) for X⊆VX \subseteq VX⊆V. For x,y∈ZVx, y \in \mathbb Z^Vx,y∈ZV, the positive support supp⁡+(x−y)={v:x(v)>y(v)}\operatorname{supp}^+(x-y) = \{v : x(v) > y(v)\}supp+(x−y)={v:x(v)>y(v)} and negative support supp⁡−(x−y)={v:x(v)<y(v)}\operatorname{supp}^-(x-y) = \{v : x(v) < y(v)\}supp−(x−y)={v:x(v)<y(v)} record where xxx exceeds, and falls short of, yyy. A nonempty set B⊆ZVB \subseteq \mathbb Z^VB⊆ZV is M-convex if it satisfies the exchange axiom (B-EXC[Z]): for all x,y∈Bx, y \in Bx,y∈B 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) has 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 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], or S[Z]S[\mathbb Z]S[Z] when integer-valued) 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) for all X,YX, YX,Y. Its base polyhedron is 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)}. The Lovász extension ρ^:RV→R∪{±∞}\hat\rho : \mathbb R^V \to \mathbb R \cup \{\pm\infty\}ρ^​:RV→R∪{±∞} linearly interpolates ρ\rhoρ off {0,1}V\{0,1\}^V{0,1}V: sorting the distinct values of p∈RVp \in \mathbb R^Vp∈RV as p^1>⋯>p^m\hat p_1 > \cdots > \hat p_mp^​1​>⋯>p^​m​ and setting Ui={v:p(v)≥p^i}U_i = \{v : p(v) \ge \hat p_i\}Ui​={v:p(v)≥p^​i​}, it is ρ^(p)=∑i=1m−1(p^i−p^i+1)ρ(Ui)+p^mρ(Um)\hat\rho(p) = \sum_{i=1}^{m-1}(\hat p_i - \hat p_{i+1})\rho(U_i) + \hat p_m \rho(U_m)ρ^​(p)=∑i=1m−1​(p^​i​−p^​i+1​)ρ(Ui​)+p^​m​ρ(Um​).

Formalization targets

Goal: discrete separation for M-convex sets

B1∩B2=∅  ⟹  ∃ p∗∈{0,1}V∪{0,−1}V,inf⁡x∈B1⟨p∗,x⟩−sup⁡x∈B2⟨p∗,x⟩≥1,B_1 \cap B_2 = \emptyset \implies \exists\, p^* \in \{0,1\}^V \cup \{0,-1\}^V,\quad \inf_{x \in B_1}\langle p^*, x\rangle - \sup_{x \in B_2}\langle p^*, x\rangle \ge 1,B1​∩B2​=∅⟹∃p∗∈{0,1}V∪{0,−1}V,x∈B1​inf​⟨p∗,x⟩−x∈B2​sup​⟨p∗,x⟩≥1,

for M-convex sets B1,B2⊆ZVB_1, B_2 \subseteq \mathbb Z^VB1​,B2​⊆ZV (Theorem 4.21). This is the weakest stable form of the result — it asserts only the existence of a combinatorially special separator, not any bound tied to ∣V∣|V|∣V∣ or a particular construction, so it is not invalidated by a sharper algorithm for finding p∗p^*p∗.

Supporting structural targets

Eleven further results build the calculus this goal rests on: the hyperplane property of M-convex sets (Prop. 4.1), an equivalent one-sided exchange axiom (Prop. 4.2), nonemptiness and the support-function identity for B(ρ)B(\rho)B(ρ) (Props. 4.4-4.5), integrality of B(ρ)B(\rho)B(ρ) for integer-valued ρ\rhoρ (Prop. 4.6), the hole-free property identifying an M-convex set with the integer points of its own convex hull (Thm. 4.12), the two-way polyhedral description of M-convex sets via induced submodular functions (Props. 4.13-4.14), the equivalence of submodularity with convexity of the Lovász extension (Thm. 4.16, due to Lovász), integrality of the intersection of M-convex sets (Thm. 4.22), and Minkowski-sum identities for base polyhedra and M-convex sets (Thm. 4.23).

Significance

The discrete separation theorem is what makes M-convexity discrete rather than merely a polyhedral fact: ordinary separation of two disjoint convex sets by a hyperplane is classical, but here the separator is forced into {0,1}V∪{0,−1}V\{0,1\}^V \cup \{0,-1\}^V{0,1}V∪{0,−1}V — a purely combinatorial object — with no loss of strength. This is the mechanism behind integrality results across combinatorial optimization (e.g., that the intersection of two integral base polyhedra is integral, Theorem 4.22, used pervasively in matroid intersection and submodular flow algorithms). The structural results (holes, Minkowski sums, the Lovász-extension convexity equivalence) are the working toolkit every later use of M-convexity in the book — proximity theorems for M-convex functions (chunks 06+), the discrete conjugacy theorem, submodular flows — draws on without restating.

None of these results are open: Murota attributes the exchange-axiom theory to the matroid and submodular-function literature it systematizes, citing Edmonds, Frank, and Lovász by name for the specific theorems. What this mission produces is a machine-checked formal statement of each result exactly as the book states it, in a shared Lean vocabulary (ExchangeAxiomB, BasePolyhedron, LovaszExtension) that the rest of the Discrete Convex Analysis series builds on; no result here has a prior formalization on the platform (see Formalization scope).

Difficulty

The separation theorem is not proved by convex separation directly — the whole point is that the naive proof (apply the ordinary hyperplane separation theorem to the convex hulls of B1,B2B_1, B_2B1​,B2​, then argue the separator can be taken {0,1}\{0,1\}{0,1}-valued) does not go through, because convex separation alone gives no control over the separator's coefficients. The book instead derives it from Edmonds's intersection theorem (Theorem 4.18, chunk 04-mconvex-sets) applied to a submodular/supermodular pair built from B1,B2B_1, B_2B1​,B2​'s associated set functions (Theorem 4.15), routed through Frank's discrete separation theorem (Theorem 4.17) — a genuine two-step reduction, not a direct argument. A second, independent difficulty sits in the supporting results: the hole-free property (Theorem 4.12) requires an explicit induction reducing an arbitrary convex combination representing an integer point to a single element of BBB, a combinatorial exchange argument with no shortcut through general polyhedral theory.

Formalization scope

Ground-set elements are a Fintype V with DecidableEq; M-convex sets are Set (V → ℤ); submodular/supermodular functions are Finset V → WithTop ℝ / WithBot ℝ; base polyhedra are Set (V → ℝ). The Lovász extension is formalized directly from the book's own sorted-values construction (SortedValues, LevelSet, Eq. (4.4)-(4.6)), not via an equivalent closed form. Since WithTop ℝ carries no Module ℝ structure, convexity for Theorem 4.16 is stated via a bespoke nonnegative-scalar action (ScalarWithTop) rather than Mathlib's ConvexOn — this changes no mathematical content, only its packaging (see MODERATION_NOTES.md). No numeric constants are hard-coded anywhere in this mission (rule 7 is vacuous). The goal's hypothesis (ExchangeAxiomB plus Nonempty on each BiB_iBi​) is exactly the book's own definition of M-convexity — no weaker substitute (e.g. requiring a specific ρ\rhoρ witness in the hypothesis rather than deriving one, or dropping the {0,1}/{0,−1}\{0,1\}/\{0,-1\}{0,1}/{0,−1} constraint on p∗p^*p∗ in favor of a generic separator) would be faithful, and both trivializations are ruled out by construction. This mission's definitions (ExchangeAxiomB, BasePolyhedron, SubmodularSetFunction, LovaszExtension) are redeclared from chunk 04-mconvex-sets 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 twelve sorrys are welcome; the hole-free property (Theorem 4.12) and the goal are the two with the most independent proof content.

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, 1970, pp. 69-87.
  • A. Frank, "An algorithm for submodular functions on graphs," Annals of Discrete Mathematics, 16 (1982), pp. 97-120.
  • L. Lovász, "Submodular functions and convexity," in Mathematical Programming: The State of the Art, Springer, 1983, pp. 235-257.
29 thms3 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
Convex OptimizationDiscrete GeometryOperations Research+1·Captain: Shuze Chen

Discrete Convex Analysis VI: Quasi M-Convex Functions and the Quasi-Proximity TheoremTextbook

Motivation

Convexity is normally defined additively — a function's value at a mixture is bounded by the mixture of its values — but many of the properties that make convexity useful in optimization (a local minimum is global, level sets are well-behaved) survive under a much weaker, purely ordinal notion: quasi-convexity, which compares function values rather than adding them. A nondecreasing rescaling of a convex function is generally not convex, but it is always quasi-convex — so a theory built only on ordinal comparisons automatically covers every such rescaling for free, at the cost of a more delicate proof architecture (since the algebraic cancellations available to additive convexity are no longer available).

Chapter 6's second half asks exactly how far this idea extends in the discrete setting: does the M-convexity exchange axiom have an ordinal, quasi-convex relaxation that still supports the same strong minimization theory — an optimality criterion, a minimizer-cut lemma, and, most significantly, a proximity theorem with the same explicit distance bound? This mission formalizes the chapter's answer: yes, and the relevant relaxed class, functions satisfying condition (SSQM≠_{\ne}=​), is large enough to include every strictly increasing rescaling of an M-convex function, a class the M-convex theory of chunk 06 alone says nothing about.

Setting

Let VVV be a finite ground set and f:ZV→R∪{+∞}f : \mathbb Z^V \to \mathbb R \cup \{+\infty\}f:ZV→R∪{+∞} with nonempty effective domain. Building on chunk 06's M-convex exchange axiom (M-EXC[Z]), this chapter introduces several ordinal relaxations. fff is weakly quasi M-convex, satisfying (QMw), if for every pair of distinct points x,y∈dom⁡fx, y \in \operatorname{dom} fx,y∈domf there exist uuu in the positive support and vvv in the negative support of x−yx - yx−y with f(x−χu+χv)≤f(x)f(x - \chi_u + \chi_v) \le f(x)f(x−χu​+χv​)≤f(x) or f(y+χu−χv)≤f(y)f(y + \chi_u - \chi_v) \le f(y)f(y+χu​−χv​)≤f(y) — an "or" where (M-EXC[Z]) demands an additive inequality. Two further conditions restrict attention to points of different function value and sharpen the conclusion to a three-way trichotomy (strictly better on one side, or exactly tied on both): (SSQM≠_{\ne}=​) quantifies universally over uuu (as in (M-EXC[Z])), while (SSQM≠,w_{\ne,w}=,w​) quantifies existentially over both uuu and vvv (as in (QMw)). The linear perturbation of fff by p:V→Rp : V \to \mathbb Rp:V→R is 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 6.78 (the quasi M-proximity theorem)

Let fff satisfy (SSQM≠_{\ne}=​), n=∣V∣n = |V|n=∣V∣, α\alphaα a positive integer. If xα∈dom⁡fx_\alpha \in \operatorname{dom} fxα​∈domf satisfies 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 u,v∈Vu, v \in Vu,v∈V, then arg⁡min⁡f≠∅\arg\min f \ne \emptysetargminf=∅ and there is x∗∈arg⁡min⁡fx^* \in \arg\min fx∗∈argminf with ∥xα−x∗∥∞≤(n−1)(α−1)\|x_\alpha - x^*\|_\infty \le (n-1)(\alpha - 1)∥xα​−x∗∥∞​≤(n−1)(α−1) — verbatim the same conclusion, and the same exact bound, as chunk 06's Theorem 6.37(1), now established for the strictly larger class satisfying (SSQM≠_{\ne}=​) rather than the M-convex exchange axiom itself.

Milestones: Theorems 6.68(2), 6.76, 6.77

Theorem 6.68(2): fff satisfies (M-EXC[Z]) if and only if every linear perturbation f[p]f[p]f[p] satisfies (QMw) — quantifying exactly how much weaker (QMw) is pointwise, and how the gap closes once quantified over every perturbation. Theorem 6.76 (the quasi M-optimality criterion): the direct analogue of chunk 06's Theorem 6.26 for the quasi-convexity classes — a purely pairwise local check still characterizes global (or, in the (QMw) case, strict unique) optimality. Theorem 6.77 (the quasi M-minimizer cut): chunk 06's Theorem 6.28 continues to hold verbatim when its M-convexity hypothesis is replaced by (SSQM≠_{\ne}=​) — the structural fact the proximity theorem's proof is built from survives the relaxation intact.

Significance

The result itself. The proximity theorem is the result algorithms actually use: a scaling algorithm for minimizing quasi-convex functions of this kind inherits exactly the same correctness guarantee, with exactly the same distance bound, as the M-convex case — this is a genuine broadening of chapter 10's algorithmic reach, not a restatement dressed in weaker hypotheses. Every strictly increasing scalar transformation of an M-convex objective (a common modeling device — re-expressing a cost in utility units, or applying a monotone risk measure) now falls under a proximity theorem, whereas prior to this chapter's relaxation such a transformation would generally destroy M-convexity itself and leave optimization theory silent on the transformed problem.

Formalizing it. No matching item exists on the platform for quasi M-convexity in any of its forms. Formalizing Theorem 6.78 requires first pinning down (SSQM≠_{\ne}=​) exactly (there are six closely related axiom variants in this section of the book, only three of which — (QMw), (SSQM≠_{\ne}=​), (SSQM≠,w_{\ne,w}=,w​) — are needed for this mission's chosen results), and this mission also captures, via Theorem 6.68(2), the precise sense in which these relaxed conditions are strictly weaker than plain M-convexity while remaining tightly connected to it.

Difficulty

The natural first instinct, given how close the quasi-convexity axioms look to (M-EXC[Z]), is to try to prove Theorem 6.78 by directly imitating chunk 06's proof of Theorem 6.37 line by line. This mostly works — the proof structure (fix a target coordinate, build a chain of strictly decreasing values via repeated exchange steps, bound the chain's length using the scaled hypothesis) survives verbatim — but every step that chunk 06's proof took by adding two instances of the exchange inequality together must be replaced by an ordinal argument, since (SSQM≠_{\ne}=​) only ever asserts a disjunction of value comparisons, never an additive inequality relating four function values simultaneously the way (M-EXC[Z])'s 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​) does. The book's proof handles this by working with strict inequalities and the trichotomy structure of (SSQM≠_{\ne}=​) directly rather than algebraic cancellation — the same overall architecture, but every arithmetic step rebuilt as a case analysis on which disjunct of (SSQM≠_{\ne}=​) fires.

Formalization scope

This mission builds directly on chunk 06's published items (CharVec, SuppPos, SuppNeg, DomZ, MExchangeAxiom, ArgMin), per the platform's textbook convention that a later chapter of the same book imports an earlier one's definitions rather than redrafting them; its own namespace DiscreteConvex.MConvexFunctions.Quasi nests under chunk 06's DiscreteConvex.MConvexFunctions accordingly. Δf(z;v,u) (Eq. (6.2)) is never reified as a separate object; every occurrence is unfolded directly into an f-value comparison, avoiding WithTop ℝ subtraction throughout, consistent with chunk 06's own convention.

A trivializing formalization of the goal would silently strengthen (SSQM≠_{\ne}=​) back to plain M-convexity (making this mission redundant with chunk 06's Theorem 6.37) or loosen the exact bound (n−1)(α−1)(n-1)(\alpha-1)(n−1)(α−1) to an unspecified function of n,αn, \alphan,α; neither is done. Six axiom variants appear in this section of the book ((QM), (SSQM), (QMw), (SSQMw_ww​), (SSQM≠_{\ne}=​), (SSQM≠,w_{\ne,w}=,w​)); only the three actually needed by this mission's four items are drafted, and Theorem 6.68's first part (an implication chain among the other three) is left out — see MODERATION_NOTES.md. Contributions building the polyhedral M-convex-function bridge (§6.11–6.12, Theorems 6.59–6.64), the level-set characterizations (Theorems 6.72, 6.74), or the scaled quasi M-minimizer cut (Theorem 6.79, the direct generalization of Theorem 6.77 drafted here) are welcome.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • M. Avriel, W. E. Diewert, S. Schaible, I. Zang, Generalized Concavity, Plenum Press, 1988.
8 thms2 active usersReviewed
CombinatoricsConvex OptimizationDiscrete Geometry+2·Captain: Shuze Chen

Discrete Convex Analysis XXI: The Exchange Axiom as Local OptimalityTextbook

Motivation

Convexity on the integer lattice cannot be defined by the classical secant-line inequality alone: a function can be midpoint-convex along every line and still admit no useful global optimality theory, because integer points off a chosen line are invisible to it. M-convex functions, introduced by Murota, resolve this by replacing the secant condition with an exchange axiom directly generalizing the basis-exchange property of matroids and the convex-hull structure of network flows: a function on the integer lattice is M-convex if, whenever two points can be improved by moving one coordinate up and a compensating coordinate down, at least one such move weakly improves the sum of the two function values. This single axiom turns out to be equivalent to several strikingly different-looking properties — invariance under a wide family of domain operations, supermodularity in the M♮ (translation-invariant) case, and, most importantly, a local-to-global optimality principle: a point is a global minimizer of an M-convex function if and only if no single coordinate exchange improves it. This mission develops the algebraic core of that theory — the exchange axiom's basic consequences, its equivalent local and dynamic reformulations, and the operations that preserve it — building toward the theorem that recasts M-convexity itself as an algorithmically meaningful local-search guarantee.

Companion mission 06-mconvex-functions-i (Discrete Convex Analysis V) covers this chapter's own primary line of development: the equivalence of M-convexity and M♮-convexity with their respective exchange axioms (Theorem 6.2), the M-optimality criterion (Theorem 6.26), a minimizer-cut lemma (Theorem 6.28), and the M-proximity theorem (Theorem 6.37, 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 that chapter leaves for a second pass: the domain structure of M- and M♮-convex functions, worked examples (quadratic forms, quasi-separable functions), the operations that preserve M-convexity, supermodularity of the M♮-convex case, the descent-direction property, and — this mission's goal — the equivalence of the exchange axiom with a dynamic sequential-improvement property.

Setting

Fix a finite ground set VVV. 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 M-convex if it satisfies the exchange axiom (M-EXC[Z]): for x,y∈dom⁡fx, y \in \operatorname{dom} fx,y∈domf and u∈supp⁡+(x−y)u \in \operatorname{supp}^+(x-y)u∈supp+(x−y) (coordinates where xxx exceeds yyy), there is v∈supp⁡−(x−y)v \in \operatorname{supp}^-(x-y)v∈supp−(x−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​).

Writing 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 (a lift to one extra coordinate), fff is M♮^\natural♮-convex if f~\tilde ff~​ is M-convex; M♮-convexity is a genuine generalization of M-convexity (every M-convex function is M♮-convex, but not conversely) and coincides with it exactly when dom⁡f\operatorname{dom} fdomf lies on a single hyperplane. The linear-weighted function f[p](x)=f(x)−⟨p,x⟩f[p](x) = f(x) - \langle p, x \ranglef[p](x)=f(x)−⟨p,x⟩ (for p∈RVp \in \mathbb R^Vp∈RV) is the standard device for testing local optimality under an arbitrary reweighting.

Formalization targets

Goal: the exchange axiom as sequential improvement

f is M-convex  ⟺  ∀p∈RV, ∀x,y∈dom⁡f, f[p](x)>f[p](y)  ⟹  f[p](x)>min⁡u∈supp⁡+(x−y) min⁡v∈supp⁡−(x−y)f[p](x−χu+χv),f \text{ is M-convex} \iff \forall p \in \mathbb R^V,\ \forall x, y \in \operatorname{dom} f,\ f[p](x) > f[p](y) \implies f[p](x) > \min_{u \in \operatorname{supp}^+(x-y)}\ \min_{v \in \operatorname{supp}^-(x-y)} f[p](x - \chi_u + \chi_v),f is M-convex⟺∀p∈RV, ∀x,y∈domf, f[p](x)>f[p](y)⟹f[p](x)>u∈supp+(x−y)min​ v∈supp−(x−y)min​f[p](x−χu​+χv​),

with the analogous statement for M♮-convexity (Theorem 6.24). This is the weakest stable form: it makes no reference to a specific algorithm, only to the existence of an improving single exchange whenever the current point is suboptimal under any linear reweighting — a property a faster algorithm could exploit without invalidating the characterization itself.

Supporting structural targets

Eleven further results build the vocabulary and toolkit this goal draws on: the domain structure of M-convex and M♮-convex functions (Propositions 6.1, 6.7), the equivalence of the exchange axiom with a local, bounded-distance version (Theorem 6.4), worked examples establishing M-convexity for quadratic forms, univariate, conservation-law, and quasi-separable functions (Propositions 6.8-6.9), the domain and range operations preserving M-convexity (Theorem 6.13, Proposition 6.14), supermodularity of the M♮-convex case (Theorem 6.19), the descent-direction property (Proposition 6.23) that Theorem 6.24 generalizes, and a discrete subgradient inequality (Proposition 6.25).

Significance

Theorem 6.24 is the bridge between the static exchange axiom (a property of function values at pairs of points) and the dynamic behavior of local-search algorithms: it says a greedy single-coordinate-exchange step, applied to any linearly reweighted version of an M-convex function, always finds a strict improvement when one exists. This is exactly the guarantee that makes steepest-descent-type algorithms for M-convex function minimization correct, and it is the theorem chapter 10's algorithmic analysis (Schrijver-type methods) relies on implicitly whenever it argues that local exchange steps make global progress. The descent-direction property (Proposition 6.23) is the special case p=0p=0p=0, isolating the core combinatorial fact before the reweighting machinery is added. The operations catalog (Theorem 6.13) is the practical toolkit that lets later chapters build complex M-convex functions (network flow costs, matroid rank functions composed with linear maps) from simple pieces without re-verifying the exchange axiom from scratch each time.

None of these results are open — they are Murota's own systematic development of the exchange- axiom theory, with worked examples drawn from classical quadratic and separable function theory. What this mission contributes is a faithful, machine-checked formal statement of each, sharing the Lean vocabulary (MExchangeAxiom, MNaturalConvex, LinearWeight) the rest of the Discrete Convex Analysis series builds on; no comparable formalization exists on the platform (see Formalization scope).

Difficulty

The forward direction of Theorem 6.24 (M-convex   ⟹  \implies⟹ sequential improvement) follows in one step from Proposition 6.23 applied to f[p]f[p]f[p], itself M-convex by Theorem 6.13(3) — routine once those two pieces are in hand. The converse is the substantial direction: it must derive the full static exchange axiom from a property that only ever exhibits some improving exchange at some linear weighting, for every pair of suboptimal points — the proof constructs an explicit adversarial weighting ppp designed so that failure of the local exchange step at that specific ppp forces the domain itself to be M-convex (via Theorem 4.3) and then forces the local exchange axiom (M-EXCloc[Z]) via a bipartite-matching argument on the coordinates that differ, finally invoking Theorem 6.4 to lift locality to the full exchange axiom. No shortcut bypasses this two-stage reduction (domain structure, then local exchange) — attempting to verify (M-EXC[Z]) directly from (M-SI[Z]) without first pinning down that dom⁡f\operatorname{dom} fdomf is M-convex fails because the exchange axiom's own statement presupposes a well-structured domain.

Formalization scope

Ground-set elements are a Fintype V with DecidableEq; functions are (V → ℤ) → WithTop ℝ. SuppPos/SuppNeg are Finset V (not Set V), matching how the (M-SI[Z])/(M♮-SI[Z]) axioms and the descent-direction property use Finset.inf, whose value on an empty index set is ⊤ — exactly the book's own stated convention for an empty minimum. No Module ℝ or ConvexOn machinery is used for WithTop ℝ-valued arithmetic; scalar actions by positive reals (PosScalarMul, Theorem 6.13(1)) and by naturals (FCheck's flow coefficients, Proposition 6.25) are built directly from the native order and AddMonoid structure. No numeric constants are hard-coded anywhere in this mission (rule 7 is vacuous); Proposition 6.8's quadratic-form conditions are stated with the book's own literal coefficients (000, and the min/≥ structure of Eq. (6.25)-(6.28)), not a special case. Theorem 6.13's parts (7) (aggregation) and (8) (integer infimal convolution) are not restated here since the book itself proves them only later via Chapter 9's network-transformation machinery — see the Difficulty note and MODERATION_NOTES.md; this is not a trivializing omission, since the six operations that are included already exercise every domain- and range-transformation technique this mission's goal needs. This mission's definitions (MExchangeAxiom, MNaturalConvex, CharVec, DomZ, SuppPos, SuppNeg, CharVecOpt) are redeclared from chunk 06-mconvex-functions-i 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's converse direction and Theorem 6.13's operations are the two with 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 exchange axiom and its equivalent local/dynamic reformulations).
38 thms3 active usersReviewed
CombinatoricsConvex OptimizationDiscrete Geometry+2·Captain: Shuze Chen

Discrete Convex Analysis XXII: Convex Extensibility of M-Convex FunctionsTextbook

Motivation

A discrete function defined only on the integer lattice cannot, by itself, be minimized by the tools of continuous optimization — gradients and convexity in the classical sense simply do not apply. Murota's theory of M-convex functions closes this gap by showing that the exchange axiom alone, a purely combinatorial condition, is enough to guarantee that a discrete function behaves exactly like a convex one: its minimizers form a well-structured (M-convex) set, it can be extended to a genuine convex function on real space without gaining any new local minima, and its behavior under a change of price vector (in the economic interpretation where the function is a cost and its argument a bundle of goods) satisfies the same gross substitutes law economists have studied since Kelso and Crawford's matching-market models. This mission develops the second half of that connection: from local optimality (established in the companion mission) to the full structural picture — minimizer sets, price-substitution laws, and the extension of M-convex functions to genuine convex functions in real variables.

Companion mission 06-mconvex-functions-i (Discrete Convex Analysis V) and sibling mission 22-ch06b-mconvexfunctions (Discrete Convex Analysis XXI) cover this chapter's optimality theory (the M-optimality criterion, the exchange axiom as sequential improvement) and its algebraic toolkit (domain operations, worked examples). This mission builds the vocabulary those results also need (redeclared here, since sibling drafts cannot yet import one another) and proves the results on minimizer structure, gross substitutability, and convex extension that this chapter's remaining sections develop: the M-convexity of minimizer sets, the gross substitutes and stepwise gross substitutes properties and their characterizing role, a minimizer-cut theorem with scaling, integral convexity of M♮-convex functions, and — this mission's goal — the theorem characterizing M-convexity entirely through the polyhedral structure of a function's convex extension.

Setting

Fix a finite ground set VVV. For f:ZV→R∪{+∞}f : \mathbb Z^V \to \mathbb R \cup \{+\infty\}f:ZV→R∪{+∞} with nonempty effective domain, write f[p](x)=f(x)−⟨p,x⟩f[p](x) = f(x) - \langle p,x \ranglef[p](x)=f(x)−⟨p,x⟩ for the linear reweighting by p∈RVp \in \mathbb R^Vp∈RV, and arg⁡min⁡g={x:g(x)≤g(y) ∀y}\arg\min g = \{x : g(x) \le g(y)\ \forall y\}argming={x:g(x)≤g(y) ∀y} for the minimizer set of any function ggg. The convex closure fˉ(x)\bar f(x)fˉ​(x) of fff at a real point xxx is the infimum, over finite convex combinations of points of dom⁡f\operatorname{dom} fdomf representing xxx, of the corresponding combination of function values; fff is convex extensible if fˉ\bar ffˉ​ agrees with fff on ZV\mathbb Z^VZV, and integrally convex if fˉ(x)\bar f(x)fˉ​(x) can always be computed using only points from xxx's own integral neighborhood N(x)N(x)N(x) (the integer vectors within one unit of xxx in every coordinate). 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 for every α∈[0,α0]\alpha \in [0,\alpha_0]α∈[0,α0​].

Formalization targets

Goal: convex extensibility characterizes M-convexity

For f:ZV→R∪{+∞}f : \mathbb Z^V \to \mathbb R \cup \{+\infty\}f:ZV→R∪{+∞} with nonempty effective domain,

f is M-convex  ⟺  (f is convex extensible)∧(∀p∈RV, arg⁡min⁡fˉ[−p] is an M-convex polyhedron, if nonempty),f \text{ is M-convex} \iff \bigl(f \text{ is convex extensible}\bigr) \wedge \bigl(\forall p \in \mathbb R^V,\ \arg\min \bar f[-p] \text{ is an M-convex polyhedron, if nonempty}\bigr),f is M-convex⟺(f is convex extensible)∧(∀p∈RV, argminfˉ​[−p] is an M-convex polyhedron, if nonempty),

with the M♮-analogue using M♮-convex polyhedra (Theorem 6.43). This is the weakest stable form: it characterizes M-convexity purely by properties of the (unique) convex closure, without reference to any specific algorithm for computing it or any bound on the polyhedron's complexity.

Supporting structural targets

Ten further results build the toolkit this goal draws on and the picture it completes: the M-convexity of minimizer sets (Proposition 6.29), the gross substitutes and stepwise gross substitutes properties and the theorems showing they characterize M-convexity and M♮-convexity among convex-extensible functions (Propositions 6.32-6.33, 6.35, Theorems 6.34, 6.36), a minimizer-cut theorem with scaling used algorithmically in Chapter 10 (Theorem 6.39), integral convexity of M♮-convex functions (Theorem 6.42), a shared-coefficient convex-combination theorem for pairs of M♮-convex functions used in Chapter 8's separation theorem (Theorem 6.44), and the polyhedral-M-convexity of an M-convex function's convex extension together with the correspondence between polyhedral M♮-convexity and the real exchange axiom (Theorems 6.45, 6.47).

Significance

Theorem 6.43 is what makes the whole edifice of M-convex function theory a genuine extension of M-convex set theory (chapters 4-5) rather than a separate parallel development: it says that knowing a function's convex extension is polyhedral, with every price-weighted minimizer set an M-convex polyhedron, is not merely a consequence of M-convexity but an exact characterization of it. This is the theorem that lets later results (the discrete conjugacy theorem of Chapter 8, the separation theorems for M♮-convex functions) move freely between the discrete and continuous pictures. The gross substitutes property (Propositions 6.32-6.36) is independently significant outside this book: it is the exact condition, discovered independently in mathematical economics (Kelso-Crawford, Gul-Stacchetti), under which competitive equilibria with indivisible goods are guaranteed to exist — Murota's theorem that gross substitutability characterizes M-convexity (among convex-extensible functions) is what unifies the economic and combinatorial literatures on this question, taken up again in Chapter 11.

None of these results are open — they are Murota's systematic account of a theory with roots in matroid theory, submodular optimization, and mathematical economics. What this mission contributes is a faithful, machine-checked formal statement of each, extending the shared Lean vocabulary (MExchangeAxiom, ConvexClosureVal, ArgMinOn) the Discrete Convex Analysis series builds on; no comparable formalization exists on the platform (see Formalization scope).

Difficulty

The forward direction of Theorem 6.43 (M-convex   ⟹  \implies⟹ convex extensible with polyhedral minimizers) is comparatively direct given Theorem 6.42 and Proposition 6.29. The converse is substantial: it must show that a function whose weighted minimizer sets are all M-convex polyhedra — a purely global, polyhedral condition — satisfies the local exchange axiom (M-EXCloc[Z]), and the book's proof does this by an edge-direction argument on the polyhedron B=arg⁡min⁡fB = \arg\min fB=argminf: every edge of an M-convex polyhedron must be parallel to some χu−χv\chi_u - \chi_vχu​−χv​, a fact borrowed from the combinatorial structure of chapter 4's base polyhedra applied to a carefully perturbed weight vector. No shortcut through convex analysis alone succeeds, because ordinary polyhedral theory says nothing about which combinatorial directions a polyhedron's edges must follow — that content comes entirely from the M-convexity of the minimizer sets, not from convexity of the closure by itself.

Formalization scope

Ground-set elements are a Fintype V with DecidableEq; functions are (V→ℤ)→WithTop ℝ (integer domain) or (V→ℝ)→WithTop ℝ (real domain, for the polyhedral theorems). The convex closure is built directly from finite convex-combination representations rather than an abstract closure operator, and integral convexity compares it against the same construction restricted to each point's integral neighborhood (Fintype.piFinset of per-coordinate Finset.Icc). Real M-convex/M♮-convex polyhedra are defined as convex hulls of M-convex/M♮-convex integer sets, reusing chapters 4-5's own characterization. The real-variable exchange axioms (Theorems 6.45, 6.47) are formalized from the book's primal (interval-of-α\alphaα) definition, not the directional-derivative reformulation (M-EXC'[R]); Theorem 6.47's own three-way equivalence is correspondingly stated with only its first two legs (see Difficulty and MODERATION_NOTES.md/HARD.md — this is a documented scope choice, not a trivializing omission, since the six results using the primal axiom already exercise the chapter's real- variable machinery in full). No numeric constants are hard-coded anywhere in this mission beyond the book's own literal coefficients in Theorem 6.39's cut bound ((n-1)(α-1)). This mission's definitions are redeclared from chunks 06-mconvex-functions-i and 22-ch06b-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's converse direction and Theorem 6.44's shared-coefficient construction carry the most independent proof content.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • A. S. Kelso Jr. and V. P. Crawford, "Job matching, coalition formation, and gross substitutes," Econometrica, 50 (1982), pp. 1483-1504.
  • F. Gul and E. Stacchetti, "Walrasian equilibrium with gross substitutes," Journal of Economic Theory, 87 (1999), pp. 95-124.
47 thms4 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
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