Solution (source code)

= Solution

Subtract adjacent values of the <cumulative distribution function>. Its value at zero is zero, so the same calculation applies also to $k=1$:
$$
\boxed{p_k=\left(1-\frac1{k+1}\right)-\left(1-\frac1k\right)=\frac1{k(k+1)},\qquad k\geq1.}
$$
The <probability mass function> sums to one by a <telescoping series>.