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!)
A118395 Expansion of e.g.f. exp(x + x^3). 8
1, 1, 1, 7, 25, 61, 481, 2731, 10417, 91225, 681121, 3493711, 33597961, 303321877, 1938378625, 20282865331, 211375647841, 1607008257841, 18157826367937, 212200671085975, 1860991143630841, 22560913203079021, 289933758771407521, 2869267483843753147, 37116733726117707025 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,4
COMMENTS
Equals row sums of triangle A118394.
LINKS
FORMULA
E.g.f.: 1 + x/(1+x)*(G(0) - 1) where G(k) = 1 + (1+x^2)/(k+1)/(1-x/(x+(1)/G(k+1) )), recursively defined continued fraction. - Sergei N. Gladkovskii, Feb 04 2013
a(n) ~ 3^(n/3-1/2) * n^(2*n/3) * exp((n/3)^(1/3)-2*n/3). - Vaclav Kotesovec, Jun 02 2013
E.g.f.: A(x) = exp(x+x^3) satisfies A' - (1+3*x^2)*A = 0. - Gheorghe Coserea, Aug 24 2015
a(n+1) = a(n) + 3*n*(n-1)*a(n-2). - Gheorghe Coserea, Aug 24 2015
a(n) = n! * Sum_{k=0..floor(n/3)} binomial(n-2*k,k)/(n-2*k)!. - Seiichi Manyama, Feb 25 2022
MAPLE
with(combstruct):seq(count(([S, {S=Set(Union(Z, Prod(Z, Z, Z)))}, labeled], size=n)), n=0..22); # Zerinvary Lajos, Mar 18 2008
MATHEMATICA
CoefficientList[Series[E^(x+x^3), {x, 0, 20}], x]* Range[0, 20]! (* Vaclav Kotesovec, Jun 02 2013 *)
T[n_, k_] := n!/(k!(n-3k)!);
a[n_] := Sum[T[n, k], {k, 0, Floor[n/3]}];
a /@ Range[0, 24] (* Jean-François Alcover, Nov 04 2020 *)
PROG
(PARI) a(n)=n!*polcoeff(exp(x+x^3+x*O(x^n)), n)
(PARI)
N=33; x='x+O('x^N);
egf=exp(x+x^3);
Vec(serlaplace(egf))
/* Joerg Arndt, Sep 15 2012 */
(PARI) a(n) = n!*sum(k=0, n\3, binomial(n-2*k, k)/(n-2*k)!); \\ Seiichi Manyama, Feb 25 2022
(Magma) [n le 3 select 1 else Self(n-1) + 3*(n-2)*(n-3)*Self(n-3): n in [1..26]]; // Vincenzo Librandi, Aug 25 2015
(Sage)
def a(n):
if (n<3): return 1
else: return a(n-1) + 3*(n-1)*(n-2)*a(n-3)
[a(n) for n in (0..25)] # G. C. Greubel, Feb 18 2021
CROSSREFS
Sequence in context: A033814 A321165 A344560 * A118396 A330044 A193375
KEYWORD
nonn
AUTHOR
Paul D. Hanna, May 07 2006
EXTENSIONS
Missing a(0)=1 prepended by Joerg Arndt, Sep 15 2012
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 29 03:51 EDT 2024. Contains 371264 sequences. (Running on oeis4.)