OFFSET
0,19
LINKS
Sean A. Irvine, Java program (github)
FORMULA
G.f.: x^2*Product_{k>1} (1+x^prime(k)). - Vladeta Jovovic, Jul 20 2003
MAPLE
with(numtheory):
b:= proc(n, i) option remember; local g;
if n=0 then [1, 0]
elif i<1 then [0, 0]
else g:= `if`(ithprime(i)>n, [0$2], b(n-ithprime(i), i-1));
b(n, i-1) +g +[0, `if`(i=1, g[1], 0)]
fi
end:
a:= n-> b(n, pi(n))[2]:
seq(a(n), n=0..80); # Sean A. Irvine, after Alois P. Heinz, Jul 29 2019
MATHEMATICA
max = 100; (* number of terms *)
CoefficientList[x^2*Product[1+x^Prime[k], {k, 2, PrimePi[max]}]+O[x]^max, x] (* Jean-François Alcover, Sep 07 2022, after Vladeta Jovovic *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Vladeta Jovovic, Jul 20 2003
a(0)-a(6) prepended by Sean A. Irvine, Jul 29 2019
STATUS
approved