login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A228037 Odd-index terms of Agoh's congruence A046094: a(n) is conjectured to be 1 iff 2n+1 is prime. 4
0, 1, 1, 1, 3, 1, 1, 5, 1, 1, 7, 1, 5, 9, 1, 1, 11, 0, 1, 13, 1, 1, 24, 1, 7, 17, 1, 0, 19, 1, 1, 21, 13, 1, 23, 1, 1, 25, 0, 1, 27, 1, 17, 29, 1, 13, 31, 0, 1, 33, 1, 1, 56, 1, 1, 37, 1, 0, 39, 0, 11, 41, 25, 1, 43, 1, 19, 45, 1, 1, 47, 0, 29, 49, 1, 1, 51, 0, 1, 53, 0, 1, 88, 1, 13, 57, 1, 25, 59, 1, 1, 61, 37, 0, 63, 1, 1, 65, 1 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,5

COMMENTS

Except for A046094(2) = 1, the even-index terms of A046094 are all zero since Bernoulli(2n+1) = 0 for n > 0.

LINKS

Alois P. Heinz, Table of n, a(n) for n = 0..10000

FORMULA

a(n) = - (2n+1)*Bernoulli(2n) mod 2n+1.

EXAMPLE

-(2*1+1)*Bernoulli(2*1) = -3*(1/6) = -1/2 == -2 == 1 mod 3, so a(1) = 1.

MAPLE

a:= n-> -(2*n+1)*bernoulli(2*n) mod (2*n+1):

seq(a(n), n=0..100);  # Alois P. Heinz, Aug 13 2013

MATHEMATICA

a[ n_ ] := Mod[ Numerator[ -(2 n + 1)* BernoulliB[ 2 n]] * PowerMod[ Denominator[(2 n + 1)* BernoulliB[ 2 n]], -1, 2 n + 1], 2 n + 1]

CROSSREFS

a(n) = A046094(2n+1).

Sequence in context: A087501 A294951 A101443 * A184726 A046230 A046229

Adjacent sequences:  A228034 A228035 A228036 * A228038 A228039 A228040

KEYWORD

nonn

AUTHOR

Jonathan Sondow, Aug 13 2013

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 October 21 01:18 EDT 2019. Contains 328291 sequences. (Running on oeis4.)