login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A306946 Expansion of e.g.f. Product_{k>=1} 1/(1 - x^(k^2)/k^2). 0
1, 1, 2, 6, 30, 150, 900, 6300, 52920, 516600, 5166000, 56826000, 689396400, 9135126000, 127891764000, 1918376460000, 32083427376000, 547888316976000, 9941031356256000, 188879595768864000, 3817981890122400000, 80266319244951840000, 1769328161437636800000 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

LINKS

Table of n, a(n) for n=0..22.

FORMULA

a(n) ~ 2 * n!. - Vaclav Kotesovec, Mar 18 2019

MATHEMATICA

nmax = 22; CoefficientList[Series[Product[1/(1 - x^k^2/k^2), {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]!

nmax = 22; CoefficientList[Series[Exp[Sum[Sum[Boole[IntegerQ[d^(1/2)]] d^(1 - k/d), {d, Divisors[k]}] x^k/k, {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]!

a[n_] := a[n] = If[n == 0, 1, (n - 1)! Sum[Sum[Boole[IntegerQ[d^(1/2)]] d^(1 - k/d), {d, Divisors[k]}] a[n - k]/(n - k)!, {k, 1, n}]]; Table[a[n], {n, 0, 22}]

CROSSREFS

Cf. A007841, A319104.

Sequence in context: A122763 A166078 A319365 * A319104 A326756 A327927

Adjacent sequences:  A306943 A306944 A306945 * A306947 A306948 A306949

KEYWORD

nonn

AUTHOR

Ilya Gutkovskiy, Mar 17 2019

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 18 16:27 EST 2020. Contains 331011 sequences. (Running on oeis4.)