A171381 Numbers k such that (3^k + 1)/2 is prime.
1, 2, 4, 16, 32, 64
Offset: 1
Examples
(3^(2^0)+1)/2 = (3^1+1)/2 = 2 is prime. (3^(2^1)+1)/2 = (3^2+1)/2 = 5 is prime. (3^(2^2)+1)/2 = (3^4+1)/2 = 41 is prime. (3^(2^3)+1)/2 = (3^8+1)/2 = 3281 is divisible by 17=1+2^4. (3^(2^4)+1)/2 = (3^16+1)/2 = 21523361 is prime. (3^(2^5)+1)/2 = (3^32+1)/2 = 926510094425921 is prime. (3^(2^6)+1)/2 = (3^64+1)/2 = 1716841910146256242328924544641 is prime. (3^(2^7)+1)/2 = (3^128+1)/2 is divisible by 257=1+2^8, so 2^7 is not a term. (3^(2^8)+1)/2 = (3^256+1)/2 is divisible by 1+2^9*24, so 2^8 is not a term. (3^(2^15)+1)/2 is divisible by 2^(2^4)+1, so 2^15 is not a term. - _Georgi Guninski_, Jun 13 2010 (3^(2^19)+1)/2 is divisible by 13631489, so 2^19 is not a term. - _Paul Zimmermann_, Jun 14 2010 (3^(2^20)+1)/2 is 5-composite so 2^20 is not a term. - _Serge Batalov_, Jun 14 2010 According to PFGW, 2^20 is not in the sequence: PFGW Version 3.3.4.20100405.x86_Stable [GWNUM 25.14] (3^1048576+1)/2 is composite: RES64: [9EE4CA1AABB9A816] (3229 sec) Base 5. Base 3 is useless here [cf. comment by _Don Reble_ - Ed.]. - _Georgi Guninski_, Jun 15 2010 2^23 is not in the sequence (listed as "Composite but no factor known" on the second Keller link). - _Serge Batalov_, Jun 18 2010 Verified 2^23 yields a composite: base 2 (PFGWv3.3.1). - _Paul Bourdelais_, Apr 25 2011
Links
- Anders Björn and Hans Riesel, Factors of generalized Fermat numbers, Math. Comp. 67 (1998), no. 221, pp. 441-446.
- I. J. Calvo, A note on factors of generalized Fermats numbers, Applied Math. Letters 13, (2000), pp. 1-5. [Gives divisibility criteria for 3^(2^m)+1 by primes of the form p=3*2^n+1 (p=7, 13, 97, 193 ...) (Theorem 2.1) and for primes of this form when they are divisors of Fermat numbers (Theorem 2.2)]
- M. F. Hasler and G. Guninski, Eliminating some further terms
- W. Keller, Factors of generalized Fermat numbers found after Bjorn & Riesel
- W. Keller, Prime factors of generalized Fermat numbers F'_m(3) = F_m(3) / 2 and complete factoring status
Programs
-
Magma
IsPrime((3^(2^15)+1) div 2); // shows that 15 is not a term - Jon E. Schoenfield, Jun 13 2010
-
PARI
is_A171381(n)=ispseudoprime(3^n\2+1) \\ M. F. Hasler, Oct 02 2012
Extensions
Edited by N. J. A. Sloane, Dec 09 2009
Incorrect terms a(7)-a(15) deleted by Jon E. Schoenfield, Jun 12 2010
The next term, if it exists, is at least 2^19. - Georgi Guninski, Jun 13 2010
k=2^21 yields a number divisible by 1+2^22*37. - M. F. Hasler, Jun 14 2010
Edited by N. J. A. Sloane, Jun 12 2010 - Jun 16 2010
Edited by M. F. Hasler, Oct 02 2012
Comments