OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
From Robert Israel, Sep 20 2016: (Start)
a(m+(n+1)*2^n) = a(m) + 4^(n+1) for 1 <= m <= (n+1)*2^n.
a(m+(n+1)*2^n) = A000695(m-(n+1)*2^n-1) + 3*4^(n+1) for (n+1)*2^n < m <= (n+2)*2^n. (End)
MAPLE
N:= 5: # to get all terms < 4^(N+1)
B0:= {0, 1}:
B1:= {3}:
for d from 1 to N do
B1:= B1 union map(`+`, B0, 3*4^d) union map(`+`, B1, 4^d);
B0:= B0 union map(`+`, B0, 4^d);
od:
sort(convert(B1, list)); # Robert Israel, Sep 20 2016
MATHEMATICA
Select[Range[400], DigitCount[#, 4, 2]==0&&DigitCount[#, 4, 3]==1&] (* Harvey P. Dale, Dec 25 2012 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
STATUS
approved