cp's OEIS Frontend

A328945 Numbers m that are neither arithmetic (A003601) nor harmonic (A001599).

%I A328945 #18 Dec 14 2023 19:31:50
%S A328945 2,4,8,9,10,12,16,18,24,25,26,32,34,36,40,48,50,52,58,63,64,72,74,75,
%T A328945 76,80,81,82,84,88,90,98,100,104,106,108,112,117,120,121,122,124,128,
%U A328945 130,136,144,146,148,152,156,160,162,170,171,172,175,176,178,180
%N A328945 Numbers m that are neither arithmetic (A003601) nor harmonic (A001599).
%C A328945 Numbers m such that neither the arithmetic mean of the divisors of m nor the harmonic mean of the divisors of m is an integer.
%C A328945 Numbers m such that neither A(m) = A000203(m)/A000005(m) nor H(m) = m * A000005(m)/A000203(m) is an integer.
%C A328945 Corresponding values of A(m): 3/2, 7/3, 15/4, 13/3, 9/2, 14/3, 31/5, 13/2, 15/2, 31/3, 21/2, 21/2, 27/2, ...
%C A328945 Corresponding values of H(m): 4/3, 12/7, 32/15, 27/13, 20/9, 18/7, 80/31, 36/13, 16/5, 75/31, 52/21, 64/21, ...
%H A328945 Robert Israel, <a href="/A328945/b328945.txt">Table of n, a(n) for n = 1..10000</a>
%p A328945 filter:= proc(n) local D,d,t;
%p A328945   D:=numtheory:-divisors(n);
%p A328945   d:= nops(D);
%p A328945   convert(D,`+`) mod d <> 0 and not ((d/add(1/t,t=D))::integer)
%p A328945 end proc:
%p A328945 select(filter, [$1..200]); # _Robert Israel_, Dec 14 2023
%t A328945 Select[Range[180], !Divisible[DivisorSigma[1, #], DivisorSigma[0, #]] && !Divisible[# * DivisorSigma[0, #], DivisorSigma[1, #]] &] (* _Amiram Eldar_, Nov 01 2019 *)
%o A328945 (Magma) [m: m in [1..10^5] | not IsIntegral(m * NumberOfDivisors(m) / SumOfDivisors(m)) and not IsIntegral(SumOfDivisors(m) / NumberOfDivisors(m))]
%o A328945 (PARI) isok(m) = my(f = factor(m), prd = sigma(f)/numdiv(f)); (denominator(prd) != 1) && (denominator(m/prd) != 1); \\ _Michel Marcus_, Nov 05 2019
%Y A328945 Cf. A000005, A000203, A001599, A003601, A007340, A328944.
%K A328945 nonn
%O A328945 1,1
%A A328945 _Jaroslav Krizek_, Oct 31 2019