A182226 Abundant numbers m (in A005101) whose sum of divisors less than their abundance A033880(m) does not exceed A033880(m).
18, 20, 70, 78, 88, 102, 104, 114, 138, 174, 186, 196, 222, 246, 258, 282, 318, 354, 366, 368, 402, 426, 438, 464, 474, 498, 534, 582, 606, 618, 642, 650, 654, 678, 762, 786, 822, 834, 836, 894, 906, 942, 978, 1002, 1014, 1038, 1074, 1086, 1146, 1158, 1182, 1194, 1266, 1338, 1362, 1374, 1398
Offset: 1
Keywords
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
Programs
-
Mathematica
aQ[n_] := (ab = DivisorSigma[1, n] - 2 n) > 0 && DivisorSum[n, # &, # < ab &] <= ab; Select[Range[1400], aQ] (* Amiram Eldar, Sep 08 2019 *)
-
PARI
for(n=1,1999,sigma(n)>2*n||next;is_A182225(n)||print1(n","))
Comments