Let $X$ and $Y$ be i.i.d. random variables with $M_X(t) = (1-p + pe^t)^n$. Define $S = X + Y$. What is the probability generating function $G_S(z)$ of $S$?