OFFSET
1,3
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.
This sequence contains all nonnegative integers. In particular, a(prime(n)!) = n.
EXAMPLE
2045 has prime indices {3,80} with combined divisors {1,2,3,4,5,8,10,16,20,40,80}, so a(2045) = 11. In fact, 2045 is the least number with this property.
MATHEMATICA
Table[Length[Union@@Divisors/@PrimePi/@First/@If[n==1, {}, FactorInteger[n]]], {n, 100}]
PROG
(PARI) a(n) = my(list=List(), f=factor(n)); for (i=1, #f~, fordiv(primepi(f[i, 1]), d, listput(list, d))); #Set(list); \\ Michel Marcus, May 02 2024
CROSSREFS
a(prime(n)) = A000005(n).
Positions of ones are A000079 except for 1.
a(n!) = A000720(n).
a(prime(n)!) = a(prime(A005179(n))) = n.
Counting prime factors instead of divisors gives A303975.
Positions of 2's are A371127.
A001221 counts distinct prime factors.
A003963 gives product of prime indices.
A355741 counts choices of a prime factor of each prime index.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 15 2024
STATUS
approved