OFFSET
1,2
COMMENTS
A000043(k) - 1 is a term for all k >= 1. - Amiram Eldar, Aug 22 2019
No more terms below 1206. - Amiram Eldar, Aug 23 2019
EXAMPLE
MATHEMATICA
f[n_] := DivisorSigma[1, EulerPhi[n]] - EulerPhi[DivisorSigma[1, n]]; pow2Q[n_] := n == 2^IntegerExponent[n, 2]; aQ[n_] := pow2Q[1 - f[2^n]]; Select[Range[130], aQ] (* Amiram Eldar, Aug 22 2019 *)
PROG
(PARI) f(n) = sigma(eulerphi(n)) - eulerphi(sigma(n)); \\ A065395
ispp2(k) = isprimepower(k, &p) && (p==2);
isok(n) = ispp2(1-f(2^n)); \\ Michel Marcus, Aug 22 2019
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Labos Elemer, Mar 03 2004
EXTENSIONS
Name and example edited by Michel Marcus, Aug 22 2019
a(17)-a(18) from Amiram Eldar, Aug 23 2019
STATUS
approved