OFFSET
1,2
COMMENTS
Here H(k) is the k-th harmonic number, Sum_{i=1..k} 1/i, and sigma(k) is the sum of the divisors of k. This function is derived from an inequality in Lagarias's paper. The condition f(k) > 0 for all k > 1 is equivalent to the Riemann hypothesis (RH). Every colossally abundant number (A004490) is here. Every superabundant number (A004394) greater than 665280 appears to be here also. This sequence is meaningless if the RH is false.
LINKS
T. D. Noe, Table of n, a(n) for n = 1..84
Jeffrey C. Lagarias, An elementary problem equivalent to the Riemann hypothesis, arXiv:math/0008177 [math.NT], 2000-2001; Amer. Math. Monthly 109 (#6, 2002), 534-543.
MATHEMATICA
(* This is just a naive recomputation of a dozen terms. *)
H = HarmonicNumber;
f[k_] := H[k] + Exp[H[k]] Log[H[k]] - DivisorSigma[1, k];
okQ[k_] := AllTrue[Range[k+1, 2k], f[#] > f[k]&];
Reap[For[k = 1, k < 10^4, k = If[k >= 60, k+60, k+1], If[okQ[k], Print[k]; Sow[k]]]][[2, 1]] (* Jean-François Alcover, Dec 11 2018 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
T. D. Noe, Apr 23 2010
STATUS
approved