Para la convergencia, ver la respuesta de Ma Ming.
Considere la posibilidad de:
$$f(n,x)=\prod_{k=1}^n(1+x^k)$$
Estamos interesados en la asymptotical comportamiento de $f(n,x)$$n\to\infty$$0\le x<1$. Aviso que el primer $m$ términos en que el producto puede contribuir a los términos de la orden de $x^m$. Cada ocurren $x^m$ corresponde a una partición de $m$ de la forma $\lambda_1>\lambda_2>\ldots>\lambda_\ell>0$, $\lambda_1+\ldots+\lambda_\ell=m$ a través de $x^m=x^{\lambda_1}\ldots x^{\lambda_\ell}$. Deje $p(m)$ el número de particiones de $m$. Esto muestra que:
$$\bigg[f(n,x)\bigg]_{x^m}=\tilde{p}(m)$$
si $n\ge m$. Así tenemos:
$$f(n,x)=\sum_{k=0}^n\tilde{p}(k)x^k+o(x^{n+1})$$
En el límite de $n\to\infty$ consigue la generación de la función de restricción de las particiones:
$$f(x)=\prod_{k=1}^\infty(1+x^k)=\sum_{k=0}^\infty\tilde{p}(k)x^k$$
Esto no tiene ninguna forma cerrada, que yo sepa, y creo que no puede ser evaluado para determinados $x>0$ (excepto numéricamente, por supuesto).
Nota: Obviamente el caso que nos interesa es $x=\frac{1}{3}$, por lo que tenemos $a_{n+1}=f(n,x)$.