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.

A286854 Numbers k such that k == 1 or -1 (mod 6) but k does not divide the numerator of Bernoulli(2*k).

This page as a plain text file.
%I A286854 #23 Feb 09 2021 01:55:10
%S A286854 55,253,275,385,605,715,935,1045,1081,1265,1375,1595,1705,1711,1771,
%T A286854 1925,2035,2255,2365,2485,2585,2695,2783,2915,3025,3245,3289,3355,
%U A286854 3403,3575,3685,3905,4015,4235,4301,4345,4565,4675,4807,4895,5005,5225,5335,5405,5555
%N A286854 Numbers k such that k == 1 or -1 (mod 6) but k does not divide the numerator of Bernoulli(2*k).
%p A286854 isa := n -> abs(mods(n, 6)) = 1 and modp(numer(bernoulli(2*n)), n) <> 0:
%p A286854 select(isa, [$1..2255]); # _Peter Luschny_, Aug 02 2017
%t A286854 Select[Range@9999,0 != Mod[Numerator@BernoulliB[2 #], #] && MemberQ[{1, 5}, Mod[#, 6]] &]
%o A286854 (PARI) isok(n) = (((n % 6) == 1) || ((n % 6) == 5)) && (numerator(bernfrac(2*n)) % n); \\ _Michel Marcus_, Aug 02 2017
%Y A286854 Cf. A000367, A286853 (1st differences).
%K A286854 nonn
%O A286854 1,1
%A A286854 _Bill Gosper_, Aug 01 2017