OFFSET
1,2
LINKS
Daniel Suteu, Table of n, a(n) for n = 1..85
Eric Weisstein's World of Mathematics, Bernoulli Number
Wikipedia, Von Staudt-Clausen theorem
MATHEMATICA
Reap[For[n = record = 0, n < 1000, n = n + 2, If[(d = Denominator[BernoulliB[n]]) > record, Sow[d]; record = d]]][[2, 1]] (* Jean-François Alcover, Nov 09 2012 *)
PROG
(PARI) b(n) = if((n==0) || (n>1 && n%2==1), 1, my(d=divisors(n)); prod(k=1, #d, if(isprime(d[k]+1), d[k]+1, 1))); \\ more efficient than denominator(bernfrac(n))
lista(n) = { my(m=0); for(k=0, n, my(d=b(k)); if(d>m, m=d; print1(d, ", "))); }
lista(1000); \\ Daniel Suteu, Dec 22 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Eric W. Weisstein, Nov 08 2004
EXTENSIONS
a(21) from Seiichi Manyama, Jan 21 2017
STATUS
approved