Esta respuesta se basa en la generación de la función de generalizada polinomios de Laguerre
\begin{align*}
L_k^{(\alpha)}(t)=\sum_{i=0}^k(-1)^k\binom{k+\alpha}{k-i}\frac{t^i}{i!}\tag{1}
\end{align*}
Los polinomios de Laguerre tienen notables propiedades combinatorias y uno de ellos es precisamente el adecuado para responder a los problemas de este tipo. Esto está muy bien presentado en el Conteo de palabras con Laguerre serie por Jair Taylor.
Nos codificar los colores (b)lue (r)ed, (g)creen y (y)ellow de las bolas con las letras
\begin{align*}
\{b,r,g,y\}
\end{align*}
y están en busca de palabras de longitud $8$ construido a partir de
\begin{align*}
b,b,r,r,g,g,y,y
\end{align*}
que tienen la propiedad de que no contienen consecutivos iguales letras. Estas palabras se llaman Carlitz palabras o Smirnov palabras.
Nos encontramos en la sección 2 de la referida papel polinomios de Laguerre $l_k(t)$ definido por su función de la generación de
\begin{align*}
\sum_{k=0}^\infty l_k(t)x^k=\exp\left({\frac{tx}{1+x}}\right)
\end{align*}
El primer par de dichos polinomios son
\begin{align*}
l_0(t)&=1\\
l_1(t)&=t\\
l_2(t)&=\frac{1}{2}t^2-t\\
l_3(t)&=\frac{1}{6}t^3-t^2+t\tag{2}
\end{align*}
Ellos son una forma específica de polinomios de Laguerre (1), es decir, $$l_k(t)=(-1)^kL_k^{(-1)}(t)$$
Teorema 2.1 en el referido documento señala: Dados los números enteros no negativos $n_1,\ldots,n_k$, el número de $k$-ary Carlitz palabras con la letra $i$ utiliza exactamente $n_i$ veces es
\begin{align*}
\int_{0}^\infty e^{-t}\left(\prod_{i=1}^kl_{n_i}(t)\right)\,dt\tag{3}
\end{align*}
Ya tenemos cuatro caracteres $b,r,g,y$ cada ocurren dos veces, hemos establecido
\begin{align*}
&n_1=n_2=n_3=n_4=2
\end{align*}
Aplicamos el teorema 2.1. y obtener el uso de (2) y (3) y con la ayuda de Wolfram Alpha
\begin{align*}
\int_{0}^\infty&e^{-t}\left(\prod_{i=1}^4l_{n_i}(t)\right)\,dt\\
&=\int_{0}^\infty e^{-t}\left(l_2(t)\right)^4\,dt\\
&=\int_{0}^\infty e^{-t}\left(\frac{1}{2}t^2-t\right)^4\,dt\\
&=\int_{0}^\infty e^{-t}\left(\frac{1}{16}t^{8}-\frac{1}{2}t^{7}+\frac{3}{2}t^{6}
-2t^5+t^4\right)\,dt\\
&=\color{blue}{864}
\end{align*}
Variación 2:
Otra variación es directamente sobre la base de un generación de la función de Smirnov palabras. (Véase el ejemplo III.24 Smirnov palabras de la Analítica de la Combinatoria de Philippe Flajolet y Robert Sedgewick para obtener más información.)
Una generación de función para el número de Smirnov palabras sobre una de cuatro letras del alfabeto $V=\{b,r,g,y\}$ está dado por
\begin{align*}
\left(1-\frac{4z}{1+z}\right)^{-1}
\end{align*}
El número de todos Smirnov palabras de longitud $8$ a través de una de cuatro letras del alfabeto es, por tanto,
\begin{align*}
[z^8]\left(1-\frac{4z}{1+z}\right)^{-1}
\end{align*}
Ya que queremos contar el número de palabras de longitud $8$ con cada personaje en $V$ se producen dos veces, hacemos un seguimiento de cada personaje. De nuevo con la ayuda de Wolfram Alpha, obtenemos
\begin{align*}
[b^2r^2g^2y^2]\left(1-\frac{b}{1+b}-\frac{r}{1+r}-\frac{g}{1+g}-\frac{y}{1+y}\right)^{-1}=\color{blue}{864}
\end{align*}