login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A266266 Primes p such that p-1 = phi(p) = k*phi(p-k) for some number 1 <= k < p-1. 1
5, 7, 11, 13, 17, 73, 257, 2593, 65537 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Corresponding values of numbers k: 2, 3, 5, 3, 2, 3, 2, 3, 2, ...
83623937 is also a term of this sequence (with k = 2).
For all primes p we have: phi(p) = k*phi(p-k) if k = p - 1.
Primes from A266267.
The first 4 known Fermat primes > 3 from A019434 are in sequence.
LINKS
EXAMPLE
17 is in the sequence because phi(17) = 16 = 2*phi(15) = 2*8.
MATHEMATICA
Select[Prime@ Range@ 500, Function[p, AnyTrue[Range[p - 2], p - 1 == # EulerPhi[p - #] &]]] (* Michael De Vlieger, Jan 09 2016, Version 10 *)
PROG
(Magma) Set(Sort([[n: k in [1..n-2] | IsPrime(n) and EulerPhi(n) eq k*EulerPhi(n-k)]: n in [1..10000]]))
(Magma) Set(Sort([5] cat [n: n in [6..100000], k in [1..5] | IsPrime(n) and EulerPhi(n) eq k*EulerPhi(n-k)]))
(PARI) listp(nn) = {forprime(p=2, nn, for (k=1, p-2, if (eulerphi(p) == k*eulerphi(p-k), print1(p, ", "); break)); ); } \\ Michel Marcus, Dec 27 2015
CROSSREFS
Sequence in context: A317250 A007529 A287956 * A246463 A314297 A108409
KEYWORD
nonn,more
AUTHOR
Jaroslav Krizek, Dec 26 2015
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 16 18:22 EDT 2024. Contains 371750 sequences. (Running on oeis4.)