Solution
= Solution
{c}
Put $p_*=p/[p+(1-p)r]$. Then
$$
G_Z(s)=\frac{p_*}{1-(1-p_*)s}.
$$
Uniqueness of the <probability generating function> therefore gives
$$
\boxed{Z\sim\operatorname{Geom}_0\!\left(\frac p{p+(1-p)r}\right).}
$$
This is <Bernoulli thinning of a geometric count>, using the convention of counting failures before success.