OFFSET
1,5
COMMENTS
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..10000
FORMULA
Fully multiplicative with a(prime(n)) = A003963(n).
EXAMPLE
94 has prime indices {1,15} with prime indices {{},{2,3}} with product a(94) = 6.
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Table[Times@@Join@@primeMS/@primeMS[n], {n, 100}]
PROG
(PARI) \\ here b(n) is A003963(n).
b(n) = { my(f=factor(n)); prod(i=1, #f~, my([p, e]=f[i, ]); primepi(p)^e) }
a(n) = { my(f=factor(n)); prod(i=1, #f~, my([p, e]=f[i, ]); b(primepi(p))^e) } \\ Andrew Howroyd, Nov 10 2025
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Gus Wiseman, Mar 25 2019
STATUS
approved
