|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
The corresponding record values are in A353244.
|
|
LINKS
|
|
|
FORMULA
|
|
|
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
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|