OFFSET
1,2
COMMENTS
If a(n) differs from 1, then a(n) is the minimal prime divisor of A064078(n);
If a(n) > 1 then a(n) is the index where n occurs first in A014664. - M. F. Hasler, Feb 21 2016
Bang's theorem (special case of Zsigmondy's theorem, see links): a(n)>1 for all n>6. - Jeppe Stig Nielsen, Aug 31 2020
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 1..1206
Dario Alejandro Alpern, Factorization using the Elliptic Curve Method
Will Edgington, Factored Mersenne Numbers [from Internet Archive Wayback Machine]
Wikipedia, Zsigmondy's theorem
PROG
(PARI) A112927(n, f=factor(2^n-1)[, 1])=!for(i=1, #f, znorder(Mod(2, f[i]))==n&&return(f[i])) \\ Use the optional 2nd arg to give a list of pseudoprimes to try when factoring of 2^n-1 is too slow. You may try factor(2^n-1, 0)[, 1]. - M. F. Hasler, Feb 21 2016
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Vladimir Shevelev, Aug 25 2008
STATUS
approved