OFFSET
1,1
COMMENTS
Consists of all primes, powers of 2, and powers of 2 times Mersenne primes (A000268). - Robert Israel, Jan 29 2018
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
A051953(n) = 2^w.
EXAMPLE
For p prime, cototient(p)=1; for x in the set {49152,57344,63488,65024,65528,65536}, x-phi(x) = 2^15 = 32768.
MAPLE
filter:= proc(n) local x; x:= n - numtheory:-phi(n); x = 2^padic:-ordp(x, 2) end proc:
select(filter, [$1..300]); # Robert Israel, Jan 29 2018
MATHEMATICA
Select[Range[200], IntegerQ[Log[2, #-EulerPhi[#]]]&] (* Harvey P. Dale, Dec 14 2011 *)
PROG
(PARI) select(x->hammingweight(x-eulerphi(x))==1, [1..200]) \\ M. F. Hasler, Nov 10 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Jan 18 2000
EXTENSIONS
Edited by M. F. Hasler, Nov 10 2016
STATUS
approved
