Solution (source code)

= Solution
{c}

Since $Y\ge1$, its finite <mean> satisfies $1\le\mu<\infty$. Every $q_k=kp_k/\mu$ is nonnegative, and
$$
\sum_{k\ge1}q_k=\frac1\mu\sum_{k\ge1}kp_k=\frac{\mathbb E Y}{\mu}=1.
$$
Thus $(q_k)$ is a \b[<probability mass function>]. It defines the <size-biased distribution> of $Y$: an outcome of size $k$ receives $k$ times its original relative weight.