Discounted Bellman solution over compact confidence sets
ProvedBanditAlgorithm.mdp_compact_confidence_discounted_bellman_solutionFix states and actions, a reward function , and for each state-action pair a nonempty compact set of probability vectors on the state space. For every discount factor there are a value function , an action map and transition rows such that
Together these say that solves the discounted Bellman optimality equation of the extended MDP, whose actions are the pairs with , and that the maximum is attained at .
This is the fixed point computed by the extended value iteration of UCRL2: the confidence sets are the balls around the empirical transition rows, and the extended MDP is the one whose transition matrices range over the confidence region. It is not a finite MDP — its action set is a continuum — so the usual finite-action argument does not apply verbatim.
What replaces finiteness is compactness. The inner maximisation is attained because a continuous image of a compact set is a compact set of reals and therefore contains its supremum. Moreover is -Lipschitz for the supremum norm, because every is a probability vector and hence . Consequently the operator is a -contraction of , and Banach's fixed point theorem supplies ; the maximisers over the finite action set and over the compact confidence set supply and . Evaluating the fixed point equation at a state where is largest gives , and at a state where it is smallest gives .
import Mathlib.Topology.MetricSpace.Contracting import Mathlib.Analysis.SpecificLimits.Basic
theorem BanditAlgorithm.mdp_compact_confidence_discounted_bellman_solution
{S A : ℕ} (hS : 0 < S) (hA : 0 < A)
(r : Fin S → Fin A → ℝ) (hr : ∀ s a, r s a ∈ Set.Icc (0 : ℝ) 1)
(C : Fin S → Fin A → Set (Fin S → ℝ)) (hCne : ∀ s a, (C s a).Nonempty)
(hCcomp : ∀ s a, IsCompact (C s a))
(hCprob : ∀ s a, ∀ p ∈ C s a, (∀ s', 0 ≤ p s') ∧ ∑ s', p s' = 1)
(γ : ℝ) (hγ0 : 0 ≤ γ) (hγ1 : γ < 1) :
∃ (V : Fin S → ℝ) (f : Fin S → Fin A) (q : Fin S → Fin S → ℝ),
(∀ s, V s ∈ Set.Icc (0 : ℝ) (1 / (1 - γ))) ∧
(∀ s a, ∀ p ∈ C s a, r s a + γ * ∑ s', p s' * V s' ≤ V s) ∧
(∀ s, q s ∈ C s (f s)) ∧
(∀ s, V s = r s (f s) + γ * ∑ s', q s s' * V s') := by
sorry