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!)
A304909 Expansion of x * (d/dx) Product_{k>=0} 1/(1 - x^(2^k)). 1
0, 1, 4, 6, 16, 20, 36, 42, 80, 90, 140, 154, 240, 260, 364, 390, 576, 612, 828, 874, 1200, 1260, 1628, 1702, 2256, 2350, 2964, 3078, 3920, 4060, 4980, 5146, 6464, 6666, 8092, 8330, 10224, 10508, 12540, 12870, 15600, 15990, 18900, 19350, 23056, 23580, 27508, 28106, 33216, 33908 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
COMMENTS
Sum of all parts of all partitions of n into powers of 2.
Convolution of the sequences A018819 and A038712.
LINKS
FORMULA
G.f.: x * (d/dx) Product_{k>=0} (1 + x^(2^k))^(k+1).
G.f.: Sum_{i>=0} 2^i*x^(2^i)/(1 - x^(2^i)) * Product_{j>=0} 1/(1 - x^(2^j)).
a(n) = n*A018819(k).
MATHEMATICA
nmax = 49; CoefficientList[Series[x D[Product[1/(1 - x^2^k), {k, 0, Floor[Log[nmax]/Log[2]] + 1}], x], {x, 0, nmax}], x]
nmax = 49; CoefficientList[Series[x D[Product[(1 + x^2^k)^(k + 1), {k, 0, Floor[Log[nmax]/Log[2]] + 1}], x], {x, 0, nmax}], x]
a[n_] := a[n] = If[n == 0, 1, Sum[Sum[d Boole[d == 2^IntegerExponent[d, 2]], {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[n a[n], {n, 0, 49}]
CROSSREFS
Sequence in context: A117988 A263636 A238411 * A302119 A120542 A132344
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, May 20 2018
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 April 26 21:53 EDT 2024. Contains 372004 sequences. (Running on oeis4.)