The OEIS is supported by the many generous donors to the OEIS Foundation.

 Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 60, we have over 367,000 sequences, and we’ve crossed 11,000 citations (which often say “discovered thanks to the OEIS”). Other ways to Give
 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A096012 Numbers k such that k^2+1 and (k+2)^2+1 are both prime; twin k^2+1 primes. 14
 2, 4, 14, 24, 54, 124, 204, 384, 464, 634, 644, 714, 1094, 1144, 1174, 1244, 1274, 1314, 1374, 1564, 1614, 1674, 1684, 1964, 2054, 2084, 2094, 2404, 2454, 2534, 2664, 2834, 2924, 3134, 3304, 3534, 3754, 3774, 4024, 4154, 4174, 4364, 4604, 4614, 4734, 4784 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 LINKS Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Seiichi Manyama) FORMULA a(k) = A108814(k) - 1. - Jeppe Stig Nielsen, Feb 26 2016 MATHEMATICA Select[Range[5000], AllTrue[{#^2+1, (#+2)^2+1}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Jul 23 2014 *) Select[Range[5000], PrimeQ[#^2 + 1] && PrimeQ[(# + 2)^2 + 1] &] (* Vincenzo Librandi, Feb 27 2016 *) PROG (Magma) [n: n in [1..5000] | IsPrime(n^2+1) and IsPrime((n+2)^2+1)]; // Vincenzo Librandi, Feb 27 2016 (PARI) isok(n) = isprime(n^2+1) && isprime((n+2)^2+1); \\ Michel Marcus, Feb 27 2016 CROSSREFS Cf. A005574, A108814. Sequence in context: A283379 A090889 A350050 * A261562 A102930 A357996 Adjacent sequences: A096009 A096010 A096011 * A096013 A096014 A096015 KEYWORD easy,nonn AUTHOR Jason Earls, Jul 20 2004 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.

Last modified December 1 23:26 EST 2023. Contains 367503 sequences. (Running on oeis4.)