In this paper we provide faster algorithms for approximately solving discounted Markov Decision Processes in multiple parameter regimes. Given a discounted Markov Decision Process (DMDP) with vertical bar S vertical bar states, vertical bar A vertical bar actions, discount factor gamma is an element of (0, 1), and rewards in the range [-M, M], we show how to compute an epsilon-optimal policy, with probability 1 - delta in time(1) (O) over tilde((vertical bar S vertical bar(2)vertical bar A vertical bar + vertical bar S parallel to A vertical bar/(1 - gamma)(3)) log (M/epsilon) log (1/delta)). This contribution reflects the first nearly linear time, nearly linearly convergent algorithm for solving DMDP's for intermediate values of gamma. We also show how to obtain improved sublinear time algorithms and provide an algorithm which computes an epsilon-optimal policy with probability 1 - delta in time (O) over tilde(vertical bar S parallel to A vertical bar M-2/(1 -gamma)(4)epsilon(2) log(1/delta)) provided we can sample from the transition function in O(1) time. Interestingly, we obtain our results by a careful modification of approximate value iteration. We show how to combine classic approximate value iteration analysis with new techniques in variance reduction. Our fastest algorithms leverage further insights to ensure that our algorithms make monotonic progress towards the optimal value. This paper is one of few instances in using sampling to obtain a linearly convergent linear programming algorithm and we hope that the analysis may be useful more broadly.