Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A166333
The largest prime that divides A027642(n) (the denominator of the Bernoulli number B_n), or 1 if A027642(n) is 1.
1
1, 2, 3, 1, 5, 1, 7, 1, 5, 1, 11, 1, 13, 1, 3, 1, 17, 1, 19, 1, 11, 1, 23, 1, 13, 1, 3, 1, 29, 1, 31, 1, 17, 1, 3, 1, 37, 1, 3, 1, 41, 1, 43, 1, 23, 1, 47, 1, 17, 1, 11, 1, 53, 1, 19, 1, 29, 1, 59, 1, 61, 1, 3, 1, 17, 1, 67, 1, 5, 1, 71, 1, 73, 1, 3, 1, 5, 1, 79, 1, 41, 1, 83, 1, 43, 1, 3, 1, 89, 1
OFFSET
0,2
COMMENTS
The largest member of the extended prime list A008578 which divides the denominator of Bernoulli(n).
Essentially A073409 padded with 1's.
LINKS
FORMULA
a(n) = A006530(A027642(n)). - Antti Karttunen, Dec 19 2018
PROG
(PARI)
A006530(n) = if(n>1, vecmax(factor(n)[, 1]), 1);
A166333(n) = A006530(denominator(bernfrac(n))); \\ Antti Karttunen, Dec 19 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul Curtz, Oct 12 2009
EXTENSIONS
Edited and extended by R. J. Mathar, Oct 21 2009
Name and comment swapped by Antti Karttunen, Dec 19 2018
STATUS
approved