OFFSET
1,2
MATHEMATICA
c = 0; a = h = 1; t = {}; Do[a = LCM[a, n]; h = h + 1/n; b = a/Denominator[h]; If[b > c, c = b; AppendTo[t, n]], {n, 10^6}]; t
PROG
(PARI) lista(nn) = {rec = 0; for (n=1, nn, new = lcm(vector(n, k, k))/denominator(sum(k=1, n, 1/k)); if (new > rec, print1(n, ", "); rec = new); ); } \\ Michel Marcus, Mar 07 2018
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Robert G. Wilson v, Sep 19 2005
EXTENSIONS
a(27)-a(28) from Amiram Eldar, Dec 18 2018
a(29)-a(31) from Chai Wah Wu, Mar 08 2021
a(32) from Chai Wah Wu, Mar 14 2021
STATUS
approved