A146345 Indices in A146326 where records occur.
1, 2, 6, 18, 31, 43, 94, 106, 151, 211, 331, 394, 526, 694, 751, 886, 919, 1114, 1324, 1726, 1759, 1831, 2011, 2311, 2326, 2671, 3019, 3691, 3754, 3931, 4174, 4951, 4999, 5119, 6211, 6406, 7606, 8254, 8719, 8779, 9244, 9619, 9739, 10399, 10651, 12919, 13126
Offset: 1
Keywords
Links
- Robert G. Wilson v and Amiram Eldar, Table of n, a(n) for n = 1..416 (terms 1..218 from Robert G. Wilson v)
Programs
-
Maple
A146326 := 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..1400)] ; RECORDS(a26)[2] ; # R. J. Mathar, Sep 06 2009
-
Mathematica
f[n_] := Length@ContinuedFraction[(1 + Sqrt[n])/2][[-1]]; mx = -1; k = 1; lst = {}; While[k < 14000, a = f@k; If[a > mx, mx = a; AppendTo[lst, k]]; k++]; lst (* Robert G. Wilson v, Apr 11 2017 *)
Extensions
19 replaced by 18, 394 inserted, 4 more terms added by R. J. Mathar, Sep 06 2009
More terms from Robert G. Wilson v, Apr 11 2017