OFFSET
1,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..5000
FORMULA
G.f.: Sum(prime(k)!*x^(1/2*prime(k)^2+1/2*prime(k))/Product(1-x^j, j = 1 .. prime(k)), k = 1 .. infinity).
MAPLE
b:= proc(n, i) option remember; `if`(n=0, [1],
`if`(n>i*(i+1)/2, [], zip((x, y)->x+y, b(n, i-1),
`if`(i>n, [], [0, b(n-i, i-1)[]]), 0)))
end:
a:= proc(n) local l; l:= b(n$2);
add(`if`(isprime(i), l[i+1]*i!, 0), i=2..nops(l)-1)
end:
seq(a(n), n=1..50); # Alois P. Heinz, Nov 20 2012
MATHEMATICA
CoefficientList[ Series[ Sum[ Prime[k]!* x^(Prime[k]^2/2 + Prime[k]/2)/Product[1 - x^j, {j, Prime[k]}], {k, 44}], {x, 0, 44}], x] (* Robert G. Wilson v, Feb 04 2005 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Vladeta Jovovic, Jan 31 2005
EXTENSIONS
More terms from Robert G. Wilson v, Feb 04 2005
STATUS
approved