OFFSET
1,1
COMMENTS
Sequence given for a(1)-a(15) in Note 6.2, p. 13 of Dixit and others.
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..500
M. W. Coffey, V. de Angelis, A. Dixit, V. H. Moll, et al., The Zagier polynomials. Part II: Arithmetic properties of coefficients, arXiv:1303.6590 [math.NT], 2013.
A. Dixit, V. H. Moll, Ch. Vignat, The Zagier modification of Bernoulli numbers and a polynomial extension. Part I, arXiv:1209.4110v1 [math.NT], Sep 18, 2012.
MAPLE
A216912 := n -> denom(add(binomial(2*n+r, 2*r)*bernoulli(r)/(2*n+r), r=0..2*n))/4;
seq(A216912(i), i=1..38); # Peter Luschny, Sep 20 2012
MATHEMATICA
a[n_] := Denominator[Sum[Binomial[2n+r, 2r]*(BernoulliB[r]/(2n+r)), {r, 0, 2n}]]/4;
Array[a, 38] (* Jean-François Alcover, Jul 14 2018, after Peter Luschny *)
PROG
(PARI) a(n) = denominator(sum(k=0, 2*n, binomial(2*n+k, 2*k)*bernfrac(k)/(2*n+k)))/4; \\ Michel Marcus, Jul 14 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Jonathan Vos Post, Sep 20 2012
EXTENSIONS
a(16)-a(38) from Peter Luschny, Sep 20 2012
STATUS
approved