cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A362284 a(n) is the least number k such that A138705(k) = n, or -1 if no such k exists.

Original entry on oeis.org

0, 1, 5, 10, 23, 6, 8, 20, 50, 18, 69, 66, 100, 45, 56, 42, 30, 40, 96, 99, 72, 234, 36, 348, 156, 200, 168, 120, 405, 390, 216, 315, 90, 198, 280, 270, 792, 210, 180, 624, 1120, 360, 576, 1188, 1134, 420, 750, 1140, 504, 1728, 660, 600, 690, 540, 630, 1380, 810
Offset: 1

Views

Author

Amiram Eldar, Apr 14 2023

Keywords

Crossrefs

Programs

  • Mathematica
    seq[len_, nmax_] := Module[{s = Table[-1, {len}], c = 0, n = 0, i}, While[c < len && n <= nmax, i = Length[ContinuedFraction[Abs[BernoulliB[2*n]]]]; If[i <= len && s[[i]] < 0, c++; s[[i]] = n]; n++]; s]; seq[60, 10^4]
  • PARI
    lista(len, nmax) = {my(s = vector(len,i,-1), c = 0, n = 0, i); while(c < len && n <= nmax, i = #contfrac(abs(bernfrac(2*n))); if(i <= len && s[i] < 0, c++; s[i] = n); n++); s;}

A362286 Record values in A138705.

Original entry on oeis.org

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

Views

Author

Amiram Eldar, Apr 14 2023

Keywords

Comments

The corresponding indices of records are in A362285.

Crossrefs

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,", "))); }

Formula

a(n) = A138705(A362285(n)).
Showing 1-2 of 2 results.