OFFSET
0,5
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..16383
FORMULA
G.f.: Product_{k>=0} (1-x^(6*4^k))/(1-x^(4^k)).
G.f. A(x) satisfies: A(x) = (1 + x + x^2 + x^3 + x^4 + x^5) * A(x^4). - Ilya Gutkovskiy, Jul 09 2019
EXAMPLE
a(17) = 3 because 17=16+1=4+4+4+4+1=4+4+4+1+1+1+1+1.
MAPLE
b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<0, 0,
add(b(n-j*4^i, i-1), j=0..min(5, n/4^i))))
end:
a:= n-> b(n, ilog[4](n)):
seq(a(n), n=0..120); # Alois P. Heinz, May 01 2018
MATHEMATICA
m = 100; A[_] = 1;
Do[A[x_] = (1+x+x^2+x^3+x^4+x^5) * A[x^4] + O[x]^m // Normal, {m}];
CoefficientList[A[x], x] (* Jean-François Alcover, Oct 06 2019, after Ilya Gutkovskiy *)
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Seiichi Manyama, May 01 2018
STATUS
approved