|
%I #27 Apr 02 2020 07:50:13
%S 1,2,6,18,31,43,94,106,151,211,331,394,526,694,751,886,919,1114,1324,
%T 1726,1759,1831,2011,2311,2326,2671,3019,3691,3754,3931,4174,4951,
%U 4999,5119,6211,6406,7606,8254,8719,8779,9244,9619,9739,10399,10651,12919,13126
%N Indices in A146326 where records occur.
%H Robert G. Wilson v and Amiram Eldar, <a href="/A146345/b146345.txt">Table of n, a(n) for n = 1..416</a> (terms 1..218 from Robert G. Wilson v)
%p 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
%t 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 *)
%Y Cf. A000290, A078370, A146326-A146345, A146348-A146360, A146363.
%K nonn
%O 1,2
%A _Artur Jasinski_, Oct 30 2008
%E 19 replaced by 18, 394 inserted, 4 more terms added by _R. J. Mathar_, Sep 06 2009
%E More terms from _Robert G. Wilson v_, Apr 11 2017
|
