OFFSET
1,2
COMMENTS
The corresponding record values are in A353244.
MATHEMATICA
m = Table[Max[ContinuedFraction[HarmonicNumber[n]]], {n, 1, 10^4}]; Map[FirstPosition[m, #][[1]] &, Union@FoldList[Max, m]]
PROG
(Python)
from itertools import count, islice
from fractions import Fraction
from sympy.ntheory.continued_fraction import continued_fraction
def A353243_gen(): # generator of terms
k, c = Fraction(), 0
for n in count(1):
k += Fraction(1, n)
if c < (m := max(continued_fraction(k))):
c = m
yield n
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Apr 08 2022
STATUS
approved