OFFSET
1,2
COMMENTS
phi(k) = sigma(core(k)) | sigma(k) so this is a subsequence of A020492. - David A. Corneth, Sep 08 2020
LINKS
David A. Corneth, Table of n, a(n) for n = 1..16559 (first 183 terms from Amiram Eldar, terms calculated from terms provided by Jud McCranie in A020492)
EXAMPLE
140 is in the sequence as phi(140) = 48 = sigma(35) = sigma(core(140)). - David A. Corneth, Sep 08 2020
MATHEMATICA
core[n_] := Times @@ Power @@@ ({#[[1]], Mod[ #[[2]], 2]} & /@ FactorInteger[n]); Select[Range[10^5], EulerPhi[#] == DivisorSigma[1, core[#]] &] (* Amiram Eldar, Jul 11 2019 after Zak Seidov at A007913 *)
PROG
(PARI) for(n=1, 10^6, if(eulerphi(n)==sigma(core(n)), print1(n, ", ")))
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Benoit Cloitre, Apr 17 2002
EXTENSIONS
More terms from Amiram Eldar, Jul 11 2019
STATUS
approved