login
A353244
Record values of A070266.
1
1, 2, 5, 12, 13, 50, 61, 68, 198, 1090, 1812, 2362, 32334, 1517757, 4150055, 8850618, 10459960, 34693154, 80385407, 277184250, 316848842, 426258508, 2150821469, 2901738400, 49091889811, 331178405563
OFFSET
1,2
COMMENTS
The corresponding indices are in A353243.
FORMULA
a(n) = A070266(A353243(n)).
MATHEMATICA
Union @ FoldList[Max, Table[Max[ContinuedFraction[HarmonicNumber[n]]], {n, 1, 10^4}]]
PROG
(Python)
from itertools import count, islice
from fractions import Fraction
from sympy.ntheory.continued_fraction import continued_fraction
def A353244_gen(): # generator of terms
k, c = Fraction(), 0
for n in count(1):
k += Fraction(1, n)
if c < (m := max(continued_fraction(k))):
yield (c := m)
A353244_list = list(islice(A353244_gen(), 10)) # Chai Wah Wu, Apr 08 2022
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Apr 08 2022
STATUS
approved