OFFSET
1,2
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
David A. Corneth, Table of n, a(n) for n = 1..10000 (first 3000 terms from Amiram Eldar)
FORMULA
Sum_{n>=1} 1/a(n) = Product_{p prime > 3} (1 + 1/(p^2*(p-1))) = (72/95)*A065483 = 1.0154153584... . - Amiram Eldar, May 28 2022
EXAMPLE
The initial terms together with their prime indices:
1: {}
125: {3,3,3}
343: {4,4,4}
625: {3,3,3,3}
1331: {5,5,5}
2197: {6,6,6}
2401: {4,4,4,4}
3125: {3,3,3,3,3}
4913: {7,7,7}
6859: {8,8,8}
12167: {9,9,9}
14641: {5,5,5,5}
15625: {3,3,3,3,3,3}
16807: {4,4,4,4,4}
24389: {10,10,10}
28561: {6,6,6,6}
29791: {11,11,11}
42875: {3,3,3,4,4,4}
MATHEMATICA
Select[Range[10000], #==1||!MemberQ[FactorInteger[#], {_?(#<5&), _}|{_, _?(#<3&)}]&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, May 16 2022
STATUS
approved