A362286 Record values in A138705.
1, 2, 3, 6, 7, 10, 17, 23, 33, 39, 42, 46, 49, 54, 55, 57, 66, 73, 78, 83, 85, 95, 100, 105, 118, 133, 157, 162, 183, 201, 220, 224, 234, 242, 262, 272, 273, 287, 309, 314, 366
Offset: 1
Programs
-
Mathematica
seq[kmax_] := Module[{s = {}, mx = 0, m}, Do[m = Length[ContinuedFraction[ Abs[BernoulliB[2*k]]]]; If[m > mx, mx = m; AppendTo[s, m]], {k, 0, kmax}]; s]; seq[1000]
-
PARI
lista(kmax) = {my(mx = 0, m); for(k = 0, kmax, m = #contfrac(abs(bernfrac(2*k))); if(m > mx, mx = m; print1(m,", "))); }
Comments