login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A246005 Least k such that ((2n+1)^k-1)/2n is prime, or 0 if no such k exists. 1
3, 3, 5, 0, 17, 5, 3, 3, 19, 3, 5, 0, 3, 5, 7, 3, 313, 13, 349, 3, 5, 19, 127, 0, 4229, 11, 17, 3, 3, 7, 5, 19, 19, 3, 3, 5, 3, 3, 5, 0, 5, 5, 7, 3, 4421, 7, 7, 17, 3, 3, 19, 3, 17, 17, 3, 23, 7, 3, 3, 0, 43, 0, 5, 5, 3, 13, 1171, 11, 163, 3, 3, 5, 3, 7, 13, 3, 3, 17, 13, 3, 7, 5, 3, 0, 181, 3, 5, 5, 19, 17, 223 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
a(92) > 10000, a(93)..a(133) = {37, 3, 17, 5, 11, 31, 577, 271, 3, 19, 13, 3, 41, 137, 3, 281, 13, 7, 239, 0, 5, 11, 3, 113, 7, 7, 5, 17, 0, 3, 17, 5, 7, 19, 5, 23, 2011, 31, 5, 5, 13}, a(134) > 10000, a(135)..a(139) = {41, 37, 5, 5, 3}, a(140) > 10000, a(141)..a(150) = {29, 5, 3, 0, 13, 3, 17, 17, 113, 193}.
LINKS
FORMULA
a(n) = A084740(2n+1).
EXAMPLE
a(23) = 127 because 2 * 23 + 1 = 47, (47^k-1)/46 is composite for k = 2, 3, ..., 126 and prime for k = 127.
PROG
(PARI) a(n) = {l=List([4, 12, 24, 40, 60, 62, 84]); for(q=1, 91, if(n==l[q], return(0))); k=1; while(k, s=((2*n+1)^prime(k)-1)/(2*n); if(ispseudoprime(s), return(prime(k))); k++)} \\ Eric Chen, Nov 14 2014
CROSSREFS
Sequence in context: A197137 A133456 A271710 * A103786 A067462 A320298
KEYWORD
nonn
AUTHOR
Eric Chen, Nov 13 2014
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 18 22:56 EDT 2024. Contains 370952 sequences. (Running on oeis4.)