OFFSET
1,2
FORMULA
P(m) = Product[1 + t^(2*A045751(n) + 1), {n, 1, m}] a(n) = Length[CoefficientList[P(n),x]].
MATHEMATICA
a = Flatten[Table[If[PrimeQ[4*n + 1] == False, n, {}], {n, 0, 50}]]; Table[Length[CoefficientList[Product[1 + t^(2*a[[n]] + 1), {n, 1, m}], t]], { m, 0, Length[a]}]
CROSSREFS
KEYWORD
nonn,uned
AUTHOR
Roger L. Bagula, May 16 2007
STATUS
approved