login
A353243
Indices of records of A070266.
1
1, 2, 3, 4, 10, 11, 18, 27, 30, 32, 43, 69, 70, 264, 1409, 3027, 7471, 8946, 10576, 12595, 14034, 22849, 37124, 70083, 107868, 469850
OFFSET
1,2
COMMENTS
The corresponding record values are in A353244.
FORMULA
A070266(a(n)) = A353244(n).
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
A353243_list = list(islice(A353243_gen(), 10)) # Chai Wah Wu, Apr 08 2022
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Apr 08 2022
STATUS
approved