login
This site is supported by donations to The OEIS Foundation.

 

Logo


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

Table of n, a(n) for n=1..45.

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.

License Agreements, Terms of Use, Privacy Policy. .

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