OFFSET
1,3
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = f(A000961(n)), where f(1) = 1, and f(p^e) = 2^(e-2) if p = 2 and e > 2, and f(p^e) = (p-1)*p^(e-1) otherwise.
MATHEMATICA
f[p_, e_] := If[p == 2 && e > 2, 2^(e-2), (p-1)*p^(e-1)]; s[n_] := If[n == 1, 1, If[PrimePowerQ[n], f @@ (FactorInteger[n][[1]]), Nothing]]; Array[s, 200] (* Amiram Eldar, Apr 05 2025 *)
PROG
(Haskell)
a207193 1 = 1
a207193 n | p == 2 && e > 2 = 2 ^ (e - 2)
| otherwise = (p - 1) * p ^ (e - 1)
where p = a025473 n; e = a025474 n
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Feb 16 2012
STATUS
approved
