Suppose I have two quantum channels. Assume they they consist of $r\in \mathbb{Z}$ applications of unitaries, $U$ and $V$ respectively. Let the error between the channels acting on some state $\rho$ be: $$\varepsilon(U,V) = U^{\dagger r} \rho U^r - V^{\dagger r} \rho V^r. $$
Ideally, I want an expression for $\epsilon(U,V)$ in terms of the difference between these unitaries and $r$. Ideally I want an expression for $\varepsilon(U,V)$ as defined above, not the norm of the expression.
My first attempt was to write that $U^r - V^r = \sum_{k=1}^{r-1} U^{r-k}(U-V)V^k$ and substitute into the expression for $\varepsilon(U,V)$. However, this introduces a lot of cross terms which make things messy, and one ends up with $O(r^2)$ many terms when it feels like there should only be $O(r)$. Is there a better way to do this?