OFFSET
1,1
FORMULA
MATHEMATICA
Select[Range[300000], # != 2^IntegerExponent[#, 2] && DivisorSigma[1, #] + 2*EulerPhi[#] == DivisorSigma[1, EulerPhi[#]] + 2*# &] (* Amiram Eldar, Dec 08 2023 *)
PROG
(PARI) test(n)=my(p=eulerphi(n)); sigma(n)+2*p==sigma(p)+2*n
for(k=1, 20, for(n=1<<k+1, 2<<k-1, if(test(n), print1(n", ")))) \\ Charles R Greathouse IV, Oct 12 2011
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Naohiro Nomoto, Oct 12 2011
EXTENSIONS
a(6)-a(9) from Charles R Greathouse IV, Oct 12 2011
a(10)-a(12) from Amiram Eldar, Dec 08 2023
STATUS
approved