A361007 a(n) = numerator of (zeta(2*n,1/4) + zeta(2*n,3/4))/Pi^(2*n) where zeta is the Hurwitz zeta function.

0, 2, 8, 64, 2176, 31744, 2830336, 178946048, 30460116992, 839461371904, 232711080902656, 39611984424992768, 955693069653508096, 1975371841521663868928, 1124025625663103358205952, 369906947004953565463576576, 278846808228005417477465964544

Offset

0,2

Comments

The function (zeta(2*n,1/4) + zeta(2*n,3/4))/Pi^(2*n) is rational for every positive integer n.

Links

Table of n, a(n) for n=0..16.

Formula

a(n) = numerator( [x^(2*n-1)] tan(2*x) ).

a(n) = numerator( (-1)^(n + 1)*4^(2*n)*(2^(2*n) - 1)*B(2*n)/(2*(2*n)!) ) where B(2*n) are Bernoulli numbers.

Example

tan(2*x) = 2*x + (8/3)*x^3 + (64/15)*x^5 + (2176/315)*x^7 + (31744/2835)*x^9 + ...

Mathematica

Table[(Zeta[2*n, 1/4] + Zeta[2*n, 3/4])/Pi^(2*n), {n, 0, 25}] //
  FunctionExpand // Numerator
Table[4^(2 k) (2^(2 k) - 1) (-1)^(k + 1) BernoulliB[2 k]/(2 (2 k)!), {k, 0, 25}] // Numerator

Prog

(PARI) my(x='x+O('x^50), v = Vec(tan(2*x)/x)); apply(numerator, vector(#v\2, k, v[2*k-1])) \\ Michel Marcus, Apr 09 2023

Crossrefs

Cf. A036279 (denominators).

Cf. A000364, A046982, A173945, A173947, A173948, A173949, A173953, A173954, A173955, A173982, A173983, A173984, A173987, A360945, A360966.

Sequence in context: A322066 A255133 A073113 * A091794 A013134 A009602

Adjacent sequences: A361004 A361005 A361006 * A361008 A361009 A361010

Keyword

nonn,frac

Author

Artur Jasinski, Mar 15 2023

Status

approved