OFFSET
1,2
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
MATHEMATICA
q[k_] := DivisorSigma[1, k] - EulerPhi[k] <= DivisorSigma[0, k]^PrimeNu[k]; Select[Range[160], q] (* Amiram Eldar, Sep 15 2024 *)
PROG
(PARI) is(n)=sigma(n)-eulerphi(n)<=numdiv(n)^omega(n) \\ Charles R Greathouse IV, Feb 14 2013
(Haskell)
a110086 n = a110086_list !! (n-1)
a110086_list = filter (\x -> a051612 x <= a110088 x) [1..]
-- Reinhard Zumkeller, Aug 05 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jul 11 2005
STATUS
approved