This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A053577 #18 Jan 29 2018 21:12:49
%S A053577 2,3,4,5,6,7,8,11,12,13,14,16,17,19,23,24,28,29,31,32,37,41,43,47,48,
%T A053577 53,56,59,61,62,64,67,71,73,79,83,89,96,97,101,103,107,109,112,113,
%U A053577 124,127,128,131,137,139,149,151,157,163,167,173,179,181,191,192,193,197
%N A053577 Cototient function n - phi(n) is a power of 2.
%C A053577 Consists of all primes, powers of 2, and powers of 2 times Mersenne primes (A000268). - _Robert Israel_, Jan 29 2018
%H A053577 Robert Israel, <a href="/A053577/b053577.txt">Table of n, a(n) for n = 1..10000</a>
%F A053577 A051953(n) = 2^w.
%e A053577 For p prime, cototient(p)=1; for x in the set {49152,57344,63488,65024,65528,65536}, x-phi(x) = 2^15 = 32768.
%p A053577 filter:= proc(n) local x; x:= n - numtheory:-phi(n); x = 2^padic:-ordp(x,2) end proc:
%p A053577 select(filter, [$1..300]); # _Robert Israel_, Jan 29 2018
%t A053577 Select[Range[200],IntegerQ[Log[2,#-EulerPhi[#]]]&] (* _Harvey P. Dale_, Dec 14 2011 *)
%o A053577 (PARI) select(x->hammingweight(x-eulerphi(x))==1,[1..200]) \\ _M. F. Hasler_, Nov 10 2016
%Y A053577 Cf. A000010, A000268, A051953, A053578.
%K A053577 nonn
%O A053577 1,1
%A A053577 _Labos Elemer_, Jan 18 2000
%E A053577 Edited by _M. F. Hasler_, Nov 10 2016