OFFSET
1,2
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..100
W. Sierpinski, Number Of Divisors And Their Sum, Monogr. Matemat. 42 (1964) chapter IV
EXAMPLE
For n=9, tau(phi(9)) = tau(6)=4 equals sigma(rad(9)) = sigma(3) = 4 which adds n=9 to the sequence.
MAPLE
with(numtheory):for n from 1 to 1500000 do : t1:= ifactors(n)[2] : t2 :=mul(t1[i][1], i=1..nops(t1)): if tau(phi(n)) = sigma(t2) then print (n): else fi: od :
MATHEMATICA
Select[Range[200000], DivisorSigma[0, EulerPhi[#]] == DivisorSigma[1, Times @@ FactorInteger[#][[All, 1]]] & ] (* Jean-François Alcover, Sep 12 2011 *)
PROG
(Magma) [1] cat [m:m in [2..200000]|#Divisors(EulerPhi(m)) eq &+Divisors(&*PrimeDivisors(m))]; // Marius A. Burtea, Jul 10 2019
(PARI) isok(n) = numdiv(eulerphi(n)) == sigma(factorback(factorint(n)[, 1])); \\ Michel Marcus, Jul 10 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Feb 23 2010
EXTENSIONS
Unspecific references removed by R. J. Mathar, Mar 26 2010
STATUS
approved