Encontrar la enésima fórmula de una fórmula recursiva $a_n=a_{n-1}+n(n-1)a_{n-2}$

Question

$$a_n=a_{n-1}+n(n-1)a_{n-2}$$
$$a_0=1, a_1=-\frac{1}{2}$$
¿Es posible encontrar una fórmula explícita para $a_n$ sólo con el uso de $a_0$ y $a_1$ ?
Sé cómo resolver este problema si $a_n=Aa_{n-1}+Ba_{n-2}$ donde $A$ y $B$ es algún número real(constante), pero en este problema no es un caso.

Answer 1

Dejemos que $\displaystyle\;b_n = \frac{a_n}{n!}$ la relación de recurrencia en cuestión puede transformarse en

$$n ( b_n - b_{n-2} ) = b_{n-1}\quad\text{ with }\quad \begin{cases} b_0 = 1\\b_1 = -\frac12\end{cases}\tag{*1}$$

Dejemos que $f(z) = \sum_{n=0}^\infty b_n z^n$ sea el OGF para $b_n$ . Si hacemos un múltiplo de $(*1)$ por $z^n$ y empezar a sumar desde $n = 2$ . Obtenemos

$$\begin{align} & z\frac{d}{dz}\left[ (1-z^2) f(z) - 1 + \frac{z}{2} \right] = z(f(z)-1)\\ \iff & (1-z^2) \frac{df(z)}{dz} - (1+2z) f(z) = -\frac{3}{2} \end{align} $$ La resolución de la EDO nos da

$$f(z) = \frac{5 - 3\sin^{-1}(z)}{2(1-z)\sqrt{1-z^2}} - \frac{3}{2(1-z)}$$

Ampliar $f(z)$ como una serie de potencias en $z$ obtenemos lo siguiente expresión fea de $a_n$ .

$$a_n = \frac{n!}{2}\left[ 5\sum_{p=0}^{\lfloor \frac{n}{2}\rfloor} \frac{\binom{2p}{p}}{4^p} - 3\sum_{s=0}^{\lfloor\frac{n-1}{2}\rfloor} \sum_{p=0}^s \frac{\binom{2p}{p}\binom{2s-2p}{s-p}}{4^s(2p+1)} - 3 \right]$$ Como doble comprobación, los primeros valores de $a_n$ según esta fórmula son

$$(2a_0,2a_1,\ldots) = 2, -1, 3, -3, 33, -27, 963, -171, 53757, 41445,4879575, 9438525, \ldots$$

y satisfacen la relación de recurrencia.

Actualización

Resulta que podemos simplificar un poco el lío.

La relación de recurrencia en $b_n$ puede reescribirse como $b_n = (n+1)(b_{n+1} - b_{n-1})$ . El lado derecho tiene la forma de una diferencia finita. Podemos utilizarla para un nivel de suma. El resultado final es que para todos los $n \ge 0$ ,

$$a_n = \frac{(n+1)!}{2}\left[ 5\frac{\binom{2r}{r}}{4^r} - 3\sum_{p=0}^s \frac{\binom{2p}{p}\binom{2s-2p}{s-p}}{4^s(2p+1)} \right] \quad\text{ with }\quad \begin{cases} r = \lfloor \frac{n+1}{2}\rfloor\\ s = \lfloor \frac{n}{2}\rfloor \end{cases} $$

Answer 2

Para la segunda pregunta, observe que

$$\begin{align*}\begin{bmatrix}a_n\\a_{n-1}\end{bmatrix}&=\begin{bmatrix}Aa_{n-1}+Ba_{n-2}\\ a_{n-1}\end{bmatrix}=\begin{bmatrix}A &B\\ 1 & 0\end{bmatrix}\begin{bmatrix}a_{n-1}\\a_{n-2}\end{bmatrix}=\begin{bmatrix}A&B\\1&0\end{bmatrix}^2\begin{bmatrix}a_{n-2}\\ a_{n-3}\end{bmatrix}\\&=\cdots=\begin{bmatrix}A&B\\1&0\end{bmatrix}^{n-2}\begin{bmatrix}a_1\\a_0\end{bmatrix}\end{align*}$$

Con un argumento similar, obtendríamos que

$$\begin{align*}\begin{bmatrix}a_n\\a_{n-1}\end{bmatrix}&=\begin{bmatrix}a_{n-1}+n(n-1)a_{n-2}\\ a_{n-1}\end{bmatrix}=\begin{bmatrix}1 &n(n-1)\\ 1 & 0\end{bmatrix}\begin{bmatrix}a_{n-1}\\a_{n-2}\end{bmatrix}\\&=\begin{bmatrix}1&n(n-1)\\1&0\end{bmatrix}\begin{bmatrix}1 &(n-1)(n-2)\\1&0\end{bmatrix}\begin{bmatrix}a_{n-2}\\ a_{n-3}\end{bmatrix}\\&=\cdots \\&=\left(\prod_{k=0}^{n-2} \begin{bmatrix}1&(n-k)(n-k-1)\\1&0\end{bmatrix}\right)\begin{bmatrix}a_1\\a_0\end{bmatrix}\end{align*}$$