Is a Good Representation Sufficient for Sample Efficient Reinforcement Learning?

Simon S. Du, Sham M. Kakade, Ruosong Wang, Lin F. Yang

Introduction

Modern reinforcement learning (RL) problems are often challenging due to the huge state space. To tackle this challenge, function approximation schemes are often employed to provide a compact representation, so that reinforcement learning can generalize across states. A common paradigm is to first use a feature extractor to transform the raw input to features (a succinct representation) and then apply a linear predictor on top of the features. Traditionally, the feature extractor is often handcrafted (Sutton & Barto 2018), while more modern methods often train a deep neural network to extract features. The hope of this paradigm is that, if there exists a good low dimensional (linear) representation, then efficient reinforcement learning is possible.

Empirically, combining various RL function approximation algorithms with neural networks for feature extraction has lead to tremendous successes on various tasks (Mnih et al. 2015; Schulman et al. 2015; Schulman et al. 2017). A major problem, however, is that these methods often require a large amount of samples to learn a good policy. For example, deep QQ-network requires millions of samples to solve certain Atari games (Mnih et al. 2015). Here, one may wonder if there are fundamental statistical limitations on such methods, and, if so, under what conditions it would be possible to efficiently learn a good policy?

In the supervised learning context, it is well-known that empirical risk minimization is a statistically efficient method when using a low-complexity hypothesis space (Shalev-Shwartz & Ben-David 2014), e.g. a hypothesis space with bounded VC dimension. For example, polynomial number of samples suffice for learning a near-optimal dd-dimensional linear classifier, even in the agnostic setting Here we only study the sample complexity and ignore the computational complexity.. In contrast, in the more challenging RL setting, we seek to understand if efficient learning is possible (say from a sample complexity perspective) when we have access to an accurate (and compact) parametric representation — e.g. our policy class contains a near-optimal policy or our hypothesis class accurately approximates the optimal value function. In particular, this work focuses on the following question:

Is a good representation sufficient for sample-efficient reinforcement learning?

This question has largely been studied only with respect to approximation error in the more classical approximate dynamic programming literature, where it is known that algorithms are stable to certain worst-case approximation errors. With regards to sample efficiency, this question is largely unexplored, where the extant body of literature mainly focuses on conditions which are sufficient for efficient reinforcement learning though there is little understanding of what are necessary conditions for efficient reinforcement learning. In reinforcement learning, there is no direct analogue of empirical risk minimization as in the supervised learning context, and it is not evident what are the statistical limits of learning based on properties of our underlying hypothesis class (which may be value-based, policy-based, or model-based).

Many recent works have provided polynomial upper bounds under various sufficient conditions, and in what follows we list a few examples. For value-based learning, the work of Wen & Van Roy 2013 showed that for deterministic systems MDPs where both reward and transition are deterministic., if the optimal QQ-function can be perfectly predicted by linear functions of the given features, then the agent can learn the optimal policy exactly with polynomial number of samples. Recent work (Jiang et al. 2017) further showed that if certain complexity measure called Bellman rank is bounded, then the agent can learn a near-optimal policy efficiently. For policy-based learning, Agarwal et al. 2019 gave polynomial upper bounds which depend on a parameter that measures the difference between the initial distribution and the distribution induced by the optimal policy.

Our Contributions. This paper gives, perhaps surprisingly, strong negative results to this question. The main results are exponential lower bounds in terms of planning horizon HH for value-based, model-based, and policy-based algorithms with given good representations Our results can be easily extend to infinite horizon MDPs with discount factors by replacing the planning horizon HH with 11−γ\frac{1}{1-\gamma}, where γ\gamma is the discount factor. We omit the discussion on discount MDPs for simplicity. . Notably, the requirements on the representation that suffice for sample efficient RL are even more stringent than the more traditional approximation viewpoint. A comprehensive summary of previous upper bounds and our lower bounds is given in Table 1, and here we briefly summarize our hardness results.

For value-based learning, we show even if QQ-functions of all policies can be approximated by linear functions of the given representation with approximation error δ=Ω(Hd)\delta=\Omega\left(\sqrt{\frac{H}{d}}\right) where dd is the dimension of the representation and HH is the planning horizon, then the agent still needs to sample exponential number of trajectories to find a near-optimal policy.

We show even if optimal policy can be perfectly predicted by a linear function of the given representation with a strictly positive margin, the agent still requires exponential number of trajectories to find a near-optimal policy.

These lower bounds hold even in deterministic systems and even if the agent knows the transition model. Note these negative results apply to the case where the QQ-function, the model, or the optimal policy can be predicted well by a linear function of the given representation. Since the class of linear functions is a strict subset of many more complicated function classes, including neural networks in particular, our negative results imply lower bounds for these more complex function classes as well. Our results highlight the following conceptual insights:

The requirements on the representation that suffice for sample efficient RL are significantly more stringent than the more traditional approximation viewpoint; our statistical lower bounds show that there are hard thresholds on the worst-case approximation quality of the representation which are not necessary from the approximation viewpoint.

Since our lower bounds apply even when the agent knows the transition model, the hardness is not due to the difficulty of exploration in the standard sense. The unknown reward function is sufficient to make the problem exponentially difficult.

Our lower bounds are not due to the agent’s inability to perform efficient supervised learning, since our assumptions do admit polynomial sample complexity upper bounds if the data distribution is fixed.

Our lower bounds are not pathological in nature and suggest that these concerns may arise in practice. In a precise sense, almost all feature extractors induce a hard MDP instance in our construction (see Section 4.4).

Instead, one interpretation is that the hardness is due to a distribution mismatch in the following sense: the agent does not know which distribution to use for minimizing a (supervised) learning error (see Kakade 2003 for discussion), and even a known transition model is not information-theoretically sufficient to reduce the sample complexity.

Furthermore, our work implies several interesting exponential separations on the sample complexity between: 1) value-based learning with perfect representation and value-based learning with a good-but-not-perfect representation, 2) value-based learning and policy-based learning, 3) policy-based learning and supervised learning and 4) reinforcement learning and imitation learning. We provide more details in Section 5.

Related Work

A summary of previous upper bounds, together with lower bounds proved in this paper, is provided in Table 1. Some key assumptions are formally stated in Section 3 and Section 4. Our lower bounds highlight that classical complexity measures in supervised learning including small approximation error and margin, and standard assumptions in reinforcement learning including optimality gap and deterministic systems, are not enough for efficient RL with function approximation. We need additional assumptions, e.g., ones used in previous upper bounds, for efficient RL.

Existing exponential lower bounds, to our knowledge, construct unstructured MDPs with an exponentially large state space and reduce a bandit problem with exponentially many arms to an MDP (Krishnamurthy et al. 2016; Sun et al. 2017). However, these lower bounds cannot apply to MDPs whose transition models, value functions, or policies can be approximated with some natural function classes, e.g., linear functions, neural networks, etc. The current paper gives the first set of lower bounds for RL with linear function approximation (and thus also hold for super classes of linear functions such as neural networks).

2 Previous Upper Bounds

We divide previous algorithms (with provable guarantees) into three classes: those that utilize uncertainty-based bonuses (e.g. UCB variants or Thompson sampling variants); approximate dynamic programming variants (which often make assumptions with respect to concentrability coefficients); and direct policy search-based methods (such as conserve policy iteration (CPI, see Kakade 2003) or policy gradient methods, which make assumptions with respect to distribution mismatch coefficients). The first class of methods include those based on witness rank, Belman rank, and the Eluder dimension, while the latter two classes of algorithms make assumptions either on concentrability coefficients or on distribution mismatch coefficients (see Agarwal et al. 2019; Scherrer 2014 for discussions).

Uncertainty bonus-based algorithms. Now we discuss existing theoretical results on value-based learning with function approximation. The most relevant work is Wen & Van Roy 2013 which showed in deterministic systems, if the optimal QQ-function is within a pre-specified function class which has bounded Eluder dimension, for which the class of linear functions is a special case, then the agent can learn the optimal policy using polynomial number of samples. This result has recently been generalized by Du et al. 2019a which can deal with stochastic reward and low variance transition but requires strictly positive optimality gap. As we listed in Table 1, it is an open problem whether the condition that the optimal QQ-function is linear itself is sufficient for efficient RL.

Li et al. 2011 proposed a QQ-learning algorithm which requires the Know-What-It-Knows oracle. However, it is in general unknown how to implement such oracle in practice. Jiang et al. 2017 proposed the concept of Bellman Rank to characterize the sample complexity of value-based learning methods and gave an algorithm that has polynomial sample complexity in terms of the Bellman Rank, though the proposed algorithm is not computationally efficient. Bellman rank is bounded for a wide range of problems, including MDP with small number of hidden states, linear MDP, LQR, etc. Later work gave computationally efficient algorithms for certain special cases (Dann et al. 2018; Du et al. 2019a; Yang & Wang 2019b; Jin et al. 2019). Recently, Witness rank, a generalization of Bellman rank to model-based methods, is studied in Sun et al. 2019.

Direct policy search-based algorithms. Stronger guarantees over approximate dynamic programming-based algrithm can be obtained with direct policy search-based methods, where instead of having a bounded concentrability coefficient, one only needs to have a bounded distribution mismatch coefficient. The latter assumption requires the agent to have access to a “good” initial state distribution (e.g. a measure which has coverage over where an optimal policy tends to visit); note that this assumption does not make restrictions over the class of MDPs. There are two classes of algorithms that fall into this category. First, there is Conservative Policy Iteration (Kakade & Langford 2002), along with Policy Search by Dynamic Programming (PSDP) (Bagnell et al. 2004), and other boosting-style of policy search-based methods Scherrer & Geist 2014; Scherrer 2014, which have guarantees in terms of bounded distribution mismatch ratio. Second, more recently, Agarwal et al. 2019 showed that policy gradient styles of algorithms also have comparable guarantees.

Recent extensions. Subsequent to this work, the work by Van Roy & Dong 2019 and Lattimore & Szepesvari 2019 made notable contributions to the misspecified linear bandit problem. In particular, both papers found that Theorem 4.1 in our paper can be extended to the misspecified linear bandit problem and gave upper bounds for this problem showing that our lower bound has tight dependency on δ\delta and dd. Lattimore & Szepesvari 2019 further gave an upper bound for the setting where the QQ-functions of all policies can be approximated by linear functions with small approximation errors and the agent can interact with the environment using a generative model. This upper bound also demonstrates that our lower bound has tight dependency on δ\delta and dd.

Preliminaries

Throughout this paper, for a given integer HH, we use [H][H] to denote the set {0,1,…,H−1}\{0,1,\ldots,H-1\}.

In this paper we prove lower bounds for deterministic systems, i.e., MDPs with deterministic transition PP, deterministic reward RR. In this setting, PP and RR can be regarded as functions instead of distributions. Since deterministic systems are special cases of general stochastic MDPs, lower bounds proved in this paper still hold for more general MDPs.

2 QQ-function and Optimality Gap

Here, ρ\rho is the smallest reward-to-go difference between the best set of actions and the rest. Recently, Du et al. 2019b gave a provably efficient QQ-learning algorithm based on this assumption and Simchowitz & Jamieson 2019 showed that with this condition, the agent only incurs logarithmic regret in the tabular setting.

3 Query Models

Here we discuss three possible query oracles interacting with the MDP.

RL: The most basic and weakest query oracle for MDP is the standard reinforcement learning query oracle where the agent can only interact with the MDP by choosing actions and observe the next state and the reward.

Generative Model: A stronger query model assumes the agent can transit to any state (Kearns & Singh 2002; Kakade 2003; Sidford et al. 2018). This query model is available in certain robotic applications where one can control the robot to reach the target state.

Known Transition: The strongest query model considered is that the agent can not only transit to any state, but also knows the whole transition function. In this model, only the reward is unknown.

In this paper, we will prove lower bounds for the strongest Known Transition query oracle. Therefore, our lower bounds also apply to RL and Generative Model query oracles.

Main Results

In this section we formally present our lower bounds. We also discuss proof ideas in Section 4.4.

Here δ\delta is the approximation error, which indicates the quality of the representation. If δ=0\delta=0, then QQ-function can be perfectly predicted by a linear function of ϕ(⋅,⋅)\phi\left(\cdot,\cdot\right). In general, δ\delta becomes smaller as we increase the dimension of ϕ\phi, since larger dimension usually has more expressive power. When the feature extractor is strong enough, previous papers (Chen & Jiang 2019; Farahmand 2011) assume that linear functions of ϕ\phi can approximate the QQ-function of any policy.

In the theoretical reinforcement learning literature, Assumption 4.2 is often called the (approximate) policy completeness assumption. This assumption is crucial in proving polynomial sample complexity guarantee for value iteration type of algorithms (Chen & Jiang 2019; Farahmand 2011).

The following theorem shows when δ=Ω(Hd)\delta=\Omega\left(\sqrt{\frac{H}{d}}\right), the agent needs to sample exponential number of trajectories to find a near-optimal policy.

There exists a family of MDPs with ∣A∣=2|\mathcal{A}|=2 and a feature extractor ϕ\phi that satisfy Assumption 4.2, such that any algorithm that returns a 1/21/2-optimal policy with probability 0.90.9 needs to sample Ω(min⁡{∣S∣,2H,exp⁡(dδ2/16)})\Omega\left(\min\{|\mathcal{S}|,2^{H},\exp(d\delta^{2}/16)\}\right) trajectories.

Note this lower bound also applies to MDPs that satisfy Assumption 4.1, since Assumption 4.2 is strictly stronger. We would like to emphasize that since linear functions is a subclass of more complicated function classes, e.g., neural networks, our lower bound also holds for these function classes. Moreover, in many scenarios, the feature extractor ϕ\phi is the last layer of a neural network. Modern neural networks are often over-parameterized, which makes dd large. In this case, dd is much larger than HH. Thus, our lower bound holds even if the representation has small approximation error. Furthermore, the assumption that ∣A∣=2|\mathcal{A}|=2 is only for simplicity. Our lower bound can be easily generalized to the case that ∣A∣>2|\mathcal{A}|>2, in which case the sample complexity lower bound is Ω(min⁡{∣S∣,∣A∣H,exp⁡(dδ2/16)})\Omega\left(\min\{|\mathcal{S}|,|\mathcal{A}|^{H},\exp(d\delta^{2}/16)\}\right).

2 Lower Bound for Model-based Learning

It has been shown in Yang & Wang 2019b; Yang & Wang 2019a; Jin et al. 2019 if ∥P(⋅∣s,a)−⟨ψ(⋅),ϕ(s,a)⟩∥1\|P\left(\cdot\mid s,a\right)-\langle\psi(\cdot),\phi\left(s,a\right)\rangle\|_{1} is bounded, then the problem admits an algorithm with polynomial sample complexity. Now we show that when δ=Ω(Hd)\delta=\Omega\left(\sqrt{\frac{H}{d}}\right) in Assumption 4.3, the agent needs exponential number of samples to find a near-optimal policy.

There exists a family of MDPs with ∣A∣=2|\mathcal{A}|=2 and a feature extractor ϕ\phi that satisfy Assumption 4.3, such that any algorithm that returns a 1/21/2-optimal policy with probability 0.90.9 needs to sample Ω(min⁡{∣S∣,2H,exp⁡(dδ2/16)})\Omega\left(\min\{|\mathcal{S}|,2^{H},\exp(d\delta^{2}/16)\}\right) trajectories.

Again, our lower bound can be easily generalized to the case that ∣A∣>2|\mathcal{A}|>2.

3 Lower Bound for Policy-based Learning

Similar to value-based learning, a natural assumption for policy-based learning is that the optimal policy is realizable Unlike value-based learning, it is hard to define completeness on the policy-based learning with function approximation, since not all policy has the arg max⁡\argmax form. , i.e., the optimal policy is linear.

Here we discuss another assumption. For learning a linear classifier in the supervised learning setting, one can reduce the sample complexity significantly if the optimal linear classifier has a margin.

Here we restrict the linear coefficients and features to have unit norm for normalization. Note that Assumption 4.5 is strictly stronger than Assumption 4.4. Now we present our result for linear policy.

There exists an absolute constant △0\triangle_{0}, such that for any △≤△0\triangle\leq\triangle_{0}, there exists a family of MDPs with ∣A∣=2|\mathcal{A}|=2 and a feature extractor ϕ\phi that satisfy Assumption 3.1 with ρ=12min⁡{H,d}\rho=\frac{1}{2\min\{H,d\}} and Assumption 4.5, such that any algorithm that returns a 1/41/4-optimal policy with probability at least 0.90.9 needs to sample Ω(min⁡{2H,2d})\Omega\left(\min\{2^{H},2^{d}\}\right) trajectories.

Again, our lower bound can be easily generalized to the case that ∣A∣>2|\mathcal{A}|>2.

Compared with Theorem 4.1, Theorem 4.3 is even more pessimistic, in the sense that even with perfect representation with benign properties (gap and margin), the agent still needs to sample exponential number of samples. It also suggests that policy-based learning could be very different from supervised learning.

4 Proof Ideas

All our lower bound are proved based on reductions from the following hard instance. In this instance, both the transition PP and the reward RR are deterministic. There are HH levels of states, which form a full binary tree of depth HH. There are 2h2^{h} states in level hh, and thus 2H−12^{H}-1 states in total. Among all the 2H−12^{H-1} states in level H−1H-1, there is only one state with reward R=1R=1, and for all other states in the MDP, the corresponding reward value R=0R=0. Intuitively, to find a 1/21/2-optimal policy for such MDPs, the agent must enumerate all possible states in level H−1H-1 to find the state with reward R=1R=1. Doing so intrinsically induces a sample complexity of Ω(2H)\Omega(2^{H}). This intuition is formalized in Theorem A.1 using Yao’s minimax principle (Yao 1977).

Lower bound for value-based and model-based learning

Lower bound for policy-based learning.

It is straightfoward to construct a set of feature vectors for the binary tree instance so that Assumption 4.4 holds, even if d=1d=1. We set ϕ(s,a)\phi(s,a) to be +1+1 if a=a1a=a_{1} and −1-1 if a=a2a=a_{2}. For each level hh, for the unique state ss in level hh with Q∗=1Q^{*}=1, we set θh\theta_{h} to be 11 if π∗(s)=a1\pi^{*}(s)=a_{1} and −1-1 if π∗(s)=a2\pi^{*}(s)=a_{2}. With this construction, Assumption 4.4 holds.

Separations

Perfect representation vs. good-but-not-perfect representation. For value-based learning in deterministic systems, Wen & Van Roy 2013 showed polynomial sample complexity upper bound when the representation can perfectly predict the QQ-function. In contrast, if the representation is only able to approximate the QQ-function, then the agent requires exponential number of trajectories. This exponential separation demonstrates a provable exponential benefit of better representation.

Value-based learning vs. policy-based learning. Note that if the optimal QQ-function can be perfectly predicted by the provided representation, then the optimal policy can also be perfectly predicted using the same representation. Since Wen & Van Roy 2013 showed polynomial sample complexity upper bound when the representation can perfectly predict the QQ-function, our lower bound on policy-based learning, which applies to perfect representations, thus demonstrates that the ability of predicting the QQ-function is much stronger than that of predicting the optimal policy.

Supervised learning vs. reinforcement learning. For policy-based learning, if the planning horizon H=1H=1, the problem becomes learning a linear classifier, for which there are polynomial sample complexity upper bounds. For policy-based learning, the agent needs to learn HH linear classifiers sequentially. Our lower bound on policy-based learning shows the sample complexity dependency on HH is exponential.

Imitation learning vs. reinforcement learning. In imitation learning (IL), the agent can observe trajectories induced by the optimal policy (expert). If the optimal policy is linear in the given representation, it can be shown that the simple behavior cloning algorithm only requires polynomial number of samples to find a near-optimal policy (Ross et al. 2011). Our Theorem 4.3 shows if the agent cannot observe expert’s behavior, then it requires exponential number of samples. Therefore, our lower bound shows there is an exponential separation between policy-based RL and IL when function approximation is used.

Acknowledgments

The authors would like to thank Yuping Luo, Wenlong Mou, Martin Wainwright, Mengdi Wang and Yifan Wu for insightful discussions. Also, the authors would also like to gratefully acknowledge Benjamin Van Roy, Shi Dong, Tor Lattimore and Csaba Szepesvári for sharing a draft of their work and their comments. Simon S. Du is supported by NSF grant DMS-1638352 and the Infosys Membership. Sham M. Kakade acknowledges funding from the Washington Research Foundation Fund for Innovation in Data-Intensive Discovery; the NSF award CCF 1740551; and the ONR award N00014-18-1-2247. Ruosong Wang is supported in part by NSF IIS1763562, AFRL CogDeCON FA875018C0014, and DARPA SAGAMORE HR00111990016. Part of this work was done while Simon S. Du was visiting Google Brain Princeton and Ruosong Wang was visiting Princeton University.

References

Appendix A Proofs of Lower Bounds

We first introduce the INDEX-QUERY problem, which will be useful in our lower bound arguments.

In the INDQn\mathsf{INDQ}_{n} problem, there is an underlying integer i∗∈[n]i^{*}\in[n]. The algorithm sequentially (and adaptively) outputs guesses i∈[n]i\in[n] and queries whether i=i∗i=i^{*}. The goal is to output i∗i^{*}, using as few queries as possible.

For a real number δ∈(0,1)\delta\in(0,1), we say a randomized algorithm A\mathcal{A} is δ\delta-correct for INDQn\mathsf{INDQ}_{n}, if for any underlying integer i∗∈[n]i^{*}\in[n], with probability at least 1−δ1-\delta, A\mathcal{A} outputs i∗i^{*}.

The following theorem states the query complexity of INDQn\mathsf{INDQ}_{n} for 0.10.1-correct algorithms, whose proof is provided in Section B.1.

Any 0.10.1-correct algorithm A\mathcal{A} for INDQn\mathsf{INDQ}_{n} requires at least 0.9n0.9n queries in the worst case.

In this section we prove Theorem 4.1. We need the following existential result, whose proof is provided in Section B.2.

∥pi∥2=1\|p_{i}\|_{2}=1 for all 0≤i≤n−10\leq i\leq n-1;

∣⟨pi,pj⟩∣≤ε\left|\langle p_{i},p_{j}\rangle\right|\leq\varepsilon for any 0≤i,j≤n−10\leq i,j\leq n-1 with i≠ji\neq j.

Now we give the construction of the hard MDP instances. We first define the transitions and the reward functions. In the hard instances, both the rewards and the transitions are deterministic. There are HH levels of states, and level h∈[H]h\in[H] contains 2h2^{h} distinct states. Thus we have ∣S∣=2H−1|\mathcal{S}|=2^{H}-1. If ∣S∣>2H−1|\mathcal{S}|>2^{H}-1 we simply add dummy states to the state space S\mathcal{S}. We use s0,s1,…,s2H−2s_{0},s_{1},\ldots,s_{2^{H}-2} to name these states. Here, s0s_{0} is the unique state in level h=0h=0, s1s_{1} and s2s_{2} are the two states in level h=1h=1, s3s_{3}, s4s_{4}, s5s_{5} and s6s_{6} are the four states in level h=2h=2, etc. There are two different actions, a1a_{1} and a2a_{2}, in the MDPs. For a state sis_{i} in level hh with h<H−1h<H-1, playing action a1a_{1} transits state sis_{i} to state s2i+1s_{2i+1} and playing action a2a_{2} transits state sis_{i} to state s2i+2s_{2i+2}, where s2i+1s_{2i+1} and s2i+2s_{2i+2} are both states in level h+1h+1. See Figure 1 for an example with H=3H=3.

In our hard instances, r(s,a)=0r(s,a)=0 for all (s,a)(s,a) pairs except for a unique state ss in level H−2H-2 and a unique action a∈{a1,a2}a\in\{a_{1},a_{2}\}. It is convenient to define r‾(s′)=r(s,a)\overline{r}(s^{\prime})=r(s,a), if playing action aa transits ss to s′s^{\prime}. For our hard instances, we have r‾(s)=1\overline{r}(s)=1 for a unique node ss in level H−1H-1 and r‾(s)=0\overline{r}(s)=0 for all other nodes.

By construction, for each level h∈[H]h\in[H], there is a unique state shs_{h} in level hh and action ah∈{a1,a2}a_{h}\in\{a_{1},a_{2}\}, such that Q∗(sh,ah)=1Q^{*}(s_{h},a_{h})=1. For all other (s,a)(s,a) pairs such that s≠shs\neq s_{h} or a≠aha\neq a_{h}, it is satisfied that Q∗(s,a)=0Q^{*}(s,a)=0. For a given level hh and policy π\pi, we take θhπ\theta_{h}^{\pi} to be Qπ(sh,ah)⋅ϕ(sh,ah)Q^{\pi}(s_{h},a_{h})\cdot\phi(s_{h},a_{h}). Now we show that ∣Qπ(s,a)−⟨θhπ,ϕ(s,a)⟩∣≤δ|Q^{\pi}(s,a)-\langle\theta_{h}^{\pi},\phi(s,a)\rangle|\leq\delta for all states ss in level hh and a∈{a1,a2}a\in\{a_{1},a_{2}\}.

In this case, we have Qπ(s,a)=0Q^{\pi}(s,a)=0 and ⟨θhπ,ϕ(s,a)⟩=0\langle\theta_{h}^{\pi},\phi(s,a)\rangle=0, since θhπ\theta_{h}^{\pi} and ϕ(s,a)\phi(s,a) do not have a common non-zero coordinate.

In this case, by the second property of P\mathcal{P} in Lemma A.1 and the fact that Qπ(sh,ah)≤1Q^{\pi}(s_{h},a_{h})\leq 1, we have ∣⟨θhπ,ϕ(s,a)⟩∣≤δ|\langle\theta_{h}^{\pi},\phi(s,a)\rangle|\leq\delta. Meanwhile, we have Qπ(s,a)=0Q^{\pi}(s,a)=0.

In this case, we have ⟨θhπ,ϕ(s,a)⟩=Qπ(sh,ah)\langle\theta_{h}^{\pi},\phi(s,a)\rangle=Q^{\pi}(s_{h},a_{h}).

Finally, we prove any algorithm that solves these MDP instances and succeeds with probability at least 0.90.9 needs to sample at least 920⋅2H\frac{9}{20}\cdot 2^{H} trajectories. We do so by providing a reduction from INDQ2H−1\mathsf{INDQ}_{2^{H-1}} to solving MDPs. Suppose we have an algorithm for solving these MDPs, we show that such an algorithm can be transformed to solve INDQ2H−1\mathsf{INDQ}_{2^{H-1}}. For a specific choice of i∗i^{*} in INDQ2H−1\mathsf{INDQ}_{2^{H-1}}, there is a corresponding MDP instance with

Notice that for all MDPs that we are considering, the transition and features are always the same. Thus, the only thing that the learner needs to learn by interacting with the environment is the reward value. Since the reward value is non-zero only for states in level H−1H-1, each time the algorithm for solving MDP samples a trajectory that ends at state sis_{i} where sis_{i} is a state in level H−1H-1, we query whether i∗=i−2H−1+1i^{*}=i-2^{H-1}+1 or not in INDQ2H−1\mathsf{INDQ}_{2^{H-1}}, and return reward value 1 if i∗=i−2H−1+1i^{*}=i-2^{H-1}+1 and 0 otherwise. If the algorithm is guaranteed to return a 1/21/2-optimal policy, then it must be able to find i∗i^{*}.

A.2 Proof of Lower Bound for Model-based Learning

We use the same construction as in the proof of Theorem 4.1. Note we just need to verify that the construction satisfies Assumption 4.3. By construction, for all h∈{1,2,…,H−1}h\in\{1,2,\ldots,H-1\}, for each state s′s^{\prime} in level hh, there exists a unique (s,a)(s,a) pair such that playing action aa transits ss to s′s^{\prime}, and we take ψ(s′)=ϕ(s,a)\psi(s^{\prime})=\phi(s,a). We also take βh=0\beta_{h}=0 for h∈{0,1,…,H−4,H−3}h\in\{0,1,\ldots,H-4,H-3\} and βH−2=ϕ(s,a)\beta_{H-2}=\phi(s,a) where (s,a)(s,a) is the unique pair with R(s,a)=1R(s,a)=1. Now, according to the design of ϕ(⋅,⋅)\phi(\cdot,\cdot) and Lemma A.1, Assumption 4.3 is satisfied. ∎

A.3 Proof of Lower Bound for Policy-based Learning

In this section, we present our hardness results for linear policy learning. We first prove a weaker lower bound which only satisfies Assumption 4.4, and then prove Theoerem 4.3.

Verifying Assumption 4.4.

Recall that for each level h∈[H]h\in[H], there is a unique state shs_{h} in level hh and action ah∈{a1,a2}a_{h}\in\{a_{1},a_{2}\}, such that Q∗(sh,ah)=1Q^{*}(s_{h},a_{h})=1. For all other (s,a)(s,a) pairs such that s≠shs\neq s_{h} or a≠aha\neq a_{h}, it is satisfied that Q∗(s,a)=0Q^{*}(s,a)=0. We simply take θh\theta_{h} to be 11 if ah=a1a_{h}=a_{1}, and take θh\theta_{h} to be −1-1 if ah=a2a_{h}=a_{2}.

Using the same lower bound argument (by reducing INDEX-QUERY to MDPs), we have the following theorem.

There exists a family of MDPs and a feature map ϕ(⋅,⋅)\phi\left(\cdot,\cdot\right) that satisfy Assumption 4.4 with d=1d=1 and Assumption 3.1 with ρ=1\rho=1, such that any algorithm that returns a 1/21/2-optimal policy with probability at least 0.90.9 needs to sample Ω(2H)\Omega\left(2^{H}\right) trajectories.

Proof of Theoerem 4.3

Now we prove Theoerem 4.3. In order to prove Theoerem 4.3, we need the following geometric lemma whose proof is provided in Section B.3.

Now we are ready to prove Theorem 4.3. In the proof we assume H=dH=d, since otherwise we can take HH and dd to be min⁡{H,d}\min\{H,d\} by decreasing the planning horizon HH or adding dummy dimensions to the feature extractor ϕ\phi.

We define a set of 2H−12^{H-1} deterministic MDPs. The transitions of these hard instances are exactly the same as those in Section A.1. The main difference is in the definition of the feature map ϕ(⋅,⋅)\phi(\cdot,\cdot) and the reward function. Again in the hard instances, r(s,a)=0r(s,a)=0 for all ss in the first H−2H-2 levels. Using the terminology in Section A.1, we have r‾(s)=0\overline{r}(s)=0 for all states in the first H−1H-1 levels. Now we define r‾(s)\overline{r}(s) for states ss in level H−1H-1. We do so by recursively defining the optimal value function V∗(⋅)V^{*}(\cdot). The initial state s0s_{0} in level 00 satisfies V∗(s0)=1/2V^{*}(s_{0})=1/2. For each state sis_{i} in the first H−2H-2 levels, we have V∗(s2i+1)=V∗(si)V^{*}(s_{2i+1})=V^{*}(s_{i}) and V∗(s2i+2)=V∗(si)−1/2HV^{*}(s_{2i+2})=V^{*}(s_{i})-1/2H. For each state sis_{i} in the level h=H−2h=H-2, we have r‾(s2i+1)=V∗(si)\overline{r}(s_{2i+1})=V^{*}(s_{i}) and r‾(s2i+2)=V∗(si)−1/2H\overline{r}(s_{2i+2})=V^{*}(s_{i})-1/2H. This implies that ρ=1/2H\rho=1/2H. In fact, this implies a stronger property that each state has a unique optimal action. See Figure 2 for an example with H=3H=3.

To define 2H−12^{H-1} different MDPs, for each state ss in level H−1H-1 of the MDP defined above, we define a new MDP by changing r‾(s)\overline{r}(s) from its original value to 11. This also affects the definition of the optimal VV function for states in the first H−1H-1 levels. In particular, for each level i∈{0,1,2,…,H−2}i\in\{0,1,2,\ldots,H-2\}, we have changed the VV value of a unique state in level ii from its original value (at most 1/21/2) to 11. By doing so we have defined 2H−12^{H-1} different MDPs. See Figure 3 for an example with H=3H=3.

We now show that all the 2H−12^{H-1} MDPs constructed above satisfy the linear policy assumption. Namely, we show that for any state ss in level H−1H-1, after changing r‾(s)\overline{r}(s) to be 1, the resulting MDP satisfies the linear policy assumption. As in Section A.1, for each level h∈[H]h\in[H], there is a unique state shs_{h} in level hh and action ah∈{a1,a2}a_{h}\in\{a_{1},a_{2}\}, such that Q∗(sh,ah)=1Q^{*}(s_{h},a_{h})=1. For all other (s,a)(s,a) pairs such that s≠shs\neq s_{h} or a≠aha\neq a_{h}, it is satisfied that Q∗(s,a)=0Q^{*}(s,a)=0. For each level hh, if ah=a1a_{h}=a_{1}, then we take (θh)H/2=1(\theta_{h})_{H/2}=1 and (θh)H=−1(\theta_{h})_{H}=-1, and all other entries in θh\theta_{h} are zeros. If ah=a2a_{h}=a_{2}, we use p‾\overline{p} to denote the vector formed by the first H/2H/2 coordinates of ϕ(sh,a2)\phi(s_{h},a_{2}). By construction, we have p‾∈P‾\overline{p}\in\overline{\mathcal{P}}. We take θh=[ωp‾;0]\theta_{h}=[\omega_{\overline{p}};0] in this case. In any case, we have ∥θh∥2≤2\|\theta_{h}\|_{2}\leq\sqrt{2}. Now for each level hh, if ah=a1a_{h}=a_{1}, then for all states ss in level hh, we have π∗(s)=a1\pi^{*}(s)=a_{1}. In this case, ⟨ϕ(s,a1),θh⟩=1\langle\phi(s,a_{1}),\theta_{h}\rangle=1 and ⟨ϕ(s,a2),θh⟩=−1\langle\phi(s,a_{2}),\theta_{h}\rangle=-1 for all states in level hh, and thus Assumption 4.5 is satisfied. If ah=a2a_{h}=a_{2}, then π∗(sh)=a2\pi^{*}(s_{h})=a_{2} and π∗(s)=a1\pi^{*}(s)=a_{1} for all states s≠shs\neq s_{h} in level hh. By construction, we have ⟨θh,ϕ(s,a1)⟩=0\langle\theta_{h},\phi(s,a_{1})\rangle=0 for all states ss in level hh, since θh\theta_{h} and ϕ(s,a1)\phi(s,a_{1}) do not have a common non-zero entry. We also have ⟨θh,ϕ(sh,a2)⟩≥2△\langle\theta_{h},\phi(s_{h},a_{2})\rangle\geq 2\triangle and ⟨θh,ϕ(s,a2)⟩≤−2△\langle\theta_{h},\phi(s,a_{2})\rangle\leq-2\triangle for all states s≠shs\neq s_{h} in level hh. Finally, we normalize all θh\theta_{h} and ϕ(s,a)\phi(s,a) so that they all have unit norm. Since ∥ϕ(s,a)∥2=2\|\phi(s,a)\|_{2}=\sqrt{2} for all (s,a)(s,a) pairs before normalization, Assumption 4.5 is still satisfied after normalization.

Finally, we prove any algorithm that solves these MDP instances and succeeds with probability at least 0.90.9 needs to sample at least Ω(2H)\Omega(2^{H}) trajectories. We do so by providing a reduction from INDQ2H−1\mathsf{INDQ}_{2^{H-1}} to solving MDPs. Suppose we have an algorithm for solving these MDPs, we show that such an algorithm can be transformed to solve INDQ2H−1\mathsf{INDQ}_{2^{H-1}}. For a specific choice of i∗i^{*} in INDQ2H−1\mathsf{INDQ}_{2^{H-1}}, there is a corresponding MDP instance with

Notice that for all MDPs that we are considering, the transition and features are always the same. Thus, the only thing that the learner needs to learn by interacting with the environment is the reward value. Since the reward value is non-zero only for states in level H−1H-1, each time the algorithm for solving MDP samples a trajectory that ends at state sis_{i} where sis_{i} is a state in level H−1H-1, we query whether i∗=i−2H−1+1i^{*}=i-2^{H-1}+1 or not in INDQ2H−1\mathsf{INDQ}_{2^{H-1}}, and return reward value 1 if i∗=i−2H−1+1i^{*}=i-2^{H-1}+1 and it original reward value otherwise. If the algorithm is guaranteed to return a 1/41/4-optimal policy, then it must be able to find i∗i^{*}.

Appendix B Technical Proofs

The proof is a straightforward application of Yao’s minimax principle Yao 1977. We provide the full proof for completeness.

Consider an input distribution where i∗i^{*} is drawn uniformly at random from [n][n]. Suppose there is a 0.10.1-correct algorithm for INDQn\mathsf{INDQ}_{n} with worst-case query complexity TT such that T<0.9nT<0.9n. By averaging, there is a deterministic algorithm A′\mathcal{A}^{\prime} with worst-case query complexity TT, such that

We may assume that the sequence of queries made by A′\mathcal{A}^{\prime} is fixed. This is because (i) A′\mathcal{A}^{\prime} is deterministic and (ii) before A′\mathcal{A}^{\prime} correctly guesses i∗i^{*}, all responses that A′\mathcal{A}^{\prime} receives are the same (i.e., all guesses are incorrect). We use S={s1,s2,…,sm}S=\{s_{1},s_{2},\ldots,s_{m}\} to denote the sequence of queries made by A′\mathcal{A}^{\prime}. Notice that mm is the worst-case query complexity of A′\mathcal{A}^{\prime}. Suppose m<0.9nm<0.9n, there exist 0.1n0.1n distinct i∈[n]i\in[n] such that A′\mathcal{A}^{\prime} will never guess ii, and will be incorrect if i∗i^{*} equals ii, which implies

B.2 Proof of Lemma A.1

We need the following tail inequality for random unit vectors, which will be useful for the proof of Lemma A.1.

In particular, when β≥6\beta\geq 6,we have

It is clear that ∥qi∥2=1\|q_{i}\|_{2}=1 for all i∈[n]i\in[n], since each qiq_{i} is drawn from the unit sphere. We now prove that for any i,j∈[n]i,j\in[n] with i≠ji\neq j, with probability at least 1−1n21-\frac{1}{n^{2}}, we have ∣⟨qi,qj⟩∣≤ε\left|\langle q_{i},q_{j}\rangle\right|\leq\varepsilon. Notice that this is sufficient to prove the lemma, since by a union bound over all the (n2)=n(n−1)/2\binom{n}{2}=n(n-1)/2 possible pairs of (i,j)(i,j), this implies that Q\mathcal{Q} satisfies the two desired properties with probability at least 1/21/2.

B.3 Proof of Lemma A.2

To prove the lemma, it suffices to show that N\mathcal{N} satisfies the property equation \refeqn:distequation~\ref{eqn:dist}. Consider a point x∈Nx\in\mathcal{N}, let AA be a hyperplane that is perpendicular to xx (notice that xx is a also a vector) and separates xx and every other points in N\mathcal{N}. We let the distance between xx and AA be the largest possible, i.e., AA contains a point in N\{x}\mathcal{N}\backslash\{x\}. Since xx is on the unit sphere and N\mathcal{N} is a ϵ\sqrt{\epsilon}-packing, we have that xx is at least ϵ\sqrt{\epsilon} away from every point on the spherical cap not containing xx, defined by the cutting plane AA. More formally, let bb be the intersection point of the line segment oxox and AA. Then

where z∈N∩Az\in\mathcal{N}\cap A. Notice that the distance between xx and the convex hull of N\{x}\mathcal{N}\backslash\{x\} is lower bounded by the distance between xx and AA, which is given by ∣bx∣|bx|. Consider the triangles defined by x,z,o,bx,z,o,b. We have bz⊥oxbz\perp ox (note that bzbz lies inside AA). By Pythagorean theorem, we have

Solve the above three equations for ∣bx∣|bx|, we have

Appendix C Exact Linear Q∗Q^{*} + Gap in Generative Model

In this section we present and prove the following theorem.

We first describe the algorithm. For each level, the agent first construct a barycentric spanner Λh≜{ϕ(sh1,ah1),…ϕ(shd,ahd)}⊂Φh≜{ϕ(s,a)}s∈Sh,a∈A\Lambda_{h}\triangleq\left\{\phi(s_{h}^{1},a_{h}^{1}),\ldots\phi(s_{h}^{d},a_{h}^{d})\right\}\subset\Phi_{h}\triangleq\left\{\phi\left(s,a\right)\right\}_{s\in\mathcal{S}_{h},a\in\mathcal{A}} (Awerbuch & Kleinberg 2008). We have the property that any ϕ(s,a)\phi(s,a) with sh∈Sh,a∈As_{h}\in\mathcal{S}_{h},a\in\mathcal{A}, we have cs,a1,…,cs,ad∈c_{s,a}^{1},\ldots,c_{s,a}^{d}\in such that ϕ(s,a)=∑i=1dcs,aiϕ(shi,ahi)\phi(s,a)=\sum_{i=1}^{d}c_{s,a}^{i}\phi(s_{h}^{i},a_{h}^{i}).

The algorithm learns the optimal policy from h=H−1,…,0h=H-1,\ldots,0. At any level hh, we assume the agent has learned the optimal policy πh′∗\pi_{h^{\prime}}^{*} at level h′=h+1,…,H−1h^{\prime}=h+1,\ldots,H-1.

Appendix D Linear QπQ^{\pi} for all π\pi in Generative Model

In this section we present and prove the following theorem.

Now we can bound the sub-optimality of π^≜π0∘⋯∘πH−1\hat{\pi}\triangleq\pi_{0}\circ\cdots\circ\pi_{H-1}: