OFFSET
1,2
PROG
(PARI) for(s=1, 100000, if(sum(i=1, length(binary(s)), if(component(binary(s), i)-1, 0, 1))^sum(i=1, length(binary(s)), if(component(binary(s), i), 0, 1))==eulerphi(s), print1(s, ", ")))
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Jun 22 2002
STATUS
approved