A353243 - OEIS

%I #10 Apr 09 2022 02:24:10

%S 1,2,3,4,10,11,18,27,30,32,43,69,70,264,1409,3027,7471,8946,10576,

%T 12595,14034,22849,37124,70083,107868,469850

%N Indices of records of A070266.

%C The corresponding record values are in A353244.

%F A070266(a(n)) = A353244(n).

%t m = Table[Max[ContinuedFraction[HarmonicNumber[n]]], {n, 1, 10^4}]; Map[FirstPosition[m, #][[1]] &, Union@FoldList[Max, m]]

%o (Python)

%o from itertools import count, islice

%o from fractions import Fraction

%o from sympy.ntheory.continued_fraction import continued_fraction

%o def A353243_gen(): # generator of terms

%o     k, c = Fraction(), 0

%o     for n in count(1):

%o         k += Fraction(1,n)

%o         if c < (m := max(continued_fraction(k))):

%o             c = m

%o             yield n

%o A353243_list = list(islice(A353243_gen(),10)) # _Chai Wah Wu_, Apr 08 2022

%Y Cf. A001008, A002805, A070266, A091532, A353244.

%K nonn,more

%O 1,2

%A _Amiram Eldar_, Apr 08 2022