login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A299468 p-adic digits of the unique simple zero of the p-adic zeta-function zeta_{(p,l)} with (p,l) = (37,32). 3
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 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,1

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.

LINKS

Bernd C. Kellner and Jonathan Sondow, Table of n, a(n) for n = 0..98

Bernd C. Kellner, The Bernoulli Number Page

Bernd C. Kellner, On irregular prime power divisors of the Bernoulli numbers, Math. Comp. 76 (2007) 405-441.

FORMULA

0 <= a(n) <= 36.

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

EXAMPLE

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

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]

CROSSREFS

Cf. A251782, A299466, A299467, A299469, A299470.

Sequence in context: A274579 A152578 A300529 * A200974 A184331 A267432

Adjacent sequences:  A299465 A299466 A299467 * A299469 A299470 A299471

KEYWORD

nonn

AUTHOR

Bernd C. Kellner and Jonathan Sondow, Mar 28 2018

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 19 12:36 EST 2020. Contains 331049 sequences. (Running on oeis4.)