login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A238445 Primes p such that f(f(p)) is prime, where f(x) = x^5-x^4-x^3-x^2-x-1. 0
3, 13, 61, 103, 193, 199, 307, 431, 569, 977, 1201, 1451, 1481, 1609, 1669, 1889, 2371, 2381, 2711, 2819, 3083, 3469, 4289, 4337, 4567, 5231, 5501, 6733, 7043, 7253, 7351, 7549, 8707, 9257, 9497, 10039, 10687, 11491, 12227, 12517, 12941, 13397 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

LINKS

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

EXAMPLE

3 is prime. 3^5-3^4-3^3-3^2-3-1 = 122 and 122^5-122^4-122^3-122^2-122-1 = 26803717321 is a prime number. Thus, 3 is a member of this sequence.

PROG

(Python)

import sympy

from sympy import isprime

def f(x):

..return x**5-x**4-x**3-x**2-x-1

{print(p) for p in range(10**5) if isprime(p) and isprime(f(f(p)))}

CROSSREFS

Cf. A125083, A237640.

Sequence in context: A112731 A106884 A232611 * A328704 A112568 A104089

Adjacent sequences:  A238442 A238443 A238444 * A238446 A238447 A238448

KEYWORD

nonn

AUTHOR

Derek Orr, Feb 26 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 January 21 05:23 EST 2020. Contains 331104 sequences. (Running on oeis4.)