login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A195509 Expansion of e.g.f. (exp(x*exp(x)) + exp(x/exp(x)))/2. 1
1, 1, 1, 4, 25, 96, 481, 3368, 20721, 141760, 1146721, 9098112, 77652169, 726208640, 6891125697, 69344336896, 738718169569, 8076031881216, 92647353941569, 1106883171037184, 13616813607795321, 174298975125127168, 2304515271134124577 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,4
LINKS
FORMULA
E.g.f.: Sum_{n>=0} x^n*cosh(n*x)/n!.
a(n) = Sum_{m=0..floor((n-1)/2)} binomial(n,2*m)*(n-2*m)^(2*m) for n >= 1. - Vladimir Kruchinin, Mar 10 2013 [Edited by Petros Hadjicostas, May 06 2020]
EXAMPLE
E.g.f.: A(x) = 1 + x + x^2/2! + 4*x^3/3! + 25*x^4/4! + 96*x^5/5! + 481*x^6/6! +...
MAPLE
a := proc(n) local m: add(binomial(n, 2*m)*(n - 2*m)^(2*m), m = 0 .. floor(1/2*n - 1/2)): end proc:
seq(a(n), n = 1..30); # Petros Hadjicostas, May 06 2020 (for n >= 1)
PROG
(PARI) {a(n)=local(X=x+x*O(x^n), A=1+X); A=(exp(X*exp(X))+exp(X/exp(X)))/2; n!*polcoeff(A, n)}
(PARI) {a(n)=n!*polcoeff(sum(m=0, n, x^m*cosh(m*x+x*O(x^n))/m!), n)}
(PARI) a(n) = sum(k=0, n\2, (n-2*k)^(2*k)*binomial(n, 2*k)); \\ Seiichi Manyama, Feb 15 2023
CROSSREFS
Sequence in context: A027764 A095669 A323967 * A027949 A226547 A264167
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Sep 19 2011
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 18:04 EDT 2024. Contains 371254 sequences. (Running on oeis4.)