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!)
A140053 Indices k such that A114850(m)+A114850(k) is prime for some m>k. 1
1, 1, 2, 5, 8, 7, 11, 6, 33, 14, 62, 57, 22, 7, 86, 61, 28, 70, 66, 134, 77, 131, 107, 58, 161, 252, 240, 52, 155, 32, 152, 322, 167, 200, 284, 258, 28, 173, 95, 563, 369, 57, 58, 126, 113, 369 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,3
COMMENTS
The associated primes b(n), which grow too quickly for many to be given as a sequence themselves, are {primes of the form A114850(a) + A114850(b)} = {primes of the form A114850(a) + A114850(b)} and begin as follows. b(1) = 437893890380859631 = 256 + 437893890380859375 = 4^4 + 15^15 = semiprime(1)^semiprime(1) + semiprime(6)^semiprime(6).
b(2) = 88817841970012523233890533447265881 = 256 + 88817841970012523233890533447265625 = 4^4 + 25^25 = semiprime(1)^semiprime(1) + semiprime(9)^semiprime(9).
b(3) = 46656 + 88817841970012523233890533447265625 = 6^6 + 24^25 = semiprime(2)^semiprime(2) + semiprime(9)^semiprime(9). This is to A068145 "Primes of the form a^a + b^b" as A001358 semiprimes is to A000040 primes; and as A114850 "(n-th semiprime)^(n-th semiprime)" is to A051674 "(n-th prime)^(n-th prime)."
M. F. Hasler gave the present definition which allows us to list merely the indices, which in the 3 examples above, are [6, 1],[9, 1],[9, 2]. The first 13 [m,k] value pairs are (as found by M. F. Hasler as an extension) are [6, 1], [9, 1], [9, 2], [19, 5], [20, 8], [25, 7], [33, 11], [38, 6], [40, 33], [59, 14], [69, 62], [76, 57], [99, 22]. Hence our sequence begins a(1) = 6, a(2) = 9, a(3) = 9. For the sequence of corresponding k values {6, 9, 9, 19, 20, ...}, see A140052.
LINKS
FORMULA
A001358(a(n))^A001358(a(n)) + A001358(A140052(n))^A001358(A140052(n)) is prime.
EXAMPLE
a(1) = 1 because semiprime(6)^semiprime(6) + semiprime(1)^semiprime(1) = 15^15 + 4^4 = 437893890380859375 + 256 = 437893890380859631 is prime.
MAPLE
? t=0; A001358=vector(100, i, until(bigomega(t++)==2, ); t); ? for(i=1, #A001358, for(j=1, i-1, ispseudoprime(A001358[i]^A001358[i]+A001358[j]^A001358[j]) | next; print1([i, j]", ")))
CROSSREFS
Sequence in context: A154127 A250206 A138371 * A346180 A103311 A019824
KEYWORD
more,nonn
AUTHOR
Jonathan Vos Post, May 03 2008
EXTENSIONS
a(14)-a(46) from Donovan Johnson, Nov 11 2008
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 23 06:04 EDT 2024. Contains 371906 sequences. (Running on oeis4.)