%I #11 Feb 23 2014 09:31:49
%S 122,340352,830519696,11479086422,266390469692,310503441398,
%T 2718130415306,14837993872846,59538248604388,889257663626476,
%U 2496623039993996,6427431330617746,7120028814392596,10777302002014868,12942591289426088,24039736320940828
%N Numbers n of the form p^5 - Phi_5(p) (for prime p) such that n^5 - Phi_5(n) is also prime.
%C All numbers are congruent to 2 mod 10, 6 mod 10, or 8 mod 10.
%C x^5 - Phi_5(x) = x^5-x^4-x^3-x^2-x-1.
%e 122 = 3^5-3^4-3^3-3^2-3^1-1 (3 is prime) and 122^5-122^4-122^3-122^2-122^1-1 = 26803717321 is prime. Thus, 122 is a member of this sequence.
%o (Python)
%o import sympy
%o from sympy import isprime
%o def poly5(x):
%o ..if isprime(x):
%o ....f = x**5-x**4-x**3-x**2-x-1
%o ....if isprime(f**5-f**4-f**3-f**2-f-1):
%o ......return True
%o ..return False
%o x = 1
%o while x < 10**5:
%o ..if poly5(x):
%o ....print(x**5-x**4-x**3-x**2-x-1)
%o ..x += 1
%Y Cf. A125083, A237527, A237528, A237639.
%K nonn
%O 1,1
%A _Derek Orr_, Feb 10 2014
|