OFFSET
1,3
COMMENTS
a(34) > 3, and depends on the full factorization of the 296-digit composite number (139^139 - 1)/138. - Tyler Busby, Jan 22 2023
Sequence continues as ?, 8, ?, 5, 8, 4, 5, ?, 8, ?, 8, 7, 6, 3, 3, ..., where ? represents uncertain terms. - Tyler Busby, Jan 22 2023
LINKS
EXAMPLE
a(3)=2, since (5^5 - 1)/(5 - 1) = 11*71.
PROG
(PARI) a(n) = my(p=prime(n)); omega((p^p-1)/(p-1)); \\ Michel Marcus, May 22 2022
(Python)
from sympy import factorint, prime
def a(n): p = prime(n); return len(factorint((p**p-1)//(p-1)))
print([a(n) for n in range(1, 12)]) # Michael S. Branicky, May 23 2022
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Luis H. Gallardo, May 20 2022
EXTENSIONS
a(18)-a(33) from Amiram Eldar, May 20 2022
STATUS
approved