OFFSET
0,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..200
FORMULA
E.g.f.: 1/3 + (2/3)*Sum_{n>=1} exp(n^2*x) / 2^n = Sum_{n>=0} a(n)*x^n/n!.
a(n) = (4/3) * Sum_{k=0..2*n} k! * Stirling2(2*n, k) for n>0 with a(0)=1.
a(n) = Sum_{k=1..[(2*n+1)/3]} 2 * (3*k-1)! * Stirling2(2*n+1, 3*k) for n>0 with a(0)=3, after Vladimir Kruchinin in A242858.
EXAMPLE
E.g.f.: E(x) = 1 + 4*x^2/2! + 100*x^4/4! + 6244*x^6/6! + 727780*x^8/8! +...
where E(x) = 1/(5 - 4*cosh(x)) = -exp(x) / (2 - 5*exp(x) + 2*exp(2*x)).
ALTERNATE GENERATING FUNCTION.
E.g.f.: A(x) = 1 + 4*x + 100*x^2/2! + 6244*x^3/3! + 727780*x^4/4! +...
where 3*A(x) = 1 + 2*exp(x)/2 + 2*exp(4*x)/2^2 + 2*exp(9*x)/2^3 + 2*exp(16*x)/2^4 + 2*exp(25*x)/2^5 + 2*exp(36*x)/2^6 + 2*exp(49*x)/2^7 +...
PROG
(PARI) /* E.g.f.: 1/(5 - 4*cosh(x)) */
{a(n) = local(X=x+O(x^(2*n+1))); (2*n)!*polcoeff( 1/(5 - 4*cosh(X)), 2*n)}
for(n=0, 20, print1(a(n), ", "))
(PARI) /* Formula for a(n): */
{Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
{a(n) = if(n==0, 1, sum(k=1, (2*n+1)\3, 2*(3*k-1)! * Stirling2(2*n+1, 3*k)))}
for(n=0, 20, print1(a(n), ", "))
(PARI) /* Formula for a(n): */
{Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
{a(n) = if(n==0, 1, (4/3)*sum(k=0, 2*n, k! * Stirling2(2*n, k) ))}
for(n=0, 20, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 19 2014
STATUS
approved