OFFSET
1,6
EXAMPLE
MATHEMATICA
Table[FactorInteger[Abs@ Numerator[BernoulliB[2 n]/(2 n)]][[-1, 1]], {n, 25}] (* Michael De Vlieger, Jan 21 2017 *)
PROG
(PARI) a(n) = if(abs(numerator(bernfrac(2*n) / (2*n))) == 1, 1, vecmax(factor(abs(numerator(bernfrac(2*n) / (2*n))))[, 1]));
for(n=1, 25, print1(a(n), ", ")) \\ Indranil Ghosh, Mar 11 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 20 2017
EXTENSIONS
a(20)-a(25) from Michael De Vlieger, Jan 21 2017
STATUS
approved