OFFSET
0,4
COMMENTS
Integer component of the numerator of a close variant of Euler's infinite prime product zeta(2n) = Product_{k>=1} (prime(k)^(2n))/(prime(k)^(2n)-1), namely with all minus signs changed into plus signs, as follows: zeta(4n)/zeta(2n) = Product_{k>=1} (prime(k)^(2n))/(prime(k)^(2n)+1). The transcendental component is Pi^(2n).
For a detailed account of the results, including proof and relation to the zeta function, see Links for the PDF file submitted as supporting material.
The reference to Apostol is to a discussion of the equivalence of 1) zeta(2s)/zeta(s) and 2) a related infinite prime product, that is, Product_{sigma>1} prime(n)^s/(prime(n)^s + 1), with s being a complex variable such that s = sigma + i*t where sigma and t are real (following Riemann), using a type of proof different from the one posted below involving zeta(4n)/zeta(2n). On this, see also Hardy and Wright cited below. - Leo Depuydt, Nov 22 2013, Nov 27 2013
The background of the sequence is now described in the link below to L. Depuydt, The Prime Sequence ... . - Leo Depuydt, Aug 22 2014
From Robert Israel, Aug 22 2014: (Start)
Numerator of (-1)^n*B(4*n)*4^n*(2*n)!/(B(2*n)*(4*n)!), where B(n) are the Bernoulli numbers (see A027641 and A027642).
Not the same as abs(A001067(2*n)): they differ first at n=17.
(End)
REFERENCES
T. M. Apostol, Introduction to Analytic Number Theory, Springer, 1976, p. 231.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Fourth Edition, Clarendon Press, 1960, p. 255.
LINKS
MAPLE
seq(numer((-1)^n*bernoulli(4*n)*4^n*(2*n)!/(bernoulli(2*n)*(4*n)!)), n=0..100); # Robert Israel, Aug 22 2014
MATHEMATICA
Numerator[Table[Zeta[4n]/(Zeta[2n] * Pi^(2n)), {n, 0, 15}]] (* T. D. Noe, Nov 18 2013 *)
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Leo Depuydt, Nov 07 2013
STATUS
approved