OFFSET
0,2
COMMENTS
Refer to A161779 for the analogous sequence based on the factorials.
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..10000 (terms 0..5000 from Alois P. Heinz)
MAPLE
p:= n-> `if`(n=0, 1, ithprime(n+1)):
b:= proc(n, i) option remember; `if`(i>n, 0,
`if`(irem(n, i, 'r')=0, p(r), 0)+
add(p(j)*b(n-i*j, i+1), j=0..n/i))
end:
a:= n-> `if`(n=0, 1, b(n, 1)):
seq(a(n), n=0..45); # Alois P. Heinz, Jul 27 2019
MATHEMATICA
p[n_] := If[n==0, 1, Prime[n+1]];
b[n_, i_] := b[n, i] = If[i>n, 0, If[Mod[n, i]==0, p[n/i], 0] + Sum[p[j] b[n - i j, i+1], {j, 0, n/i}]];
a[n_] := If[n==0, 1, b[n, 1]];
a /@ Range[0, 45] (* Jean-François Alcover, Nov 20 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Gary W. Adamson, Jun 24 2009
EXTENSIONS
Definition and comment corrected by Omar E. Pol, Aug 18 2011
Correct offset and a(13)-a(36) from Alois P. Heinz, Jul 27 2019
STATUS
approved