OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..6
MAPLE
b:= proc(n) b(n):= `if`(n=0, 1, b(n-1)*ithprime(n)) end:
g:= proc(n) option remember; `if`(n=0, 1, add(g(n-j)*add(
`if`(d::odd, d, 0), d=numtheory[divisors](j)), j=1..n)/n)
end:
a:= n-> g(b(n)):
seq(a(n), n=0..5);
MATHEMATICA
$RecursionLimit = 2^13;
b[n_] := b[n] = If[n == 0, 1, b[n - 1]*Prime[n]];
g[n_] := g[n] = If[n == 0, 1, Sum[g[n - j]*Sum[
If[OddQ[d], d, 0], {d, Divisors[j]}], {j, 1, n}]/n];
a[n_] := g[b[n]];
Table[a[n], {n, 0, 5}] (* Jean-François Alcover, May 02 2022, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Apr 06 2021
STATUS
approved