OFFSET
0,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..363
Z. Li, Z. Li and Y. Cao, Enumeration of symplectic and orthogonal injective partial transformations, Discrete Math., 306 (2006), 1781-1787. (The function s_n.)
FORMULA
E.g.f.: (1/(1-4x))*exp(x/(1-4x)).
a(n) = (8*n-3)*a(n-1) - 16*(n-1)^2*a(n-2). - Vaclav Kotesovec, Oct 09 2013
a(n) ~ n^(n+1/4) * exp(sqrt(n)-n-1/8) * 4^n * (1 + 37/(96*sqrt(n))). - Vaclav Kotesovec, Oct 09 2013
Sum_{n>=0} a(n) * x^n / (n!)^2 = exp(4*x) * BesselI(0,2*sqrt(x)). - Ilya Gutkovskiy, Jul 17 2020
MAPLE
seq(sum('binomial(k, i)^2*i!*4^i', 'i'=0..k), k=0..30);
MATHEMATICA
f[n_] := Sum[k!*4^k*Binomial[n, k]^2, {k, 0, n}]; Table[ f[n], {n, 0, 16}] (* or *)
Range[0, 16]! CoefficientList[ Series[1/(1 - 4x)*Exp[x/(1 - 4x)], {x, 0, 16}], x] (* Robert G. Wilson v, Mar 16 2005 *)
PROG
(PARI) a(n)=my(t=1); sum(i=1, n, t*=i; binomial(n, i)^2*t<<(2*i), 1) \\ Charles R Greathouse IV, Oct 23 2023
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Miklos Kristof, Mar 16 2005
EXTENSIONS
More terms from Robert G. Wilson v, Mar 16 2005
STATUS
approved