Sample Complexity Bounds for Linear Constrained MDPs with a Generative Model
Xingtu Liu, Lin F. Yang, Sharan Vaswani
TL;DR
This work addresses sample-efficient learning for constrained MDPs with linear function approximation under a generative model. It develops a model-based primal-dual framework that reduces CMDPs to unconstrained linear MDPs via a three-oracle setup (data collection, MDP solver, policy evaluation) and proves a meta-theorem bounding overall sample complexity. The authors instantiate the framework with Least-Squares MDVI (LS-MDVI) and Least-Squares Policy Evaluation (LS-PE), deriving near-optimal bounds: $\tilde{O}\left(\frac{d^2}{(1-\gamma)^4\varepsilon^2}\right)$ for relaxed feasibility and $\tilde{O}\left(\frac{d^2}{(1-\gamma)^6\varepsilon^2\zeta^2}\right)$ for strict feasibility, plus a matching lower bound up to horizon factors. They also show that tabular CMDPs inherit near-optimal complexities, and discuss the framework’s flexibility (e.g., GSS), variance-weighted LS limitations, online/offline comparisons, and implications for constrained RL in large or continuous spaces.
Abstract
We consider infinite-horizon $γ$-discounted (linear) constrained Markov decision processes (CMDPs) where the objective is to find a policy that maximizes the expected cumulative reward subject to expected cumulative constraints. Given access to a generative model, we propose to solve CMDPs with a primal-dual framework that can leverage any black-box unconstrained MDP solver. For linear CMDPs with feature dimension $d$, we instantiate the framework by using mirror descent value iteration (\texttt{MDVI})~\citep{kitamura2023regularization} an example MDP solver. We provide sample complexity bounds for the resulting CMDP algorithm in two cases: (i) relaxed feasibility, where small constraint violations are allowed, and (ii) strict feasibility, where the output policy is required to exactly satisfy the constraint. For (i), we prove that the algorithm can return an $ε$-optimal policy with high probability by using $\tilde{O}\left(\frac{d^2}{(1-γ)^4ε^2}\right)$ samples. For (ii), we show that the algorithm requires $\tilde{O}\left(\frac{d^2}{(1-γ)^6ε^2ζ^2}\right)$ samples, where $ζ$ is the problem-dependent Slater constant that characterizes the size of the feasible region. Furthermore, we prove a lower-bound of $Ω\left(\frac{d^2}{(1-γ)^5ε^2ζ^2}\right)$ for the strict feasibility setting. We note that our upper bounds under both settings exhibit a near-optimal dependence on $d$, $ε$, and $ζ$. Finally, we instantiate our framework for tabular CMDPs and show that it can be used to recover near-optimal sample complexities in this setting.
