%I A000978 M2413 N0956 #158 Oct 04 2026 13:21:45 %S A000978 3,5,7,11,13,17,19,23,31,43,61,79,101,127,167,191,199,313,347,701, %T A000978 1709,2617,3539,5807,10501,10691,11279,12391,14479,42737,83339,95369, %U A000978 117239,127031,138937,141079,267017,269987,374321,986191,4031399 %N A000978 Wagstaff numbers: numbers k such that (2^k + 1)/3 is prime. %C A000978 It is easy to see that the definition implies that k must be an odd prime. - _N. J. A. Sloane_, Oct 06 2006 %C A000978 The terms from a(32) on only give probable primes as of 2018. Caldwell lists the largest certified primes. - _Jens Kruse Andersen_, Jan 10 2018 %C A000978 Prime numbers of the form 1+Sum_{i=1..m} 2^(2i-1). - _Artur Jasinski_, Feb 09 2007 %C A000978 There is a new conjecture stating that a Wagstaff number is prime under the following condition (based on DiGraph cycles under the LLT): Let p be a prime integer > 3, N(p) = 2^p+1 and W(p) = N(p)/3, S(0) = 3/2 (or 1/4) and S(i+1) = S(i)^2 - 2 (mod N(p)). Then W(p) is prime iff S(p-1) == S(0) (mod W(p)). - _Tony Reix_, Sep 03 2007 %C A000978 As a member of the DUR team (Diepeveen, Underwood, Reix), and thanks to the LLR tool built by Jean Penne, I've found a new and big Wagstaff PRP: (2^4031399+1)/3 is Vrba-Reix PRP! This Wagstaff number has 1,213,572 digits and today is the 3rd biggest PRP ever found. I've done a second verification on a Nehalem core with the PFGW tool. - _Tony Reix_, Feb 20 2010 %C A000978 13347311 and 13372531 were found to be terms of this sequence (maybe not the next ones) by _Ryan Propper_ in September 2013. - _Max Alekseyev_, Oct 07 2013 %C A000978 The next term is larger than 10 million. - _Gord Palameta_, Mar 22 2019 %C A000978 Ryan Propper found another likely term, 15135397, though it only corresponds to a probable prime. - _Charles R Greathouse IV_, Jul 01 2021 %C A000978 Conjecture: a(1) = 3; and then, for n > 0, a(n+1) is the least odd k such that lpf((2^k + 1)/3) > lpf((2^a(n) + 1)/3), where lpf(w) is the least prime factor of w. - _Thomas Ordowski_, Sep 13 2026 %D A000978 J. Brillhart et al., Factorizations of b^n +- 1. Contemporary Mathematics, Vol. 22, Amer. Math. Soc., Providence, RI, 2nd edition, 1985; and later supplements. %D A000978 N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence). %D A000978 N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence). %D A000978 S. S. Wagstaff, Jr., personal communication. %H A000978 John D. Brillhart, Derrick H. Lehmer, John L. Selfridge, Bryant Tuckerman, and Samuel S. Wagstaff, Factorizations of b^n +- 1, Contemp. Math. 22, Amer. Math. Soc., Providence, RI, 3rd edition, 2002. %H A000978 Chris Caldwell's The Top Twenty, Wagstaff. %H A000978 Chris Caldwell, New Mersenne Conjecture. %H A000978 Alexey Dolotov, Three Brillhart-Lehmer-Selfridge primality proofs for Wagstaff numbers, arXiv:2605.18555 [math.NT], 2026. %H A000978 Harvey Dubner, Generalized repunit primes, Math. Comp., 61 (1993), 927-930. [Annotated scanned copy] %H A000978 Harvey Dubner and Torbjörn Granlund, Primes of the Form (b^n+1)/(b+1), J. Int. Seq. 3 (2000), Art. P00.2.7. %H A000978 Editor's Note, Table of Wagstaff primes sent by D. H. Lehmer (annotated and scanned copy) %H A000978 J. E. Foster, Problem 174, A solution in primes, Math. Mag., 27 (1954), 156-157. %H A000978 Richard K. Guy, Letter to N. J. A. Sloane, Aug 1986 %H A000978 Richard K. Guy, Letter to N. J. A. Sloane, 1987 %H A000978 Henri Lifchitz, Mersenne and Fermat primes field. %H A000978 Henri Lifchitz and Renaud Lifchitz, PRP Top Records. %H A000978 Henri Lifchitz and Renaud Lifchitz, PRP Records. %H A000978 Gord Palameta, There are no new Wagstaff primes with exponent below 10 million, mersenneforum.org %H A000978 Ryan Propper et al., New Wagstaff PRP exponents, mersenneforum.org %H A000978 PRP top list: PRP top [From _Tony Reix_, Feb 20 2010] %H A000978 Tony Reix, Yahoo Primeform Group Message 10184 dd. Feb 20, 2010, reconstruction in html. %H A000978 Tony Reix, Some Maths about the Vrba-Reix PRP test [From _Tony Reix_, Feb 20 2010] %H A000978 Djurre G. Sikkema, Probable primality testing for Wagstaff prime, Bachelor's project mathematics, Univ. Groningen (Netherlands 2024). See p. 32. %H A000978 Samuel S. Wagstaff, Jr., The Cunningham Project. %H A000978 Samuel S. Wagstaff, Jr., Two Mersenne Prime Conjectures, J. Int. Seq. 28 (2025), Art. 25.7.2. See p. 2. %H A000978 Eric Weisstein's World of Mathematics, Integer Sequence Primes. %H A000978 Eric Weisstein's World of Mathematics, Repunit. %H A000978 Eric Weisstein's World of Mathematics, Repunit Prime. %H A000978 Eric Weisstein's World of Mathematics, Wagstaff Prime. %H A000978 Wikipedia, Wagstaff prime %H A000978 Robert G. Wilson v, Letter to N. J. A. Sloane, circa 1991. %F A000978 a(n) = A107036(n) for n>1. - _Alexander Adamchuk_, Feb 10 2007 %t A000978 Select[Range[5000], PrimeQ[(2^# + 1)/3] &] (* _Michael De Vlieger_, Jan 10 2018 *) %t A000978 (* Alternative: *) %t A000978 Select[Prime[Range[2,500]],PrimeQ[(2^#+1)/3]&] (* _Harvey P. Dale_, Jun 13 2022 *) %o A000978 (PARI) forprime(p=2,5000,if(ispseudoprime(2^p\/3),print1(p", "))) \\ _Charles R Greathouse IV_, Jul 15 2011 %o A000978 (Haskell) %o A000978 a000978 n = a000978_list !! (n-1) %o A000978 a000978_list = filter ((== 1) . a010051 . a001045) a065091_list %o A000978 -- _Reinhard Zumkeller_, Mar 24 2013 %o A000978 (Python) %o A000978 from gmpy2 import divexact %o A000978 from sympy import prime, isprime %o A000978 A000978 = [p for p in (prime(n) for n in range(2,10**2)) if isprime(divexact(2**p+1,3))] # _Chai Wah Wu_, Sep 04 2014 %Y A000978 Cf. A107036 (indices of prime Jacobsthal numbers). %Y A000978 Cf. A000979, A124400, A124401, A127955, A127956, A127957, A127958, A127936, A010051, A065091, A001045. %K A000978 nonn,hard,nice,more,changed %O A000978 1,1 %A A000978 _N. J. A. Sloane_, _Robert G. Wilson v_ %E A000978 a(30) from Kamil Duszenko (kdusz(AT)wp.pl), Feb 03 2003; a(30) was proved prime by Francois Morain with FastECPP. - _Tony Reix_, Sep 03 2007 %E A000978 a(31)-a(39) from _Robert G. Wilson v_, Apr 11 2005 %E A000978 a(40) from Vincent Diepeveen (diep(AT)xs4all.nl) added by _Alexander Adamchuk_, Jun 19 2008 %E A000978 a(41) from _Tony Reix_, Feb 20 2010