OFFSET
1,3
LINKS
Daniel Suteu, Table of n, a(n) for n = 1..90
Eric Weisstein's World of Mathematics, Bernoulli Number
Wikipedia, Von Staudt-Clausen theorem
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(k, ", "))); }
lista(100000); \\ Daniel Suteu, Dec 23 2018
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Eric W. Weisstein, Nov 08 2004
EXTENSIONS
Corrected and extended by Daniel Suteu, Dec 23 2018
STATUS
approved