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.

A245089 The unique integer r with |r| < prime(n)/2 such that B_{prime(n)-2}(1/3) == r (mod prime(n)), where B_m(x) denotes the Bernoulli polynomial of degree m.

Original entry on oeis.org

-2, -1, 4, -6, 8, -6, -10, -5, 3, -16, 4, 6, 3, 6, -11, -29, 2, 7, 21, 4, -16, -23, -5, 43, 14, 3, -32, 26, 13, -23, 64, 52, -30, -74, -17, -33, 37, -82, -68, 55, -78, 96, 79, 22, -81, -26, -7, 70, -38, 9, 3, -118, 128, -123, -67, -69, -78, -138, 30, -60, -19, 88, -26, 110, 27, 63, -82, 138
Offset: 3

Views

Author

Zhi-Wei Sun, Jul 11 2014

Keywords

Comments

Conjecture: a(n) = 0 infinitely often. In other words, there are infinitely many primes p > 3 such that B_{p-2}(1/3) == 0 (mod p).
This seems reasonable in view of the standard heuristic arguments. Our computation shows that if a(n) = 0 then n > 2600 and hence prime(n) > 23000.
Zhi-Wei Sun made many conjectures on congruences involving B_{p-2}(1/3), see the Sci. China Math. paper and arXiv:1407.0967.
The first value of n with a(n) = 0 is 18392. For the prime p = prime(18392) = 205129, we have B_{p-2}(1/3) == 20060*p (mod p^2). - Zhi-Wei Sun, Dec 13 2014

Examples

			a(3) = -2 since B_{prime(3)-2}(1/3) = B_3(1/3) = 1/27 == -2 (mod prime(3)=5).
		

Crossrefs

Programs

  • Mathematica
    rMod[m_,n_]:=Mod[Numerator[m]*PowerMod[Denominator[m],-1,n],n,-n/2]
    a[n_]:=rMod[BernoulliB[Prime[n]-2,1/3],Prime[n]]
    Table[a[n],{n,3,70}]