|
|
A100195
|
|
Numbers n such that the denominator of BernoulliB(n) is a record.
|
|
5
|
|
|
0, 1, 2, 4, 6, 10, 12, 30, 36, 60, 72, 108, 120, 144, 180, 240, 360, 420, 540, 840, 1008, 1080, 1200, 1260, 1620, 1680, 2016, 2160, 2520, 3360, 3780, 5040, 6480, 7560, 8400, 10080, 12600, 15120, 25200, 30240, 42840, 45360, 55440, 60480, 75600, 85680, 100800
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,3
|
|
LINKS
|
|
|
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, ", "))); }
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|