OFFSET
1,1
COMMENTS
The sequence D = (denominator(S(k, 2n+1)))_{k>0} is periodic for any n>0 . i.e. for n=2, D has period {1, 4, 12, 6, 6, 12, 12, 3, 1, 4, 4, 6, 6, 12, 12, 3, 3, 4, 4, 2, 6, 12, 12, 3, 3, 12, 4, 2, 2, 12, 12, 3, 3, 12, 12, 2, 2, 4, 12, 3, 3, 12, 12, 6, 2, 4, 4, 3, 3, 12, 12, 6, 6, 4, 4, 1, 3, 12, 12, 6, 6, 12, 4, 1, 1, 12, 12, 6, 6, 12, 12, 1} of length 72 and reaches 12 as maximum value, hence a(2)=12
FORMULA
a(n) = denominator(B(2*n)*(n+1/2)) where B() is Bernoulli number. - Michael Somos, Sep 05 2020
PROG
(PARI) H(n, r)=sum(i=1, n, 1/i^r); S(n, s)=sum(k=1, n, k^s*H(k, 2))-H(n, 2)*H(n, -s); a(n)=vecmax(vector(100, i, denominator(S(i, 2*n+1))))
(PARI) {a(n) = if(n<1, 0, denominator(bernfrac(2*n)*(n+1/2)))}; /* Michael Somos, Sep 05 2020 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Aug 10 2004
STATUS
approved