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!)
A294394 E.g.f.: exp(Sum_{n>=1} A000593(n) * x^n). 4
1, 1, 3, 31, 145, 1641, 17731, 194503, 2676801, 40644145, 667689571, 11514903951, 227665389073, 4578990563161, 100913115588195, 2372334731747191, 57930324367791361, 1509398686720812513, 41341036374519788611, 1184009909077133031295 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

LINKS

Seiichi Manyama, Table of n, a(n) for n = 0..434

FORMULA

a(0) = 1 and a(n) = (n-1)! * Sum_{k=1..n} k*A000593(k)*a(n-k)/(n-k)! for n > 0.

a(n) ~ Pi^(1/3) * exp((3*Pi)^(2/3) * n^(2/3) / 2^(4/3) - 1/24 - n) * n^(n - 1/6) / (2^(1/6) * 3^(2/3)). - Vaclav Kotesovec, Sep 07 2018

MATHEMATICA

a[n_] := a[n] = If[n == 0, 1, Sum[k*Sum[-(-1)^d*k/d, {d, Divisors[k]}]*a[n - k], {k, 1, n}]/n]; Table[n!*a[n], {n, 0, 20}] (* Vaclav Kotesovec, Sep 07 2018 *)

PROG

(PARI) N=66; x='x+O('x^N); Vec(serlaplace(exp(sum(k=1, N, sumdiv(k, d, d*(d%2))*x^k))))

CROSSREFS

E.g.f.: exp(Sum_{n>=1} (Sum_{d|n and d is odd} d^k) * x^n): A294392 (k=0), this sequence (k=1), A294395 (k=2).

Cf. A000009, A000593, A301799, A301800.

Sequence in context: A069630 A069615 A087389 * A064055 A227049 A078604

Adjacent sequences:  A294391 A294392 A294393 * A294395 A294396 A294397

KEYWORD

nonn

AUTHOR

Seiichi Manyama, Oct 30 2017

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 28 06:43 EST 2022. Contains 350654 sequences. (Running on oeis4.)