OFFSET
1,2
COMMENTS
LINKS
Clark Kimberling, Table of n, a(n) for n = 1..12
EXAMPLE
Let r = sqrt(1/2). Then
1 + 1/2 < r < 1 + 1/2 + 1/3, so a(1) = 3.
1 + 1/2 + 1/3 -1/24 < r, so a(2) = 24.
1 + 1/2 + 1/3 -1/24 + 1/815 > r, so a(3) = 815.
MATHEMATICA
$MaxExtraPrecision = Infinity;
nn = 12; f[n_] := 1/n; r = Sqrt[1/2]; s = 0; b[1] = NestWhile[# + 1 &, 1, ! (s += f[#]) > r &]; u[1] = Sum[f[n], {n, 1, b[1]}]; c[1] = Floor[1/(u[1] - r)]; v[1] = u[1] - 1/c[1]; n = 1; While[n < nn/2, n++; b[n] = Floor[1/(r - v[n - 1])]; u[n] = v[n - 1] + 1/b[n]; c[n] = Floor[1/(u[n] - r)]; v[n] = u[n] - 1/c[n]]; a = Riffle[Table[b[i], {i, 1, nn/2}], Table[c[i], {i, 1, nn/2}]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, May 24 2013
STATUS
approved