OFFSET
0,2
LINKS
Eric Weisstein's World of Mathematics, Primorial.
FORMULA
G.f.: Product_{k>=1} 1 / (1 - x^k)^prime(k)#.
MAPLE
p:= proc(n) option remember; `if`(n<1, 1, p(n-1)*ithprime(n)) end:
a:= proc(n) option remember; `if`(n=0, 1, add(a(n-j)*
add(d*p(d), d=numtheory[divisors](j)), j=1..n)/n)
end:
seq(a(n), n=0..20); # Alois P. Heinz, Jan 25 2025
MATHEMATICA
nmax = 20; CoefficientList[Series[Product[1/(1 - x^k)^Product[Prime[j], {j, k}], {k, 1, nmax}], {x, 0, nmax}], x]
primorial[n_] := Product[Prime[j], {j, 1, n}]; a[n_] := a[n] = If[n == 0, 1, Sum[Sum[d primorial[d], {d, Divisors[j]}] a[n - j], {j, 1, n}]/n]; Table[a[n], {n, 0, 20}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 25 2025
STATUS
approved