cp's OEIS Frontend

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.

%I A246515 #11 Sep 14 2014 00:25:14
%S A246515 1,2,4,9,10,11,12,13,14,15,18,19,22,23,25,26,27,30,37,38,40,41,42,44,
%T A246515 45,48,51,52,53,54,55,56,58,59,64,67,68,70,71,72,73,74,79,82,83,84,86,
%U A246515 88,89,93
%N 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.
%t A246515 iend=150;
%t A246515 aim=CoefficientList[Sum[x^(i^3),{i,0,iend^(1/3)+1}],x];
%t A246515 c[1]=1;
%t A246515 index=1;
%t A246515 For[i=1,i<iend,i++,
%t A246515 ser=Mod[CoefficientList[Series[1/( Product[(1-x^c[j]),{j,1,index}] ),{x,0,iend+1}],x],2];
%t A246515 If[ser[[i+2]]!=aim[[i+2]],c[index+1]=i+1;index++]
%t A246515 ];
%t A246515 Table[c[i],{i,1,50}]
%K A246515 nonn
%O A246515 1,2
%A A246515 _David S. Newman_, Aug 28 2014