This site is supported by donations to The OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A237366 Primes p such that f(f(p)) is prime where f(x) = x^2+x+1. 0
 7, 19, 31, 67, 127, 181, 223, 241, 331, 367, 409, 463, 487, 673, 709, 751, 811, 823, 883, 997, 1117, 1231, 1321, 1489, 1549, 1861, 1933, 2083, 2179, 2287, 2473, 2551, 2707, 2803, 2851, 2857, 2917, 2971, 3067, 3361, 3499, 3559, 3691, 3847, 3931 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 LINKS EXAMPLE 31 is prime and (31^2+31+1)^2+(31^2+31+1)+1 = 987043 is prime. Thus, 31 is a member of this sequence. PROG (Python) import sympy from sympy import isprime {print(n) for n in range(10**4) if isprime(n) and isprime((n**2+n+1)**2+(n**2+n+1)+1)} (PARI) s=[]; forprime(p=2, 4000, if(isprime(p^4+2*p^3+4*p^2+3*p+3), s=concat(s, p))); s \\ Colin Barker, Feb 07 2014 CROSSREFS Cf. A237360, A002383. Sequence in context: A212492 A234310 A141338 * A216531 A216564 A145042 Adjacent sequences:  A237363 A237364 A237365 * A237367 A237368 A237369 KEYWORD nonn AUTHOR Derek Orr, Feb 06 2014 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified May 26 19:44 EDT 2019. Contains 323597 sequences. (Running on oeis4.)