Otro enfoque que vale la pena mencionar es el de la exponencial en la generación de funciones (egfs). Deje $x$ enumerar las vocales, $y$ enumerar las consonantes y $z$ enumerar total de cartas. A continuación, el egf factores para las vocales se dividen de la siguiente manera
III
$$1+xz+\frac{1}{2!}(xz)^2+\frac{1}{3!}(xz)^3$$
EE
$$1+xz+\frac{1}{2!}(xz)^2$$
A, O
$$1+xz$$
and consonants:
FF, NN, TT
$$1+yz+\frac{1}{2!}(yz)^2$$
D, R
$$1+yz$$
The product of these is the egf for words
$$\begin{multline}f(x,y,z)=\left(1+xz+\frac{1}{2!}(xz)^2+\frac{1}{3!}(xz)^3\right)\times\\\left(1+xz+\frac{1}{2!}(xz)^2\right)\left(1+\vphantom{\frac{1}{1}}xz\right)^2\left(1+yz+\frac{1}{2!}(yz)^2\right)^3\times\\\left(1+\vphantom{\frac{1}{1}}yz\right)^2\end{multline}$$
and we want the $x^3y^4z^7/7!$ coefficient of $f(x,y,z)$. El uso de la salvia de entrada
y,z=var('y','z')
show(factorial(7)*expand((1+z*x+(z*x)^2/2+(z*x)^3/6)*(1+z*x+(z*x)^2/2)*(1+z*x)^2*(1+z*y+(z*y)^2/2)^3*(1+z*y)^2).coefficient(z^7).coefficient(x^3).coefficient(y^4))
da la respuesta:
$$532\,770\tag{Answer}$$