OFFSET
1,1
COMMENTS
This sequence does not contain only prime powers. Least term that has a prime factor which is not of the form m^2 + 1 is 35831808 = 2^14 * 3^7. The next one is 102503232 = 2^6 * 3^6 * 13^3. There are infinitely many such numbers.
EXAMPLE
125 = 5^3 is a term because phi(5^3) = 10^2 and cototient(5^3) = 5^2.
MATHEMATICA
Select[Range[10^6], Times @@ Boole@ Map[Or[# == 1, GCD @@ FactorInteger[#][[All, 2]] > 1] &, {#, EulerPhi@ #, # - EulerPhi@ #}] > 0 &] (* Michael De Vlieger, Jan 14 2017 *)
PROG
(PARI) is(n) = ispower(eulerphi(n)) && ispower(n-eulerphi(n)) && ispower(n);
CROSSREFS
KEYWORD
nonn
AUTHOR
Altug Alkan, Jan 13 2017
STATUS
approved