OFFSET
1,5
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..10000
EXAMPLE
We have 91 = prime(4) * prime(6), so a(91) = pi(4) * pi(6) = 6.
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Table[Times@@PrimePi/@primeMS[n], {n, 100}]
PROG
(PARI) a(n) = { my(f=factor(n)); prod(i=1, #f~, my([p, e]=f[i, ]); primepi(primepi(p))^e) } \\ Andrew Howroyd, Nov 10 2025
CROSSREFS
KEYWORD
nonn,mult,easy
AUTHOR
Gus Wiseman, Oct 25 2022
STATUS
approved
