OFFSET
1,3
COMMENTS
LINKS
FORMULA
Completely multiplicative with a(prime(t)) = A228732(t). - Andrew Howroyd, Aug 01 2018
MATHEMATICA
r[n_] := FactorInteger[n][[1, 1]];
s[n_] := n/r[n];
A[n_] := A[n] = If[n==1, {1, 1}, If[PrimeOmega[n]==1, {A[PrimePi[n]][[2]], A[PrimePi[n]] // Total}, A[r[n]] * A[s[n]]]];
a[n_] := A[n][[1]];
a /@ Range[1, 80] (* Jean-François Alcover, Sep 20 2019 *)
PROG
(Haskell) see A184165.
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Emeric Deutsch and Reinhard Zumkeller, Sep 01 2013
STATUS
approved