OFFSET
1,6
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
FORMULA
EXAMPLE
The prime indices of 6500 are {1,1,3,3,3,6} with fixed points at positions {1,3,6}, so a(6500) = 3.
MAPLE
f:= proc(n) local F, J, t;
F:= sort(ifactors(n)[2], (s, t) -> s[1]<t[1]);
J:= map(t -> numtheory:-pi(t[1])$t[2], F);
nops(select(t -> J[t]=t, [$1..nops(J)]));
end proc:
map(f, [$1..200]); # Robert Israel, Apr 11 2023
MATHEMATICA
pq[y_]:=Length[Select[Range[Length[y]], #==y[[#]]&]];
Table[pq[Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]], {n, 100}]
PROG
(PARI) A352822(n) = { my(f=factor(n), i=0, c=0); for(k=1, #f~, while(f[k, 2], f[k, 2]--; i++; c += (i==primepi(f[k, 1])))); (c); }; \\ Antti Karttunen, Apr 11 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Apr 05 2022
EXTENSIONS
Data section extended up to 105 terms by Antti Karttunen, Apr 11 2022
STATUS
approved