General-profile Equation (5.3) rate bookkeeping
Provedmme_stothers_general_profile_rate_below_marginal_multinomialEquation (5.3) rate bookkeeping at an arbitrary integral ten-class profile.
Fix an integral witness for the Davie--Stothers ten Table-1 symmetry classes: a vector
of strictly positive natural numbers (base), its weighted total
the induced normalized profile , and the nine integral marginal counts obtained by applying the Equation (5.2) matrix to the unnormalized . Because is linear, is exactly the paper's marginal , and summing the nine rows of gives , so is the address length at scale .
Then there is a constant such that for all large ,
In words: at any integral profile, the nine-letter marginal multinomial carries the entire
scalar part of the global rate of Equation (5.3) that is not already accounted for by the ten
Table-1 constituent values , up to a single uniform loss which is
irrelevant in the laser limit. The two sides are an identity up to Stirling factors: the
entropy factor of globalRate is exactly the exponential growth rate of the
multinomial coefficient, and the pair inside globalRate cancels on the
diagonal .
This is the profile-parametric form of the published fixed-witness statement
mme_stothers_fixed_profile_rate_below_marginal_multinomial, which is recovered verbatim by
taking ,
and
.
No numerical property of that witness is used: the proof needs only positivity of the ten
counts, and it is therefore reusable at every profile appearing in the general Theorem 5.3
optimization (and at the different profiles of the DWZ and More-Asymmetry fourth-power tables).
import Mathlib.Analysis.SpecialFunctions.Exp import Mathlib.Analysis.SpecialFunctions.Sqrt import Mathlib.Data.Nat.Choose.Multinomial import Definitions.Def_mme_stothers_fourth_data open MME BigOperators Filter set_option autoImplicit false
theorem mme_stothers_general_profile_rate_below_marginal_multinomial
(tau : ℝ) (a : Fin 10 → ℝ) (base : Fin 10 → ℕ) (marg : Fin 9 → ℕ) (D : ℕ)
(hbase : ∀ r, 0 < base r)
(hD : D = ∑ r : Fin 10, MME.StothersFourth.classMultiplicity r * base r)
(ha : ∀ i, a i = (base i : ℝ) / (D : ℝ))
(hmarg : ∀ j, (marg j : ℝ) =
MME.StothersFourth.Q (fun i ↦ (base i : ℝ)) j) :
∃ C : ℝ, 0 ≤ C ∧
∀ᶠ m : ℕ in Filter.atTop,
(MME.StothersFourth.globalRate 6 tau a a) ^ (3 * D * m) *
Real.exp (-C * Real.sqrt (((3 * D * m + 1 : ℕ) : ℝ))) ≤
(Nat.multinomial Finset.univ (fun j : Fin 9 ↦ marg j * m) : ℝ) *
(∏ r : Fin 10,
(MME.StothersFourth.classValue 6 tau r) ^
(MME.StothersFourth.classMultiplicity r * (base r * m))) := by
sorry