Loading home page…
Each mission turns a result from a paper or textbook into small Lean 4 statements anyone can tackle.
Campaigns group missions around a shared mathematical goal. Each one tracks a quantity, such as an upper or lower bound. Have a good candidate in mind? Ping us on Slack, Zulip, or WeChat.
Each mission turns a result from a paper or textbook into small Lean 4 statements anyone can tackle.
Campaigns group missions around a shared mathematical goal. Each one tracks a quantity, such as an upper or lower bound. Have a good candidate in mind? Ping us on Slack, Zulip, or WeChat.
Quantum Max-Cut asks for the largest eigenvalue of a graph's quantum spin Hamiltonian. For a graph with nonnegative rational edge weights, write
This campaign tracks expected approximation ratios for classical polynomial-time energy estimation. One uniform randomized classical algorithm must work for every input graph. For each fixed , let be the algorithm's internal fair random bits, and let be its nonnegative rational output. A ratio requires
The expectation averages the algorithm's internal randomness on a fixed graph. Its running time is bounded by a polynomial on every run. Higher ratios are better.
See the Quantum MaxCut reference for the problem definition and research context.
How much larger can block sensitivity be than sensitivity for a total Boolean function?
For a total Boolean function , a block of coordinates is sensitive at if flipping every bit in that block changes . Let count sensitive singletons and be the maximum number of pairwise disjoint sensitive blocks. Define and .
This campaign seeks larger exponents for which there exists an infinite family of total Boolean functions , indexed by , satisfying
Here denotes a quantity tending to zero as . Higher is better. The classical exponent is , and Huang's bound gives a ceiling of .
OpenAI's construction achieves an exponent strictly greater than .
One sufficient way to establish an exponent is to find a single total Boolean function with and . Input and output normalization followed by self-composition then yields the required asymptotic family.
Turn proposed improvements to integer multiplication into complete Lean proofs, and push the exponent saving further.
Harvey and van der Hoeven established an algorithm in 2021. This campaign builds on that foundation, the OpenAI manuscript, and subsequent community constructions to pursue a strict asymptotic improvement.
For two -bit integers, the target is
A positive beats asymptotically; larger is better. Every entry must exhibit one deterministic multitape Turing machine, with a fixed finite alphabet and tape count, that computes the exact product at every positive input length and meets the eventual worst-case time bound. The tracked number measures an asymptotic exponent saving.
Classical algorithms solve 3SUM in time. In a 2026 breakthrough, Alman and Vassilevska Williams gave a deterministic algorithm, refuting the integer 3SUM hypothesis. How low can the exponent go?
Building on existing Lean formalizations, this campaign tracks upper bounds for 3SUM on polynomially bounded integers, using a word RAM with -bit words, and pursues smaller exponents.
Classical algorithms solve all-pairs shortest paths in time. In a 2026 breakthrough, Alman and Vassilevska Williams refuted the APSP conjecture with a deterministic algorithm. How low can the exponent go?
Building on existing Lean formalizations, this campaign tracks upper bounds for exact APSP and pursues smaller exponents.
The irrationality measure of π quantifies how closely rational numbers can approximate it. This campaign seeks formal proofs of sharper upper bounds, starting with Mahler’s bound of 42.
The sharp Hlawka inequality for Schatten -norms is a cousin of the triangle inequality: it relates the norms of three matrices to the norms of their pairwise sums and their total sum. For complex diagonal matrices, an exact formula for the best possible comparison constant has been proved in Lean for every real . We conjecture that the same formula holds for all .
What is the smallest cutoff for which this formula holds for every real ?
References:
Is every odd number a sum of primes? This campaign tracks formalized proofs of the smallest that suffices.
Schnirelmann (1930) showed some finite works. Vinogradov (1937) showed that three is enough for all sufficiently large odd numbers. Tao (2012) proved unconditionally. Helfgott (2013) proved that every odd number greater than is a sum of three primes, though the proof is still unrefereed. Ideally, we can formalize this statement here. Note that three is optimal: is neither prime nor + prime.
Schoolbook matrix multiplication takes operations. The exponent is the infimum of all such that two matrices can be multiplied in arithmetic operations; trivially , and is conjectured but open.
Strassen gave the first nontrivial bound, , in 1969, and introduced the laser method in 1986 to reach . Coppersmith and Winograd's 1990 bound of stood for two decades. Every subsequent improvement comes from analyzing higher tensor powers of their construction with refined laser-method variants. That line reached in 2025, and the current record is , from August 2026. See Computational complexity of matrix multiplication for the full table. Can we formalize these results and even improve on them?
A retailer marking down a seasonal line, an airline raising fares as seats sell, a manufacturer pricing while restocking: each sets prices over time against a finite and changing inventory. Chapter 5 of Talluri and van Ryzin's The Theory and Practice of Revenue Management (2004) collects the structural theory of that problem. Without replenishment, the Bernoulli-arrival model of Gallego and van Ryzin (1994) gives a marginal value of capacity that falls with inventory and with time, hence prices that jump up at each sale and drift down between sales, and the deterministic fluid model bounds it from above. With replenishment, the model of Federgruen and Heching (1999) has a jointly concave and supermodular continuation value, from which the base-stock, posted-price policy follows: below a base-stock level order up to it and post a fixed price, above it order nothing and discount, more the higher the inventory. This mission formalizes those results with Proposition 5.3 as the goal.
Bernoulli demand (Sect. 5.2.2.2). One customer arrives per period with a random
willingness to pay; the firm's decision is the demand rate , the probability of a
sale at the inverse-demand price , with revenue rate
(revenueRate). The value function (5.12) is
, ,
(bernoulliValue), and (bernoulliDelta). The
deterministic model (5.1) maximizes over rates with
(deterministicValue).
Pricing with replenishment (Sect. 5.3.2). Inventory may be negative (backorders). In
period with inventory the firm orders up to at unit cost , chooses a
rate , sells the random demand
(additive, multiplicative or mixed noise), and pays the convex cost on ending inventory.
The value function (5.20) is
with the continuation value
(ReplPricing.value, contValue). Assumption 7.2, marginal revenue decreasing, is the
concavity of .
For every period, is jointly concave on , is
concave on , and has increasing differences in , the supermodularity
that the book states through its partial derivatives: replenishment_concave_supermodular.
Proposition 5.2, the marginal value of capacity in the Bernoulli model decreases in and in ; the upper bound of Sect. 5.2.2.3, the optimal deterministic revenue dominates the optimal expected stochastic revenue; and the base-stock, posted-price structure of Sect. 5.3.2.1, derived from Proposition 5.3: below the unconstrained optimum order up to it and use , above it order nothing and use a rate at least that is nondecreasing in the inventory.
Proposition 5.1 and Lemma 5-5.A.1 (the continuous-demand model without replenishment) are not targets; see the formalization scope. The deterministic sections (efficient prices, discrete price sets), the asymptotic optimality of the deterministic heuristic, the infinite-horizon stationary problem and the multiproduct and finite-population models carry no numbered results.
Proposition 5.3 is the structural core of joint pricing and inventory control: joint concavity makes the period problem a concave program, and supermodularity is what turns its solution into a policy, the base-stock, posted-price rule that Federgruen and Heching showed optimal and that later work on pricing with inventory builds on. Proposition 5.2 is the reason optimal dynamic prices in the stochastic single-item model rise at every sale and fall while inventory sits, the behaviour of Figure 5.5, and the deterministic upper bound is what justifies the fluid model as a benchmark and a heuristic, the pattern quantified in Table 5.6. None of these has a machine-checked proof; the replenishment result in particular needs the interplay of concavity, expectation and partial maximization on all of .
Proposition 5.2 is an induction whose step compares suprema over the rate interval, with the boundary condition breaking the recursion at and requiring . The deterministic bound is an induction on periods that uses the concavity of the deterministic value in the inventory (a concave program's value) to absorb the two branches of the Bernoulli recursion. The goal needs: integrability and continuity of under bounded noise; that the expectation of a concave function of an affine map is jointly concave, and its increasing differences from those of the concave integrand; that a partial supremum of a jointly concave function over the convex feasible set is concave in ; and the boundedness of the objective so that every supremum is a real number. The base-stock item is the segment argument that moves an unconstrained maximizer onto the boundary and a monotone comparative-statics argument on the supermodular objective, with maxima attained by continuity on the compact rate interval.
Periods are natural numbers with value t the value with periods to go, the
maxima are suprema, and the book's ranges are hypotheses. The demand is affine in the rate,
which is the additive and multiplicative models the book names; with a merely convex demand
(Assumption 5.1) the joint concavity of Proposition 5.3 fails when is not
monotone, and the noise has bounded support, strengthening Assumption 7.6. The
partial-derivative statements (iii)-(iv) are in difference form. Proposition 5.1 is not
formalized: its model (5.11) evaluates at negative inventories the model
does not define while truncating revenue at , and its Lemma 5-5.A.1 (joint concavity of
) is false as stated, its Hessian argument mistaking an indefinite matrix for a negative
definite one; a counterexample is in the mission's check. The deterministic model restricts
rates to , the rates the Bernoulli model can realize.
Every model in the earlier chapters of this series treats demand as the uncertain quantity and supply as given. Chapter 9 of Snyder and Shen's Fundamentals of Supply Chain Theory (2019) reverses the roles: demand is deterministic and the supplier fails. A disruption is a binary event, modeled as a two-state Markov process between up and down periods, during which nothing can be ordered. The chapter's thesis, from Snyder and Shen (2006), is that supply uncertainty is a mirror image of demand uncertainty: the optimal base-stock level has the same critical-fractile form as the newsvendor solution but with the fractile taken over the disruption-length distribution (Tomlin 2006), and consolidation, which pools demand risk, now does nothing to expected cost and multiplies its variance, the risk-diversification effect of Schmitt, Sun, Snyder and Shen (2015). The chapter closes with the reliable fixed-charge location problem of Snyder and Daskin (2005). This mission formalizes the chapter's theorems on disruptions, with the risk-diversification theorem as its goal.
A supplier that is up is disrupted next period with probability ; one that is down
recovers with probability . The disruption chain records when the supplier is up
and in the -th consecutive period of a disruption; its stationary distribution is
and (disruptionPmf),
with distribution function (disruptionCdf).
A single location faces demand per period, pays per unit held and per unit
backordered per period, and follows a base-stock policy: it orders up to in every up
period and nothing in down periods. In the -th period of a disruption it has
units on hand or backordered, so its cost is
(periodCost), and the expected cost per period is
(meanCost), with variance over the disruption state (varCost). The critical fractile
and , the smallest with , determine the
optimal level.
In the reliable fixed-charge location problem (RFLP), sites fail independently with
probability and each customer is assigned to a chain of facilities: its level- facility
serves it when the closer facilities are disrupted, until it is assigned to an emergency
facility that never fails and charges the penalty . The objective (9.61) is fixed
cost plus expected transportation cost, with coefficients
(rflpPsi, rflpCost) under the constraints (9.62)-(9.67) (RFLPFeasible).
For identical locations and the centralized location formed by merging them (demand ):
that is, an optimal single-location level scales to the optimal centralized level ,
the centralized expected cost at is times the single-location cost, and its variance is
times the single-location variance. This is risk_diversification.
Lemma 9.2, the stationary distribution and distribution function of the disruption chain; Lemma 9.4, that the optimal base-stock level is a multiple of ; Theorem 9.5, , as the least minimizer of ; and Theorem 9.10, that in every optimal RFLP solution consecutive backup assignments are ordered by cost. Theorem 9.3, the optimality of base-stock policies, is cited by the book to Song and Zipkin without a model of the policy space and is not a target; Proposition 9.1 and the multisupplier results of Sect. 9.4 are left for a later mission, as discussed below.
Theorem 9.5 is the supply-side newsvendor formula: it says exactly how much inventory buys protection against disruptions of a given length, and it underlies the disruption models used in practice for raw-material buffers. Theorem 9.9 is the chapter's central insight and the reason supply and demand uncertainty call for opposite strategies: pooling reduces expected cost under demand uncertainty but only redistributes risk under supply uncertainty, concentrating it. Its three identities are what a risk-averse planner needs to compare the two designs by a mean-variance criterion. Theorem 9.10 is what lets the RFLP be formulated without ordering constraints and solved by Lagrangian relaxation like the UFLP.
None of these results has a machine-checked proof. The book proves Theorem 9.5 and the identities behind Theorem 9.9, sketches Lemma 9.4, and leaves Lemma 9.2 and Theorem 9.10 as exercises. The formal treatment of the piecewise-linear cost and its finite differences is reusable for the yield-uncertainty and multi-supplier models of the same chapter.
The cost is an infinite series whose terms grow linearly in against a geometric weight, so every statement about it begins with summability, and the finite-difference identity requires exchanging a difference with a sum. Theorem 9.5 then needs the convexity and piecewise linearity of to pass from a sign condition on slopes at multiples of to a global minimum over all real , and the identification of the least minimizer needs the slopes to be strictly negative below . The obvious idea, treating the problem as a discrete newsvendor over multiples of , is only half of the argument: it does not by itself exclude non-multiple minimizers, which is what Lemma 9.4 asserts.
Lemma 9.2 is elementary but the stationary equations involve a series over all down states, and the proof must establish summability before manipulating it. Theorem 9.9's scaling identities are termwise, but the optimality transfer in part 1 requires the scaling to preserve minimizers, which follows from the cost identity holding for every .
Theorem 9.10 is an exchange argument on a binary program with layered constraints. The delicate case is the emergency facility: swapping it into a lower level is infeasible, and the correct move is to promote it and drop the later assignment, which changes the constraints for every higher level; the argument must show feasibility of the modified solution level by level.
The disruption distribution is given by Lemma 9.2's formula rather than defined as the stationary distribution, and Lemma 9.2 shows it satisfies the stationary equations; the theorems assume and , under which every series is a geometric series times a polynomial and is summable. The quantity enters Theorem 9.5 as a natural number characterized by and for , which exists since ; the conclusion asserts both optimality and leastness of among all real levels.
Theorem 9.9 is stated as the scaling of the single-location functions; the decentralized totals
and are the mean and variance of a sum of independent copies, which the book
asserts rather than derives, and are not modeled separately. In the RFLP, levels are indexed by
Fin m, the emergency facility is a designated index whose data satisfy the book's
conventions through the hypotheses, and demands are positive with , both needed:
with or backup assignments are free and any order is optimal.
The EOQ with disruptions (Proposition 9.1) is a renewal-reward derivation without a formal model of the renewal process in the book, and the multisupplier newsvendor of Sect. 9.4 (Lemma 9.6, Theorems 9.7 and 9.8) rests on differentiability conditions the book defers to Dada et al. and on a lemma it leaves as an exercise; both are natural extensions on the same definitions rather than targets here.
Choosing where to open distribution centers is the strategic decision that fixes a supply chain's shape for years. The basic model, the uncapacitated fixed-charge location problem (UFLP) of Balinski (1965), trades the fixed cost of opening sites against the cost of transporting demand from open sites to customers. It is NP-hard, yet routinely solved to optimality, and the reason is a fact about its relaxations: the LP relaxation is unusually tight, and Lagrangian relaxation, which Cornuejols, Fisher and Nemhauser (1977) brought to location problems, gives the same bound with a subproblem solvable by inspection. Chapter 8 of Snyder and Shen's Fundamentals of Supply Chain Theory (2019) develops the UFLP, its Lagrangian relaxation and Erlenkotter's (1978) DUALOC dual-ascent method, then the p-median problem with Hakimi's (1965) node-optimality theorem and the covering models. This mission formalizes the chapter's numbered results, with the equality of the Lagrangian and LP bounds as its goal.
Customers have demands ; candidate sites have fixed costs ; shipping
one unit from to costs . A solution opens sites () and assigns
demand fractions (, , ); its cost is
(uflpCost, UFLPFeasible). The optimal value is
(uflpOpt); relaxing to gives the LP relaxation with
value (uflpLP).
Lagrangian relaxation removes the assignment constraints and charges per unit of
violation. For fixed multipliers the subproblem (UFLP-LR) minimizes
over ,
binary, (lagrObjective, LagrFeasible), with value (zLR). It
separates by site: the benefit of opening is
(benefit), and is opened iff . The best bound is the Lagrangian dual
value (zLRbest).
DUALOC works with the condensed dual of the LP relaxation, whose variables satisfy
with . A dual solution and
a site set form a primal-dual pair (PDP) when these constraints are tight on and
every customer has a site in with ; the primal solution opens and
assigns each customer to its nearest open site (NearestIn, primalX, primalY).
For the p-median problem on a network, the customers are the nodes, is the shortest-path
distance between nodes, and a facility may sit at position along an edge of length
, at distance from node (NetPoint, netDist).
The p-center problem minimizes the largest distance from a customer to its nearest of
open sites; pCenterValue c p is its optimal value.
For every instance with at least one candidate site,
This is lagrangian_equals_lp.
Theorem 8.1, the closed form with its optimal solution; the bounds (8.16) and (8.19) ; Theorem 8.3, the variable-fixing tests; Lemma 8.4, the DUALOC duality gap ; Lemma 8.6, the characterization of complementary slackness violations; Theorem 8.7, Hakimi's theorem that some nodes are optimal among all -point sets; Lemma 8.8, the equivalence between the -center value being at most and a set cover of radius with at most sites. Proposition 8.5, which concerns the output of a specific procedure, is not a target.
Corollary 8.2 explains the behavior of every Lagrangian location code: the bound cannot beat the LP bound, so its value lies in the ease of the subproblem and in extensions to nonlinear location models (the location model with risk pooling of Chapter 12) where no LP is available. Theorem 8.1 is the subproblem solution those codes use; Theorem 8.3 is the variable-fixing device of Daskin, Snyder and others that shrinks branch-and-bound trees. Lemmas 8.4 and 8.6 are the analytical core of DUALOC, the method that made large UFLP instances solvable in the 1970s. Hakimi's theorem is the reason the -median problem is a discrete problem at all, and Lemma 8.8 is the reason -center problems are solved by bisection over covering problems rather than by their weak MIP formulation.
None of these results has a machine-checked proof. The book proves Theorems 8.1 and 8.3 and Lemma 8.6, cites Corollary 8.2 to Appendix D and Theorem 8.7 to Hakimi, and leaves Lemmas 8.4 and 8.8 as exercises. The formal treatment of the integrality property and of Lagrangian duality for a linear objective over a product of boxes is reusable for the p-median and capacitated variants the chapter goes on to discuss.
The goal is an LP duality statement in disguise, and the obvious idea, that is the dual function of the LP relaxation, is exactly what needs proof. Two facts must be established: that for fixed the subproblem over binary has the same value as over , because after the optimal is substituted the objective is linear in ; and that the supremum over of the resulting concave piecewise-linear function equals the LP minimum. The second is strong duality for a linear program, which Mathlib does not provide ready-made; it has to be obtained either through a Farkas-type argument or by exhibiting, for the LP optimum, a multiplier vector that attains it, which for this problem can be read off the LP dual. The book proves none of this; it invokes Lemma D.3.
The bounds (8.16) and (8.19) are easier but not free: the infima and suprema defining , and must be shown attained, which needs finiteness of the binary choices and compactness of the assignment polytope. Theorem 8.3 depends on the value of the subproblem with one variable forced, which is Theorem 8.1 applied to a modified instance. Hakimi's theorem needs a concavity argument in each point's position and a bookkeeping step, since moving several points to nodes may merge them and the result must still have exactly nodes. Lemma 8.8 is combinatorial and short once the -center value is identified with a minimum over -subsets.
Customers and sites are Fin n and Fin m; demands, costs and fixed costs are arbitrary
reals, as the book's formulations are, and the theorems that need it assume . Optimal
values are infima or suprema of the sets of attainable objective values, all of which are
nonempty and bounded under the stated hypotheses. The Lagrangian dual value is a supremum over
all real multiplier vectors, so Corollary 8.2 asserts in particular that the supremum equals the
attained LP value.
The DUALOC statements take the nearest-facility assignment as a function a with the
defining property, so ties are resolved by the hypothesis, and the complementary slackness
violation is written exactly as (8.51) with substituted. Hakimi's theorem is stated
for a family of points with repetition allowed, which is stronger than for a set. It assumes
what the book's network supplies: the node distances satisfy the triangle inequality, since they
are shortest-path distances, and every edge carrying a point is at least as long as the distance
between its endpoints. The argument needs both: they make the ends of an edge coincide with its
nodes, and without them a point on a short fictitious edge can beat every node. The
set covering value in Lemma 8.8 is expressed through the existence of a cover with at most
sites rather than as a natural-number infimum, whose value on infeasible instances would
falsify the equivalence.
The definition module is shared by all ten items. The Lagrangian relaxation of the -median problem (Sect. 8.3.2.2), the continuous knapsack subproblem of the capacitated problem, and Proposition 8.5 on the dual-ascent procedure are natural extensions on the same definitions.
A retailer facing uncertain single-period demand and buying from a supplier at a wholesale price orders less than the two of them together would want. The reason is not irrationality but incentives: the retailer bears the whole cost of unsold stock while the supplier collects a margin on every unit ordered, so each party marks up its own cost and the combined markup, Spengler's (1950) double marginalization, depresses the order. Pasternack (1985) showed that a buyback credit for unsold units, priced correctly, realigns the retailer with the chain, and Cachon (2003) surveys the contracts that followed. Chapter 14 of Snyder and Shen's Fundamentals of Supply Chain Theory (2019) develops this as a Stackelberg game on the newsvendor model: the supplier sets contract terms, the retailer sets the order quantity. This mission formalizes the chapter's seven theorems, with the buyback allocation theorem, which shows that buyback both coordinates the chain and can divide its profit in any proportion, as the goal.
Demand is a random variable with law on and mean . The retail price is
; the supplier's and retailer's per-unit costs are and with ;
lost sales cost the two parties goodwill penalties and with ; unsold
units salvage for (ContractData). With the expected
sales (expSales) and the expected leftover (expLeftover), a transfer
payment from retailer to supplier determines the two profits (retailerProfit,
supplierProfit),
whose sum (chainProfit) is independent of the
contract. The chain-optimal quantity maximizes ; the retailer's and supplier's
optimal quantities , maximize their own profits. A contract coordinates the
chain when , stated here as equality of the sets of maximizers, and is a
coordinating contract type when some choice of its parameters does this with both profits
positive.
The four contracts are their transfer payments. Wholesale price: . Buyback:
, the supplier crediting per unsold unit, with
and of (14.22). Revenue sharing: the retailer keeps a fraction of sales
and salvage revenue, , with of
(14.34). Quantity flexibility: the supplier reimburses the retailer's loss on
unsold units up to , , with
of (14.46). For buyback, is the
retailer's share of the chain profit (buybackShare), and (buybackB1,
buybackB2) are the credits at which one party earns everything.
Under buyback with at the chain-optimal , the retailer's profit is decreasing and the supplier's increasing in , with and
the supplier losing money for , both earning positive profit for , and
the retailer losing money for . This is buyback_allocation.
Equation (14.8), maximizes iff ; Theorem 14.1, the wholesale price contract coordinates iff , at which the supplier's profit is negative; Theorem 14.2, whenever ; Theorem 14.3, for IGFR demand the supplier's induced profit is unimodal; the identities (14.27) and (14.28), and its complement under buyback; Theorem 14.4, buyback with coordinates; Theorem 14.6, revenue sharing with coordinates; Theorem 14.7, quantity flexibility with makes optimal for the retailer.
The chapter's theorems are the analytical basis of contract design in newsvendor supply chains. Theorem 14.1 and 14.2 make double marginalization precise: coordination by price alone is possible only at a price the supplier rejects, and any acceptable price makes the retailer under-order. Theorem 14.3 is what lets the supplier optimize the wholesale price at all, and is the reason the IGFR class of Lariviere and Porteus (2001) is standard in the field. Theorems 14.4 to 14.6 show that buyback and revenue sharing coordinate and, through the share , that the chain profit can be split arbitrarily, so a coordinating contract can be made acceptable to both parties. Theorem 14.7 shows the limits: quantity flexibility coordinates the retailer but not necessarily the supplier.
None of these results has a machine-checked proof. The book proves Theorems 14.1, 14.3, 14.4, 14.5, 14.6 and 14.7 and leaves 14.2 as an exercise. The formalization of the fractile characterization and of the affine profit identities is reusable for the many contract types (sales rebates, quantity discounts) the chapter cites but does not analyze.
The wholesale price results are first-order conditions on concave functions, and the difficulty is entirely in the analysis: has derivative for every when is continuous, which is a differentiation under the integral that has to be carried out for a Lipschitz integrand and an arbitrary law, and the maximizers of the concave profit must then be identified with the solutions of the fractile equation, including existence by the intermediate value theorem. Theorem 14.1's "only if" direction requires that a coincidence of maximizer sets pins the fractile and hence the price, which is where strict monotonicity of enters.
Theorem 14.3 is the delicate one. The obvious approach, concavity of , fails:
the book stresses the function is not concave in general. The proof is a sign-change argument
on the derivative (14.16), which is times a bracket that IGFR makes decreasing,
minus a constant; one must show the derivative is positive near , eventually negative, and
crosses zero exactly once, and that the last needs both strictly decreasing and the
bracket positive at the crossing. Working with a density in Mathlib means relating the
withDensity law to the distribution function throughout.
The buyback, revenue sharing and quantity flexibility theorems are algebra once the right identity is found: is an affine function of with slope . The formalization must handle the boundary cases the book glosses over, where or and one party is indifferent among all quantities, so "the same maximizes" holds only in one direction. Theorem 14.5 also needs , a Jensen-type bound , for , and the ordering needs , which the book's proof does not isolate.
Demand is an arbitrary probability law on with finite mean, not assumed
nonnegative; the book's own examples use normal demand. Optimal quantities are maximizers over
all of (IsMaxOn … Set.univ), and coordination is the coincidence of maximizer
sets, with the degenerate endpoints of the parameter ranges stated as one-directional. Where the
book uses a density, continuity of the distribution function (NullSingletonClass) is assumed
instead, except in Theorem 14.3, where the density is explicit, continuous on
(so the exponential law is included), positive on , and IGFR on .
Strict monotonicity of on is never assumed, since it would rule out every
nonnegative demand law. Theorem 14.1 assumes and a positive chain
optimum, , both implicit in the book's proof.
The transfer payments are functions of and the profits are defined for every , so the
theorems compare values of one family of functions; the book's is with
Mathlib's cdf, and the quantity flexibility integral is an interval integral. Theorem 14.5
takes , and as hypotheses. Without positive goodwill costs
its strict inequalities and fail (at , ; at
, ).
The definition module is shared by all eleven items. The equivalence (14.42)-(14.43) of revenue sharing and buyback, the supplier's stationarity under quantity flexibility (Problem 14.11) and the allocation results for revenue sharing (14.40)-(14.41) are natural extensions on the same definitions.
Chapter 9 of Axsäter's Inventory Control opens with the observation that in a multi-echelon system it is not optimal to choose batch quantities installation by installation: the batch at one site is the demand pattern of the next site upstream. Even with constant customer demand the exact optimum can be complicated, and the book's Example 9.4 shows a four-stage serial system whose optimal batch at one stage alternates between two values over time. Roundy (1985, 1986) showed that this complexity can be avoided at a guaranteed price: restrict every batch quantity to be a power of two times a common basic quantity, nested from stage to stage, and the best such policy costs at most 2 % more than the optimum. The book presents the result for a serial system and remarks that the same approach handles assembly and distribution systems and, in Sect. 7.3.1.2, joint replenishments. It is the capstone of Chapter 9 and the multi-echelon payoff of the powers-of-two analysis of Chapter 7.
A serial system has installations; installation produces item from one unit of item , item is obtained from an outside supplier, and item 1 faces a constant, continuous final demand . Lead-times are zero, shortages are not allowed, production is instantaneous, and each batch quantity is constant over time. Installation has an ordering cost per batch and an echelon holding cost per unit and time unit, charged on the echelon stock (the stock at installation and everything downstream), so that the cost per time unit is the sum of single-item costs of the chapter 4 form,
The book first works out the two-level case, where the optimum has for a positive integer , the cost (9.9) is the EOQ cost with modified parameters and , and the best is found from by a rounding rule.
For stages, Roundy's constraints (9.16) require with nonnegative integers , so that every solution is nested. The relaxed constraints (9.18) require only ; they are implied by (9.16), the relaxed problem is convex with linear constraints, and its solution is computed by aggregating consecutive stages whose cost ratios decrease. Roundy's solution rounds to for a basic quantity , chosen as in Proposition 7.2.
For , , and any minimizer of over positive batch quantities satisfying (9.18), there exist and integers such that satisfies (9.16) and
The two-level results of Sect. 9.2.1: the equivalence of the installation and echelon cost forms (9.6) and (9.9); the optimal and cost (9.10)-(9.11) for a given ; the closed form and convexity of (9.12) and the real minimizer (9.13); and the integer rounding rule with its corollary that forces . For stages: existence of the relaxed optimum; the aggregation lemma, that forces ; that (9.16) implies (9.18), so the relaxed optimum bounds every powers-of-two policy from below; and that rounding to the nearest power of two times is monotone and lands within a factor .
The result itself. Roundy's theorem replaces an intractable lot-sizing problem by a closed-form computation with a provable 2 % guarantee, and the policies it produces are the nested, periodic schedules that production planning wants anyway. In Example 9.4 the rounding with is already within 0.7 % of the relaxed bound. The two-level analysis has its own use: the modified parameters and are what the Blackburn-Millen heuristic of Sect. 9.3.2 feeds to the Wagner-Whitin algorithm under time-varying demand, and the condition tells when a two-stage system collapses to one stage.
Formalizing it. The book proves the bound in a paragraph that leans on three earlier results:
Proposition 7.2 for the rounding, the Lagrangean relaxation for the lower bound, and the
aggregation algorithm for the structure of . Formalizing it makes explicit
what the paragraph glosses: that the multipliers vanish off tight constraints, that equal
quantities round to equal quantities, and that the comparison class is the class of nested
constant-batch policies. The published Proposition 7.2 (pot_two_percent) and the mean value
are cited as reference items. Nothing here is open; no statement has a
machine-checked proof yet.
The obvious argument, rounding item by item and invoking Proposition 7.2 for each, does not work: Proposition 7.2 bounds the rounded cost relative to each item's unconstrained optimum, and is generally not that optimum. The proof must pass through the Lagrangean relaxation (9.19)-(9.20), under which is the unconstrained optimum for the modified holding costs , apply Proposition 7.2 there, and transfer the bound back using complementary slackness: the multiplier is positive only where , and equal quantities round to equal quantities, so the correction terms vanish for both and its rounding. A solver therefore needs the KKT conditions for this convex program, or an equivalent direct argument through the aggregation structure (within an aggregate all quantities are equal and their sum is an EOQ problem). The two-level statements and the rounding lemma are elementary.
Stages are indexed by Fin N; the cost is serialCost A e d Q = ∑ i, eoqCost (A i) d (e i) (Q i)
with eoqCost imported from the chapter 4 definitions, and the two-level cost is
eoqCost (A1 + A2/k) d (e1 + k e2) Q1 by definition, so the published eoq_optimal and
eoq_cost_at_eoq apply to it directly. SerialNested and SerialPowerOfTwo are the constraints
(9.18) and (9.16) on consecutive indices; for both hold vacuously and the goal is
Proposition 7.2 for one item. potRound q Q is .
All theorems assume and positive batch quantities. The book says "nonnegative ordering and echelon holding costs"; strict positivity is needed for the relaxed problem to have a solution at all ( or sends the optimal to or ), and it is what the book's own examples satisfy. The relaxed optimum enters the goal as a hypothesis, with existence stated separately; uniqueness is not used. The constant is the exact .
Two readings are excluded. The bound is not against an arbitrary nested , for which it would be false, but against the relaxed optimum; and the comparison class is stated as nested constant-batch policies, the class the relaxed problem bounds directly, not the book's larger class of time-varying policies, which would need a dynamic model. The definitions are reusable for assembly systems and for the joint replenishment problem of Sect. 7.3.1.2, whose Roundy analysis is the same with a fictive item 0 of zero holding cost; contributions in that direction are welcome.
Every model in Chapters 5 and 6 of Axsäter's Inventory Control assumes at the outset that the ordering policy is of or type. Section 6.2 asks whether better policies exist and answers, for the case in which there are no ordering costs but every order must be a multiple of a fixed batch quantity , that none do: Proposition 6.1, "an policy is optimal", with a proof the book attributes to Chen (2000). For it is the optimality of an order-up-to- policy in the absence of ordering costs, and for continuous or Poisson demand it transfers to policies, which are then the same thing. The proposition is the one place in the book where a policy is compared against every feasible alternative rather than against other members of its own family, and its proof is short enough to be given in full, which makes it the natural capstone of Chapter 6.
Demand is compound Poisson: customers arrive according to a Poisson process with rate and each demands an integral number of units, with sizes independent and identically distributed with law on the positive integers, independent of the arrival process. The book's standing assumption that not all demands are multiples of some integer larger than one is kept. Writing for the -th arrival time, for the number of arrivals by time and , the total demand by time is .
The replenishment lead-time is constant and , the demand over a lead-time, has law . A holding cost and a shortage cost per unit and time unit are charged. There are no ordering costs, but all orders must be multiples of a given batch quantity and can only be triggered by customer demands. A policy is therefore any rule that decides at each demand epoch how many batches to order; with initial position the inventory position evolves as (Eq. 6.22), and .
The standard argument of Sect. 5.3.2 gives the cost rate at time as
the expected holding-plus-shortage cost rate of an inventory level (Eq. 6.20), which is convex in with as . The band cost is and denotes an integer minimizing . The policy orders, as soon as the position is at or below , the smallest number of batches that brings it above ; its position lives in the band from the first order on. The performance measure is the long-run average cost rate as .
Almost surely, (1) for every policy and every ,
and (2) the policy attains it:
Lemma 6.1, that is convex and minimized at the representative of in the band; the closed form of the position, until the first order and the band representative of afterwards; the uniform occupation of the band by the reduced process , the book's "the steady state distribution can be shown to be uniform", and its consequence that the long-run average of is ; Proposition 5.1 in the same ergodic form for the policy; and the two halves of the goal as separate statements.
The result itself. Proposition 6.1 is what licenses the two-parameter policies on which the rest of the book's single-echelon theory is built, and it does so for the practically important case of batch ordering (pallets, containers, production lots). Its proof also explains why the policy works: the only quantity a policy controls is the residue class of the inventory position modulo , which no policy can influence, and the position within that class, which the policy always sets to the cheapest possible value. The book extends the same reasoning to other cost structures and to periodic review.
Formalizing it. The proposition is a theorem about the class of all policies, and the book's proof is pathwise: Lemma 6.1 compares any policy with the reduced process instant by instant, and an ergodic statement about the reduced process does the rest. Formalizing it therefore forces the policy class, the demand process and the long-run average to be written down exactly, which the book never does. Nothing here is open; no statement has a machine-checked proof yet.
The pointwise comparison is elementary once Lemma 6.1 is available, and Lemma 6.1 is discrete convexity. The difficulty is entirely in the ergodic statement: that the reduced position (the band representative of ) spends a fraction of the time at each point of the band, almost surely. In discrete time this is the convergence of occupation frequencies for an irreducible random walk on with step law modulo , where irreducibility is exactly the aperiodicity assumption on , and the book's double-stochasticity argument (Eq. 5.33-5.34) identifies the uniform law as stationary. Passing to continuous time adds the exponential holding times: the time average is the arrival average weighted by i.i.d. holding times independent of the walk, which a strong law of large numbers turns back into the discrete statement. Mathlib has the strong law and the exponential law but no ergodic theorem for finite Markov chains, so that is the groundwork a solver must build. The naive route through the stationary distribution of the embedded chain at demand epochs is not enough on its own: for pure Poisson demand that chain is periodic, and the book itself notes it.
A CompoundPoissonDemand on a probability space packages the rate, the size law with
and the aperiodicity condition, and two sequences of random variables, the gaps and the sizes,
with their laws (expMeasure lam, the given pmf), independence within each sequence, and
independence between the sequences. Arrival times are partial sums of the gaps, count t is
the supremum of , and cumDemand n is the partial sum of the sizes. A policy
is a function m : ℕ → Ω → ℕ with no measurability requirement; ipPath and ipAt are the
position after the -th demand and at time ; rqIP and rqOrders are the policy,
with a theorem identifying rqOrders as a policy in the general sense. avgCost is
with for . The cost function
is sPolicyCost from the previous mission, applied to a DiscreteDemand that the hypotheses
tie to the process as the law of the demand in .
Conventions: , , ; is any minimizer of (it exists by
the divergence of , proved in the previous mission); is arbitrary. count and the
interval integral take junk values on the null set where arrivals do not tend to infinity, which
the almost-sure conclusions absorb. "" is stated as "for every ,
eventually ", avoiding a liminf on that could be junk.
Two readings that would trivialize the goal are excluded: the lower bound is over every rule, not over stationary or measurable ones, and the achievability half is a genuine limit, not a bound. The demand model and the occupation-frequency theorems are reusable for Proposition 10.1 and the batch-ordering models of Sect. 10.5; contributions establishing the ergodic theorem for irreducible chains on a finite cyclic group are welcome and would close most of this mission.
An policy has two parameters. Chapter 4 of Axsäter's Inventory Control chooses the batch quantity from a deterministic model, and Chapter 5 chooses the reorder point from a stochastic one with held fixed; the book presents this two-step practice as an adequate approximation. Section 6.1 asks what is lost by it and shows how to optimize both parameters jointly in one stochastic model. For discrete demand the answer, due to Federgruen and Zheng (1992), is an algorithm of remarkable simplicity: increase one unit at a time, keep the reorder point optimal along the way by a one-line rule, and stop at the first for which the cost goes up. The claim that this stopping rule finds the global optimum over all pairs is the capstone of Sect. 6.1.1.1 and the goal of this mission. The same idea is reused in the book for policies (Sect. 6.1.1.2) and, in continuous form, for normally distributed demand (Sect. 6.1.2).
An item is controlled by a continuous review policy with integral reorder point and batch quantity . Demand is discrete and stationary; the lead-time demand takes values in the nonnegative integers with probabilities and has a finite mean . The average demand per unit of time is . Costs are a holding cost and a shortage cost , both per unit and time unit, and an ordering cost per batch.
The building block is the cost of an policy that keeps the inventory position at a fixed integer . By the standard argument of Sect. 5.3.2 the inventory level a lead-time later is , and the average holding and shortage cost rate is (Eq. 6.3)
Under the policy the inventory position is uniformly distributed on (Proposition 5.1), so the total average cost rate is (Eq. 6.4)
and (Eq. 6.5), attained at an optimal reorder point . The ready rate links the cost to service: raising the reorder point by one unit changes the cost by (Eq. 5.60).
Let be the smallest batch quantity with and let be an optimal reorder point for . Then
The increment identity ; convexity of on together with as ; Eq. (5.60) and the convexity of in ; Eq. (5.61), that the largest with is optimal for its ; existence of an optimal for every ; the recursion (6.6)-(6.7), chosen by comparing with , and ; the equivalence and the monotonicity of that minimum in ; and the existence of some at which the costs stop decreasing.
The result itself. Joint optimization typically enlarges the batch and lowers the reorder point relative to the two-step procedure, and the book's Example 6.1 puts the resulting cost saving at a few percent. The Federgruen-Zheng procedure makes the exact joint optimum for discrete demand as cheap to compute as the two-step approximation, because each step of the recursion evaluates at two points. It is the exact benchmark against which the book's approximate techniques for normal demand (Sects. 6.1.2 and 6.1.3) are judged, and its structural core, that the optimal window of consecutive inventory positions grows one neighbour at a time, is the discrete-convexity fact behind the whole of Sect. 6.1.
Formalizing it. The mathematics is settled. What the mission produces is a Lean development in which the steps the book marks "evident" and "obvious" are separate statements: that the recursion preserves optimality, that the marginal cost of enlarging the batch is monotone, and that a minimum over exists at all. None of the statements has a machine-checked proof yet.
The obvious first idea, to argue that is convex in and stop at its first increase, does not work as stated: is a minimum over of a ratio and is not convex in general. What is true, and what the proof uses, is that the marginal cost is nondecreasing, because it equals the smaller of the two -values adjacent to the optimal window. Establishing the recursion (6.6) is where discrete convexity is needed: one must show that the best window of consecutive positions is obtained from the best window of by adding a neighbour, which fails for non-convex . The window sums are themselves convex in with increments , and the argument compares a competing window with the optimal one through the end terms. The remaining steps are finite algebra and an induction on from .
A DiscreteDemand is a function with ,
and a summable first moment; is a finite sum over , empty for
. The reorder point ranges over and the batch quantity over ,
with assumed in every statement because Lean's division by is . Convexity of a
function on is stated as nondecreasing increments, and divergence as
Tendsto g (cocompact ℤ) atTop.
enters the goal as a function CQ together with the hypothesis that CQ Q is the least
value of ; the stopping index is characterized by the two conditions
that define "the smallest with ", and by optimality at .
That these objects exist is the content of two separate items, so the goal is not vacuous: a
minimum over exists for every because , and the costs cannot decrease
forever because the average of the smallest values of tends to infinity. The
positivity of and is the book's setting and is assumed where it appears.
The uniform inventory position of Proposition 5.1, which needs the book's assumption that not all demands are multiples of an integer larger than one, is taken as given in the cost formula (6.4), as the book does; the proposition itself is the subject of the next mission of this series. The definitions are reusable for the optimization of Sect. 6.1.1.2 and for the multi-echelon batch-ordering results of Sect. 10.5; contributions formalizing the Zheng-Federgruen algorithm on top of them are welcome.
Chapter 4 of Axsäter's Inventory Control fixes the batch quantity from a deterministic model, and Chapter 5 asks the question that deterministic models cannot answer: with demand random and a replenishment lead-time , when should the next batch be ordered? Under a continuous review policy the answer is a single number, the reorder point , and Sections 5.3 through 5.9 are one sustained computation of what a given buys. The chapter's own capstone is Eq. (5.67): if backorders are charged at per unit and time unit and stock at , the cost-minimizing reorder point is exactly the one whose fill rate is . The book calls the relationship "even more striking" than its discrete counterpart, and it is used in practice in both directions: a backorder cost prescribes a service level, and a chosen service level reveals the backorder cost a planner is implicitly assuming (Eq. 5.68).
An order for a fixed batch quantity is triggered whenever the inventory position (stock on hand plus outstanding orders minus backorders) falls to the reorder point , and it arrives time units later. Demand is continuous and normally distributed; the demand over a lead-time has mean and standard deviation . Two modelling facts from the book are taken as the definition of the steady state. First (Sect. 5.3.1), the inventory position is uniformly distributed on ; the book proves this for compound Poisson demand (Proposition 5.1) and adopts it as an accurate approximation for continuous demand. Second (Sect. 5.3.2, Eq. 5.35), the inventory level a lead-time later is the inventory position now minus the demand in between,
with the two terms independent. The law of is therefore the image of the product of a uniform and a normal law under subtraction, and everything in the chapter is a functional of it:
The closed forms run through the standard normal loss function of Eq. (5.40), published with the newsboy mission and reused here, and through its integral, the second loss function of Eq. (5.64). Both are tabulated in the book's Appendix 2, and both recur in Chapters 6, 9 and 10.
For and any , a reorder point minimizes over if and only if
The biconditional carries both halves of the book's sentence: the stationary point is the optimum ("the optimal is obtained for ") and the optimum is stationary ("in the optimal solution we have ").
In the order the chapter builds them: Eq. (5.41), , with decreasing and convex; Eq. (5.39), the distribution function by conditioning on the inventory position; Eq. (5.42), its closed form ; Eq. (5.43), the density; Eq. (5.52), the fill rate ; Eq. (5.55), the expected backorders covered by one batch, and Eq. (5.54), that is the same fill rate; integrability of the cost rate and the mean ; Eq. (5.63), ; Eq. (5.64), the closed form of and ; Eq. (5.65), the cost ; Eq. (5.66), ; and the convexity of in .
The result itself. The reorder point is the one parameter of an policy that stochastic demand actually decides, and Eq. (5.67) says that deciding it by cost and deciding it by service level are the same decision, with an explicit dictionary between the two. That is why the book can present service-level constraints (Sect. 5.7) and shortage costs (Sect. 5.9) as interchangeable ways of specifying the same thing, and why it warns, immediately after Eq. (5.67), that the equivalence is only valid when is given: with an ordering cost and a joint optimization of and it fails, which is Chapter 6's problem.
The intermediate formulas have independent standing. Eq. (5.42) is the single expression from which every service measure of the chapter is computed, and Eq. (5.65) is the cost function that Chapter 6 extends by an ordering cost, Eq. (6.10), and optimizes iteratively. The second loss function returns in the periodic-review fill rate of Eq. (5.86) and in the two-echelon batch-ordering model of Sect. 10.5.
Formalizing it. Nothing here is open; the value is that the steady-state model becomes an explicit measure, so that formulas the book obtains by manipulating integrals whose existence it never questions become theorems about that measure. Integrability of the cost rate is a target of its own for exactly that reason. None of the statements has a machine-checked proof yet.
The obvious route is the book's, and it is not the hard part: once is known in closed form, every later identity is calculus on and . The work is upstream of that. The distribution function (5.39) is a conditioning argument over the product measure, a Fubini step in which the inner probability is a Gaussian tail; the density (5.43) is a derivative of a parameter-dependent integral; and Eq. (5.63) exchanges the order of two integrals over an unbounded region, which needs integrability of itself. The derivative (5.66) is then obtained from the closed form (5.65), not by differentiating under an expectation, which is what makes the goal reachable: is a smooth function of with an explicitly increasing derivative, and Eq. (5.67) follows from strict convexity together with the fact that the fill rate is a continuous, strictly increasing function of ranging over . A solver who starts from the expectation and tries to differentiate it directly will meet the kink of at and a dominated-convergence argument; the integrated route avoids both.
The inventory position is rqPosition R Q, Lebesgue measure conditioned on ; the
lead-time demand is newsboyDemand m s, Mathlib's gaussianReal m (s^2).toNNReal from the
newsboy mission, with and ; and rqLevel R Q m s is the pushforward of
their product under . Two lemmas in the definition file record that both
are probability measures when . The distribution function, ready rate and cost are the
measure of , the measure of , and a Bochner integral against this law.
Every statement assumes and . At the conditioned measure is the zero measure and every integral is , so a statement without the hypothesis would be true and empty; at the demand is a point mass and has jumps. The mean is unrestricted, as none of the formulas depends on its sign, and reorder points may be negative, which the book explicitly allows in Sect. 5.8. Costs and are positive in every statement that involves them, and is stated for because its closed form holds there.
Two trivializing readings are ruled out. The cost is not defined as a formula in but as an expectation, so the closed form (5.65) has content; and because Lean's Bochner integral of a non-integrable function is , integrability of the cost rate is stated as a theorem rather than assumed, otherwise a zero cost would make every reorder point optimal. The goal quantifies optimality over all real competing reorder points, not a neighbourhood.
A complete development needs Gaussian tail integrals, differentiation of parameter-dependent integrals, Fubini on a product of a bounded interval with the line, and the strict monotonicity of the Gaussian distribution function. The loss functions and and the inventory-level law are reusable across the rest of the series; contributions that establish the same identities for an arbitrary continuous lead-time demand with a finite mean, where Eqs. (5.39), (5.63) and (5.66) hold verbatim, are welcome.
Formalization status, 10 October 2026: all eight milestones are complete. Both corrected addenda are Proved. The subgroup upper Hilbert estimate has an accepted proof with eight Proved inputs and no Open dependencies (submission 07198c4c-ea10-4c9a-89b8-b4a5319c40b5). Together with the accepted sharp contact-rank and ambient lower Hilbert bounds, this discharges every dependency of the corrected converse. The corrected 26-part full-paper goal remains Open, with exactly two frontier leaves: quadratic parameter-stalk germs and the local complex model for the closure action. The common-shift forward conclusion, quadratic first-input addition hypothesis in the converse, all constants and degree thresholds, and all 24 unaffected goal clauses are unchanged. The separately published Chardin theorem in arbitrary characteristic remains Open; it is no longer required by this mission proof.
An auxiliary polynomial in a transcendence proof is constructed to vanish to high order at many points. A multiplicity estimate limits how often that can happen without a geometric reason: a positive-dimensional algebraic subgroup can make the apparent vanishing conditions dependent. Such estimates are used to turn analytic approximations into algebraic-independence conclusions. The theory applies to products of commutative algebraic groups with both archimedean and nonarchimedean analytic directions.
The mission's objective is a faithful formalization of the complete 1986 paper, Patrice Philippon's Lemmes de zéros dans les groupes algébriques commutatifs, including every theorem, corollary, proposition, lemma, definition, and mathematical supporting claim. All results remain required whether or not a current application uses them. Necessary source corrections are explicit and their counterexamples remain part of the completion goal. The author's 1987 corrections are applied transparently, and the addendum's additional results are recorded separately (original, errata and addenda).
Let K be the complex field or the completed algebraic closure of the field of ℓ-adic numbers for a prime ℓ, as in the source. Let G be a product of finitely many commutative algebraic groups G₁,…,Gᵣ over K, with each factor embedded as a quasi-projective variety in a projective space. Write n for the sum of their dimensions. A point of the product embedding has one block of homogeneous coordinates for each factor.
A nonzero multihomogeneous polynomial P has a degree Dᵢ in its i-th block of coordinates. Its zeros define a hypersurface of the ambient product of projective spaces. An analytic subgroup A is locally parametrized by an analytic homomorphism from a finite-dimensional additive K-space. The order of P along A at a group point g is the vanishing order of P after composing local projective coordinates with the translated parametrization. This definition must be independent of the choice of nonzero local coordinate representatives.
Let Σ be a finite set of group points containing the identity. Its n-fold sumset Σ(n) consists of all sums of n, not necessarily distinct, elements of Σ; Σ(0) is the singleton identity. For a connected algebraic subgroup H, the integer s is the analytic codimension of A∩H in A. The expression |(Σ+H)/H| counts distinct H-cosets meeting Σ. Codimension concerns the analytic tangent dimension; it does not assert that the point-set intersection is finite.
The source's Hilbert degree form ℋ(V;D₁,…,Dᵣ) is (dim V)! times the highest homogeneous part of the multigraded Hilbert–Samuel polynomial of the projective closure of V, evaluated at the degrees. It must be constructed from the coordinate ring, rather than supplied as an arbitrary numerical function. For one projective factor it is deg(V)D^(dim V). These definitions are fixed in §§2–3 of the original paper.
The completion target is all results of the paper, with an aggregate goal that requires their individual formal statements. Theorem 2.1 is one milestone within that target. For each fixed family of embedded group factors, it chooses positive integers cᵢ, each depending only on its corresponding embedding. These constants precede the analytic subgroup, the finite sampling set, the polynomial degrees, the polynomial, and the contact parameter T. If P has order at least nT+1 along A at every point of Σ(n), there is a connected algebraic subgroup H with
The same H is contained in a translate of the zero locus of P on G and is incompletely defined by equations of multidegrees at most (c₁D₁,…,cᵣDᵣ): it is an irreducible component of the common zero locus in G of equations with those degree bounds. Both geometric conclusions are part of the target. A formalization retaining only the displayed numerical inequality would omit part of the original result.
The 1986 paper has 13 numbered results. Section 2 contains Theorem 2.1 and Corollaries 2.2–2.3, including the one-dimensional analytic result and the result for disjoint group factors. Section 3 contains Lemmas 3.1–3.2, Proposition 3.3 and Lemma 3.4. Section 4 contains Propositions 4.3–4.4, Lemmas 4.5–4.6 and Proposition 4.7. Section 5 contains Lemma 5.1. Every clause of these statements belongs to the mission.
Definitions 3.5, 4.1 and 4.2, the unnumbered setup, the internal Facts A–E, the counterexample after Proposition 3.3, and the mathematical claims in remarks also require coverage. Source numbers and page references identify the correspondence between the prose and Lean declarations. The two 1987 addendum targets are explicitly corrected variants. The forward statement assumes ambient connectedness and allows one common shift , so for every sampled , with the same connected subgroup satisfying incomplete definition and the full coset-count inequality. Its embedding-dependent constants are chosen before the sample and polynomial. This is not an exact transcription of the printed strengthening. The converse assumes ambient connectedness and explicit local addition formulas of degree at most two in the first input for every embedded factor. Its subgroup may remain disconnected; the degree threshold, constant , and polynomial nonzero somewhere on are unchanged. All original sampling and contact assumptions remain. No desired estimate is encoded as a geometric assumption.
The completion goal is PhilipponMultiplicity.paper_results_corrected_erratum, a conjunction of 26 concrete propositions. This is a collection goal requiring the source statements and the additional mathematical claims in the coverage record. It is separate from the paper's original Theorem 2.1, which remains individually named and reusable.
The mission has eight milestones; seven are currently Proved:
| Source | Milestone | Status |
|---|---|---|
| Theorem 2.1 | General multiplicity estimate, with every geometric and numerical conclusion. | Proved |
| Corollary 2.2 | One-dimensional analytic-subgroup consequence. | Proved |
| Corollary 2.3 | Consequence for disjoint group factors and sampling grids. | Proved |
| Proposition 3.3 | Multigraded intersection bounds, including the multiplicity-sensitive bound. | Proved |
| Proposition 4.7 | Binomial lower bound for contact multiplicity. | Proved |
| Lemma 5.1 | Stabilizer construction and the geometric counting estimate. | Proved |
| 1987 addendum, p. 398 | Corrected common-shift vanishing in a connected ambient group. | Proved |
| 1987 addendum, p. 398 | Connected converse with first-input addition degree at most two. | Open |
The full statement package has 27 compiled theorem statements: all thirteen numbered results, both addenda, eleven supporting or correction targets, and the aggregate goal. The fourteen original admission-free definition bundles and one erratum supplement supply their actual geometric and algebraic objects. All 41 original statement and definition items have independent blind readbacks. The eight milestones above retain their existing identities.
The supporting clauses include Hilbert-polynomial existence, primary components and Facts A–E; the geometric interpretation of mixed degree; the actual counterexample after Proposition 3.3; translation operators and their comparison with intrinsic ideals; contact invariance; translation-invariance and embedding remarks; the component/stabilizer construction; counting estimates; and the exact Masser–Wüstholz Theorem I consequence claimed on p.361. Proof-internal recursive ideals and tangent/exponential arguments belong to their corresponding theorem proofs.
Source corrections are visible for review. Lemma 3.1 and Corollary 2.2 require positive equation degrees; explicit zero-degree counterexamples are required goal clauses. Both addenda require positive ambient dimension; the goal also requires their dimension-zero counterexamples. The three-generator ideal printed on p.370 is nonradical, contrary to the printed word “prime”; its valid degree-four versus length-six counterexample is preserved, together with an explicit nilpotent witness. Connectedness is stated in the translation-invariance and Lange reembedding remarks. These are mathematical corrections documented during the source comparison, beyond the author's 1987 errata. They are not silent changes to the source.
Existing completed proofs and all 24 unaffected goal clauses are retained. The corrected forward addendum is Proved. The repaired converse and the replacement collection goal remain Open; two general geometry lemmas also remain on the frontier. The existing proved local-algebra references remain reusable ingredients.
The output is a reusable development of the paper's multigraded commutative algebra, translation and differential operators, geometric multiplicity theory, and zero estimates. All source results are required for completion. A downstream Weierstrass application can consume a specialization of Theorem 2.1; its needs do not determine the scope or completion of this mission (example application).
These are established mathematical results whose proofs are being formalized. The complete original Theorem 2.1 has an accepted proof. The selected Weierstrass application uses individually named Philippon results, with its application-specific model and subgroup bridges proved in the Senthil mission. The remaining general results stay required even though that application is complete.
Vanishing conditions need not be independent: many can occur on the same component or on translates with a nontrivial stabilizer. Counting coefficients of P therefore does not bound their total multiplicity. The development needs geometric degree for actual components, multiplicities measured by lengths of localized quotient rings, and uniform control of translations in fixed projective embeddings.
Proposition 3.3 supplies multigraded intersection bounds, Proposition 4.7 controls contact multiplicity, and Lemma 5.1 connects the stabilizer to the geometric counting estimate. These three results and Theorem 2.1 now have checked proofs. The remaining corollaries, general geometric and analytic claims, examples and addenda are independent deliverables. Application-specific contact results do not discharge the remaining general statements.
The declarations use the namespace PhilipponMultiplicity. Source statements retain their conclusions, constants, quantifier order, and complex or ℓ-adic scope, with the visible boundary and wording corrections listed above. Intermediate specializations are labelled as such and do not discharge a more general source result. Natural-number and zero-degree conventions have been compared with the original scans; the discovered failures and corrected hypotheses are recorded explicitly for human review. Corrections and inferred conventions must be documented rather than silently changing the source.
The required definitions include embedded commutative algebraic groups and their connected subgroups; products of projective spaces and multihomogeneous coordinate rings; analytic local homomorphisms and intrinsic contact order; multigraded Hilbert polynomials and their degree forms; local component lengths; and finite coset counts. No model may assume the desired multiplicity inequality, hide it as a structure field, or replace geometric degree by an unconstrained function.
The 26 individual target statements, replacement collection goal, original definition bundles, and erratum supplement compile in Lean 4.33.1 with Mathlib revision 0df444a360eaa60ab8c11dca51a86af692955474. Source comparisons and independent readbacks document the fixed statements and their explicit corrections. The goal remains the concrete 26-part full-paper collection. The completed proofs supply reusable Hilbert theory, primary-component multiplicities, polynomial differential operators, bounded translation atlases and the Section 5 construction. Contributions to the remaining targets, including results unused by Senthil, complete the original scope.
The previous aggregate, original sampled-translates addendum, and original converse remain available with their prior submissions and disproofs. The replacement statements changed the mission goal explicitly before proofs were supplied. The corrected common-shift forward statement has since received an accepted proof.
Exact methods for the resource-constrained project scheduling problem — branch-and-bound over activity lists or over start-time assignments, and the lower-bound computations inside them — live or die by how much of the search space can be discarded before it is enumerated. The standard tool is constraint propagation: deducing, from the precedence, resource and time-window data, new precedence relations that every feasible schedule must satisfy, and tighter time windows for the activities. Brucker and Knust's Section 3.6 (doi:10.1007/978-3-642-23929-8) presents the family of interval consistency tests — input, output, input-or-output and their negations — that constraint-programming schedulers apply at every node of the search, following Carlier and Pinson (An algorithm for solving the job-shop problem, Management Science 35, 1989, doi:10.1287/mnsc.35.2.164) and Baptiste, Le Pape and Nuijten (Constraint-Based Scheduling, Kluwer, 2001, doi:10.1007/978-1-4615-1479-4). Every test is an instance of one theorem, Theorem 3.7, and its cumulative-resource analogue, Theorem 3.8. Those two theorems, and the tests as their corollaries, are this mission.
The instance is the RCPSP of mission I: activities with integer processing times , renewable resources with capacities and demands , and precedence arcs; a schedule is an integer start-time vector , feasible when it meets the precedences and never exceeds a capacity. Section 3.6 adds three things.
Relations. A conjunction holds in when . Two activities are parallel, , when they overlap for at least one time unit, and a disjunction is the negation of that: or . The instance carries a set of conjunctions and a set of disjunctions that every feasible schedule must satisfy; initially is the precedence relation and the pairs whose combined demand exceeds some capacity, and propagation adds to them.
Disjunctive sets. A set of at least two activities is disjunctive when any two of its members are related by a disjunction or a conjunction, so no two are ever processed together: the activities of a unit-capacity resource, the jobs of a single machine, the operations of one job in a shop. Its total processing time is .
Time windows. Each activity has a head and a deadline , and a feasible schedule has and . An activity starts first in a set when no activity of starts earlier, and ends last when none completes later.
For a cumulative resource the work of activity is and .
Let be a disjunctive set, , and proper subsets of with . If
then in every feasible schedule an activity from starts first in or an activity from ends last in .
If some nonempty has , no feasible schedule exists.
For nonempty and : if then for all ; symmetrically, if then for all .
For , : if then starts first in or ends last in , and when .
For a cumulative resource , and proper subsets of : if then an activity from starts first in or an activity from ends last in .
Theorem 3.7 is the single statement behind a whole toolbox. Every interval consistency test in the literature — the input and output tests that fix a new conjunction, the input-or-output test, the negation tests that only shrink a window — is the theorem with a particular choice of and , and the book's Section 3.6.4 derives them one by one. A propagation engine that applies these tests to a fixpoint is what makes branch-and-bound for the job shop and the RCPSP practical; Carlier and Pinson's solution of the 10×10 job-shop instance is the historical demonstration. Theorem 3.8 extends the same reasoning from disjunctive to cumulative resources by replacing "no overlap" with "at most units per time unit", the energetic-reasoning viewpoint that the rest of Section 3.6.5 develops.
The results are elementary and proved; formalizing them fixes, once, what "feasible" means in the presence of the relation sets and and time windows, on top of the RCPSP model of mission I. That layer is reusable: the start-start distance matrix of Section 3.6.2 and the symmetric-triple rules of Section 3.6.3 are statements about the same schedules and the same relations. Nothing here is on the platform or in Mathlib.
The obvious argument for Theorem 3.7 is the correct one, and its difficulty is in the bookkeeping. If no activity of starts first and none of ends last, the first starter is some and the last finisher some , and every activity of is processed inside . Since the activities of a disjunctive set are pairwise non-overlapping, their total length fits in that interval, contradicting (3.121). The formal work is the packing lemma: pairwise disjoint integer intervals inside an interval of length have total length at most , which requires ordering the activities by start time and an induction that Mathlib does not supply.
The subtle point is the restriction in (3.121). The book allows it because "an activity which starts first cannot complete also last" when there are at least two activities with positive durations. The formal statement reads "starts first" and "ends last" with , which makes the theorem true without a positivity hypothesis: when the restriction empties the index set, , the conclusion holds because cannot be both the unique first starter and the unique last finisher of a disjunctive set with two or more members. A solver should expect to handle that corner separately.
For the tests the extra step is turning "starts first" into a conjunction , which uses the disjunction between and together with positive processing times: with an activity could start at the same instant as without violating the disjunction, so the tests carry the positivity hypothesis that Theorem 3.7 itself does not need. Theorem 3.8 replaces the packing lemma by a work-counting lemma: over an interval of length a resource of capacity supplies at most units, and every activity of consumes of them.
Schedules are integer start-time vectors on Fin n, as in missions I and II, and a feasible
schedule of this mission is one that is FeasibleSchedule for the RCPSP instance (mission I),
respects the arcs of (RespectsArcs, mission II), satisfies the disjunctions of and lies
within the time windows; these four hypotheses are the book's "feasible schedule" in Section 3.6
and are carried on every statement, although the arguments use only the last two, or, for
Theorem 3.8, the resource constraint and the windows. The sets and are parameters, not
derived from the instance, since propagation enlarges them.
Every inequality "" is stated as the family of inequalities over the same index pairs. This is equivalent, avoids natural-number subtraction, and gives an empty index set the value the convention would: the hypothesis is then vacuous. "Starts first" and "ends last" use . Proper-subset hypotheses , are the book's; the first infeasibility test needs nonempty, and the input-or-output test needs .
A trivializing reading is ruled out on the disjunctive side by the nonemptiness hypotheses (an empty would make the infeasibility test's family vacuous and its conclusion false) and on the cumulative side by the observation that Theorem 3.8 with asserts infeasibility, which is the book's intended reading. Welcome contributions beyond the milestones: the input-negation and output-negation tests, the window-tightening rules of Section 3.6.4, and the SSD-matrix results of Section 3.6.2.
When several processors share a workload, the first question is how long the workload takes if it is spread out as well as possible. If a job may be interrupted and resumed later, possibly on another processor — preemption, the setting of operating systems, of communication links and of any resource that can be time-shared — the answer is a closed formula, and it is one of the oldest results in scheduling: McNaughton's wrap-around rule of 1959 (Scheduling with deadlines and loss functions, Management Science 6, doi:10.1287/mnsc.6.1.1) shows that on identical machines the optimal makespan is the larger of the longest job and the average load. Horvath, Lam and Sethi (A level algorithm for preemptive scheduling, Journal of the ACM 24, 1977, doi:10.1145/321992.321995) extended this to machines of different speeds, and Gonzalez and Sahni (Preemptive scheduling of uniform processor systems, Journal of the ACM 25, 1978, doi:10.1145/322047.322055) gave the fast algorithm with few preemptions. Brucker's Chapter 5 (doi:10.1007/978-3-540-69516-5) presents the level-algorithm version, and this mission formalizes the statements it proves about schedules.
There are jobs with processing requirements and uniform machines with speeds : running job on machine for a period of length performs units of its requirement, so the whole job would take time units there. Identical machines are the case .
A preemptive schedule is a finite list of pieces, each a job, a machine, a start time and a stop time. It is feasible for the data when every piece lies in , no two pieces on the same machine overlap, no two pieces of the same job overlap (a job is on at most one machine at any instant), and every job receives total work exactly over its pieces. Its makespan is the largest stop time; the completion time of job is the largest stop time of one of its pieces. A schedule is nonpreemptive when every job consists of a single piece.
Following Section 5.1.2 the data are sorted, and , with , and one writes , . For a set of jobs, if and otherwise: the largest combined speed that jobs can use at one instant.
The optimal makespan of is the bound (5.5):
in the sense that some feasible preemptive schedule has makespan exactly and no feasible preemptive schedule has a smaller makespan.
Every feasible preemptive schedule has makespan at least .
On identical machines, is a lower bound on the makespan and is attained by some feasible preemptive schedule.
The jobs can be scheduled preemptively within if and only if for every set of jobs.
For with nonnegative weights there is an optimal schedule without preemption.
Theorem 5.8 turns an optimization over an infinite family of schedules into a formula in the data, and the formula is tight in both directions: each of its terms is a resource bound that some schedule meets exactly. That is what makes preemptive makespan minimization one of the few parallel-machine problems that is solvable at all — its nonpreemptive counterpart is NP-hard (p. 124) — and it is why the preemptive relaxation appears as a bound inside branch-and-bound methods for the nonpreemptive problems.
Condition (5.8) is the form in which the result is reused. It is a Hall-type condition, one inequality per set of jobs, and it is exactly what Section 5.1.2 needs to prove Theorem 5.9, which decides by a maximum flow in an expanded network. Theorem 5.7 is the complementary statement for the other classical objective: for total weighted completion time preemption buys nothing, so the nonpreemptive solutions of Section 5.1.1 are optimal in the larger class too.
On status: every statement here is classical and proved, and the formalization adds a checked
model of preemptive schedules. Mathlib has no scheduling material, and the platform's
SchedulingAlgorithms series so far models only single-machine sequences (missions I, II, IV)
and two-machine permutation flow shops (mission III), none of which allow a job to be split. The
piece-list model of this mission is the first reusable object for preemptive and parallel-machine
problems, and the later sections of Chapter 5 — , — are stated in it.
The obvious first idea for the goal is to run McNaughton's rule with the speeds ignored. It fails on uniform machines: filling machines one after another does not respect the constraint that a long job on a slow machine is not done when a short job on a fast one is. The correct idea is the level algorithm — always process the jobs of highest remaining requirement on the fastest free machines, sharing machines among tied jobs — and the difficulty is in the analysis rather than the idea. The proof of Theorem 5.8 has to show that the schedule it produces ends exactly at one of the terms of : either no machine idles before the end, giving , or the machines finish in speed order with the first jobs busy from time , giving . Making that case analysis rigorous requires tracking that the order of remaining requirements is preserved over time (the invariant (5.6)) and that ties are broken consistently.
A second, formal difficulty is that the level algorithm's output is defined by continuous-time events (the next completion, the next time two levels coincide), so producing an explicit finite list of pieces with the required properties is itself a construction. Any proof must build a concrete schedule; "the infimum of makespans equals " is not the goal.
For the lower bound the trap is the opposite: it is tempting to argue only with total capacity , which gives but not . The latter needs the rule that a job is on at most one machine at a time, so that jobs run at combined speed at most ; a model that let a job be split across machines simultaneously would make the theorem false, and the definition of feasibility rules it out explicitly.
A schedule is a List of Pieces over jobs Fin n and machines Fin m, with real start and
stop times. Feasibility is the three-part condition of the Setting, disjointness of two pieces
meaning one stops no later than the other starts. Work is measured with the machine's speed, so
the same definitions cover identical machines as the constant speed . Pieces of length zero
and unsorted lists are allowed; both are harmless.
The sorted orders are hypotheses Antitone p and Antitone s, the speeds and requirements are
positive, and, where the book assumes it, . The book's normalization is
not assumed: every statement here is invariant under scaling all speeds, and the book uses the
normalization only for a running-time estimate. The bound is defined as the maximum of an
explicit nonempty finite set, so no supremum of an empty or unbounded set occurs; takes a
proof that so that is meaningful.
Two things are deliberately not stated. The level algorithm itself is not transcribed: Theorem 5.8 is stated as the existence of an optimal schedule of makespan , which is what its proof establishes. And Theorem 5.9, the flow characterization for , is left for a later mission, since it needs release times and the expanded network on top of this model.
A trivializing reading is excluded by the existential form of the goal and of Theorem 5.7: each asserts that an optimal schedule exists, not merely that any optimal schedule has a property. A proof of the goal has to construct a schedule; a proof of Theorem 5.7 has to construct a nonpreemptive one that beats every preemptive competitor. Contributions welcome beyond the milestones: a general lemma that a feasible schedule can be normalized to sorted, positive-length pieces, and a proof that (5.8) for the sets is equivalent to .
Cryptography before 1949 was a catalogue of ciphers and of the tricks that broke them. C. E. Shannon's Communication Theory of Secrecy Systems (Bell System Technical Journal 28(4):656–715, 1949) replaced the catalogue with a probabilistic model: a cipher is a family of invertible maps from messages to cryptograms, indexed by a key drawn from a known distribution, and the cryptanalyst's state of knowledge after an interception is the a posteriori distribution over messages. Part II of that paper asks when interception conveys nothing at all. The answer — perfect secrecy — is the origin of the one-time pad's security proof and of the modern habit of defining security as the indistinguishability of a posteriori from a priori beliefs.
This mission formalizes §10, Perfect Secrecy: the definition, the necessary and sufficient condition (Shannon's Theorem 6), the counting bound that the key set be at least as large as the message set, the cyclic system that attains the bound, the Latin-square description of the systems that attain it, and the entropy form of the constraint.
A finite secrecy system consists of three finite sets — the messages , the keys , the cryptograms — together with an enciphering map for each key . Each is non-singular, i.e. injective, so a receiver who knows the key deciphers unambiguously. The key is drawn from an a priori distribution , , and the message from an a priori distribution , , independently of the key.
Three derived quantities carry the theory. The key weight
is the probability that the cryptogram is given that the message is . The cryptogram probability is , the probability of obtaining from any cause. The a posteriori probability of the message after interception is , defined for cryptograms with .
A system has perfect secrecy when, for every a priori message distribution and every cryptogram that can occur, for all : interception leaves the cryptanalyst's probabilities unchanged. Quantifying over all a priori message distributions is Shannon's requirement that the equality hold "independently of the values of ".
Finally, and are the entropies of the message and key choices.
Equivalently, does not depend on : the total probability of the keys carrying to a given is the same as that of the keys carrying to the same . The goal fixes nothing about the sizes of the three sets, and no structure on the key distribution beyond normalization.
Theorem 6 is the structural fact behind every later statement in the section: the counting bound, the Latin-square description, and the entropy constraint are all read off from the independence of from . The counting bound and the entropy constraint are the precise sense in which unconditional secrecy is expensive — the key must be at least as uncertain as the message — and they are the reason practical cryptography moved to computational assumptions. The cyclic system supplies the matching construction, so the bound is sharp.
These results are classical and have textbook proofs; what this mission produces is a machine-checked development of them from a single explicit model of a finite secrecy system, including the model itself. Mathlib has extensive measure-theoretic probability and some information theory, but no formalization of Shannon's secrecy systems, of perfect secrecy, or of the results of §10. The shared definition layer — ciphers with injective enciphering maps, key weights, a posteriori probabilities, the perfect-secrecy predicate, finite entropy — is reusable for the rest of the paper: the equivocation results of §11–§13, unicity distance, and the ideal-system theorems of §16 all rest on the same objects.
The two directions of Theorem 6 are not symmetric. Sufficiency is a computation. Necessity has to extract, from a statement about a posteriori probabilities that only constrains cryptograms of positive probability and only for one distribution at a time, a statement about key weights alone; the useful instantiation is a distribution giving every message positive mass, and the cryptograms of probability zero must be handled separately rather than ignored.
The counting and Latin-square milestones are finite combinatorics, but the obvious argument — "for each message and cryptogram pick the key joining them" — needs the keys picked to have positive probability, and the model permits keys of probability zero; this is exactly where an argument that looks complete on paper leaves a gap in Lean. The entropy milestone is the hardest of the five: passes through the conditional entropies , and Mathlib's information-theoretic API is stated for measures rather than for finite weighted sums, so either a bridge or a self-contained development of finite conditional entropy is required.
Messages, keys and cryptograms are arbitrary finite types with decidable equality on cryptograms; probabilities are real-valued functions with explicit non-negativity and normalization hypotheses rather than PMF, so that all sums are finite sums. A cipher bundles the enciphering maps, their injectivity, and the key distribution. The a posteriori probability is a real quotient, so it evaluates to zero when the cryptogram has probability zero; every statement that mentions it carries the hypothesis that the cryptogram has positive probability. Entropy uses the natural logarithm, with vanishing at ; all inequalities are therefore in nats and are unaffected by the choice of base.
One trivializing formalization is ruled out explicitly: perfect secrecy is quantified over all a priori message distributions, not over one fixed distribution, and it is not vacuous — the cyclic milestone exhibits a family of ciphers satisfying it for every .
Contributions welcome beyond the milestones: finite conditional entropy and its chain rule as reusable lemmas, the equivalence of the present perfect-secrecy predicate with the statistical-independence formulation of message and cryptogram, and the extension of the model towards §11's equivocation , which is the next result of the paper to formalize.
Lagrangian mechanics is the standard reformulation of Newtonian dynamics used throughout physics: a mechanical system is described not by forces but by a single scalar function, the Lagrangian, and the equations of motion are the stationarity conditions of an integral of that function along a path. The organizing result of the theory is Noether's theorem: every continuous symmetry of the Lagrangian produces a quantity that is constant along every solution of the equations of motion. Conservation of momentum, of angular momentum, and of energy all arise this way, and the same statement — in field-theoretic form — underlies conservation laws in electrodynamics, general relativity, and the Standard Model.
This mission formalizes the Lagrangian chapter of David Tong's Classical Dynamics (University of Cambridge, Part II Mathematical Tripos), Section 2, in Lean 4. The finite-dimensional statement it targets is the one a physics student meets first, and it is the statement on which the field-theoretic versions are modelled.
The configuration space is ; a path is a map , written , and its velocity is the coordinatewise time derivative . A Lagrangian is a function
smooth in all arguments; Tong writes it .
The action of a path over is
A path is a motion (it solves Lagrange's equations, Tong eq. 2.43) when
both partial derivatives being evaluated at . The quantity is the generalised momentum conjugate to (eq. 2.44), and
is the Hamiltonian (eq. 2.47). A function of time is a constant of the motion when it takes the same value at all times (eq. 2.46).
A one-parameter family of paths is a map , , whose member is the path under study. It is a continuous symmetry of (eq. 2.52) when
the dot denoting .
If is a continuous symmetry of and the path solves Lagrange's equations, then
takes the same value at every time .
Noether's theorem is the statement that converts a structural property of a system — invariance under translations, rotations, or time shifts — into a number that does not change during the evolution. It is what licenses the use of conserved quantities as first integrals when solving mechanical problems, and it is the finite-dimensional prototype of the conservation laws that classify symmetries in field theory.
Formalizing it produces a small but reusable Lean development: a definition layer for
Lagrangians and their partial derivatives on , Lagrange's equations as a
predicate on paths, the action functional, and the notion of a constant of the motion.
Mathlib currently has the calculus these statements need — deriv, fderiv,
ContDiff, interval integrals and differentiation under the integral sign — but no
Lagrangian-mechanics layer built on top of it, so the definition layer of this mission
is itself a contribution. The results here are classical and completely proved in the
mathematical literature; what is missing is machine-checked proofs.
The physics derivations are three lines each, and all three lines are the ones that Lean does not grant for free. Tong's proof of Noether's theorem differentiates with respect to by the chain rule, exchanges with on , and then substitutes Lagrange's equations; in Lean, each of these steps requires the corresponding differentiability side conditions to be produced from the smoothness hypotheses, and the exchange of the two partial derivatives requires a symmetry-of-second-derivatives argument. The principle-of-least-action target additionally needs differentiation under the integral sign and the fundamental lemma of the calculus of variations (a continuous function integrating to zero against every smooth bump supported in vanishes there), neither of which is available in the exact shape needed.
The naive route — defining everything through fderiv on the product space and hoping
simp unfolds the chain rule — does not close these goals, because the partial
derivatives are taken in single coordinate slots and must be reassembled into a total
derivative along the path.
Conventions fixed by the Lean development and worth knowing before starting:
Fin n → ℝ for a fixed n : ℕ, and time is all of ℝ;
no interval restriction is imposed except in the action functional.derivs taken in a single slot, using
Function.update to vary one coordinate. Where a derivative fails to exist, deriv
returns 0; the smoothness hypotheses on L, on paths, and on families are stated
so that this junk value is never the content of a statement.HasDerivAt, so they assert the existence of
d/dt (∂L/∂q̇ᵢ) as well as its value.∀ t₁ t₂, F t₁ = F t₂, i.e. constancy on all of ℝ,
not merely a vanishing derivative.s-derivative vanishes for all s and t,
not only at s = 0.n = 0 is included in every statement; there it is trivially
true. It is not, however, a trivializing formalization of any target: for n ≥ 1 the
hypotheses are satisfied by standard systems (a free particle, a particle in a central
potential), so none of the statements is vacuous.Contributions welcome: proofs of any supporting target, reusable lemmas relating the
coordinate-slot derivs used here to fderiv on the product space, and worked
instances (translation invariance giving momentum conservation, rotational invariance
giving angular momentum) built on the same definition layer.
Periodic signals are the basic objects of audio processing, electrical engineering and the analytic study of partial differential equations, and the standard first-course tool for handling them is the Fourier series: a decomposition of a -periodic signal into sinusoids of frequencies . The identity that makes the decomposition quantitatively useful is Parseval's theorem: the average power of the signal equals the sum of the squared moduli of its Fourier coefficients, so energy can be accounted for frequency by frequency.
This mission formalizes the chain of results leading to that identity exactly as it is presented in the open-access undergraduate textbook Análise de Fourier — Um Livro Colaborativo (Esequia Sauter and Fabio Souto de Azevedo, orgs., UFRGS/REAMAT, 2022): the orthogonality relations of the trigonometric system (Teorema 3.2.1), the coefficient formulas (Teorema 3.2.2), the passage to the exponential form (§4.2), and the Parseval identity itself (Teorema 5.1.1).
Fix a period and write
For a real-valued integrable on the book defines the trigonometric Fourier coefficients (eq. 3.23)
and for a complex-valued the exponential coefficients (§4.2, eq. 4.8)
which for real satisfy . The average power of a -periodic signal is (Definição 5.1.1)
A trigonometric polynomial of degree (Definição 3.2.1) is
For and continuous and -periodic,
The book states the identity for a periodic function "representable by a Fourier series"; the mission fixes continuity as the concrete sufficient regularity hypothesis, which keeps the statement non-vacuous and true, and leaves the sum as an unordered sum over so that no summation order is smuggled in.
Parseval's identity converts an integral estimate into a coefficientwise one: it is the bridge between the time domain and the frequency domain, underlies the RMS computations of the book's own electrical-engineering examples (Exemplo 5.1.2), and yields closed-form evaluations of numerical series (the book's Exemplo 5.1.3 computes a series of reciprocal odd squares this way).
Mathlib already contains the abstract Hilbert-space form of the identity for the Fourier
basis on the circle AddCircle T. What this mission adds is the concrete, classroom-level
layer that an engineering or physics reader recognizes: coefficients written as interval
integrals of a function on with an explicit period , the real
form, the exponential form and the dictionary between them, and the orthogonality
relations as free-standing computations. That layer is reusable by any later mission in
this series (transformada de Fourier, equação do calor, equação da onda).
The obvious route — expand , substitute the Fourier series, and integrate term by term — is exactly the book's proof, and it is the step that a formal development cannot take for granted: interchanging the integral with an infinite sum requires a convergence theorem that the book explicitly places outside its scope (the remark after Definição 3.2.3), and pointwise convergence of the Fourier series of a merely continuous function is false in general. A complete proof must run through convergence rather than pointwise convergence, so the central task is to connect the interval-integral coefficients on with the Hilbert-space machinery for the circle, including the change of variables and the normalization of Haar measure.
The development commits to the following conventions, all fixed in the definition bundle so that every statement shares them:
The goal is not trivializable: the hypotheses (a positive period, a continuous periodic function) are satisfiable, and both sides of the identity are non-degenerate already for .
Useful infrastructure beyond the mission: the dictionary between the interval-integral
coefficients on and Mathlib's fourierCoeff on AddCircle T, and the
orthogonality relations as standalone integral computations.
Density functional theory (DFT) is the computational method behind a large fraction of present-day electronic-structure calculations in condensed-matter physics, quantum chemistry and materials science. Its appeal is a drastic reduction in dimension: instead of an -electron wavefunction living on , one works with the one-particle density on . That reduction is not a modelling approximation; it is licensed by two theorems of Hohenberg and Kohn (1964), for which Walter Kohn shared the 1998 Nobel Prize in Chemistry. Every statement in this mission comes from the two Hohenberg–Kohn theorems and the identities the source text uses to set them up.
Timeline. The Thomas–Fermi model (1927) introduced density-only energy expressions but is known not to bind molecules (Teller, 1962). Hohenberg and Kohn (1964) proved that for non-degenerate ground states the external potential is fixed by the ground-state density, and that the exact energy functional is minimized at that density. Kohn and Sham (1965) turned this into a practical scheme by reproducing the density with non-interacting orbitals. Levy (1979) replaced the Hohenberg–Kohn functional, defined only on ground-state densities, by a constrained-search functional defined on all -representable densities.
Fix an electron number. The objects are:
Two functionals act on these. The universal part is the expectation value of the kinetic-energy and electron–electron-interaction operators, which do not depend on the system. The external part is the expectation value of the external potential, and the identity that makes DFT possible is that it depends on only through its density:
The total energy is . A state is a ground state of when it minimizes , and a nondegenerate ground state when it is the strict minimizer. A density is -representable when it is the density of a nondegenerate ground state of some potential.
with equality exactly when , the ground-state density of . Part (1) is the source's Corollary 1 — the density determines every ground-state property; part (2) is its Theorem 2.
The supporting statements are the density-determines-the-state half of Theorem 1, the potential-is-determined-up-to-a-constant half of Theorem 1, the well-definedness of the universal functional , the variational principle, the Levy constrained-search identity, the identity on , and the normalization of a Kohn–Sham density built from orthonormal orbitals.
The Hohenberg–Kohn theorems are what justifies computing with rather than : without Theorem 1 the density would not carry enough information to determine the system, and without Theorem 2 there would be no variational principle to minimize. They are also the reason the exchange–correlation functional is a well-posed object of study: it is the unknown part of a functional that the theorems guarantee exists.
Formalizing them contributes a reusable variational core — ground states, universal functionals, -representability, constrained search — plus two genuinely analytic statements about : the density-pairing identity for the external energy and the "a.e. nonvanishing eigenfunction forces the potential difference to be constant" lemma. The latter two are where the mathematical work lies, and both are stated here over the concrete configuration space , not in the abstract model.
The reductio that proves Theorem 1 is short, and it is short in Lean too. The difficulty of DFT as mathematics sits elsewhere, and this mission isolates it:
sInf has to be handled.The variational layer is stated over an abstract model carrying the data actually used by the Hohenberg–Kohn argument: wavefunctions, potentials, densities, a universal functional , a density map, and an external-energy pairing. No Hilbert-space structure is assumed there, because the argument uses none; what it uses is that does not depend on and that the external energy sees only through . The two analytic milestones are stated concretely instead, on over , with .
Conventions fixed in Lean: energies are real; "nondegenerate ground state" means strict
minimizer, so degenerate ground states are outside the scope of every statement that
mentions it, exactly as in the source's caveat that the original theorems hold for
non-degenerate ground states; -representability quantifies over nondegenerate ground
states only; the Hohenberg–Kohn functional takes the junk value off the
-representable set, so statements about it always carry the representability hypothesis;
the Levy functional is an sInf, which is on an empty fibre.
The statements are not trivially satisfiable: the model has no hypothesis forcing the wavefunction set to be empty or the ground states to be unique, part (1) of the goal quantifies over every real-valued observable, and part (2) asserts an equality case, not merely an inequality. Contributions of the two analytic milestones are the most valuable: they need marginal-density machinery over product spaces that is reusable for any many-body formalization.
The discrete Fourier transform (DFT) is the basic spectral tool of digital signal processing, numerical analysis and computational number theory: it turns a finite sample of a signal into its frequency content, converts convolution into pointwise multiplication, and underlies fast integer and polynomial multiplication. Computed straight from its defining sum, a transform of length costs arithmetic operations. The fast Fourier transform (FFT) is any method that computes the same values in operations, and the radix-2 Cooley–Tukey algorithm is the most widely used of them.
The history is long and worth stating precisely. Gauss used an equivalent divide-and-conquer scheme in unpublished work of 1805 on the orbits of Pallas and Juno, without any complexity analysis. Danielson and Lanczos (1942) published the doubling step — the even/odd splitting identity that the recursion rests on — for x-ray crystallography. Good (1958) gave the prime-factor algorithm for coprime factorizations. Cooley and Tukey (1965) published the general composite-size algorithm together with the analysis that made the method famous. What is not known is a matching lower bound: whether a DFT of length can be computed in arithmetic operations is open.
Fix a length and a sequence of complex samples , modelled as a function . The DFT of length evaluated at frequency index is
and the inverse transform of a spectrum is
Only the first samples enter the sum; the frequency index is an arbitrary natural number, and the transform is -periodic in it.
The radix-2 Cooley–Tukey algorithm is the recursion on the exponent , with :
where is the transform of the even-indexed subsequence and that of the odd-indexed one. The factor is the twiddle factor of the merge stage.
The arithmetic cost of that recursion is modelled by two counters. At the merge stage of a transform of points there are butterflies; each butterfly performs one twiddle-factor multiplication and two complex additions, and the two recursive calls are on points each:
The goal fixes no operation count and no rate: it asserts only that the recursion computes exactly the transform it is supposed to compute, at every frequency index, for every input. That is the statement an implementation's correctness rests on, and it stays true whatever the cost analysis of the algorithm turns out to be.
The milestones are the intermediate statements the source separates out: -periodicity of the transform in the frequency index; the Danielson–Lanczos even/odd splitting
Fourier inversion for ; the conjugate symmetry of the spectrum of a real input; and the closed forms and of the two cost recursions, with .
Fourier inversion and the splitting identity are the two facts that make the transform usable at all: the first says no information is lost, the second is the single algebraic step from which every Cooley–Tukey variant — mixed radix, split radix, iterative bit-reversed implementations — is assembled. The operation counts are what justify the claim that the algorithm is rather than : for they are the difference between roughly and roughly complex multiplications.
None of this is open mathematics; all of it is classical. What this mission produces is a self-contained machine-checked development of the radix-2 algorithm against an explicit definition of the DFT: a recursive algorithm, a proof that it agrees with the defining sum at every index, and a proof of its arithmetic cost. Mathlib carries Fourier theory for characters of finite abelian groups and for , but the mission is stated against the elementary sum above so that the statements can be read, and audited, without that background.
The obvious induction on does not close by itself. Unfolding the recursion gives the even and odd sub-transforms evaluated at the reduced index , while the splitting identity produces them at . Closing the gap needs periodicity of the length- transform in its frequency index — which is why periodicity is a milestone rather than a step inside the main proof. The second recurring obstacle is bookkeeping with the exponential: and are the only analytic inputs, but they must be applied under a finite sum with natural-number casts, where -subtraction and division by a cast that may be zero are easy to get wrong.
Sequences are total functions , not vectors of a fixed length: the length appears only as the range of the defining sum, so the transform of length ignores for . The frequency index ranges over all of , and the goal is asserted for every such index, not only for — the periodicity of both sides makes this the stronger, and cleaner, statement.
Degenerate inputs are inside the statements, not excluded by hypothesis. The length- transform is the empty sum ; division by the cast follows Lean's convention , which is why inversion carries the hypothesis and the splitting identity carries . The real-input symmetry is stated for , with -subtraction in .
There is no trivializing reading: the FFT recursion is defined without ever mentioning the DFT, the cost counters are defined by their own recursions, and every hypothesis in the list above is satisfiable.
A complete development needs: the geometric-sum evaluation of for an -th root of unity (for inversion), the even/odd reindexing of a sum over , and the periodicity computations above. All of these are reusable beyond this mission; the splitting identity in particular is the entry point for a mixed-radix or split-radix sequel.
The NRL Plasma Formulary is a standard desk reference of the plasma-physics community: a compilation of the formulas, constants and unit conversions used in daily practice. Its opening section, "Numerical and Algebraic" (p. 3 of the 2013 edition), collects the few purely mathematical identities the rest of the handbook leans on. Two of them are exact summation formulas rather than approximations, and the first is the Rothe–Hagen identity, quoted there as valid "for all complex , , except when singular" and attributed to H. W. Gould's work on binomial coefficient summations.
Unlike the handbook's numerical entries, this identity is a theorem with a precise hypothesis set, and it is exactly the kind of entry a reader takes on trust. It generalizes the Vandermonde convolution, it is the coefficient identity underlying the generalized binomial series, and it specializes to Abel's binomial theorem. Formalizing it turns one line of a reference handbook into a machine-checked statement and produces, as a by-product, a reusable Lean development of binomial coefficients with an arbitrary complex upper index.
For a complex number and a natural number , the generalized binomial coefficient is the falling factorial divided by a factorial,
with the empty-product convention . It is a polynomial in of degree , and it agrees with the usual binomial coefficient when is a natural number.
Fix complex parameters , , and, for , consider the Rothe factor
which is defined whenever . The factor in the denominator cancels against the leading factor of the falling factorial, so extends to a polynomial in and : and, for ,
Both forms occur in the literature; the mission carries both and asks for the comparison between them, because the quotient form is the one printed in the handbook while the polynomial form is the one that carries no side condition.
For complex and , provided and for every , and ,
This is the handbook's line, with its "except when singular" proviso made explicit as the three non-vanishing hypotheses.
in terms of the polynomial form above. This version holds for all complex , with no exceptional locus, and implies the printed form wherever the latter's denominators are non-zero.
The first is Pascal's rule for a complex upper index; the second is the Vandermonde convolution over , which is the case of the goal.
The identity is the convolution law of the generalized binomial series: the formal power series solving satisfies , so the goal is the statement read off coefficientwise (Graham–Knuth–Patashnik, Concrete Mathematics, §5.4). Consequences include Abel's binomial theorem, the Lagrange-inversion count of -ary trees, and ballot-type identities in lattice-path enumeration.
Mathlib provides the Vandermonde convolution for natural-number arguments (Nat.add_choose_eq), the ascending and descending Pochhammer polynomials, and Ring.choose for binomial rings; it does not contain the Rothe–Hagen identity in any form. The four statements of this mission are therefore new formal content, and the complex-upper-index binomial API they force is reusable well beyond the mission.
The identity is classical and has been proved many times since Rothe (1793) and Hagen (1891), with the modern treatment in Gould's papers; nothing here is open mathematics. What is missing is a machine-checked proof.
The obvious attack — induction on using Pascal's rule — does not close as stated: the summand is not Pascal-stable, since shifting by moves for every at once, and the induction hypothesis is about a different family. The standard proofs instead treat both sides as polynomials in and for fixed and , verify the identity on an infinite set of points where a combinatorial reading is available, and conclude by the identity theorem for polynomials; or they extract coefficients from the generalized binomial series via Lagrange inversion. Either route needs infrastructure: a two-variable polynomial-identity argument over , or a formal-power-series compositional inverse.
A second, more prosaic difficulty is the singular locus. The printed identity divides by for every , and in Lean division by zero returns zero rather than failing, so a formalization that drops the non-vanishing hypotheses states a different — and in general false — claim.
All statements are over . The generalized binomial coefficient is defined as an explicit product over Finset.range k divided by (k ! : ℂ), so holds definitionally and no Nat.choose coercion enters. Sums run over Finset.range (n+1), with the complementary index written as the truncated natural subtraction n - k; inside that range this is the ordinary , so the truncation convention is never exercised.
The proviso "except when singular" is formalized as three explicit hypotheses — for all , for all , and — rather than by relying on Lean's junk value for division by zero. These hypotheses are satisfiable (for instance , ), so the goal is not vacuous, and they constrain only denominators, so they do not trivialize the sum.
The singularity-free milestone commits to a total function rotheA defined by cases on , with value at ; the comparison milestone pins that function to the quotient form under the non-vanishing hypothesis, so the mission cannot be satisfied by proving facts about a differently normalized object.
Contributions welcome beyond the milestones: complex-upper-index binomial API (symmetry, negation , polynomiality in the upper index), and any formal-power-series development supporting Lagrange inversion.
This mission formalizes the group theory of J. A. Wolf's Growth of finitely generated solvable groups and curvature of Riemannian manifolds, J. Differential Geometry 2 (1968) 421–446 (doi:10.4310/jdg/1214428658), namely its §3 and §4. Wolf opens with the object of study: "If a group is generated by a finite subset , then one has the 'growth function' , where is the number of distinct elements of expressible as words of length on " (p. 421). Section 3 proves that a group with a finitely generated nilpotent subgroup of finite index "is of polynomial growth, and in fact "; section 4 proves that "a polycyclic group, either has a finitely generated nilpotent subgroup of finite index and thus is of polynomial growth, or has no such subgroup and is of exponential growth" (p. 421).
Applying Milnor's companion note, which is the complete mission Milnor: growth of finitely generated solvable groups on this platform, Wolf concludes "that a finitely generated solvable group, either is polycyclic and has a nilpotent subgroup of finite index and is thus of polynomial growth, or has no nilpotent subgroup of finite index and is of exponential growth" (p. 421). That is the Milnor–Wolf theorem, the goal here, and Wolf's Theorem 4.8. It prompted the question Wolf asks on p. 422, "whether every finitely generated group , which is not of exponential growth, necessarily has a nilpotent subgroup of finite index"; Grigorchuk's groups of intermediate growth (1984) answered it in the negative, and Gromov (1981) proved that polynomial growth alone does force a nilpotent subgroup of finite index.
Growth. For a finite subset of a group , the growth function counts
the elements expressible as words of length on , a word
having length (p. 426). The published MilnorWolf.growthFunction S m takes
this as the size of the ball Chou.wordBall S m, the set of products of at most factors from
. Wolf's two kinds of bound (p. 421) are polynomial growth of degree ,
, and exponential growth, with ; the
published predicates are MilnorWolf.HasPolynomialGrowthOfDegreeLE and
Chou.HasExponentialGrowth, the latter with , which is what Theorem 4.3 actually produces.
Both quantify over some finite generating set, and Wolf shows (p. 427, p. 434) that neither depends
on which one.
Growth exponents. For a finitely generated nilpotent with lower central series
, each factor
is finitely generated abelian with free part of rank ; Wolf's (3.3)
sets and . These are
MilnorWolf.growthExponentOne and MilnorWolf.growthExponentTwo.
Polycyclic groups. A solvable group is polycyclic if it satisfies condition (1) of
Proposition 4.1: "There is a normal series
with every quotient
finite or infinite cyclic" (p. 433), where a normal series, as in Kurosh, has each term normal in
the preceding one. This is MilnorWolf.IsPolycyclic. Solvable is Mathlib's
Group.IsSolvable.
The goal is Wolf's Theorem 4.8 (p. 438): "Let be a finitely generated solvable group. If has a nilpotent subgroup of finite index, then is polycyclic and of polynomial growth of degree . If does not have a nilpotent subgroup of finite index, then is of exponential growth."
The milestones are the results it rests on, in the paper's order. Section 3 climbs from two rescaling lemmas (3.4, 3.5) and the growth function of a free abelian group (3.6) through generating sets adapted to the lower central series (3.7) to the polynomial bounds for a finitely generated nilpotent group (3.2), then passes them to a group with a finite-index subgroup (3.11). Section 4 characterises polycyclic groups (4.1, whose one hard input is Hirsch's theorem, a milestone of its own) and proves the polycyclic dichotomy (4.3). Milnor's Theorem, already proved on this platform, appears as a reference milestone: it is what carries 4.3 from polycyclic to solvable, and Wolf notes on p. 437 that "J. Milnor extended the scope of Theorem 4.3 by proving [10] that a finitely generated nonpolycyclic solvable group must be of exponential growth."
Theorem 4.8 is the first growth dichotomy: within finitely generated solvable groups there is nothing between polynomial and exponential, and the polynomial side is exactly the almost nilpotent one. It is the result Gromov's theorem generalises in one direction and Grigorchuk's groups bound in the other, and the reason a group of intermediate growth can be neither solvable nor elementary amenable, which is how those groups were first placed outside both classes. Chou's Theorem 3.2, in the mission Chou: elementary amenable groups, extends exactly this statement to elementary amenable groups and cites both halves.
Section 3 also has independent value. The polynomial bounds with the explicit exponents and from the lower central series are sharper than the usual textbook statement that a finitely generated nilpotent group has polynomial growth, and the pieces along the way are reusable: the growth function of in closed form, and the fact that polynomial growth of a given degree is inherited from a finite-index subgroup.
The hard milestone is the second half of Theorem 4.3, that a polycyclic group with no nilpotent subgroup of finite index has exponential growth. Wolf proves it through Proposition 4.4, which passes to a Malcev completion and reads the eigenvalues of the induced Lie algebra automorphisms; that machinery is not in Mathlib and is not part of this mission. A route that avoids it: a polycyclic group acts on a free abelian factor of its chain by integer matrices, so either some element has an eigenvalue off the unit circle, and then the separation argument of Wolf's p. 437 applies verbatim with that factor in place of the Lie algebra, or every eigenvalue of every element lies on the unit circle, and then Kronecker's theorem, which Mathlib has, makes them all roots of unity, so a finite-index subgroup acts unipotently and is nilpotent. Wolf records in the "Added in Proof" (p. 445) that "B. Hartley recently sent me a manuscript consisting of an alternate proof that a polycyclic group without nilpotent subgroups of finite index is of exponential growth."
Theorem 3.2 is long rather than deep: the lower bound (3.8) and the upper bound (3.9), (3.10) are two inductions along the lower central series, each carrying explicit word-length bookkeeping, and both need the adapted generating sets of Lemma 3.7, in the corrected form stated here: as printed, the lemma fails when a lower-central factor has torsion. Mathlib has no polycyclic groups, so Proposition 4.1 starts from nothing; it also lacks Hirsch's theorem and P. Hall's finite presentability results, which is why Hirsch is stated here as its own milestone rather than assumed.
No new definitions: every notion is taken from the published bundles Chou_Growth and
MilnorWolf_Growth, so the statements here are directly comparable with those of the Chou and
Milnor missions. Growth is measured on closed balls, Chou.wordBall S m being the products of at
most letters from and its cardinality, finite because is a
Finset. "Not of exponential growth" is the negation of an existential, so it is a statement about
every finite generating set.
Conditions (8) to (11) of Proposition 4.1, Proposition 4.4 and the Lie-theoretic apparatus of §4, and the Riemannian geometry of §2, §5 and §6, are outside this mission; Proposition 4.1 is stated here with its first seven conditions, which are the ones §4 uses. Hirsch's theorem is cited, not proved, by Wolf, and appears as a milestone under its own source (Proc. London Math. Soc. 44 (1938) 53–60) rather than as an omission. Cheeger's and Hartley's sharpenings, recorded in Wolf's "Added in Proof", are not included.
No hypothesis is vacuous: the free abelian groups witness Proposition 3.6 and Theorem 3.2, a polycyclic group with no nilpotent subgroup of finite index exists (the semidirect product for a hyperbolic , which is Wolf's own example of the exponential case), and the lamplighter group is finitely generated solvable and not polycyclic, so Theorem 4.8's second alternative is not empty either. Contributions welcome: proofs of any milestone, and the general library results they need, such as polycyclic groups and Hirsch's theorem.
The Schrödinger equation is the equation of motion of non-relativistic quantum mechanics: it determines how the state of an isolated quantum system evolves in time. Its time-independent form is an eigenvalue equation whose eigenvalues are the energies the system can have. The single most elementary illustration that this eigenvalue problem forces energies to be discrete is the particle in a one-dimensional box, also called the infinite potential well: a particle of mass confined to an interval by walls of infinite potential energy. It is the first worked example in essentially every introduction to quantum mechanics, and it is the model in which "quantization" can be derived in a few lines from a boundary-value problem for an ordinary differential equation.
This mission is the first in a series formalizing the mathematical content of the Schrödinger equation. It targets exactly the particle-in-a-box computation: that the admissible energies are for positive integers , and nothing else.
Fix three positive real constants: the reduced Planck constant , the mass , and the width of the box. A stationary state of energy is a function such that
together with the regularity needed to state it — is continuous on and twice differentiable at every interior point — and the Dirichlet boundary conditions imposed by the infinite walls,
Inside the box the potential vanishes, which is why only the kinetic term appears; outside the box the potential is infinite and the wave function is identically zero, which is what forces the two boundary values. A stationary state is non-trivial if for at least one ; the trivial solution is excluded because a physical state must be normalizable to unit probability.
Write for , and for the corresponding unnormalized eigenfunction.
This is the weakest statement that carries the physical content: it fixes no eigenfunction, no normalization and no phase, only the set of admissible energies.
The last two give the converse direction — every really occurs — and the normalization constant that turns into a unit vector of .
The result itself is the standard textbook derivation: confinement plus the time-independent Schrödinger equation yields a discrete spectrum, with a strictly positive ground-state energy (no state of zero energy exists) and level spacing growing like . As a piece of mathematics it is the computation of the Dirichlet spectrum of the one-dimensional Laplacian on an interval, which is the model case for spectral theory on domains.
What the mission adds is a machine-checked version of both halves of that derivation, in a form that later missions in this series can import: a reusable predicate for "stationary state with Dirichlet boundary conditions", the classification of the solutions of with one boundary zero, and the quantization step. Mathlib contains the ordinary-differential-equation uniqueness results and the trigonometric API needed, but does not contain this boundary-value classification.
The physics derivation writes down the general solution and is done. Formally, that step is the whole problem: one must prove that every twice-differentiable solution of on an interval has this form, rather than assume it. The usual route is a uniqueness theorem for the associated first-order system applied to the difference of and an explicit sine solution, which requires converting the second-order scalar equation into a system and supplying a Lipschitz estimate. Two further traps: the sign case must be excluded separately (hyperbolic solutions with two zeros are identically zero), and the equation is assumed only on the open interval, so the values at the walls are reached by continuity rather than by evaluating the equation there.
Wave functions are complex-valued functions of a real variable, ; energies, , , are real. Derivatives are the ordinary derivative of a complex-valued function of a real variable, the second derivative being the derivative of the derivative. The equation is imposed pointwise on the open interval only, with continuity on the closed interval and differentiability at interior points; nothing is assumed about outside , so the junk values of the derivative outside the box never enter a statement. The constants are hypothesized strictly positive, so no division by zero occurs in .
The goal is not vacuous: the supporting target exhibiting shows that its hypotheses are satisfiable for every , so the mission cannot be closed by proving that no non-trivial stationary state exists. Conversely the non-triviality hypothesis cannot be dropped: solves the boundary value problem for every real .
Contributions welcome: the ODE classification lemma (solutions of on an interval), orthogonality of the in , and the finite-well and harmonic-oscillator analogues, which are planned as later missions in this series.
A nucleus whose proton number or neutron number equals one of — the magic numbers — is markedly more tightly bound than the semi-empirical mass formula predicts, and markedly more stable against decay. The empirical evidence for closed nuclear shells at and protons and at neutrons was assembled by Maria Goeppert Mayer in 1948 (Phys. Rev. 74, 235), and explained in 1949 by adding spin–orbit coupling to the shell picture (Phys. Rev. 75, 1969); the shell model that grew out of this work earned Goeppert Mayer and Hans Jensen a share of the 1963 Nobel Prize in Physics.
The textbook first step toward these numbers is purely combinatorial: approximate the nuclear force by an isotropic three-dimensional quantum harmonic oscillator, count the degenerate states at each energy level, double the count for spin, and accumulate. This mission formalizes exactly that step — including its well-known failure — from the Derivation section of the Wikipedia article Magic number (physics).
For the isotropic three-dimensional harmonic oscillator the energy depends only on the sum of the three quantum numbers, so the states at level are the triples of nonnegative integers summing to — the weak -compositions of . Write
Each state accommodates two nucleons of the same kind (two spin orientations), so the shell at level holds
nucleons, and filling every shell of level below takes
nucleons. A shell closure — hence a predicted magic number — is a value . Finally, denotes the seven most widely recognized magic numbers.
The mission's goal collects the content of the derivation:
together with the comparison against : the values lie in ; the magic number is never of the form ; and is produced by the count although .
The positive half of the statement explains three of the seven magic numbers from nothing but degeneracy counting: , and are exactly the first three oscillator shell closures. The negative half is the reason the shell model needs more physics: the oscillator count misses altogether and predicts , which is not magic. Historically this gap is what spin–orbit coupling was introduced to close, and it is the precise sense in which "the correct formula for predicted magic numbers fundamentally relies on spin-orbit coupling".
As a formalization target the mission is small and self-contained: it fixes a reusable Lean model of oscillator shell filling (levels, degeneracies, spin doubling, cumulative capacity) and the stars-and-bars count behind it, against which later, more physical refinements — spin–orbit corrections, the empirical magic-number list — can be stated.
The mathematics is elementary; the work is in the counting and in the negative claims. The degeneracy is stars and bars for three boxes, most directly obtained by fibering over the first coordinate and summing . Asserting that no gives is not a finite check: it needs the closed form together with monotonicity of to rule out every above the range where is computed, since .
Everything is over and finite. The level set is a Finset (ℕ × ℕ × ℕ) cut out of the box by the equation , so its elements are ordered triples: permutations of the same multiset are distinct states, as the physics requires. Spin enters only as the literal factor in , not as an extra coordinate. The cumulative count runs over levels , so and ; the -th predicted magic number is , not . Binomial coefficients are Nat.choose, and the closed forms are stated in the division-free shape to avoid truncated natural-number division. The set of magic numbers is a hard-coded Finset ℕ transcribing the seven numbers listed by the source — it is a datum, not a derived object, and no claim in this mission asserts anything about nuclear stability itself.
The statements about and are what keep the mission from being a pure identity check: quantifies over all , and is a membership fact about the transcribed list.
A convex polytope has both geometric coordinates and a finite pattern of faces. Moving its vertices can change distances and angles while preserving which vertices belong to which faces. Enumeration by combinatorial type asks for those incidence patterns, with geometrically different realizations of the same pattern counted together. Grünbaum’s enumeration theorem establishes that the complete collection can be determined algorithmically when the dimension and number of vertices are prescribed Convex Polytopes, §5.5, p.91.
A d-polytope here is a nonempty convex hull of finitely many points in real d-dimensional coordinate space, with affine span equal to the whole space. A face is obtained by maximizing a linear functional over the polytope; the collection of faces also includes the empty face and the polytope itself. Inclusion orders this collection. Two polytopes have the same combinatorial type when their full face collections admit an order isomorphism. This notion retains the incidence structure and forgets metric measurements.
To describe a type with finite data, label the k vertices by the integers from zero through k−1. Record the set of vertex labels belonging to each nonempty proper face. This family of subsets is the polytope’s scheme. A finite list of finite lists of natural numbers encodes such a family. The order of labels and repeated occurrences of a label do not change the represented subset. A scheme is realized only when its recorded subsets are exactly the vertex sets of the nonempty proper faces: it cannot add faces, omit faces, or introduce labels outside the prescribed range.
The target is a single total computable function
For every dimension d and vertex count k, each scheme in E(d,k) must have a realizing d-polytope with exactly k vertices. Every d-polytope with k vertices must realize some scheme in that output. Finally, if polytopes realizing two output positions have isomorphic full face posets, those positions must be equal. These requirements say that the output contains precisely one representative of every combinatorial type. The existential choice of one function precedes both numerical inputs, so the algorithm must work uniformly for all dimensions and vertex counts.
The result supplies an effective finite classification at each prescribed size. It is stronger than the observation that only finitely many incidence families can be written down: those families need not all arise from real convex polytopes. It also addresses termination, since a total algorithm must return its entire finite answer for every input. The theorem is a known mathematical result in the cited textbook. The present formal target asks for a proof of its stated algorithmic conclusion; no machine-checked proof of that conclusion is asserted here.
The principal difficulty is the connection between finite incidence data and geometric realizability. Combinatorial consistency alone does not supply real coordinates for a convex polytope. The source distinguishes this realizability question from the enumeration conclusion and identifies decidability over the real numbers as relevant to it §5.5, p.91. The mission retains the complete enumeration conclusion, including realizability of each answer, coverage of every type, and absence of duplicate types.
The formal representation uses real coordinates without a rationality restriction and permits nonsimplicial polytopes. Vertices are labelled injectively and exhaustively. The number of vertices is expressed as the number of nonempty zero-dimensional exposed faces. Nonemptiness separates the empty face, whose book dimension is −1, from the natural-valued dimension used to count faces. For finite convex hulls the face collections are finite, so this cardinality has its usual meaning.
Natural-number dimensions include dimension zero. The unique point has one vertex and no nonempty proper faces, and therefore uses an empty scheme. This is an explicit extension of the source’s positive-dimensional scheme convention. Unrealizable dimension and vertex-count pairs require an empty output. Finite-polytope faces, their ordered collection, vertex counts, and finite scheme realizations are the concrete objects needed to state the goal; total computability applies to the complete finite output rather than to individual tests alone.
Reference: Branko Grünbaum, Convex Polytopes, second edition, Springer, 2003, §5.5, Theorem 2 (5.5.2), printed p.91, PDF p.117; definitions in §§2.4 and 3.1. Source text.
A vector whose coordinates are all zero or one records a collection of yes-or-no choices. Taking the convex hull of a collection of these vectors gives a 0/1-polytope, also called a binary polytope. Such objects connect discrete choices with geometry: points describe feasible combinations, while supporting inequalities describe restrictions that every feasible combination satisfies. Grünbaum's discussion in §4.9 emphasizes the role of facets of these polytopes in combinatorial optimization and cutting-plane descriptions.
The question here concerns how many facets can occur when the dimension grows. Restricting every generating point to binary coordinates gives a finite, highly structured collection of possible points. Nevertheless, this restriction still allows polytopes with very many distinct boundary faces. The theorem of Bárány and Pór establishes superexponential growth in the largest possible number of facets. The mission concerns the form of this growth asserted in Grünbaum's 2003 notes, printed page 69a, rather than a prescribed numerical constant.
Fix a natural number d. The ambient space is the real coordinate space R^d. A convex hull consists of all convex combinations of the generating points, meaning weighted averages with nonnegative weights whose sum is one. A polytope is the convex hull of a finite set of points. It is full-dimensional when its affine span is all of R^d; this rules out a polytope lying in a proper affine subspace.
For a binary polytope, every generating point belongs to {0,1}^d. A subset of this cube is automatically finite. Nonemptiness and full dimension are imposed separately in the target. Consequently, the dimension in the facet bound is the actual dimension of the polytope as well as the dimension of its coordinate space.
An exposed face is the set of all points of the polytope that maximize a given linear functional. A facet is a face of affine dimension d−1. Facets are counted as geometric sets. Two different inequalities that expose the same face do not contribute two facets. Write f_(d−1)(P) for this number.
The target is the following existence statement:
The constant c and threshold D are chosen before the dimension d. The polytope P may depend on d. The assertion therefore supplies a witness in every sufficiently large dimension. It does not merely assert the existence of one complicated polytope or of witnesses along an unspecified subsequence.
The logarithm is natural. Replacing it by another fixed base greater than one changes the admissible constant c, while preserving the shape of the theorem. The statement leaves both c and D unspecified. The Bárány–Pór paper, Theorem 1.1, supplies the asymptotic context for the book's formulation.
For any fixed c greater than one, the expression c^(d ln d) eventually exceeds A^d for every fixed A greater than one. Thus a single exponential base cannot bound the facet counts of all binary polytopes across dimensions. The binary-coordinate restriction alone does not yield that kind of uniform bound.
The underlying mathematical result is established in the literature. The formalization task is to prove its stated existence conclusion with the precise geometric definitions above. The conclusion concerns actual facets of actual polytopes, so a family of redundant inequalities or a list with repetitions would not satisfy the counting requirement.
A large supply of binary points does not by itself identify the supporting hyperplanes of their convex hull. Counting points and counting facets are different tasks. Moreover, the theorem requires the same growth constant across all sufficiently large dimensions. Verifying individual examples, even examples with many facets, leaves that uniform asymptotic requirement unresolved.
The book describes certain random polytopes as witnesses. Its stated conclusion gives no probability distribution or numerical probability bound. The target records the resulting extremal existence assertion; it imposes no extra probabilistic hypothesis on the witness.
The coordinate model is Fin d → ℝ. IsDPolytope requires a nonempty finite convex hull with full affine span. IsZeroOnePolytope specifies the hull of binary-coordinate points. faceCount P k counts nonempty exposed faces with affine dimension k. These concrete notions also apply to other questions about finite-dimensional polytopes.
Face counts use natural cardinality. On the domain of finite polytopes there are only finitely many faces, so this is the ordinary finite count. Requiring D at least two ensures that d−1 represents the facet dimension without a low-dimensional subtraction convention and that the logarithm's argument is positive. Full dimension excludes the whole polytope from the facet count, and nonemptiness excludes the empty face.
A convex combination expresses a point as a weighted average of other points. Such representations connect geometric sets with finite lists of coordinates and coefficients. A point in a convex hull may initially be described using many generators, even when the ambient space has small dimension. Carathéodory's theorem gives a bound depending only on that dimension. Grünbaum presents this result as a basic theorem of convexity in Convex Polytopes, §2.3, Theorem 5, printed page 15 (the supplied source PDF, page 33).
Fix a natural number . Real -space consists of vectors with real coordinates. Let be any subset of this space. A convex combination of points of is a finite sum of those points multiplied by nonnegative real coefficients, with the coefficients adding to one. The convex hull is the set of all such combinations; equivalently, it is the smallest convex set containing .
The generating set can be finite or infinite. It need not be closed, bounded, or convex. Membership of a point in its convex hull is the sole geometric assumption. The theorem concerns exact equality of vectors, so the representation is not an approximation or a limiting expression.
For every , find points and real numbers satisfying
This is the complete target of Theorem 2.3.5. The witnesses may depend on , , and . No single selection of points must work for every point of the hull. Repeated points are permitted, and some coefficients may be zero. Thus the fixed list of positions can express a combination that uses fewer distinct points.
The result gives a finite certificate for membership in a convex hull whose length is controlled by dimension rather than by the size of the generating set. In the plane it gives three positions, and in three-dimensional space it gives four. Infinite generating sets remain within the same statement: once a particular hull point is chosen, only finitely many generators are needed for its certificate.
The mission asks for a formal proof of this known mathematical result in its fixed-length formulation. Its content is the existence of points and normalized coefficients together. Supplying coefficients without ensuring that their points lie in , or supplying a finite representation without the dimension bound, would leave part of the target unestablished.
The definition of the convex hull guarantees a finite representation but does not directly fix its length at . The dimension bound must hold while preserving nonnegativity, normalization, membership in the original generating set, and exact equality with the selected point. None of these constraints can be replaced by a condition on nearby points or on the closure of .
The coordinate model is Fin d → ℝ. Both the point list and the coefficient list are indexed by Fin (d + 1), corresponding to the source indices from zero through . The convex hull, finite sums, and real scalar multiplication have their standard mathematical meanings. No additional geometric definition is needed.
Dimension zero has one position in the representation. If is empty, its convex hull is empty, so there is no hull point to represent. A singleton generating set and sets contained in a proper affine subspace are allowed. Neither strict positivity of the weights nor affine independence of the selected points is required. The target does not assert uniqueness of a representation.
Branko Grünbaum, Convex Polytopes, second edition, Springer, 2003, Chapter 2, §2.3, Theorem 5, printed page 15; supplied source.pdf, page 33.
A polytope can be represented by cutting a higher-dimensional simplex with an affine flat. The geometry of the cut determines the resulting polytope. Perles's prescribed-point theorem adds a constraint to this representation: the cutting flat must pass through a point chosen in advance inside the simplex. Theorem 5.1.10 in Branko Grünbaum's Convex Polytopes, second edition (2003), shows that this requirement can always be met when the simplex has enough facets. The source is §5.1, printed page 74, with the section convention introduced on printed page 71.
A polytope is the convex hull of finitely many points. In this mission it is nonempty and full-dimensional in its ambient real coordinate space. Thus a d-polytope P in R^d contains enough points to span that space affinely. A face is the set of points on which a supporting linear functional attains its maximum; a facet is a face of dimension d−1. Different supporting functionals can describe the same face, so the facet allowance counts geometric faces rather than inequalities used to describe them.
A k-simplex is the convex hull of k+1 affinely independent points. Write its vertices as v₀,…,vₖ in R^k and its convex hull as T. Affine independence means there is no nontrivial affine relation among these vertices. The simplex therefore has dimension k, and its interior is taken in the whole space R^k. No restriction is placed on its side lengths or angles. A d-flat is a translate of a d-dimensional linear subspace. A section of T by such a flat L is the entire intersection T ∩ L.
The target is the following existence assertion. For every d-polytope P with at most k+1 facets, every k-simplex T in R^k, and every p in the interior of T, there is a d-flat L such that
An affine equivalence here preserves affine combinations and is invertible on the affine spans. It can change lengths and angles. In the stated coordinates, it is represented by an injective real affine map A from R^d onto L satisfying
All input data, including the interior point p, are universally quantified before the flat and map are chosen. The flat can depend on these data. The source writes the facet allowance as f and the simplex dimension as f−1; the notation here uses f=k+1.
The prescribed-point conclusion gives control beyond the existence of some simplex section representing P. It permits the simplex and an interior point to be fixed while the cutting flat is selected to recover the given polytope. The conclusion concerns its full affine geometry: the image of every point of P lies in the section, and every point of the section belongs to that image. This is the additional representational constraint discussed immediately before Theorem 10 on printed page 74.
The difficulty lies in meeting these requirements simultaneously. A flat through the chosen point can have an intersection of the wrong affine shape. A section with the correct affine shape need not contain the prescribed point. The theorem requires both properties for every allowed simplex and point, while retaining the original facet bound. Neither a special regular simplex nor a selected interior point captures this quantifier structure.
The formal statement uses real coordinate spaces indexed by finite sets, finite convex hulls, affine spans, exposed faces, and affine maps. Nonempty faces are counted by their affine dimension. In dimension zero the facet allowance is automatic: under the convention counting the empty face it contributes at most one facet, which fits the allowance k+1. This avoids interpreting natural subtraction d−1 as a facet dimension at d=0. The zero-dimensional simplex and zero-dimensional polytope remain within the statement.
This is a known geometric theorem whose formal proof remains to be supplied. The target contains one proof placeholder. The definitions specify the underlying geometry concretely. The mission consists of the single prescribed-section theorem; its definitions support that statement, and no separate supporting theorem is included as a milestone. The source is Grünbaum, Convex Polytopes, second edition, Springer, 2003, §5.1, Theorem 10, printed page 74 (source.pdf page 100); see also §5.1, printed page 71 (PDF page 97), for the section convention.
A convex polytope has geometric coordinates and a combinatorial structure: its faces, ordered by inclusion. Different coordinates can describe the same combinatorial type. Gale transforms connect these viewpoints by encoding affine dependencies among the vertices as another point configuration, often in a smaller-dimensional space.
Let P and Q be full-dimensional d-polytopes in real coordinate space, each with n vertices. Choose injective lists V and W containing all their vertices and fix a permutation θ of the n indices. The question is whether this prescribed vertex correspondence extends to an isomorphism between the full face posets of P and Q. The empty face and the whole polytope are included in these posets.
An affine dependency of V is a list of real coefficients a whose sum is zero and whose weighted sum of the vertices is zero. These coefficient lists form a vector space of dimension n−d−1. Choose any basis and put its vectors into the columns of a matrix. The n rows of this matrix form a Gale transform G of V. A Gale transform H of W is constructed in the same way, with its own choice of basis. The rows may repeat and may be zero; there is no general-position assumption.
Theorem 5.4.5 states that the prescribed correspondence extends to an isomorphism of face posets exactly when it preserves every relative-interior test on the Gale configurations. For every subset J of vertex indices, the origin belongs to the relative interior of the convex hull of the rows G(J) if and only if it belongs to the relative interior of the convex hull of H(θ(J)). Relative interior means interior within the affine span of the selected convex hull, not interior in the whole coordinate space.
The equivalence concerns every subset of indices and both directions of implication. It retains the prescribed permutation rather than merely asking whether the polytopes have some combinatorial equivalence. Indexing the Gale points also keeps track of which original vertices they represent when several Gale rows coincide. Taking a set image does not change the convex hull of a chosen subconfiguration.
The empty subset is included: its convex hull has empty relative interior, so both membership tests are false. For a simplex, n=d+1 and the Gale space has dimension zero. Its nonempty row subconfigurations consist of the zero vector, and the relative-interior formulation still applies. Zero-dimensional polytopes are also retained.
This mission is the equivalence criterion itself, with the concrete definitions of a full-dimensional polytope, its face poset, and a Gale transform. The grouping keeps transform construction within the vocabulary of the criterion. The neighboring results on affine and projective equivalence have different conclusions and are outside this goal.
Source: Branko Grünbaum, Convex Polytopes, second edition (2003), §5.4, Theorem 5, printed page 89 (PDF page 115). The affine-dependence construction appears on printed pages 85–86 (PDF pages 111–112).