Ya que todos los coeficientes son positivos, esto equivale a solo verificar que cada número es alcanzable por al menos una combinación de los términos. $a_k$ en el de arriba tiene la combinatoria de interpretación como el número de posibles soluciones a $(c_1,c_2,\dots,c_{15})$ donde cada una de las $c_i\in\{0,1,3,4\}$$c_1+c_2+\dots+c_{15}=k$.
El uso de esta interpretación, cada una de las $a_0,a_1,\dots,a_{15}$, obviamente, son positivas, ya que para $k\in\{0,1,\dots,15\}$ uno se puede topar con al menos uno de acuerdo como $(\underbrace{1,1,\dots,1,}_k~\underbrace{0,0,\dots,0}_{15-k})$. El mismo argumento funciona para mostrar cada una de las $a_{45},a_{46},\dots,a_{60}$ utilizando en su lugar $4$'s y $3$'s en el lugar de la $0$'s y $1$'s.
Cada una de las $a_{16},\dots,a_{27}$ puede ser visto para ser positivo, ya que para $k\in\{16,\dots,27\}$ podemos usar $(4,4,4,4,\underbrace{1,1,1,\dots,1}_{k-16}~,\underbrace{0,0,\dots,0}_{27-k})$
Podemos también construir al menos un acuerdo para $28,29$ $30$ así como también el uso de unos cuantos más, $4$'s en el principio. Finalmente, mediante la sustitución de $0$'s $4$'s y viceversa, y $1$'s $3$'s y viceversa, esto le da a los arreglos para cada una de las $30,\dots 45$ así sin ningún esfuerzo adicional.
Al final todo esto demuestra que cada coeficiente de $a_0,a_1,\dots,a_{60}$ son enteros positivos y así todos los sesenta y un posible aparecen en la plena expansión de $(1+x+x^3+x^4)^{15}$ y de forma equivalente, a todos los sesenta y uno de los posibles términos aparecen en la expansión de $(x+\frac{1}{x}+x^2+\frac{1}{x^2})^{15}$
Mucho más interesante y difícil hubiera sido si algunos de los signos que fueron alterados, como en $(x-\frac{1}{x}+x^2-\frac{1}{x^2})^{15}$, donde por ejemplo el coeficiente del término de $x^0$ sería igual a cero.