OFFSET
1,1
COMMENTS
FORMULA
a(n) = denominator( ((2*abs(bernoulli(4*n))*((2*n)!)^2)/(bernoulli(2*n)^2*(4*n)!))*(4^n + 1)/(4^n - 1) ).
a(n) = denominator( (zeta(4*n)/zeta(2*n)^2) * (2^(4*n) + 2^(2*n)) / (2^(4*n) - 2^(2*n)) ).
a(n) = denominator(Product_{p prime > 2} (p^(2*n) - 1) / (p^(2*n) + 1)).
EXAMPLE
For n = 2: (zeta(8)/zeta(4)^2) * (4^2 + 1)/(4^2 - 1) = ((Pi^8/9450) / (Pi^4/90)^2) * (17/15) = (6/7) * (17/15) = 102/105 = 34/35. The denominator is 35.
For n = 3: (zeta(12)/zeta(6)^2) * (4^3 + 1)/(4^3 - 1) = ((691*Pi^12/638512875) / (Pi^6/945)^2) * (65/63) = (691/715) * (65/63) = 691/693. The denominator is 693.
PROG
(PARI) {a(n) = denominator((2*abs(bernfrac(4*n))*(2*n)!^2) / (bernfrac(2*n)^2*(4*n)!) * (4^n+1)/(4^n-1))}
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Dimitris Valianatos, Apr 16 2026
EXTENSIONS
Incorrect a(14)-a(15) removed by Sean A. Irvine, May 05 2026
a(14) from Dimitris Valianatos, May 06 2026
STATUS
approved
