Prove2Me
Navigate
MissionsFormalpediaBlogsUsersMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

One-step concentration of the configuration (Lem. 9, Gast et al. 2024)

Proved
MarkovEntanglement.rmab_one_step_concentration

by tianyipeng · Aug 28, 2026 · Mathlib c5ea003 (Lean v4.30.0)

Consider one step of the NNN-agent restless-bandit chain under an index policy, started from an arbitrary joint state sss, and let φ\varphiφ be the explicit mean-field map of the configuration, taken at the exact activation fraction ⌊αN⌋/N\lfloor \alpha N \rfloor / N⌊αN⌋/N of the NNN-agent system. Writing ϵ[1]=m[1]−φ(m[0])\epsilon[1] = m[1] - \varphi(m[0])ϵ[1]=m[1]−φ(m[0]) for the one-step deviation of the configuration from its mean-field image,

E[∥ϵ[1]∥1  ∣  m[0]]  ≤  ∣S∣N.\mathbb{E}\big[\|\epsilon[1]\|_1 \;\big|\; m[0]\big] \;\le\; \sqrt{\frac{|S|}{N}}.E[∥ϵ[1]∥1​​m[0]]≤N∣S∣​​.

The mechanism is the usual 1/N1/\sqrt{N}1/N​ of an average of NNN indicators. Conditionally on the current joint state, the index policy activates a deterministic number of agents in each local state — only which agents is random, and the agents in a state are exchangeable — so the next local states are an independent family: for each local state xxx, a fixed number of draws from P1(x,⋅)P_1(x, \cdot)P1​(x,⋅) and a fixed number from P0(x,⋅)P_0(x,\cdot)P0​(x,⋅). The next configuration is their average, its mean is exactly φ(m[0])\varphi(m[0])φ(m[0]) by definition of the mean-field map, and each coordinate has variance at most 1/N1/N1/N times its mean. Bounding each E∣ϵx∣\mathbb{E}|\epsilon_x|E∣ϵx​∣ by the standard deviation and then summing over xxx with Cauchy–Schwarz produces ∣S∣/N\sqrt{|S|/N}∣S∣/N​.

This is the input that fixes the rate in Theorem 7: everything downstream amplifies it by constants that do not depend on NNN.

Preamble
import Mathlib
import Definitions.Def_markov_entanglement_meanfield

open scoped BigOperators
open MarkovEntanglement
Formal statement
namespace MarkovEntanglement

variable {S : Type*} [Fintype S] [DecidableEq S]

/-- Lemma 9 (One-step Concentration; Lemma 1 of Gast, Gaujal and Yan 2024).  Starting from any
joint state, the configuration after one step deviates from its mean-field image by at most
`√(|S| / N)` in expected `ℓ¹` norm:

`E[‖m[1] − φ(m[0])‖₁ ∣ m[0]] ≤ √(|S| / N)`.

The expectation is over one step of the `N`-agent chain, the mean-field map is the explicit
one at the exact activation fraction `⌊αN⌋ / N` of the `N`-agent system, and the rate is the
usual `1/√N` of an average of `N` conditionally independent indicators. -/
theorem rmab_one_step_concentration
    (P0 P1 : Matrix S S ℝ) (hP0 : IsTransitionMatrix P0) (hP1 : IsTransitionMatrix P1)
    (ν : S → ℝ) (α : ℝ) (N : ℕ) (hN : 0 < N)
    (π : (Fin N → S) → (Fin N → Bool) → ℝ)
    (hπ : IsIndexPolicy ν ⌊α * (N : ℝ)⌋₊ π)
    (s : Fin N → S) :
    ∑ s' : Fin N → S, rmabStep P0 P1 π s s' *
        l1Norm (fun x => configuration s' x
          - meanFieldMap P0 P1 ν ((⌊α * (N : ℝ)⌋₊ : ℝ) / (N : ℝ)) (configuration s) x)
      ≤ Real.sqrt ((Fintype.card S : ℝ) / (N : ℝ)) := by
  sorry

/-! ### M5 — Lemma 10, multi-step concentration (Gast, Gaujal and Yan) -/

end MarkovEntanglement
Source
Nicolas Gast, Bruno Gaujal and Chen Yan, 'Reoptimization nearly solves weakly coupled Markov decision processes' (2024), Lemma 1; cited as Lemma 9 in Chen and Peng, arXiv:2506.02385v3, Appendix I, p. 43
Read-back

What the Lean code literally says, in plain math · claude-opus-5

\newcommand{\Fin}{\mathrm{Fin}}For every type SSS carrying a finiteness structure and decidable equality (both implicit, so SSS ranges over all finite types, the empty one included), for all matrices P0,P1∈RS×SP_0,P_1 \in \mathbb{R}^{S\times S}P0​,P1​∈RS×S, every function ν:S→R\nu : S \to \mathbb{R}ν:S→R, every real number α\alphaα, every natural number NNN, every function π:S{0,…,N−1}→{false,true}{0,…,N−1}→R\pi : S^{\{0,\dots,N-1\}} \to \{\mathrm{false},\mathrm{true}\}^{\{0,\dots,N-1\}} \to \mathbb{R}π:S{0,…,N−1}→{false,true}{0,…,N−1}→R, and every tuple s=(s0,…,sN−1)∈SNs = (s_0,\dots,s_{N-1}) \in S^Ns=(s0​,…,sN−1​)∈SN, the statement asserts the displayed inequality below under the following hypotheses. Write M:=⌊αN⌋NM := \lfloor \alpha N \rfloor_{\mathbb{N}}M:=⌊αN⌋N​ for the natural-number floor of the real number α⋅N\alpha \cdot Nα⋅N — this is the greatest natural number ≤αN\le \alpha N≤αN, and it is 000 whenever αN<0\alpha N < 0αN<0; nothing in the statement constrains α\alphaα to lie in [0,1][0,1][0,1], so MMM may be 000 or may exceed NNN. For a tuple u∈SNu \in S^Nu∈SN and x∈Sx \in Sx∈S put nu(x):=#{i:ui=x}n_u(x) := \#\{i : u_i = x\}nu​(x):=#{i:ui​=x} (a natural number) and cu(x):=nu(x)/Nc_u(x) := n_u(x)/Ncu​(x):=nu​(x)/N (a real quotient; N≠0N \neq 0N=0 is guaranteed by hypothesis (c) below). The hypotheses are:

  • (a) P0P_0P0​ has all entries ≥0\ge 0≥0 and each of its rows sums to 111: ∑j∈SP0(i,j)=1\sum_{j \in S} P_0(i,j) = 1∑j∈S​P0​(i,j)=1 for every iii.
  • (b) P1P_1P1​ likewise has all entries ≥0\ge 0≥0 and every row summing to 111.
  • (c) 0<N0 < N0<N.
  • (d) π\piπ is an index policy for the priority function ν\nuν at budget MMM, meaning all three of:
    • π(u,a)≥0\pi(u,a) \ge 0π(u,a)≥0 for all u∈SNu \in S^Nu∈SN and all action profiles a∈{false,true}Na \in \{\mathrm{false},\mathrm{true}\}^Na∈{false,true}N, and ∑a∈{false,true}Nπ(u,a)=1\sum_{a \in \{\mathrm{false},\mathrm{true}\}^N} \pi(u,a) = 1∑a∈{false,true}N​π(u,a)=1 for every uuu;
    • for all uuu and aaa, if π(u,a)≠0\pi(u,a) \neq 0π(u,a)=0 then aaa is exactly budgeted: #{i:ai=true}=M\#\{i : a_i = \mathrm{true}\} = M#{i:ai​=true}=M;
    • for every u∈SNu \in S^Nu∈SN and every index iii, the marginal probability that agent iii is activated, ∑a∈{false,true}N[ ai=true ] π(u,a)\sum_{a \in \{\mathrm{false},\mathrm{true}\}^N} [\,a_i = \mathrm{true}\,]\,\pi(u,a)∑a∈{false,true}N​[ai​=true]π(u,a), equals qνM(u,ui)q_\nu^M(u, u_i)qνM​(u,ui​), where for x∈Sx \in Sx∈S
qνM(u,x):={0if nu(x)=0,min⁡ ⁣(nu(x),  M−˙hν(u,x))nu(x)otherwise,hν(u,x):=∑y : ν(x)<ν(y)nu(y),q_\nu^M(u,x) := \begin{cases} 0 & \text{if } n_u(x) = 0,\\[4pt] \dfrac{\min\!\big(n_u(x),\; M \dot- h_\nu(u,x)\big)}{n_u(x)} & \text{otherwise,}\end{cases} \qquad h_\nu(u,x) := \sum_{y \,:\, \nu(x) < \nu(y)} n_u(y),qνM​(u,x):=⎩⎨⎧​0nu​(x)min(nu​(x),M−˙​hν​(u,x))​​if nu​(x)=0,otherwise,​hν​(u,x):=y:ν(x)<ν(y)∑​nu​(y),
with $\dot-$ denoting **truncated natural-number subtraction** (so $M \dot- h_\nu(u,x) = 0$ whenever $h_\nu(u,x) \ge M$), the numerator being cast to $\mathbb{R}$ after the truncated subtraction and the minimum, and the comparison $\nu(x) < \nu(y)$ being strict (ties in $\nu$ contribute nothing to $h_\nu$; $\nu$ is not assumed injective).

Under these hypotheses the conclusion is the single inequality

∑s′∈SN(∑a∈{false,true}Nπ(s,a)∏j=0N−1Kaj ⁣(sj,sj′))⋅(∑x∈S∣ cs′(x)  −  Φ(cs)(x) ∣)  ≤  ∣S∣N,\sum_{s' \in S^N} \left(\sum_{a \in \{\mathrm{false},\mathrm{true}\}^N} \pi(s,a)\prod_{j=0}^{N-1} K_{a_j}\!\left(s_j, s'_j\right)\right)\cdot\left(\sum_{x \in S}\Big|\,c_{s'}(x) \;-\; \Phi\big(c_s\big)(x)\,\Big|\right) \;\le\; \sqrt{\frac{|S|}{N}},s′∈SN∑​​a∈{false,true}N∑​π(s,a)j=0∏N−1​Kaj​​(sj​,sj′​)​⋅(x∈S∑​​cs′​(x)−Φ(cs​)(x)​)≤N∣S∣​​,

where Ktrue:=P1K_{\mathrm{true}} := P_1Ktrue​:=P1​ and Kfalse:=P0K_{\mathrm{false}} := P_0Kfalse​:=P0​; the outer sum runs over all ∣S∣N|S|^N∣S∣N tuples s′s's′; ∣S∣|S|∣S∣ is the cardinality of SSS; ⋅\sqrt{\cdot}⋅​ is the total real square root; and Φ\PhiΦ is the mean-field map at activation fraction α′:=M/N=⌊αN⌋N/N\alpha' := M/N = \lfloor \alpha N\rfloor_{\mathbb{N}}/Nα′:=M/N=⌊αN⌋N​/N (a nonnegative real, not α\alphaα itself), defined for an arbitrary m:S→Rm : S \to \mathbb{R}m:S→R by

Φ(m)(y):=∑x∈S[(m(x)−f(m,x)) P0(x,y)  +  f(m,x) P1(x,y)],f(m,x):=min⁡ ⁣(m(x),  max⁡(0,  α′−∑y : ν(x)<ν(y)m(y))).\Phi(m)(y) := \sum_{x \in S}\Big[\big(m(x) - f(m,x)\big)\,P_0(x,y) \;+\; f(m,x)\,P_1(x,y)\Big],\qquad f(m,x) := \min\!\Big(m(x),\; \max\big(0,\; \alpha' - \textstyle\sum_{y \,:\, \nu(x)<\nu(y)} m(y)\big)\Big).Φ(m)(y):=x∈S∑​[(m(x)−f(m,x))P0​(x,y)+f(m,x)P1​(x,y)],f(m,x):=min(m(x),max(0,α′−∑y:ν(x)<ν(y)​m(y))).

Note that Φ\PhiΦ is applied only to the particular argument csc_scs​, the empirical configuration of the given starting tuple sss; the left-hand side is a single weighted average of ℓ1\ell^1ℓ1 distances taken over one step from that fixed sss, not a supremum over starting tuples, not an iterate, and not a statement about any stationary or limiting object. The quantities in the bracketed weight are exactly the sum over joint action profiles of π(s,a)\pi(s,a)π(s,a) times the product over the NNN agents of the corresponding one-agent transition probability; the statement does not separately assume this weight is nonnegative or sums to 111 over s′s's′ — that is left to follow, or not, from (a), (b), (d). Degenerate instances silently included: ∣S∣=1|S| = 1∣S∣=1; SSS empty, in which case for N≥1N \ge 1N≥1 there is no tuple sss at all and the claim is vacuous for that SSS; M=0M = 0M=0 (which forces π(u,a)=0\pi(u,a) = 0π(u,a)=0 unless aaa is all-false\mathrm{false}false, and makes qνM≡0q_\nu^M \equiv 0qνM​≡0 and f≡0f \equiv 0f≡0); M>NM > NM>N, in which case no action profile can satisfy #{i:ai=true}=M\#\{i : a_i = \mathrm{true}\} = M#{i:ai​=true}=M, so the second and first clauses of (d) are jointly unsatisfiable and hypothesis (d) is impossible to meet for such α,N\alpha, Nα,N; and α<0\alpha < 0α<0, which collapses to the M=0M = 0M=0 case. The inequality is non-strict (≤\le≤), and both the truncated M−˙hν(u,x)M \dot- h_\nu(u,x)M−˙​hν​(u,x) in the discrete activation probability and the max⁡(0,⋅)\max(0,\cdot)max(0,⋅) in the continuum f(m,x)f(m,x)f(m,x) floor their arguments at zero.

Human review
  • Endorsed by Shuze Chen · Aug 30, 2026

  • Endorsed by tianyipeng · Aug 30, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactJoin Slack© 2026 Prove2Me