OFFSET
1,2
EXAMPLE
f(sqrt(2)) = f([1,2,2,...]) = [2,3,3,...] = (1 + sqrt(13)/2, so a(2) = 13.
MATHEMATICA
$MaxExtraPrecision = Infinity;
c[x_] := c[x] = FromContinuedFraction[ContinuedFraction[x] + 1]
Table[c[Sqrt[n]], {n, 1, 30}]
f[y_] := Cases[y, x_^(1/2 | -1/2) :> x, Infinity];
t = Table[f[c[Sqrt[n]]], {n, 1, 80}]; Flatten[t /. {} -> 1] (*A229957*)
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Oct 04 2013
STATUS
approved