cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-5 of 5 results.

A241885 Write the coefficient of x^n/n! in the expansion of (x/(exp(x)-1))^(1/2) as f(n)/g(n); sequence gives f(n).

Original entry on oeis.org

1, -1, 1, 1, -3, -19, 79, 275, -2339, -11813, 14217, 95265, -4634445, -193814931, 131301607, 1315505395, -3890947599, -136146236611, 46949081169401, 124889801445461, -10635113572583999, -158812278992229461, 56918172351554857, 8484151253958927197
Offset: 0

Views

Author

Jitender Singh, May 01 2014

Keywords

Comments

For g(n) see A242225(n).
The old definition was "Numerator of (B_n)^(1/2) in the Cauchy type product (sometimes known as binomial transform) where B_n is the n-th Bernoulli number".
The Nørlund polynomials N(a, n, x) with parameter a = 1/2 evaluated at x = 0 give the rational values. - Peter Luschny, Feb 18 2024

Examples

			For n=1, B_1=-1/2 and B_1^(1/2)=-1/4 so a(1)=-1.
For n=6, B_6=1/6 and B_6^(1/2)=79/86016 so a(6)=79.
1/1, -1/4, 1/48, 1/64, -3/1280, -19/3072, 79/86016, 275/49152, -2339/2949120, -11813/1310720, 14217/11534336 = A241885 / A242225.
		

Crossrefs

Cf. A242225 (denominators), A126156, A242233.

Programs

  • Maple
    g := proc(f, n) option remember; local g0, m; g0 := sqrt(f(0));
    if n=0 then g0 else if n=1 then 0 else add(binomial(n, m)*g(f,m)*g(f,n-m), m=1..n-1) fi; (f(n)-%)/(2*g0) fi end:
    a := n -> numer(g(bernoulli, n));
    seq(a(n), n = 0..23); # Peter Luschny, May 07 2014
  • Mathematica
    a := 1
    g[0] := Sqrt[f[0]]
    f[k_] := BernoulliB[k]
    g[1] := f[1]/(2 g[0]^1);
    g[k_] := (f[k] - Sum[Binomial[k, m] g[m] g[k - m], {m, 1, k - 1}])/(2 g[0])
    Table[Factor[g[k]], {k, 0, 15}] // TableForm
    (* Alternative: *)
    Table[Numerator@NorlundB[n, 1/2, 0], {n, 0, 23}]  (* Peter Luschny, Feb 18 2024 *)
  • PARI
    a(n)=numerator(sum(k=0,n,binomial(-1/2,k)*binomial(n+1/2,n-k)*stirling(n+k,k,2)/binomial(n+k,k))) \\ Tani Akinari, Oct 08 2024

Formula

Theorem: a(n)/A242225(n) = n!*A222411(n)/(A222412(n)*(-1)^n/(1-2*n)) = n!*A350194(n)/(A350154(n)*(2*n+1)). - David Broadhurst, Apr 23 2022 (see Link).
For any arithmetic function f and a positive integer k > 1, define the k-th root of f to be the arithmetic function g such that g*g*...*g(k times)=f and is determined by the following recursive formula:
g(0) = f(0)^(1/m);
g(1) = f(1)/(m*g(0)^(m-1));
g(k) = 1/(m*g(0)^(m-1))*(f(k) - Sum_{k_1+...+k_m=k,k_i= 2.
This formula is applicable for any rational root of an arithmetic function with respect to the Cauchy type product.
E.g.f: sqrt(x/(exp(x)-1)); take numerators. - Peter Luschny, May 08 2014
a(n) = numerator(Sum_{k=0..n} binomial(-1/2,k)*binomial(n+1/2,n-k)*Stirling2(n+k,k)/binomial(n+k,k)). - Tani Akinari, Oct 08 2024

Extensions

Simpler definition from N. J. A. Sloane, Apr 24 2022 at the suggestion of David Broadhurst

A222411 Numerators in Taylor series expansion of (x/(exp(x) - 1))^(3/2)*exp(x/2).

Original entry on oeis.org

1, -1, -1, 5, 7, -19, -869, 715, 2339, -200821, -12863, 2117, 7106149, -64604977, -131301607, 7629931291, 174053933, -19449462373, -46949081169401, 355455588729389, 10635113572583999, -6511303438681407901, -349640201588122693, 9112944418860287
Offset: 0

Views

Author

N. J. A. Sloane, Feb 28 2013

Keywords

Examples

			The first few fractions are 1, -1/4, -1/32, 5/384, 7/10240, -19/40960, -869/61931520, 715/49545216, ... = A222411/A222412. - _Petros Hadjicostas_, May 14 2020
		

Crossrefs

Cf. A222412 (denominators).

Programs

  • Maple
    gf:= (x/(exp(x)-1))^(3/2)*exp(x/2):
    a:= n-> numer(coeff(series(gf, x, n+3), x, n)):
    seq(a(n), n=0..25);  # Alois P. Heinz, Mar 02 2013
  • Mathematica
    Series[(x/(Exp[x]-1))^(3/2)*Exp[x/2], {x, 0, 25}] // CoefficientList[#, x]& // Numerator (* Jean-François Alcover, Mar 18 2014 *)

Formula

Theorem: A241885(n)/A242225(n) = n!*A222411(n)/(A222412(n)*(-1)^n/(1-2*n)) = n!*A350194(n)/(A350154(n)*(2*n+1)). - David Broadhurst, Apr 23 2022 (see Link).

A222412 Denominators in Taylor series expansion of (x/(exp(x) - 1))^(3/2)*exp(x/2).

Original entry on oeis.org

1, 4, 32, 384, 10240, 40960, 61931520, 49545216, 7927234560, 475634073600, 1993133260800, 177167400960, 48753634065776640, 195014536263106560, 39002907252621312000, 842462796656620339200, 2204424056667635712000, 79359266040034885632000
Offset: 0

Views

Author

N. J. A. Sloane, Feb 28 2013

Keywords

Examples

			The first few fractions are 1, -1/4, -1/32, 5/384, 7/10240, -19/40960, -869/61931520, 715/49545216, ... = A222411/A222412. - _Petros Hadjicostas_, May 14 2020
		

Crossrefs

Programs

  • Maple
    gf:= (x/(exp(x)-1))^(3/2)*exp(x/2):
    a:= n-> denom(coeff(series(gf, x, n+3), x, n)):
    seq(a(n), n=0..25);  # Alois P. Heinz, Mar 02 2013
  • Mathematica
    Series[(x/(Exp[x]-1))^(3/2)*Exp[x/2], {x, 0, 25}] // CoefficientList[#, x]& // Denominator (* Jean-François Alcover, Mar 18 2014 *)

Formula

Theorem: A241885(n)/A242225(n) = n!*A222411(n)/(A222412(n)*(-1)^n/(1-2*n)) = n!*A350194(n)/(A350154(n)*(2*n+1)). - David Broadhurst, Apr 23 2022 (see Link).

A242225 Write the coefficient of x^n/n! in the expansion of (x/(exp(x)-1))^(1/2) as f(n)/g(n); sequence gives g(n).

Original entry on oeis.org

1, 4, 48, 64, 1280, 3072, 86016, 49152, 2949120, 1310720, 11534336, 4194304, 1526726656, 2348810240, 12079595520, 3221225472, 73014444032, 51539607552, 137095356088320, 5772436045824, 3809807790243840, 725677674332160, 2023101395107840, 3166593487994880
Offset: 0

Views

Author

Jitender Singh, May 08 2014

Keywords

Comments

For f(n) see A241885(n).
The old definition was "Denominator of (B_n)^(1/2) in the Cauchy type product (sometimes known as binomial transform) where B_n is the n-th Bernoulli number".
The Nørlund polynomials N(a, n, x) with parameter a = 1/2 evaluated at x = 0 give the rational values. - Peter Luschny, Feb 18 2024

Examples

			For n=1, B_1=-1/2 and B_1^(1/2)=-1/4 so a(1)=4.
For n=6, B_6=1/6 and B_6^(1/2)=79/86016 so a(6)=86016.
		

Crossrefs

Programs

  • Maple
    g := proc(f, n) option remember; local g0, m; g0 := sqrt(f(0));
    if n=0 then g0 else if n=1 then 0 else add(binomial(n, m)*g(f, m)*g(f, n-m), m=1..n-1) fi; (f(n)-%)/(2*g0) fi end:
    a := n -> denom(g(bernoulli, n));
    seq(a(n), n=0..23);
  • Mathematica
    a := 1
    g[0] := Sqrt[f[0]]
    f[k_] := BernoulliB[k]
    g[1] := f[1]/(2 g[0]^1);
    g[k_] := (f[k] - Sum[Binomial[k, m] g[m] g[k - m], {m, 1, k - 1}])/(2 g[0])
    Table[Denominator[Factor[g[k]]], {k, 0, 15}] // TableForm
    (* Alternative: *)
    Table[Denominator@NorlundB[n, 1/2, 0], {n, 0, 23}] (* Peter Luschny, Feb 18 2024 *)

Formula

Theorem: A241885(n)/A242225(n) = n!*A222411(n)/(A222412(n)*(-1)^n/(1-2*n)) = n!*A350194(n)/(A350154(n)*(2*n+1)). - David Broadhurst, Apr 23 2022 (see Link).
For any arithmetic function f and a positive integer k>1, define the k-th root of f to be the arithmetic function g such that g*g*...*g(k times)=f and is determined by the following recursive formula:
g(0)= f(0)^{1/m};
g(1)= f(1)/(mg(0)^(m-1));
g(k)= 1/(m g(0)^{m-1})*(f(k)-sum_{k_1+...+k_m=k,k_i=2.
This formula is applicable for any rational root of an arithmetic function with respect to the Cauchy type product.

Extensions

Simpler definition from N. J. A. Sloane, Apr 24 2022 at the suggestion of David Broadhurst.

A350194 Numerators of a power series characterizing how powers of the cosine function converge to the Gaussian function.

Original entry on oeis.org

1, -1, 1, 1, -1, -19, 79, 11, -2339, -11813, 677, 2117, -308963, -64604977, 131301607, 263101079, -5614643, -1768132943, 46949081169401, 9606907803497, -10635113572583999, -158812278992229461, 8131167478793551, 9112944418860287, -40395223967437706149
Offset: 0

Views

Author

Robert B Fowler, Dec 19 2021

Keywords

Comments

See A350154 for the denominators of this sequence of rational coefficients, as well as relevant comments, formulae, and examples.

Crossrefs

Cf. A350154.

Formula

Theorem: A241885(n)/A242225(n) = n!*A222411(n)/(A222412(n)*(-1)^n/(1-2*n)) = n!*A350194(n)/(A350154(n)*(2*n+1)). - David Broadhurst, Apr 23 2022 (see Link).
Showing 1-5 of 5 results.