login
a(n) = 1 + Product_{d|n, d < n} a(d).
3

%I #25 Apr 14 2021 04:18:53

%S 2,3,3,7,3,19,3,43,7,19,3,2395,3,19,19,1807,3,2395,3,2395,19,19,3,

%T 246546091,7,19,43,2395,3,370387,3,3263443,19,19,19,96124306951,3,19,

%U 19,246546091,3,370387,3,2395,2395,19,3,109838449356687381331,7,2395,19,2395,3,246546091,19

%N a(n) = 1 + Product_{d|n, d < n} a(d).

%H Seiichi Manyama, <a href="/A343390/b343390.txt">Table of n, a(n) for n = 1..959</a>

%F a(p^k) = A000058(k) for p prime.

%p a:= proc(n) option remember;

%p 1+mul(a(d), d=numtheory[divisors](n) minus {n})

%p end:

%p seq(a(n), n=1..55); # _Alois P. Heinz_, Apr 13 2021

%t a[n_] := a[n] = 1 + Product[If[d < n, a[d], 1], {d, Divisors[n]}]; Table[a[n], {n, 55}]

%Y Cf. A000058, A067824, A068334.

%K nonn

%O 1,1

%A _Ilya Gutkovskiy_, Apr 13 2021