OFFSET
1,1
COMMENTS
An infinite sequence.
LINKS
Bartosz Sobolewski and Maciej Ulas, Values of binary partition function represented by a sum of three squares, arXiv:2211.16622 [math.NT], 2023.
FORMULA
Each term is equal to 2*b(m) or 2*b(m)+1 for some m, where b(m) = A363769(m).
EXAMPLE
a(1)=20 because b(20)=60 is not a sum of three squares and for i=1, ..., 19, the numbers b(i), i=1,...,19 are sums of three squares, where b(i) is the number of binary partitions of n.
MATHEMATICA
bin[n_] :=
bin[n] =
If[n == 0, 1,
If[Mod[n, 2] == 0, bin[n - 1] + bin[n/2],
If[Mod[n, 2] == 1, bin[n - 1]]]];
B := {}; Do[
If[Mod[bin[n]/4^IntegerExponent[bin[n], 4], 8] == 7,
AppendTo[B, n]], {n, 1000}];
B
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Maciej Ulas, Jun 21 2023
STATUS
approved