Solution
= Solution
Set $f_n=n$. A birth changes $f$ by $1$, an ordinary pair death by $-2$, and singleton death by $-1$. Hence
$$
\dot\mu=b-2d\sum_{n\geq2}np_n-Dp_1
=b-2d\mu+(2d-D)p_1.
$$
When $D=2d$ this closes to
$$
\boxed{\dot\mu=b-2d\mu},
$$
so the stationary mean is $\mu_*=b/(2d)$.