OFFSET
1,1
COMMENTS
The first 4 known Fermat primes > 3 from A019434 are in the sequence.
LINKS
Michel Marcus, Table of n, a(n) for n = 1..70
EXAMPLE
257 is a term because phi(257+1) = phi(258) = 84 = phi(phi(257-2)+1) = phi(phi(255)+1) = phi(128+1) = phi(129).
MATHEMATICA
Select[Prime@ Range[10^6], EulerPhi[# + 1] == EulerPhi[EulerPhi[# - 2] + 1] &] (* Michael De Vlieger, Apr 11 2016 *)
PROG
(Magma) [n: n in [2..10^8] | IsPrime(n) and EulerPhi(n+1) eq EulerPhi(EulerPhi(n-2)+1)]
(PARI) is(n) = eulerphi(n+1)==eulerphi(eulerphi(n-2)+1)
forprime(p=3, , if(is(p), print1(p, ", "))) \\ Felix Fröhlich, Apr 23 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Apr 11 2016
EXTENSIONS
Definition adjusted and more terms added by Felix Fröhlich, Apr 23 2016
STATUS
approved