A111008 a(n) = A000367(n)/A141590(n).

1, 1, 1, 1, 1, 5, 1, 7, 1, 1, 1, 11, 1, 13, 7, 5, 1, 17, 1, 19, 1, 1, 11, 23, 1, 25, 13, 1, 7, 29, 1, 31, 1, 11, 629, 35, 1, 37, 19, 13, 1, 41, 1, 43, 11, 5, 23, 47, 1, 49, 1, 1003, 481, 53, 1, 5, 7, 19, 29, 59, 1, 61, 2077, 103, 1, 65, 11, 67, 17, 23, 259, 71, 1, 73, 37, 25, 2489, 77

Offset

0,6

Links

Jeppe Stig Nielsen, Table of n, a(n) for n = 0..10000

Maple

A120082 := proc(n) local b; if n = 0 then b := 1 ; elif n = 1 then b := -1/4 ; elif type(n, 'odd') then b := 0; else b := bernoulli(n)/(n+1)! ; fi; numer(b) ; end:
A141590 := proc(n) A120082(2*n) ; end: A000367 := proc(n) numer(bernoulli(2*n)) ; end:
A111008 := proc(n) A000367(n)/A141590(n) ; end: seq(A111008(n), n=0..120) ; # R. J. Mathar, Sep 03 2009

Prog

(PARI) upto(N)=bernvec(N); forstep(n=0, 2*N, 2, print1(gcd(numerator(bernfrac(n)), (n+1)!), ", ")) \\ Jeppe Stig Nielsen, Jun 22 2023

Crossrefs

See A141517.

Cf. A300711.

Sequence in context: A348500 A308090 A300711 * A065330 A140215 A366284

Adjacent sequences: A111005 A111006 A111007 * A111009 A111010 A111011

Keyword

nonn

Author

Paul Curtz, Aug 25 2008

Extensions

Edited and extended by R. J. Mathar, Sep 03 2009

Status

approved