Constrained Markov Decision Process (CMDP)

• CMDP is a generalization of MDP that incorporates side constraints (e.g., safety, resource, fairness) to optimize cumulative rewards while enforcing state/action limits.
• It employs a specialized backward induction and LP-based formulation that computes randomized, nonstationary policies ensuring constraint satisfaction at every decision step.
• Simulation results in multi-agent settings demonstrate that CMDP policies maintain strict feasibility with near-optimal rewards compared to unconstrained MDP solutions.

Constrained Markov Decision Process (CMDP) is a generalization of the Markov Decision Process (MDP) framework that introduces side constraints—typically state and/or action constraints—which must be satisfied in addition to optimizing a cumulative objective. CMDPs are used to formalize safety, resource, or fairness requirements in sequential stochastic control and reinforcement learning. In contrast to unconstrained MDPs, the presence of constraints—often expressed as expectations over trajectories or as linear inequalities involving state occupancies—necessitates sophisticated mathematical and algorithmic approaches for policy synthesis and analysis.

1. Mathematical Formulation of CMDPs

A finite-horizon CMDP is specified by a tuple (𝒮, 𝒜, P, r, B, d, N) where:

• 𝒮: finite set of states
• 𝒜: finite set of actions
• P: time-dependent state transition matrices {Pₜ}
• rₜ: stagewise reward vectors
• B: state constraint matrix (e.g., capacity or safety constraint coefficients)
• d: constraint bound vector
• N: planning horizon

At each time t, the state constraint is given by

$B x_t \leq d$

where $x_t$ is the (possibly random) state occupancy vector at time t induced by the policy and system dynamics.

The canonical CMDP objective is

$$\max_{\pi} \mathbb{E}_\pi\left[\sum_{t=1}^N r_t(s_t, a_t)\right], \quad \text{subject to } B x_t \le d,\,\forall t$$

Policy $\pi$ maps histories up to t to distributions over actions; feasibility is enforced at each time along the system's evolution.

2. Randomization, Nonstationarity, and the Convex Set of Policies

Unlike unconstrained MDPs, where deterministic stationary policies are optimal, CMDPs with state or action constraints may require randomization and nonstationarity to satisfy hard constraints for all possible realizations of the system trajectory. For finite-horizon CMDPs with state constraints $B x_t \leq d$, the optimal policies must be computed over the convex set:

$$\mathcal{C} = \{ Q \in \mathbb{R}^{n \times p} : Q 1 = 1,\, Q \geq 0,\, M(Q) x \le d,\ \forall x \}$$

Here, $Q$ is a stochastic decision matrix specifying action selection, and $M(Q)$ is the transition operator induced by $Q$ (see precise definitions in (Chamie et al., 2015), eqn. set for $\mathcal{C}(x)$). This convexification enables use of LP-based methods, as the (generally nonconvex) original set of feasible Markov randomized policies is intractable.

3. Backward Induction and LP-Based Synthesis

Solving a finite-horizon CMDP with hard state constraints requires a nonstandard dynamic programming (DP) approach, since the value function does not admit a closed-form as in the unconstrained case. The main method [(Chamie et al., 2015), Algorithm 3] consists of the following backward recursion at each time t:

• Given terminal value $x_t$0.
• For t = N-1 down to 1:

$x_t$1

where $x_t$2 is the set of admissible state distributions (e.g., $x_t$3).

At each stage, the policy is obtained by solving a max–min optimization: the inner minimization finds the worst-case performance over all possible current state vectors, reflecting the “hedging” required for constraint satisfaction, while the outer maximization selects the best feasible randomized control.

4. Linear Programming Duality and Policy Computation

The computational core of this approach is reformulating the inner min–max as a (primal–dual) linear program. The minimization over $x_t$4 with affine objective and polyhedral constraints can be cast into a dual maximizing over $x_t$5 and a free variable $x_t$6:

$x_t$7

All additional constraints (such as $x_t$8, normalizations, etc.) are reduced to linear or convex constraints in the decision variables. This enables efficient computation of the nonstationary, randomized policies at scale—an essential step for practical CMDP deployments.

5. Projection Heuristic and Relation to Unconstrained MDPs

Since the optimal unconstrained MDP policy typically violates state constraints, the paper introduces a projection-based heuristic: among all LP-generated $x_t$9 in the feasible set, select the one (in a norm such as Frobenius or $$\max_{\pi} \mathbb{E}_\pi\left[\sum_{t=1}^N r_t(s_t, a_t)\right], \quad \text{subject to } B x_t \le d,\,\forall t$$0) closest to the unconstrained ($$\max_{\pi} \mathbb{E}_\pi\left[\sum_{t=1}^N r_t(s_t, a_t)\right], \quad \text{subject to } B x_t \le d,\,\forall t$$1) deterministic policy. This ensures that if $$\max_{\pi} \mathbb{E}_\pi\left[\sum_{t=1}^N r_t(s_t, a_t)\right], \quad \text{subject to } B x_t \le d,\,\forall t$$2 is feasible, it will be recovered; otherwise, the policy with maximal proximity and minimal loss in reward is used, while maintaining the guarantee $$\max_{\pi} \mathbb{E}_\pi\left[\sum_{t=1}^N r_t(s_t, a_t)\right], \quad \text{subject to } B x_t \le d,\,\forall t$$3.

6. Simulation Results and Empirical Findings

A key illustration is a multi-agent swarm navigation problem. Each agent transitions on a $$\max_{\pi} \mathbb{E}_\pi\left[\sum_{t=1}^N r_t(s_t, a_t)\right], \quad \text{subject to } B x_t \le d,\,\forall t$$4 grid with stochastic actions; a per-bin density constraint $$\max_{\pi} \mathbb{E}_\pi\left[\sum_{t=1}^N r_t(s_t, a_t)\right], \quad \text{subject to } B x_t \le d,\,\forall t$$5 enforces capacity/safety. The naive unconstrained MDP solution leads to over-concentration in high-reward bins (i.e., constraint violation). In contrast, the CMDP policy synthesized via LP backward induction always satisfies bin capacities, and the total expected reward is provably no less than $$\max_{\pi} \mathbb{E}_\pi\left[\sum_{t=1}^N r_t(s_t, a_t)\right], \quad \text{subject to } B x_t \le d,\,\forall t$$6. Empirical results indicate that the projected policy typically achieves reward levels close to the unconstrained optimum, but with strict feasibility.

Policy Type	Reward Achieved	Constraint Satisfaction?
Unconstrained MDP	Highest	Possibly violated
CMDP Synthesized	Slightly lower	Always satisfied (all $$\max_{\pi} \mathbb{E}_\pi\left[\sum_{t=1}^N r_t(s_t, a_t)\right], \quad \text{subject to } B x_t \le d,\,\forall t$$7)
Projected CMDP	Close to MDP	Always satisfied (all $$\max_{\pi} \mathbb{E}_\pi\left[\sum_{t=1}^N r_t(s_t, a_t)\right], \quad \text{subject to } B x_t \le d,\,\forall t$$8)

7. Impact and Computational Considerations

This framework is the first to provide an efficient, finite-horizon algorithm with optimality guarantees for CMDPs with state constraints (Chamie et al., 2015). The methodology extends to large-scale, multi-agent, and distributed systems where explicit state constraint satisfaction (e.g., collision avoidance, density regulation) is paramount. The approach is computationally tractable for moderate $$\max_{\pi} \mathbb{E}_\pi\left[\sum_{t=1}^N r_t(s_t, a_t)\right], \quad \text{subject to } B x_t \le d,\,\forall t$$9 and $\pi$0 due to convexity and the reduction to LPs.

It is important to note that the method is independent of the initial state distribution, and all policies can be pre-computed offline for deployment. Furthermore, by recasting the inner minimization as a dual LP, the approach remains practically implementable even when nonstationarity and randomization are necessary for constraint satisfaction.

8. Theoretical Guarantees

The central result is the provision of a computable lower bound $\pi$1 on achievable reward for any initial distribution $\pi$2. The constructed policy sequence ensures feasibility at every step, and the results generalize to systems where constraints are central to system safety, reliability, or resource allocation.

This methodology fundamentally extends the scope of MDP optimization to realistic systems in which stringent state-space constraints—reflecting physical limits or safety considerations—cannot be ignored, offering a blueprint for synthesis, analysis, and real-world deployment.