A246515 a(n) such that the series 1/((1-x^a(1))(1-x^a(2))...)=1+x+x^8+x^27+...+x^(j^3)+... when read mod 2.

1, 2, 4, 9, 10, 11, 12, 13, 14, 15, 18, 19, 22, 23, 25, 26, 27, 30, 37, 38, 40, 41, 42, 44, 45, 48, 51, 52, 53, 54, 55, 56, 58, 59, 64, 67, 68, 70, 71, 72, 73, 74, 79, 82, 83, 84, 86, 88, 89, 93

Offset

1,2

Links

Table of n, a(n) for n=1..50.

Mathematica

iend=150;
aim=CoefficientList[Sum[x^(i^3), {i, 0, iend^(1/3)+1}], x];
c[1]=1;
index=1;
For[i=1, i<iend, i++,
ser=Mod[CoefficientList[Series[1/( Product[(1-x^c[j]), {j, 1, index}] ), {x, 0, iend+1}], x], 2];
If[ser[[i+2]]!=aim[[i+2]], c[index+1]=i+1; index++]
];
Table[c[i], {i, 1, 50}]

Crossrefs

Sequence in context: A279035 A076661 A258710 * A275658 A267414 A072583

Adjacent sequences: A246512 A246513 A246514 * A246516 A246517 A246518

Keyword

nonn

Author

David S. Newman, Aug 28 2014

Status

approved