nn=10;邮编:A279484=Rest[CoefficientList[系列[产品[(1-x^(k^3)),{k,nn}],{x,0,nn^3}],x]];选择[范围[nn^3],邮编:A279484[[#]]==0&]
nmax=1000;nn=Floor[nmax^(1/3)]+1;poly=ConstantArray[0,nn^3+1];poly[[1]]=1;poly[[3]]=0;Do[Do[poly[[j+1]]-=poly[[j-k^3+1]],{j,nn^3,k^3,-1}];,{k,2,nn}];邮编:A279484=取[poly,{2,nmax+1}];选择[Range[nmax],邮编:A279484[[#]]==0&]
|