%I M2896 N1161 #73 Aug 23 2024 12:17:06
%S 3,11,43,683,2731,43691,174763,2796203,715827883,2932031007403,
%T 768614336404564651,201487636602438195784363,
%U 845100400152152934331135470251,56713727820156410577229101238628035243,62357403192785191176690552862561408838653121833643
%N Wagstaff primes: primes of form (2^p + 1)/3.
%C Also, the primes with prime indices in the Jacobsthal sequence A001045.
%C Indices n such that (2^n + 1)/3 is prime are listed in A000978. - _Alexander Adamchuk_, Oct 03 2006
%C Primes in A126614. - _Omar E. Pol_, Nov 05 2013
%D N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H T. D. Noe, <a href="/A000979/b000979.txt">Table of n, a(n) for n=1..20</a>
%H P. Berrizbeitia, F. Luca, and R. Melham, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL13/Luca/luca33.html">On a compositeness test for (2^p+1)/3</a>, JIS 13 (2010) 10.1.7
%H C. Caldwell's The Top Twenty, <a href="https://t5k.org/top20/page.php?id=67">Wagstaff</a>.
%H Ernest G. Hibbs, <a href="https://www.proquest.com/openview/4012f0286b785cd732c78eb0fc6fce80">Component Interactions of the Prime Numbers</a>, Ph. D. Thesis, Capitol Technology Univ. (2022), see p. 33.
%H Editor's Note, <a href="http://www.jstor.org/stable/3029104">Table of Wagstaff primes sent by D. H. Lehmer</a>, Math. Mag., 27 (1954), 156-157.
%H Editor's Note, <a href="/A000979/a000979.pdf">Table of Wagstaff primes sent by D. H. Lehmer</a> (annotated and scanned copy)
%H Djurre G. Sikkema, <a href="https://fse.studenttheses.ub.rug.nl/33908">Probable primality testing for Wagstaff prime</a>, Bachelor's project mathematics, Univ. Groningen (Netherlands 2024). See p. 32.
%H S. S. Wagstaff, Jr., <a href="http://homes.cerias.purdue.edu/~ssw/cun/index.html">The Cunningham Project</a>.
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Wagstaff_prime">Wagstaff prime</a>
%t Select[ Array[(2^# + 1)/3 &, 190], PrimeQ] (* _Vladimir Joseph Stephan Orlovsky_, Apr 03 2010 *)
%o (Haskell)
%o a000979 n = a000979_list !! (n-1)
%o a000979_list = filter ((== 1) . a010051) a007583_list
%o -- _Reinhard Zumkeller_, Mar 24 2013
%o (Python)
%o from gmpy2 import divexact
%o from sympy import prime, isprime
%o A000979 = [p for p in (divexact(2**prime(n)+1,3) for n in range(2,10**2)) if isprime(p)] # _Chai Wah Wu_, Sep 04 2014
%o (PARI) forprime(p=2,10000,if(ispseudoprime(2^p\/3),print1(2^p\/3,","))) \\ _Edward Jiang_, Sep 05 2014
%Y Cf. A000978, A049883, A001045, A127962.
%Y Cf. A010051; subsequence of A007583.
%K nonn,changed
%O 1,1
%A _N. J. A. Sloane_
|