The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A360945 a(n) = numerator of (Zeta(2*n+1,1/4) - Zeta(2*n+1,3/4))/Pi^(2*n+1) where Zeta is the Hurwitz zeta function. 2
1, 2, 10, 244, 554, 202084, 2162212, 1594887848, 7756604858, 9619518701764, 59259390118004, 554790995145103208, 954740563911205348, 32696580074344991138888, 105453443486621462355224, 7064702291984369672858925136, 4176926860695042104392112698 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
COMMENTS
The function (Zeta(2*n+1,1/4) - Zeta(2*n+1,3/4))/Pi^(2*n+1) is rational for every positive integer n.
For denominators see A360966.
(Zeta(2*n+1,1/4) + Zeta(2*n+1,3/4))/Zeta(2*n+1) = 4*16^n - 2*4^n; see A193475.
For numerators of the function (Zeta(2*n,1/4) + Zeta(2*n,3/4))/Pi^(2*n) see A361007.
For denominators of the function (Zeta(2*n,1/4) + Zeta(2*n,3/4))/Pi^(2*n) see A036279.
(Zeta(2*n,1/4) - Zeta(2*n,3/4))/beta(2*n) = 16^n (see A001025) where beta is the Dirichlet beta function.
From the above formulas we can express Zeta(k,1/4) and Zeta(k,3/4) for every positive integer k.
LINKS
FORMULA
a(n) = A046982(2*n).
(Zeta(2*n + 1, 1/4) - Zeta(2*n + 1, 3/4))/(Pi^(2*n + 1)) = A000364(n)*(2*n + 1)*2^(2*n)/(2*n + 1)!.
EXAMPLE
a(0) = 1 because lim_{n->0} (Zeta(2*n+1,1/4) - Zeta(2*n+1,3/4))/Pi^(2*n+1) = 1.
a(3) = 244 because (Zeta(2*3+1,1/4) - Zeta(2*3+1,3/4))/Pi^(2*3+1) = 244/45.
MATHEMATICA
Table[(Zeta[2*n + 1, 1/4] - Zeta[2*n + 1, 3/4]) / Pi^(2*n + 1), {n, 1, 25}] // FunctionExpand // Numerator (* Vaclav Kotesovec, Feb 27 2023 *)
t[0, 1] = 1; t[0, _] = 0;
t[n_, k_] := t[n, k] = (k-1) t[n-1, k-1] + (k+1) t[n-1, k+1];
a[n_] := Sum[t[2n, k]/(2n)!, {k, 0, 2n+1}] // Numerator;
Table[a[n], {n, 0, 16}] (* Jean-François Alcover, Mar 15 2023 *)
a[n_] := SeriesCoefficient[Tan[x+Pi/4], {x, 0, 2n}] // Numerator;
Table[a[n], {n, 0, 16}] (* Jean-François Alcover, Apr 15 2023 *)
PROG
(PARI) a(n) = numerator(abs(eulerfrac(2*n))*(2*n + 1)*2^(2*n)/(2*n + 1)!); \\ Michel Marcus, Apr 11 2023
CROSSREFS
Sequence in context: A346222 A289948 A282567 * A308756 A225371 A088310
KEYWORD
nonn,frac
AUTHOR
Artur Jasinski, Feb 26 2023
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 24 11:42 EDT 2024. Contains 372773 sequences. (Running on oeis4.)