A353243 Indices of records of A070266.
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
Programs
-
Mathematica
m = Table[Max[ContinuedFraction[HarmonicNumber[n]]], {n, 1, 10^4}]; Map[FirstPosition[m, #][[1]] &, Union@FoldList[Max, m]]
-
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
Comments