OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..416
FORMULA
MAPLE
146326 := proc(n) if not issqr(n) then numtheory[cfrac]( (1+sqrt(n))/2, 'periodic', 'quotients') ; nops(%[2]) ; else 0 ; fi; end: read("transforms") ; a26 := [seq(A146326(n), n=1..1000)] ; RECORDS(a26)[1] ; # R. J. Mathar, Sep 06 2009
MATHEMATICA
f[n_] := Length @ ContinuedFraction[(1 + Sqrt[n])/2][[-1]]; fmax = -1; seq = {}; Do[f1 = f[n]; If[f1 > fmax, fmax = f1; AppendTo[seq, f1]], {n, 1, 10^4}]; seq (* Amiram Eldar, Apr 02 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 30 2008
EXTENSIONS
42 inserted by R. J. Mathar, Sep 06 2009
a(1) inserted and more terms added by Amiram Eldar, Apr 02 2020
STATUS
approved