I would like to understand better the scaling of the following summation as a function of $r$ and $p > 1$:

$$ S_r(p) := \sum_{m=1}^{r} \left( \sum_{k=r-m+1}^{r} \left( \frac{k^q}{\sum_{k'=1}^{r} (k')^q} \right)^{1/p} \right)^{-1} \:, \:\: \frac{1}{p} + \frac{1}{q} = 1 \:. $$

In particular, I am interested in upper bounds on $S_r(p)$. Let us consider the two extremes of $p$. When $p$ approaches $1$, we have that $S_r(1) = r$. On the other hand, when $p$ approaches $+\infty$, we have that $S_r(\infty) = H_r = \Theta(\log{r})$, where $H_r$ is the $r$-th Harmonic number. I would like to understand how to analyze the interpolation of $S_r(p)$ between $r$ and $\log{r}$. The intermediate $p=2$ case is also tractable, and is given by $S_r(2) = H_{2r} \sqrt{\frac{2r(1+r)}{3(1+2r)}} = \Theta(\sqrt{r} \log{r} )$.

I did some numerical calculations here. In particular, in the notebook there is a picture of the curve for a few values of $r$.