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-5 of 5 results.

A251782 Least even integer k such that numerator(B_k) == 0 (mod 37^n).

Original entry on oeis.org

32, 284, 37580, 1072544, 55777784, 325656968, 42764158652, 2444284077476, 46872402575720, 4093248733492712, 167845040875289732, 4841789050865438960, 235423026877046134208, 7818983737604766777920, 95503904455394036720840, 6908622244227620311285724, 114945213060615779807957456
Offset: 1

Views

Author

Keywords

Comments

37 is the first irregular prime. The corresponding entry for the second irregular prime 59 is A299466, and for the third irregular prime 67 is A299467.
The p-adic digits of the unique simple zero of the p-adic zeta function zeta_{(p,l)} with (p,l)=(37,32) were used to compute the sequence (see the Mathematica program below). This corresponds with Table A.2 in Kellner (2007). The sequence is increasing, but some consecutive entries are identical, e.g., entries 18 / 19 and 80 / 81. This is caused only by those p-adic digits that are zero.

Examples

			a(3) = 37580 because the numerator of B_37580 is divisible by 37^3 and there is no even integer less than 37580 for which this is the case.
		

Crossrefs

Programs

  • Mathematica
    p = 37; l = 32; LD = {7, 28, 21, 30, 4, 17, 26, 13, 32, 35, 27, 36, 32, 10, 21, 9, 11, 0, 1, 13, 6, 8, 10, 11, 10, 11, 32, 13, 30, 10, 6, 8, 2, 12, 1, 8, 2, 5, 3, 10, 19, 8, 4, 7, 19, 27, 33, 29, 29, 11, 2, 23, 8, 34, 5, 8, 35, 35, 13, 31, 29, 6, 7, 22, 13, 29, 7, 15, 22, 20, 19, 29, 2, 14, 2, 2, 31, 11, 4, 0, 27, 8, 10, 23, 17, 35, 15, 32, 22, 14, 7, 18, 8, 3, 27, 35, 33, 31, 6}; CalcIndex[L_, p_, l_, n_] := l + (p - 1) Sum[L[[i + 1]] p^i , {i, 0, n - 2}]; Table[CalcIndex[LD, p, l, n], {n, 1, Length[LD] + 1}] // TableForm

Formula

Numerator(B_{a(n)}) == 0 (mod 37^n).

Extensions

Edited for consistency with A299466 and A299467 by Bernd C. Kellner and Jonathan Sondow, Feb 20 2018

A299467 Least even integer k such that numerator(B_k) == 0 (mod 67^n).

Original entry on oeis.org

58, 3292, 153640, 12597148, 846312184, 52715297638, 320040068824, 370475739904372, 23170872799129498, 532379740455157312, 111861518490094080436, 1314934469494256636776, 291496130251698265225984, 7852328398132458266800348, 1925603427201316655808983674
Offset: 1

Views

Author

Keywords

Comments

67 is the third irregular prime. The corresponding entry for the first irregular prime 37 is A251782, and for the second irregular prime 59 is A299466.
The p-adic digits of the unique simple zero of the p-adic zeta function zeta_{(p,l)} with (p,l)=(67,58) were used to compute the sequence (see the Mathematica program below). This corresponds with Table A.2 in Kellner (2007). The sequence is increasing, but some consecutive entries are identical, e.g., entries 22 / 23 and 84 / 85. This is caused only by those p-adic digits that are zero.

Examples

			a(3) = 153640 because the numerator of B_153640 is divisible by 67^3 and there is no even integer less than 153640 for which this is the case.
		

Crossrefs

Programs

  • Mathematica
    p = 67; l = 58; LD = {49, 34, 42, 42, 39, 3, 62, 57, 19, 62, 10, 36, 14, 53, 57, 16, 60, 22, 41, 21, 25, 0, 56, 21, 24, 52, 33, 28, 51, 34, 60, 8, 47, 39, 42, 33, 14, 66, 50, 48, 45, 28, 61, 50, 27, 8, 30, 59, 32, 15, 3, 1, 54, 12, 30, 20, 14, 12, 10, 49, 33, 49, 54, 13, 26, 42, 8, 58, 12, 63, 19, 16, 48, 15, 2, 13, 1, 23, 2, 44, 64, 25, 40, 0, 16, 58, 44, 31, 62, 47, 61, 46, 9, 2, 50, 1, 62, 34, 31}; CalcIndex[L_, p_, l_, n_] := l + (p - 1) Sum[L[[i + 1]] p^i , {i, 0, n - 2}]; Table[CalcIndex[LD, p, l, n], {n, 1, Length[LD] + 1}] // TableForm

Formula

Numerator(B_{a(n)}) == 0 (mod 67^n).

A299468 p-adic digits of the unique simple zero of the p-adic zeta-function zeta_{(p,l)} with (p,l) = (37,32).

Original entry on oeis.org

7, 28, 21, 30, 4, 17, 26, 13, 32, 35, 27, 36, 32, 10, 21, 9, 11, 0, 1, 13, 6, 8, 10, 11, 10, 11, 32, 13, 30, 10, 6, 8, 2, 12, 1, 8, 2, 5, 3, 10, 19, 8, 4, 7, 19, 27, 33, 29, 29, 11, 2, 23, 8, 34, 5, 8, 35, 35, 13, 31, 29, 6, 7, 22, 13, 29, 7, 15, 22, 20, 19, 29, 2, 14, 2, 2, 31, 11, 4, 0, 27, 8, 10, 23, 17, 35, 15, 32, 22, 14, 7, 18, 8, 3, 27, 35, 33, 31, 6
Offset: 0

Views

Author

Keywords

Comments

The p-adic digits are used to compute A251782(n) = least even integer k such that numerator(B_k) == 0 (mod 37^n) (see 2nd formula below and the program in A251782).
The algorithm used in the Mathematica program below is from Kellner 2007, Prop. 5.3, p. 428.
The corresponding sequences for (p,l) = (59,44) and (p,l) = (67,58) are A299469 and A299470, respectively.

Examples

			The zero is given by a(0) + a(1)*p + a(2)*p^2 + ... with p = 37.
		

Crossrefs

Programs

  • Mathematica
    n = 99; p = 37; l = 32;
    ModR[x_, m_] := Mod[Mod[Numerator[x], m] PowerMod[Denominator[x], -1, m], m];
    B[n_] := -(1 - p^(n - 1)) BernoulliB[n]/n;
    T[r_, k_, x_] := Sum[(-1)^(j + k) Binomial[j, k] Binomial[x, j], {j, k, r}];
    zt = Table[ModR[B[l + (p - 1) k]/p, p^n], {k, 0, n}];
    Z[n_] := zt[[n + 1]]; d = Mod[Z[0] - Z[1], p]; t = 0; L = {};
    For[r = 1, r <= n, r++, x = Mod[Sum[Z[k] T[r, k, t], {k, 0, r}], p^r];
      s = ModR[x/(d*p^(r - 1)), p]; AppendTo[L, s]; t += s*p^(r - 1)];
    Print[L]

Formula

0 <= a(n) <= 36.
l + (p - 1)*Sum_{i=0..n-2} a(i)*p^i = A251782(n) with (p,l) = (37,32).

A299469 p-adic digits of the unique simple zero of the p-adic zeta-function zeta_{(p,l)} with (p,l) = (59,44).

Original entry on oeis.org

15, 25, 40, 36, 18, 11, 17, 28, 58, 9, 51, 13, 25, 41, 44, 17, 43, 35, 21, 10, 21, 38, 9, 12, 40, 43, 45, 30, 41, 0, 3, 25, 34, 49, 45, 9, 19, 48, 57, 11, 13, 29, 28, 44, 41, 37, 33, 29, 43, 8, 57, 12, 48, 15, 15, 53, 57, 16, 51, 16, 54, 30, 9, 26, 8, 49, 22, 58, 11, 42, 28, 36, 33, 45, 24, 32, 18, 12, 29, 45, 40, 27, 19, 40, 41, 11, 42, 49, 35, 41, 57, 54, 33, 0, 34, 34, 49, 6, 31
Offset: 0

Views

Author

Keywords

Comments

The p-adic digits are used to compute A299466(n) = least even integer k such that numerator(B_k) == 0 (mod 59^n) (see 2nd formula below and the program in A299466).
The algorithm used in the Mathematica program below is from Kellner 2007, Prop. 5.3, p. 428.
The corresponding sequences for (p,l) = (37,32) and (p,l) = (67,58) are A299468 and A299470, respectively.

Examples

			The zero is given by a(0) + a(1)*p + a(2)*p^2 + ... with p = 59.
		

Crossrefs

Programs

  • Mathematica
    n = 99; p = 59; l = 44;
    ModR[x_, m_] := Mod[Mod[Numerator[x], m] PowerMod[Denominator[x], -1, m], m];
    B[n_] := -(1 - p^(n - 1)) BernoulliB[n]/n;
    T[r_, k_, x_] := Sum[(-1)^(j + k) Binomial[j, k] Binomial[x, j], {j, k, r}];
    zt = Table[ModR[B[l + (p - 1) k]/p, p^n], {k, 0, n}];
    Z[n_] := zt[[n + 1]]; d = Mod[Z[0] - Z[1], p]; t = 0; L = {};
    For[r = 1, r <= n, r++, x = Mod[Sum[Z[k] T[r, k, t], {k, 0, r}], p^r];
      s = ModR[x/(d*p^(r - 1)), p]; AppendTo[L, s]; t += s*p^(r - 1)];
    Print[L]

Formula

0 <= a(n) <= 58.
l + (p - 1)*Sum_{i=0..n-2} a(i)*p^i = A299466(n) with (p,l) = (59,44).

A299470 p-adic digits of the unique simple zero of the p-adic zeta-function zeta_{(p,l)} with (p,l) = (67,58).

Original entry on oeis.org

49, 34, 42, 42, 39, 3, 62, 57, 19, 62, 10, 36, 14, 53, 57, 16, 60, 22, 41, 21, 25, 0, 56, 21, 24, 52, 33, 28, 51, 34, 60, 8, 47, 39, 42, 33, 14, 66, 50, 48, 45, 28, 61, 50, 27, 8, 30, 59, 32, 15, 3, 1, 54, 12, 30, 20, 14, 12, 10, 49, 33, 49, 54, 13, 26, 42, 8, 58, 12, 63, 19, 16, 48, 15, 2, 13, 1, 23, 2, 44, 64, 25, 40, 0, 16, 58, 44, 31, 62, 47, 61, 46, 9, 2, 50, 1, 62, 34, 31
Offset: 0

Views

Author

Keywords

Comments

The p-adic digits are used to compute A299467(n) = least even integer k such that numerator(B_k) == 0 (mod 67^n) (see 2nd formula below and the program in A299467).
The algorithm used in the Mathematica program below is from Kellner 2007, Prop. 5.3, p. 428.
The corresponding sequences for (p,l) = (37,32) and (p,l) = (59,44) are A299468 and A299469, respectively.

Examples

			The zero is given by a(0) + a(1)*p + a(2)*p^2 + ... with p = 67.
		

Crossrefs

Programs

  • Mathematica
    n = 99; p = 67; l = 58;
    ModR[x_, m_] := Mod[Mod[Numerator[x], m] PowerMod[Denominator[x], -1, m], m];
    B[n_] := -(1 - p^(n - 1)) BernoulliB[n]/n;
    T[r_, k_, x_] := Sum[(-1)^(j + k) Binomial[j, k] Binomial[x, j], {j, k, r}];
    zt = Table[ModR[B[l + (p - 1) k]/p, p^n], {k, 0, n}];
    Z[n_] := zt[[n + 1]]; d = Mod[Z[0] - Z[1], p]; t = 0; L = {};
    For[r = 1, r <= n, r++, x = Mod[Sum[Z[k] T[r, k, t], {k, 0, r}], p^r];
      s = ModR[x/(d*p^(r - 1)), p]; AppendTo[L, s]; t += s*p^(r - 1)];
    Print[L]

Formula

0 <= a(n) <= 66.
l + (p - 1)*Sum_{i=0..n-2} a(i)*p^i = A299467(n) with (p,l) = (67,58).
Showing 1-5 of 5 results.