|
%I #13 Aug 04 2019 05:49:58
%S 1,2,3,4,6,8,10,12,16,18,20,24,30,36,48,60,72,84,90,108,120,144,168,
%T 180,240,336,360,420,480,504,630,720,840,1008,1080,1260,1440,1680,
%U 2160,2520,3360,3780,3960,4200,4320,4620,5040,7560,9240,10080,12600,13860,15120
%N Positions of records in the sequence of harmonic means, i.e., in the sequence of rationals A099377(.)/A099378(.).
%H Amiram Eldar, <a href="/A179971/b179971.txt">Table of n, a(n) for n = 1..200</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/OresConjecture.html">Ore's Conjecture</a>
%e The sequence of harmonic means starts 1 < 4/3 < 3/2 < 12/7, increasing from the first to the fourth, which adds 1 to 4 to the sequence.
%e The fifth harmonic mean is 5/3, smaller than 12/7 and not a record, so 5 is not in the sequence.
%p hm := proc(n) option remember; n* numtheory[tau](n)/numtheory[sigma](n) ; end proc:
%p A179971 := proc(n) option remember; if n = 1 then 1; else for k from procname(n-1)+1 do if hm(k) > hm(procname(n-1)) then return k; end if; end do; end if; end proc:
%p seq(A179971(n),n=1..40) ; # _R. J. Mathar_, Aug 06 2010
%t f[n_] := f[n] = DivisorSigma[0, n]/Plus @@ (1/Divisors@n); k = 1; mx = 0; lst = {}; While[k < 18480, a = f@k; If[a > mx, mx = a; AppendTo[lst, k]]; k++ ]; lst
%Y Cf. A099377, A099378.
%K nonn
%O 1,2
%A _Robert G. Wilson v_, Aug 04 2010
%E Definition rephrased by _R. J. Mathar_, Aug 06 2010
|
