A342182 Sum_{n>=0} a(n) * x^n / (n!)^2 = 1 / (1 - x * BesselI(0,2*sqrt(x))).

1, 1, 8, 117, 3184, 134025, 8141436, 672837277, 72634878016, 9923765772177, 1673881314096700, 341631408064928421, 82978986493779894288, 23653894531273155603961, 7819996460332550715977588, 2967815528758036870644773925, 1281517958938232539844046259456

Offset

0,3

Links

Table of n, a(n) for n=0..16.

Formula

Sum_{n>=0} a(n) * x^n / (n!)^2 = 1 / (1 - Sum_{n>=1} x^n / ((n-1)!)^2).

a(0) = 1; a(n) = Sum_{k=0..n-1} (binomial(n,k) * (n-k))^2 * a(k).

a(n) ~ n!^2 / ((1 + r^(3/2)*BesselI(1, 2*sqrt(r))) * r^n), where r = 0.592860029867912878114616561736048937618032595935338954527835... is the root of the equation r*BesselI(0, 2*sqrt(r)) = 1. - Vaclav Kotesovec, May 04 2024

Mathematica

nmax = 16; CoefficientList[Series[1/(1 - x BesselI[0, 2 Sqrt[x]]), {x, 0, nmax}], x] Range[0, nmax]!^2
a[0] = 1; a[n_] := a[n] = Sum[(Binomial[n, k] (n - k))^2 a[k], {k, 0, n - 1}]; Table[a[n], {n, 0, 16}]

Prog

(PARI) a(n) = {n!^2*polcoef(1 / (1 - sum(k=1, n, x^k / ((k-1)!)^2) + O(x*x^n)), n)} \\ Andrew Howroyd, Mar 04 2021

Crossrefs

Cf. A006153, A101514, A336228, A337591.

Sequence in context: A046914 A364984 A341197 * A194494 A230911 A267752

Adjacent sequences: A342179 A342180 A342181 * A342183 A342184 A342185

Keyword

nonn

Author

Ilya Gutkovskiy, Mar 04 2021

Status

approved