OFFSET
1,2
FORMULA
{ k : d(k) + phi(k) - 1 > k/2 }.
EXAMPLE
30 is not in the sequence because d(30) + phi(30) - 1 = 8 + 8 - 1 = 15. There are as many divisors and coprimes as there are numbers j <= 30 that neither divide nor are coprime to 30.
50 is not here because d(50) + phi(50) - 1 = 6 + 20 - 1 = 25. There are as many divisors and coprimes as there are numbers j < 50 that neither divide nor are coprime to 50.
146 is here because d(146) + phi(146) - 1 = 4 + 72 - 1 = 75; 146/2 = 73, and 75 > 73.
61455 is here because d(61455) + phi(61455) - 1 = 16 + 30720 - 1 = 30735; 61455/2 = 30727 + 1/2, and 30735 > 61455/2.
MATHEMATICA
Do[r=EulerPhi[n]; d=DivisorSigma[0, n]; u=n-r-d+1; If[Greater[n-u, n/2], Print[n, {d, r, u}]], {n, 1, 100}]
(* Second program: *)
Select[Range[120], DivisorSigma[0, #] + EulerPhi[#] - 1 > #/2 &] (* Michael De Vlieger, Aug 22 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, May 07 2003
EXTENSIONS
Data corrected and entry edited by Michael De Vlieger, Aug 22 2023
STATUS
approved