OFFSET
1,2
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
KEYWORD
nonn
AUTHOR
David S. Newman, Aug 28 2014
STATUS
approved