OFFSET
0,10
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000
FORMULA
G.f.: Product_{i>0} (1+x^i)/(1+x^prime(i)).
MAPLE
b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
b(n, i-1)+ `if`(i>n or isprime(i), 0, b(n-i, i-1))))
end:
a:= n-> b(n$2):
seq(a(n), n=0..70); # Alois P. Heinz, May 29 2013
MATHEMATICA
CoefficientList[ Series[ Product[(1 + x^n)/(1 + x^Prime[n]), {n, 70}], {x, 0, 67}], x] (* Robert G. Wilson v, Aug 02 2004 *)
PROG
(Haskell)
a096258 = p a018252_list where
p _ 0 = 1
p (k:ks) m = if m < k then 0 else p ks (m - k) + p ks m
-- Reinhard Zumkeller, Jan 15 2012
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Vladeta Jovovic, Jul 31 2004
EXTENSIONS
More terms from Robert G. Wilson v, Aug 02 2004
STATUS
approved