High-Dimensional Probability I: Approximate Carathéodory's TheoremTextbook
Motivation
Many arguments in high-dimensional geometry, statistics and computer science need to approximate a point of a convex set by an average of a handful of extreme points, rather than represent it exactly. The classical Carathéodory theorem (1907) answers the exact question: every point of the convex hull of a set is a convex combination of at most points of . That bound is tight and grows with the dimension , which makes it useless whenever is large — exactly the regime of interest in high-dimensional probability.
B. Maurey's empirical method — an unpublished 1980–81 result reported by G. Pisier, "Remarques sur un résultat non publié de B. Maurey," Séminaire d'Analyse Fonctionnelle 1980–1981 — and later applied by B. Carl to bound covering numbers of operators between Banach spaces (Inequalities of Bernstein-Jackson-type and the degree of compactness of operators in Banach spaces, Ann. Inst. Fourier 35(3), 1985, 79–118), replaces the exact question with an approximate one and removes the dimension dependence entirely: to approximate to accuracy , the number of points needed depends only on , never on . Vershynin's High-Dimensional Probability opens with this result as its "Appetizer," using it to illustrate the book's central theme — that randomness is a tool for constructing deterministic combinatorial objects — before any probabilistic machinery has been introduced.
Setting
A convex combination of finitely many points is a sum with and . The convex hull of a set is the set of all convex combinations of all finite collections of points of . The diameter of is , the Euclidean norm throughout.
The classical Carathéodory theorem states that every is a convex combination of at most points of — with generally unavoidable, attained by a simplex. The question this mission answers is different: given that we are willing to approximate rather than represent it exactly, and willing to use only combinations with equal coefficients (an average of points, with repetition allowed), how large must be as a function of the desired accuracy?
Formalization targets
Goal — Theorem 0.0.2, Approximate Carathéodory's theorem
The quantifiers are exactly this order: for every bounded , every point of its convex hull, and every target , such an averaging set exists. This is the weakest stable statement carrying the theorem's content — the number of points does not depend on the dimension , and the coefficients are forced to be uniform — and it is the form the corollary below invokes directly.
Milestone — Corollary 0.0.4, Covering polytopes by balls
This is a direct application of the goal to computational geometry's covering problem: how many balls of radius are needed to cover a polytope, and where should they be centered.
Significance
The result itself. The approximate Carathéodory theorem is the prototype of a dimension-free approximation result: whenever a set is bounded, a fixed number of points (depending only on the target accuracy, not the ambient dimension) suffices to approximate any point of its convex hull. This is what makes possible dimension-independent covering-number bounds such as Corollary 0.0.4, which in turn are the starting point for the book's later treatment of entropy, packing and generic chaining (Chapters 4, 7–8). The technique generalizes far beyond : it underlies covering-number bounds for operators between Banach spaces (Carl's original application) and is a recurring device in learning theory for bounding the size of an -net of a hypothesis class.
Formalizing it. Both results are elementary and already fully proved in the literature; no open mathematical content remains. What this mission contributes is a machine-checked, faithful Lean statement of Maurey's construction and its corollary, phrased over Mathlib's existing convex-hull and metric-diameter machinery, so that later missions in this series (concentration inequalities, Johnson–Lindenstrauss, chaining) can build on a verified base case of "probability constructs a deterministic covering," and so that the empirical method itself becomes a reusable, linked component on the platform. The proof of the goal (via the probabilistic argument sketched by the book: interpret a convex combination as a probability distribution, average i.i.d. samples, and bound the variance) is left open for solvers.
Difficulty
The identity that makes the proof work — averaging independent copies of a random vector concentrates around its mean at rate in mean-square — is a two-line computation once the convex combination is reinterpreted probabilistically. The step that is easy to miss is this reinterpretation itself: nothing in the statement mentions probability, so the "obvious" attack of manipulating the convex-combination weights directly, or trying to construct by some explicit combinatorial recipe, does not see a path to a bound independent of . The probabilistic argument produces the points non-constructively, via an averaging/existence argument (the expected squared distance is small, so some realization achieves it) rather than an explicit formula — a solver has to introduce a probability space and a random vector that does not appear anywhere in the formal statement to be proved.
Formalization scope
Both results are stated over EuclideanSpace ℝ (Fin n) for an explicit dimension n : ℕ, so
‖·‖ is the Euclidean norm and Mathlib's Metric.diam is used directly for
. The convex hull is Mathlib's convexHull ℝ T; by
Mathlib's convexHull_eq, this already coincides with the book's own definition of a convex
combination of finitely many points of , so no bespoke convex-combination definition is
introduced — this mission needs no supporting definitions of its own. In the corollary, "a
polytope with vertices" is formalized, following the book's own proof, as
for a finite vertex set with , rather than via a
separate Polytope structure (which Mathlib does not provide and the book's argument does not
need). The covering bound is an exponent, not a product with
— matching the book's own proof, which counts the ordered -tuples of vertices with
repetition, ; the typeset "" in
the corollary statement is the same juxtaposition-as-exponent notation the proof uses for
"" one line earlier.
A trivializing formalization is ruled out explicitly: the goal must hold for every integer and every , not merely some convenient choice — e.g. together with trivially satisfies the inequality but proves nothing about the theorem's actual content, that a fixed, dimension-independent works uniformly over all points of the hull. The formal statement quantifies , then , then , and only then asserts existence of the , exactly in that order.
Classical Carathéodory (Theorem 0.0.1, stated for context in the source but not used by
either formalized result's proof) is not drafted here: Mathlib already proves the corresponding
statement via affine independence
(Caratheodory.eq_pos_convex_span_of_mem_convexHull, Analysis/Convex/Caratheodory.lean),
from which the book's " points" bound follows via
AffineIndependent.card_le_finrank_succ. It is not added as a kind: reference milestone
because no corresponding theorem is yet published on the Prove2Me platform to point at (checked
2026-09-17: GET /theorems?q=Caratheodory returns only unrelated tropical-convexity results),
and re-drafting existing Mathlib content as a new platform theorem would duplicate rather than
reuse it.
Selected references
- R. Vershynin, High-Dimensional Probability: An Introduction with Applications in Data Science, Cambridge University Press, 2018, DOI 10.1017/9781108231596, Appetizer (pp. 1–5).
- G. Pisier, "Remarques sur un résultat non publié de B. Maurey," Séminaire d'Analyse Fonctionnelle (Maurey–Schwartz), 1980–1981, exposé no. 5. numdam.org/item/SAF_1980-1981____A5_0
- B. Carl, "Inequalities of Bernstein-Jackson-type and the degree of compactness of operators in Banach spaces," Annales de l'Institut Fourier, 35(3), 1985, 79–118. numdam.org/item/AIF_1985__35_3_79_0