A328945 Numbers m that are neither arithmetic (A003601) nor harmonic (A001599).
2, 4, 8, 9, 10, 12, 16, 18, 24, 25, 26, 32, 34, 36, 40, 48, 50, 52, 58, 63, 64, 72, 74, 75, 76, 80, 81, 82, 84, 88, 90, 98, 100, 104, 106, 108, 112, 117, 120, 121, 122, 124, 128, 130, 136, 144, 146, 148, 152, 156, 160, 162, 170, 171, 172, 175, 176, 178, 180
Offset: 1
Keywords
Links
- Robert Israel, Table of n, a(n) for n = 1..10000
Programs
-
Magma
[m: m in [1..10^5] | not IsIntegral(m * NumberOfDivisors(m) / SumOfDivisors(m)) and not IsIntegral(SumOfDivisors(m) / NumberOfDivisors(m))]
-
Maple
filter:= proc(n) local D,d,t; D:=numtheory:-divisors(n); d:= nops(D); convert(D,`+`) mod d <> 0 and not ((d/add(1/t,t=D))::integer) end proc: select(filter, [$1..200]); # Robert Israel, Dec 14 2023
-
Mathematica
Select[Range[180], !Divisible[DivisorSigma[1, #], DivisorSigma[0, #]] && !Divisible[# * DivisorSigma[0, #], DivisorSigma[1, #]] &] (* Amiram Eldar, Nov 01 2019 *)
-
PARI
isok(m) = my(f = factor(m), prd = sigma(f)/numdiv(f)); (denominator(prd) != 1) && (denominator(m/prd) != 1); \\ Michel Marcus, Nov 05 2019
Comments