login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

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