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!)
A279883 Primes of the form (prime(j)-1)^(prime(j)+1) + 1. 1
2, 17, 185302018885184100000000000000000000000000000001 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Prime terms from A279884.
If a(4) exists, it must be bigger than (prime(2200)-1)^(prime(2200)+1) + 1 = 19422^19424 + 1.
Corresponding pairs of numbers (j, prime(j)): (1, 2); (2, 3); (11, 31); ...
It is extremely unlikely that a(4) exists. The term a(1)=1^3+1 is special. For other terms, note that b^k+1 can only be prime if k is a power of 2, so say k=2^p. Otherwise, if k has an odd factor, b^k+1 is algebraically factorizable. Therefore terms a(2) and later are of form (2^p-2)^(2^p)+1 with the additional obstruction that 2^p-1 must be a prime (namely a Mersenne prime), meaning that p (the exponent of the exponent) must be prime. With generalized Fermat numbers b^(2^p)+1 already primality tested to a high limits, the first undecided possibility for a(4) is 2147483646^2147483648+1 (has j=105097565 and p=31). - Jeppe Stig Nielsen, Sep 15 2022
LINKS
PROG
(Magma) [(NthPrime(n)-1)^(NthPrime(n)+1)+1: n in[1..200] | IsPrime((NthPrime(n)-1)^(NthPrime(n)+1)+1)]
(PARI) print1(1^3+1, ", "); forprime(p=2, 19, if(isprime(2^p-1), a=(2^p-2)^(2^p)+1; ispseudoprime(a)&&print1(a, ", "))) \\ Jeppe Stig Nielsen, Sep 15 2022
CROSSREFS
Sequence in context: A112969 A208208 A290189 * A266166 A077452 A353309
KEYWORD
nonn,bref
AUTHOR
Jaroslav Krizek, Dec 21 2016
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 April 16 14:17 EDT 2024. Contains 371740 sequences. (Running on oeis4.)