OFFSET
1,1
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..959
FORMULA
a(p^k) = A000058(k) for p prime.
MAPLE
a:= proc(n) option remember;
1+mul(a(d), d=numtheory[divisors](n) minus {n})
end:
seq(a(n), n=1..55); # Alois P. Heinz, Apr 13 2021
MATHEMATICA
a[n_] := a[n] = 1 + Product[If[d < n, a[d], 1], {d, Divisors[n]}]; Table[a[n], {n, 55}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Apr 13 2021
STATUS
approved