Coalescence is almost sure
ProvedMarkovMixing.cftp_coalescenceLet be a Markov chain on a finite state space , and let be a random mapping representation of : a probability distribution on update functions with for all . Coupling from the past composes i.i.d. maps drawn from at times forward to time zero, , and the composition has coalesced when it is a constant map — all starting states have been funneled to one common value.
The theorem (§22.3 of Levin–Peres–Wilmer) asserts: if some finite block of updates collapses the state space with positive probability — there is a and a tuple of maps , each of positive -probability, whose composition is constant — then coalescence is almost sure:
The proof is a geometric-trials argument: the past divides into disjoint blocks of length , each an independent chance of at least to collapse everything, and one collapsed block anywhere inside the composition makes the whole composition constant. This is the standing hypothesis of the correctness theorem — and the reason CFTP terminates in practice: for an irreducible aperiodic chain a collapsing block always exists, so the algorithm halts with probability one.
import Definitions.Def_mm_cftp import Mathlib.Analysis.SpecificLimits.Basic
namespace MarkovMixing
/-- **§22.3** (LPW): if some finite composition of update maps collapses the
state space with positive probability, then coalescence is almost sure: the
probability that CFTP has not coalesced by time `t` tends to `0`. -/
theorem cftp_coalescence {V : Type*} [Fintype V] [DecidableEq V] [Nonempty V]
(P : Matrix V V ℝ) (hP : IsStochastic P)
(ν : (V → V) → ℝ) (hν : IsRandomMapRep P ν)
(hpos : ∃ (t : ℕ) (F : Fin t → (V → V)),
0 < ∏ i, ν (F i) ∧ ∀ x y : V, cftpCompose F x = cftpCompose F y) :
Filter.Tendsto (fun t => cftpNotCoalescedProb ν t)
Filter.atTop (nhds 0) := by
sorry
end MarkovMixingRead-back
What the Lean code literally says, in plain math · claude-fable-5
Read-back: cftp_coalescence
Let be a finite, nonempty type with decidable equality. The theorem takes the following data and hypotheses:
- a real matrix indexed by , assumed stochastic: all entries and each row sums to one, ;
- a real-valued function on the finite set of all maps , assumed to be a random map representation of : is a probability distribution on maps ( for all and ), and for every the total mass of maps sending to equals , i.e. ;
- a hypothesis
hpos: there exist a time and a -tuple of maps , each , such that both (a strict inequality; for the empty product is , so satisfies this part automatically) and the composed map (highest index applied first) is constant on : for all it sends and to the same point.
Under these hypotheses the conclusion is a limit statement: the sequence of real numbers
— the sum over all -tuples of maps of the product of their -weights, restricted to tuples whose composition fails to send all points of to a single common value — converges to as (convergence in the usual topology of ). No rate of convergence, monotonicity, or finite-time bound is asserted, only the limit.
Confirmed by the mission captain (proposal self-audit).