Integrability of composite of Riemann integrable functions given condition This is a question from a past qualifying exam, which I am studying for. The question has been asked before here, and has an answer, but the answer uses Lebesgue's criterion for Riemann integrability, which is disallowed on the exam. Is there a more elementary way to solve this question?

Let $f: [0,1] \to \mathbb{R}$ and $g: [0,1] \to [0,1]$ be two Riemann integrable functions. Assume that $|g(x) - g(y)| \geq \alpha |x-y|$ for any $x,y \in [0,1]$ and some fixed $\alpha \in (0,1)$. Show that $f \circ g$ is Riemann integrable.

Some thoughts have been bounding the intervals in which $f$ has a large oscillation by its integrability, and trying use the condition on $g$ to control the growth of these interval lengths. However, I am unsure how to apply the Riemann integrability of $g$. 
 A: I have not worked through all of the details, but here is a sketch of an idea, too long for a comment. I will put question mark at the part I have not thought through yet. 
Set $I:=[0,1]$ and choose an integer $k$ such that $\frac{1}{k}<\frac{\epsilon}{2}.$ The set $D_k=\{x\in I:\text{osc}_x\ f\ge1/k\}$ has measure zero so it has a countable covering by open sets $J_j = (a_j, b_j)$ whose total length is less than $\frac{\epsilon}{2}.$ Now, for each $x\in I\setminus D_k$ there is an open interval $x\in I_x\subseteq I\setminus D_k$ such that $\sup\ f-\inf\ f<1/k$ on $I_x$ (because $\text{osc}_x\ f<1/k$). Then, the $J_j$ and $I_x$ form an open cover of $I$. Let $\lambda$ be the Lebesgue number of this cover and take any partition $Q=\{y_i\}$ of $I$ such that $|Q|<\lambda$ and 
$[y_i,y_{i+1}]\subseteq \text{im}\ g$. ???
Let $M_i,m_i$ be the maxima, resp. minima of $f$ on $[y_i,y_{i+1}].$
Then let $x_i=g^{-1}(y_i)$.  Since $g$ is injective, the $x_i$ form a partition $P$ of $I$ and
$U(f\circ g,P)-L(f\circ g,P)=\sum_i(M_i-m_i)|g^{-1}(y_{i+1})-g^{-1}(y_i)|\le$
$\frac{1}{\alpha}\sum_i(M_i-m_i)(y_{i+1}-y_i).$ 
By construction, $[y_i,y_{i+1}]$ is either in one of the $J_j$ or one of the $I_x.$ Now split this sum up into those subintervals of $P$ that lie in one of the $J_j$ and those that lie in one of the  $I_x$. The set-up in first paragraph shows that the sum is less than $\epsilon.$
