A119766 Numbers n such that numerator of Bernoulli(n)/n is (apart from sign) 1 or a prime.

1, 2, 4, 6, 8, 10, 12, 14, 16, 18, 26, 34, 36, 38, 42, 74, 114, 118, 396, 674, 1870, 4306, 22808

Offset

1,2

Comments

In 1911 Ramanujan believed that the numerator of Bernoulli(n)/n for n even was (apart from sign) always either 1 or a prime. This is false.

References

S. Ramanujan, Some properties of Bernoulli's numbers, J. Indian Math. Soc., 3 (1911), 219-234.

Links

Table of n, a(n) for n=1..23.

K. Ono, Honoring a gift from Kumbakonam, Notices Amer. Math. Soc., 53 (2006), 640-651.

Example

As an example, Bernoulli(20)/20 = -174611/6600, but 174611 = 283*617. - Robert G. Wilson v, Jun 22 2006

Maple

A119766 := proc(nmax) local numr; for n from 2 to nmax by 2 do numr := abs(numer(bernoulli(n)/n)) ; if numr = 1 or isprime(numr) then print(n) ; fi ; od ; end : A119766(2000) ; # R. J. Mathar, Jun 21 2006

Mathematica

OldPrimeQ[n_] := Abs[n]==1 || PrimeQ[Abs[n]]; Select[Range[2000], OldPrimeQ[Numerator[BernoulliB[ # ]/# ]] &] (* T. D. Noe, Jun 20 2006 *)

Crossrefs

Cf. A112548, A001067.

Sequence in context: A309239 A152966 A066122 * A030143 A343335 A110725

Adjacent sequences: A119763 A119764 A119765 * A119767 A119768 A119769

Keyword

nonn

Author

N. J. A. Sloane, Jun 19 2006

Extensions

a(21) and a(22) from T. D. Noe, Jun 20 2006

Status

approved