login
A111008
a(n) = A000367(n)/A141590(n).
4
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
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
KEYWORD
nonn
AUTHOR
Paul Curtz, Aug 25 2008
EXTENSIONS
Edited and extended by R. J. Mathar, Sep 03 2009
STATUS
approved