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-10 of 20 results. Next

A056981 a(n) = A002596(n)^2.

Original entry on oeis.org

1, 1, 1, 1, 25, 49, 441, 1089, 184041, 511225, 5909761, 17631601, 863948449, 2704312009, 34493775625, 111759833025, 93990019574025, 312541206957225, 4201942893536025, 14258670483605625, 780804795682244025
Offset: 0

Views

Author

Keywords

Comments

Numerator of h^(2n) in the Kummer-Gauss series for the perimeter of an ellipse.

Crossrefs

Programs

  • Mathematica
    Numerator[Table[Binomial[1/2, n]^2, {n, 50}]] (* Artur Jasinski, Aug 30 2016 *)

Extensions

New description from N. J. A. Sloane, Feb 18 2004

A046161 a(n) = denominator of binomial(2n,n)/4^n.

Original entry on oeis.org

1, 2, 8, 16, 128, 256, 1024, 2048, 32768, 65536, 262144, 524288, 4194304, 8388608, 33554432, 67108864, 2147483648, 4294967296, 17179869184, 34359738368, 274877906944, 549755813888, 2199023255552, 4398046511104, 70368744177664, 140737488355328, 562949953421312
Offset: 0

Views

Author

Keywords

Comments

Also denominator of e(0,n) (see Maple line). - N. J. A. Sloane, Feb 16 2002
Denominator of coefficient of x^n in (1+x)^(k/2) or (1-x)^(k/2) for any odd integer k. - Michael Somos, Sep 15 2004
Numerator of binomial(2n,n)/4^n = A001790(n).
Denominators in expansion of sqrt(c(x)), c(x) the g.f. of A000108. - Paul Barry, Jul 12 2005
Denominator of 2^m*Gamma(m+3/4)/(Gamma(3/4)*Gamma(m+1)). - Stephen Crowley, Mar 19 2007
Denominator in expansion of Jacobi_P(n,1/2,1/2,x). - Paul Barry, Feb 13 2008
This sequence equals the denominators of the coefficients of the series expansions of (1-x)^((-1-2*n)/2) for all integer values of n; see A161198 for detailed information. - Johannes W. Meijer, Jun 08 2009
Numerators of binomial transform of 1, -1/3, 1/5, -1/7, 1/9, ... (Madhava-Gregory-Leibniz series for Pi/4): 1, 2/3, 8/15, 16/35, 128/315, 256/693, .... First differences are -1/3, -2/15, -8/105, -16/315, -128/3465, -256/9009, ... which contain the same numerators, negated. The second differences are 1/5, 2/35, 8/315, 16/1155, 128/15015, ... again with the same numerators. Second column: 2/3, -2/15, 2/35, -2/63, 2/99; see A000466(n+1) = A005563(2n+1). Third column: 8*(1/15, -1/105, 1/315, -1/693, ...), see A061550. See A173294 and A173296. - Paul Curtz, Feb 16 2010
0, 1, 5/3, 11/5, 93/35, 193/63, 793/231, ... = (0 followed by A120778(n))/A001790(n) is the binomial transform of 0, 1, -1/3, 1/5, -1/7, 1/9, ... . See A173755 and formula below. - Paul Curtz, Mar 13 2013
Numerator of power series of arcsin(x)/sqrt(1-x^2), centered at x=0. - John Molokach, Aug 02 2013
Denominators of coefficients in the Taylor series expansion of Sum_{n>=0} exp((-1)^n * Euler(2*n)*x^n/(2*n)), see A280442 for the numerators. - Johannes W. Meijer, Jan 05 2017
Denominators of Pochhammer(n+1, -1/2)/sqrt(Pi). - Adam Hugill, Sep 11 2022
a(n) is the denominator of the mean value of cos(x)^(2*n) from x = 0 to 2*Pi. - Stefano Spezia, May 16 2023
4^n/binomial(2n,n) is the expected value of the number of socks that are randomly drawn out of a drawer of n different pairs of socks, when one sock is drawn out at a time until a matching pair is found (King and King, 2005). - Amiram Eldar, Jul 02 2023
a(n) is the denominator of (1/Pi) * Integral_{x=-oo..+oo} sech(x)^(2*n+1) dx. The corresponding numerator is A001790(n). - Mohammed Yaseen, Jul 29 2023
a(n) is the numerator of Integral_{x=0..Pi/2} sin(x)^(2*n+1) dx. The corresponding denominator is A001803(n). - Mohammed Yaseen, Sep 22 2023

Examples

			sqrt(1+x) = 1 + (1/2)*x - (1/8)*x^2 + (1/16)*x^3 - (5/128)*x^4 + (7/256)*x^5 - (21/1024)*x^6 + (33/2048)*x^7 + ...
binomial(2n,n)/4^n => 1, 1/2, 3/8, 5/16, 35/128, 63/256, 231/1024, 429/2048, 6435/32768, ...
The sequence e(0,n) begins 1, 3/2, 21/8, 77/16, 1155/128, 4389/256, 33649/1024, 129789/2048, 4023459/32768, ...
		

References

  • W. Feller, An Introduction to Probability Theory and Its Applications, Vol. 1, 2nd ed. New York: Wiley, 1968; Chap. III, Eq. 4.1.
  • B. D. Hughes, Random Walks and Random Environments, Oxford 1995, vol. 1, p. 513, Eq. (7.282).
  • Eli Maor, e: The Story of a Number. Princeton, New Jersey: Princeton University Press (1994), p. 72.
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 6, equations 6:14:5 - 6:14:9 at pages 50-51.

Crossrefs

Cf. A161198 triangle related to the series expansions of (1-x)^((-1-2*n)/2) for all values of n.

Programs

  • Magma
    [Denominator(Binomial(2*n,n)/4^n): n in [0..30]]; // Vincenzo Librandi, Jul 18 2015
    
  • Maple
    e := proc(l,m) local k; add(2^(k-2*m)*binomial(2*m-2*k,m-k)* binomial(m+k, m) *binomial(k,l), k=l..m); end: seq(denom(e(0,n)), n = 0..24);
    Z[0]:=0: for k to 30 do Z[k]:=simplify(1/(2-z*Z[k-1])) od: g:=sum((Z[j]-Z[j-1]), j=1..30): gser:=series(g, z=0, 27): seq(denom(coeff(gser, z, n)), n=-1..23); # Zerinvary Lajos, May 21 2008
    A046161 := proc(n) option remember: if n = 0 then 1 else 2^A001511(n) * procname(n-1) fi: end: A001511 := proc(n): padic[ordp](2*n, 2) end: seq(A046161(n), n = 0..24); # Johannes W. Meijer, Nov 04 2012
    A046161 := n -> 4^n/2^add(i,i=convert(n, base, 2)):
    seq(A046161(n), n=0..24); # Peter Luschny, Apr 08 2014
  • Mathematica
    a[n_, m_] := Binomial[n - m/2 + 1, n - m + 1] - Binomial[n - m/2, n - m + 1]; s[n_] := Sum[ a[n, k], {k, 0, n}]; Table [Denominator[s[n]], {n, 0, 26}] (* Michele Dondi (bik.mido(AT)tiscalinet.it), Jul 11 2002 *)
    Denominator[Table[Binomial[2n,n]/4^n,{n,0,30}]] (* Harvey P. Dale, Oct 29 2012 *)
    Table[Denominator@LegendreP[2n,0],{n,0,24}] (* Andres Cicuttin, Jan 22 2018 *)
  • Maxima
    a(n) := denom(binomial(-1/2,n));
    makelist(a(n),n,0,24); /* Peter Luschny, Nov 21 2012 */
    
  • PARI
    a(n)=if(n<0,0,denominator(binomial(2*n,n)/4^n)) /* Michael Somos, Sep 15 2004 */
    
  • PARI
    a(n)=my(s=n);while(n>>=1,s+=n);2^s \\ Charles R Greathouse IV, Apr 07 2012
    
  • PARI
    a(n)=denominator(I^-n*pollegendre(n,I/2)) \\ Charles R Greathouse IV, Mar 18 2017
    
  • Python
    def A046161(n): return 1<<(n<<1)-n.bit_count() # Chai Wah Wu, Nov 15 2022
  • Sage
    def A046161(n):
        A005187 = lambda n: A005187(n//2) + n if n > 0 else 0
        return 2^A005187(n)
    [A046161(n) for n in (0..24)]  # Peter Luschny, Nov 16 2012
    

Formula

a(n) = 2^(2*n - 1 - A048881(n-1)), if n > 0.
a(n) = 2^A005187(n).
a(n) = 4^n/2^A000120(n). - Michael Somos, Sep 15 2004
a(n) = 2^A001511(n)*a(n-1) with a(0) = 1. - Johannes W. Meijer, Nov 04 2012
a(n) = denominator(binomial(-1/2,n)). - Peter Luschny, Nov 21 2012
a(n) = (0 followed by A120778(n)) + A001790(n). - Paul Curtz, Mar 13 2013
a(n) = 2^n*A060818(n). - Johannes W. Meijer, Jan 05 2017
a(n)/A001790(n) ~ sqrt(Pi*n) (King and King, 2005). - Amiram Eldar, Jul 02 2023

A001790 Numerators in expansion of 1/sqrt(1-x).

Original entry on oeis.org

1, 1, 3, 5, 35, 63, 231, 429, 6435, 12155, 46189, 88179, 676039, 1300075, 5014575, 9694845, 300540195, 583401555, 2268783825, 4418157975, 34461632205, 67282234305, 263012370465, 514589420475, 8061900920775, 15801325804719
Offset: 0

Views

Author

Keywords

Comments

Also numerator of e(n-1,n-1) (see Maple line).
Leading coefficient of normalized Legendre polynomial.
Common denominator of expansions of powers of x in terms of Legendre polynomials P_n(x).
Also the numerator of binomial(2*n,n)/2^n. - T. D. Noe, Nov 29 2005
This sequence gives the numerators of the Maclaurin series of the Lorentz factor (see Wikipedia link) of 1/sqrt(1-b^2) = dt/dtau where b=u/c is the velocity in terms of the speed of light c, u is the velocity as observed in the reference frame where time t is measured and tau is the proper time. - Stephen Crowley, Apr 03 2007
Truncations of rational expressions like those given by the numerator operator are artifacts in integer formulas and have many disadvantages. A pure integer formula follows. Let n$ denote the swinging factorial and sigma(n) = number of '1's in the base-2 representation of floor(n/2). Then a(n) = (2*n)$ / sigma(2*n) = A056040(2*n) / A060632(2*n+1). Simply said: this sequence is the odd part of the swinging factorial at even indices. - Peter Luschny, Aug 01 2009
It appears that a(n) = A060818(n)*A001147(n)/A000142(n). - James R. Buddenhagen, Jan 20 2010
The convolution of sequence binomial(2*n,n)/4^n with itself is the constant sequence with all terms = 1.
a(n) equals the denominator of Hypergeometric2F1[1/2, n, 1 + n, -1] (see Mathematica code below). - John M. Campbell, Jul 04 2011
a(n) = numerator of (1/Pi)*Integral_{x=-oo..+oo} 1/(x^2-2*x+2)^n dx. - Leonid Bedratyuk, Nov 17 2012
a(n) = numerator of the mean value of cos(x)^(2*n) from x = 0 to 2*Pi. - Jean-François Alcover, Mar 21 2013
Constant terms for normalized Legendre polynomials. - Tom Copeland, Feb 04 2016
From Ralf Steiner, Apr 07 2017: (Start)
By analytic continuation to the entire complex plane there exist regularized values for divergent sums:
a(n)/A060818(n) = (-2)^n*sqrt(Pi)/(Gamma(1/2 - n)*Gamma(1 + n)).
Sum_{k>=0} a(k)/A060818(k) = -i.
Sum_{k>=0} (-1)^k*a(k)/A060818(k) = 1/sqrt(3).
Sum_{k>=0} (-1)^(k+1)*a(k)/A060818(k) = -1/sqrt(3).
a(n)/A046161(n) = (-1)^n*sqrt(Pi)/(Gamma(1/2 - n)*Gamma(1 + n)).
Sum_{k>=0} (-1)^k*a(k)/A046161(k) = 1/sqrt(2).
Sum_{k>=0} (-1)^(k+1)*a(k)/A046161(k) = -1/sqrt(2). (End)
a(n) = numerator of (1/Pi)*Integral_{x=-oo..+oo} 1/(x^2+1)^n dx. (n=1 is the Cauchy distribution.) - Harry Garst, May 26 2017
Let R(n, d) = (Product_{j prime to d} Pochhammer(j / d, n)) / n!. Then the numerators of R(n, 2) give this sequence and the denominators are A046161. For d = 3 see A273194/A344402. - Peter Luschny, May 20 2021
Using WolframAlpha, it appears a(n) gives the numerator in the residues of f(z) = 2z choose z at odd negative half integers. E.g., the residues of f(z) at z = -1/2, -3/2, -5/2 are 1/(2*Pi), 1/(16*Pi), and 3/(256*Pi) respectively. - Nicholas Juricic, Mar 31 2022
a(n) is the numerator of (1/Pi) * Integral_{x=-oo..+oo} sech(x)^(2*n+1) dx. The corresponding denominator is A046161. - Mohammed Yaseen, Jul 29 2023
a(n) is the numerator of (1/Pi) * Integral_{x=0..Pi/2} sin(x)^(2*n) dx. The corresponding denominator is A101926(n). - Mohammed Yaseen, Sep 19 2023

Examples

			1, 1, 3/2, 5/2, 35/8, 63/8, 231/16, 429/16, 6435/128, 12155/128, 46189/256, ...
binomial(2*n,n)/4^n => 1, 1/2, 3/8, 5/16, 35/128, 63/256, 231/1024, 429/2048, 6435/32768, ...
		

References

  • P. J. Davis, Interpolation and Approximation, Dover Publications, 1975, p. 372.
  • W. Feller, An Introduction to Probability Theory and Its Applications, Vol. 1, 2nd ed. New York: Wiley, 1968; Chap. III, Eq. 4.1.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 6, equation 6:14:6 at page 51.
  • J. V. Uspensky and M. A. Heaslet, Elementary Number Theory, McGraw-Hill, NY, 1939, p. 102.

Crossrefs

Cf. A060818 (denominator of binomial(2*n,n)/2^n), A061549 (denominators).
Cf. A123854 (denominators).
Cf. A161198 (triangle of coefficients for (1-x)^((-1-2*n)/2)).
Cf. A163590 (odd part of the swinging factorial).
Cf. A001405.
First column and diagonal 1 of triangle A100258.
Bisection of A036069.
Bisections give A061548 and A063079.
Inverse Moebius transform of A180403/A046161.
Numerators of [x^n]( (1-x)^(p/2) ): A161202 (p=5), A161200 (p=3), A002596 (p=1), this sequence (p=-1), A001803 (p=-3), A161199 (p=-5), A161201 (p=-7).

Programs

  • Magma
    A001790:= func< n | Numerator((n+1)*Catalan(n)/4^n) >;
    [A001790(n): n in [0..40]]; // G. C. Greubel, Sep 23 2024
  • Maple
    e := proc(l,m) local k; add(2^(k-2*m)*binomial(2*m-2*k,m-k)*binomial(m+k,m)*binomial(k,l),k=l..m); end;
    # From Peter Luschny, Aug 01 2009: (Start)
    swing := proc(n) option remember; if n = 0 then 1 elif irem(n, 2) = 1 then swing(n-1)*n else 4*swing(n-1)/n fi end:
    sigma := n -> 2^(add(i,i=convert(iquo(n,2),base,2))):
    a := n -> swing(2*n)/sigma(2*n); # (End)
    A001790 := proc(n) binomial(2*n, n)/4^n ; numer(%) ; end proc : # R. J. Mathar, Jan 18 2013
  • Mathematica
    Numerator[ CoefficientList[ Series[1/Sqrt[(1 - x)], {x, 0, 25}], x]]
    Table[Denominator[Hypergeometric2F1[1/2, n, 1 + n, -1]], {n, 0, 34}]   (* John M. Campbell, Jul 04 2011 *)
    Numerator[Table[(-2)^n*Sqrt[Pi]/(Gamma[1/2 - n]*Gamma[1 + n]),{n,0,20}]] (* Ralf Steiner, Apr 07 2017 *)
    Numerator[Table[Binomial[2n,n]/2^n, {n, 0, 25}]] (* Vaclav Kotesovec, Apr 07 2017 *)
    Table[Numerator@LegendreP[2 n, 0]*(-1)^n, {n, 0, 25}] (* Andres Cicuttin, Jan 22 2018 *)
    A = {1}; Do[A = Append[A, 2^IntegerExponent[n, 2]*(2*n - 1)*A[[n]]/n], {n, 1, 25}]; Print[A] (* John Lawrence, Jul 17 2020 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( pollegendre(n), n) * 2^valuation((n\2*2)!, 2))};
    
  • PARI
    a(n)=binomial(2*n,n)>>hammingweight(n); \\ Gleb Koloskov, Sep 26 2021
    
  • Sage
    # uses[A000120]
    @CachedFunction
    def swing(n):
        if n == 0: return 1
        return swing(n-1)*n if is_odd(n) else 4*swing(n-1)/n
    A001790 = lambda n: swing(2*n)/2^A000120(2*n)
    [A001790(n) for n in (0..25)]  # Peter Luschny, Nov 19 2012
    

Formula

a(n) = numerator( binomial(2*n,n)/4^n ) (cf. A046161).
a(n) = A000984(n)/A001316(n) where A001316(n) is the highest power of 2 dividing C(2*n, n) = A000984(n). - Benoit Cloitre, Jan 27 2002
a(n) = denominator of (2^n/binomial(2*n,n)). - Artur Jasinski, Nov 26 2011
a(n) = numerator(L(n)), with rational L(n):=binomial(2*n,n)/2^n. L(n) is the leading coefficient of the Legendre polynomial P_n(x).
L(n) = (2*n-1)!!/n! with the double factorials (2*n-1)!! = A001147(n), n >= 0.
Numerator in (1-2t)^(-1/2) = 1 + t + (3/2)t^2 + (5/2)t^3 + (35/8)t^4 + (63/8)t^5 + (231/16)t^6 + (429/16)t^7 + ... = 1 + t + 3*t^2/2! + 15*t^3/3! + 105*t^4/4! + 945*t^5/5! + ... = e.g.f. for double factorials A001147 (cf. A094638). - Tom Copeland, Dec 04 2013
From Ralf Steiner, Apr 08 2017: (Start)
a(n)/A061549(n) = (-1/4)^n*sqrt(Pi)/(Gamma(1/2 - n)*Gamma(1 + n)).
Sum_{k>=0} a(k)/A061549(k) = 2/sqrt(3).
Sum_{k>=0} (-1)^k*a(k)/A061549(k) = 2/sqrt(5).
Sum_{k>=0} (-1)^(k+1)*a(k)/A061549(k) = -2/sqrt(5).
a(n)/A123854(n) = (-1/2)^n*sqrt(Pi)/(gamma(1/2 - n)*gamma(1 + n)).
Sum_{k>=0} a(k)/A123854(k) = sqrt(2).
Sum_{k>=0} (-1)^k*a(k)/A123854(k) = sqrt(2/3).
Sum_{k>=0} (-1)^(k+1)*a(k)/A123854(k) = -sqrt(2/3). (End)
a(n) = 2^A007814(n)*(2*n-1)*a(n-1)/n. - John Lawrence, Jul 17 2020
Sum_{k>=0} A086117(k+3)/a(k+2) = Pi. - Antonio Graciá Llorente, Aug 31 2024
a(n) = A001803(n)/(2*n+1). - G. C. Greubel, Sep 23 2024

A001803 Numerators in expansion of (1 - x)^(-3/2).

Original entry on oeis.org

1, 3, 15, 35, 315, 693, 3003, 6435, 109395, 230945, 969969, 2028117, 16900975, 35102025, 145422675, 300540195, 9917826435, 20419054425, 83945001525, 172308161025, 1412926920405, 2893136075115, 11835556670925
Offset: 0

Views

Author

Keywords

Comments

a(n) is the denominator of the integral from 0 to Pi of (sin(x))^(2*n+1). - James R. Buddenhagen, Aug 17 2008
a(n) is the denominator of (2n)!!/(2*n + 1)!! = 2^(2*n)*n!*n!/(2*n + 1)! (see Andersson). - N. J. A. Sloane, Jun 27 2011
a(n) = (2*n + 1)*A001790(n). A046161(n)/a(n) = 1, 2/3, 8/15, 16/35, 128/315, 256/693, ... is binomial transform of Madhava-Gregory-Leibniz series for Pi/4 (i.e., 1 - 1/3 + 1/5 - 1/7 + ... ). See A173384 and A173396. - Paul Curtz, Feb 21 2010
a(n) is the denominator of Integral_{x=-oo..oo} sech(x)^(2*n+2) dx. The corresponding numerator is A101926(n). - Mohammed Yaseen, Jul 25 2023

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 798.
  • G. Prévost, Tables de Fonctions Sphériques. Gauthier-Villars, Paris, 1933, pp. 156-157.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 6, equation 6:14:9 at page 51.

Crossrefs

The denominator is given in A046161.
Largest odd divisors of A001800, A002011, A002457, A005430, A033876, A086228.
Bisection of A004731, A004735, A086116.
Second column of triangle A100258.
Cf. A002596 (numerators in expansion of (1-x)^(1/2)).
Cf. A161198 (triangle related to the series expansions of (1-x)^((-1-2*n)/2)).
A163590 is the odd part of the swinging factorial, A001790 at even indices. - Peter Luschny, Aug 01 2009

Programs

  • Julia
    A001803(n) = sum(<<(A001790(k), A005187(n) - A005187(k)) for k in 0:n) # Peter Luschny, Oct 03 2019
    
  • Magma
    A001803:= func< n | Numerator(Binomial(n+2,2)*Catalan(n+1)/4^n) >;
    [A001803(n): n in [0..30]]; // G. C. Greubel, Apr 27 2025
    
  • Maple
    swing := proc(n) option remember; if n = 0 then 1 elif irem(n, 2) = 1 then swing(n-1)*n else 4*swing(n-1)/n fi end:
    sigma := n -> 2^(add(i,i= convert(iquo(n,2),base,2))):
    a := n -> swing(2*n+1)/sigma(2*n+1); # Peter Luschny, Aug 01 2009
    A001803 := proc(n) (2*n+1)*binomial(2*n,n)/4^n ; numer(%) ; end proc: # R. J. Mathar, Jul 06 2011
    a := n -> denom(Pi*binomial(n, -1/2)): seq(a(n), n = 0..22); # Peter Luschny, Dec 06 2024
  • Mathematica
    Numerator/@CoefficientList[Series[(1-x)^(-3/2),{x,0,25}],x]  (* Harvey P. Dale, Feb 19 2011 *)
    Table[Denominator[Beta[1, n + 1, 1/2]], {n, 0, 22}] (* Gerry Martens, Nov 13 2016 *)
  • PARI
    a(n) = numerator((2*n+1)*binomial(2*n,n)/(4^n)); \\ Altug Alkan, Sep 06 2018
    
  • SageMath
    def A001803(n): return numerator((n+1)*binomial(2*n+2,n+1)/2^(2*n+1))
    print([A001803(n) for n in range(31)]) # G. C. Greubel, Apr 27 2025

Formula

a(n) = (2*n + 1)! /(n!^2*2^A000120(n)) = (n + 1)*binomial(2*n+2,n+1)/2^(A000120(n)+1). - Ralf Stephan, Mar 10 2004
From Johannes W. Meijer, Jun 08 2009: (Start)
a(n) = numerator( (2*n+1)*binomial(2*n,n)/(4^n) ).
(1 - x)^(-3/2) = Sum_{n>=0} ((2*n+1)*binomial(2*n,n)/4^n)*x^n. (End)
Truncations of rational expressions like those given by the numerator or denominator operators are artifacts in integer formulas and have many disadvantages. A pure integer formula follows. Let n$ denote the swinging factorial and sigma(n) = number of '1's in the base-2 representation of floor(n/2). Then a(n) = (2*n+1)$ / sigma(2*n+1) = A056040(2*n+1) / A060632(2*n+2). Simply said: This sequence gives the odd part of the swinging factorial at odd indices. - Peter Luschny, Aug 01 2009
a(n) = denominator(Pi*binomial(n, -1/2)). - Peter Luschny, Dec 06 2024

A161198 Triangle of polynomial coefficients related to the series expansions of (1-x)^((-1-2*n)/2).

Original entry on oeis.org

1, 1, 2, 3, 8, 4, 15, 46, 36, 8, 105, 352, 344, 128, 16, 945, 3378, 3800, 1840, 400, 32, 10395, 39048, 48556, 27840, 8080, 1152, 64, 135135, 528414, 709324, 459032, 160720, 31136, 3136, 128
Offset: 0

Views

Author

Johannes W. Meijer, Jun 08 2009, Jul 22 2011

Keywords

Comments

The series expansion of (1-x)^((-1-2*n)/2) = sum(b(p)*x^p, p=0..infinity) for n = 0, 1, 2, .. can be described with b(p) = (F(p,n)/ (2*n-1)!!)*(binomial(2*p,p)/4^(p)) with F(x,n) = 2^n * product( x+(2*k-1)/2, k=1..n). The roots of the F(x,n) polynomials can be found at p = (1-2*k)/2 with k from 1 to n for n = 0, 1, 2, .. . The coefficients of the F(x,n) polynomials lead to the triangle given above. The triangle row sums lead to A001147.
Quite surprisingly we discovered that sum(b(p)*x^p, p=0..infinity) = (1-x)^(-1-2*n)/2, for n = -1, -2, .. . We assume that if m = n+1 then the value returned for product(f(k), k = m..n) is 1 and if m> n+1 then 1/product(f(k), k=n+1..m-1) is the value returned. Furthermore (1-2*n)!! = (-1)^(n+1)/(2*n-3)!! for n = 1, 2, 3 .. . This leads to b(p) = ((-1-2*n)!!/ G(p,n))*(binomial(2*p,p) /4^(p)) for n = -1, -2, .. . For the G(p,n) polynomials we found that G(p,n) = F(-p,-n). The roots of the G(p,n) polynomials can be found at p=(2*k-1)/2 with k from 1 to (-n) for n = -1, -2, .. . The coefficients of the G(p,n) polynomials lead to a second triangle that stands with its head on top of the first one. It is remarkable that the row sums lead once again to A001147.
These two triangles together look like an hourglass so we propose to call the F(p,n) and the G(p,n) polynomials the hourglass polynomials.
Triangle T(n,k), read by rows, given by (1, 2, 3, 4, 5, 6, 7, 8, 9, ...) DELTA (2, 0, 2, 0, 2, 0, 2, 0, 2, ...) where DELTA is the operator defined in A084938. Philippe Deléham, May 14 2015.

Examples

			From _Gary W. Adamson_, Jul 19 2011: (Start)
The first few rows of matrix M are:
  1, 2,  0,  0, 0, ...
  1, 3,  2,  0, 0, ...
  1, 4,  5,  2, 0, ...
  1, 5,  9,  7, 2, ...
  1, 6, 14, 16, 9, ... (End)
The first few G(p,n) polynomials are:
  G(p,-3) = 15 - 46*p + 36*p^2 - 8*p^3
  G(p,-2) = 3 - 8*p + 4*p^2
  G(p,-1) = 1 - 2*p
The first few F(p,n) polynomials are:
  F(p,0) = 1
  F(p,1) = 1 + 2*p
  F(p,2) = 3 + 8*p + 4*p^2
  F(p,3) = 15 + 46*p + 36*p^2 + 8*p^3
The first few rows of the upper and lower hourglass triangles are:
  [15, -46, 36, -8]
  [3, -8, 4]
  [1, -2]
  [1]
  [1, 2]
  [3, 8, 4]
  [15, 46, 36, 8]
		

Crossrefs

Cf. A001790 [(1-x)^(-1/2)], A001803 [(1-x)^(-3/2)], A161199 [(1-x)^(-5/2)] and A161201 [(1-x)^(-7/2)].
Cf. A002596 [(1-x)^(1/2)], A161200 [(1-x)^(3/2)] and A161202 [(1-x)^(5/2)].
A046161 gives the denominators of the series expansions of all (1-x)^((-1-2*n)/2).
A028338 is a scaled triangle version, A039757 is a scaled signed triangle version and A109692 is a transposed scaled triangle version.
A001147 is the first left hand column and equals the row sums.
A004041 is the second left hand column divided by 2, A028339 is the third left hand column divided by 4, A028340 is the fourth left hand column divided by 8, A028341 is the fifth left hand column divided by 16.
A000012, A000290, A024196, A024197 and A024198 are the first (n-m=0), second (n-m=1), third (n-m=2), fourth (n-m=3) and fifth (n-m=4) right hand columns divided by 2^m.
A074599 * A025549 is not always equals the second left hand column.
Cf. A029635. [Gary W. Adamson, Jul 19 2011]

Programs

  • Maple
    nmax:=7; for n from 0 to nmax do a(n,n):=2^n: a(n,0):=doublefactorial(2*n-1) od: for n from 2 to nmax do for m from 1 to n-1 do a(n,m) := 2*a(n-1,m-1)+(2*n-1)*a(n-1,m) od: od: seq(seq(a(n,k), k=0..n), n=0..nmax);
    nmax:=7: M := Matrix(1..nmax+1,1..nmax+1): A029635 := proc(n,k): binomial(n,k) + binomial(n-1,k-1) end: for i from 1 to nmax do for j from 1 to i+1 do M[i,j] := A029635(i,j-1) od: od: for n from 0 to nmax do B := M^n: for m from 0 to n do a(n,m):= B[1,m+1] od: od: seq(seq(a(n,m), m=0..n), n=0..nmax);
    A161198 := proc(n,k) option remember; if k > n or k < 0 then 0 elif n = 0 and k = 0 then 1 else 2*A161198(n-1, k-1) + (2*n-1)*A161198(n-1, k) fi end:
    seq(print(seq(A161198(n,k), k = 0..n)), n = 0..6);  # Peter Luschny, May 09 2013
  • Mathematica
    nmax = 7; a[n_, 0] := (2*n-1)!!; a[n_, n_] := 2^n; a[n_, m_] := a[n, m] = 2*a[n-1, m-1]+(2*n-1)*a[n-1, m]; Table[a[n, m], {n, 0, nmax}, {m, 0, n}] // Flatten (* Jean-François Alcover, Feb 25 2014, after Maple *)
  • PARI
    for(n=0,9, print(Vec(Ser( 2^n*prod( k=1,n, x+(2*k-1)/2 ),,n+1))))  \\ M. F. Hasler, Jul 23 2011
    
  • Sage
    @CachedFunction
    def A161198(n,k):
        if k > n or k < 0 : return 0
        if n == 0 and k == 0: return 1
        return 2*A161198(n-1,k-1)+(2*n-1)*A161198(n-1,k)
    for n in (0..6): [A161198(n,k) for k in (0..n)]  # Peter Luschny, May 09 2013

Formula

a(n,m) := coeff(2^(n)*product((x+(2*k-1)/2),k=1..n), x, m) for n = 0, 1, .. ; m = 0, 1, .. .
a(n, m) = 2*a(n-1,m-1)+(2*n-1)*a(n-1,m) with a(n, n) = 2^n and a(n, 0) = (2*n-1)!!.
a(n,m) = the (m+1)-th term in the top row of M^n, where M is an infinite square production matrix; M[i,j] = A029635(i,j-1) = binomial(i, j-1) + binomial(i-1, j-2) with A029635 the (1.2)-Pascal triangle, see the examples and second Maple program. [Gary W. Adamson, Jul 19 2011]
T(n,k) = 2^k * A028338(n,k). - Philippe Deléham, May 14 2015

A098597 Numerator of Catalan(n)/2^(2n+1). Also, numerators of (2n-1)!!/(n+1)!. Odd part of the n-th Catalan number.

Original entry on oeis.org

1, 1, 1, 5, 7, 21, 33, 429, 715, 2431, 4199, 29393, 52003, 185725, 334305, 9694845, 17678835, 64822395, 119409675, 883631595, 1641030105, 6116566755, 11435320455, 171529806825, 322476036831, 1215486600363, 2295919134019, 17383387729001, 32968493968795, 125280277081421
Offset: 0

Views

Author

Michael Somos, Sep 15 2004

Keywords

Comments

Also numerators of g.f. c(x/2) = (1-sqrt(1-2x))/x where c(x) = g.f. of A000108. - Paul Barry, Sep 04 2007
Also numerator of x(n)=Sum(x(k)*x(n-k-1):0<=kA086117(n). - Reinhard Zumkeller, Feb 06 2008
Also numerator of (1/Pi)*int(x^n*sqrt((1-x)/x), x=0..1). - Groux Roland, Mar 17 2011
The negative of this sequence appears in the A-sequence of the Riordan triangle A084930 as numerators 4, -2, -seq(a(n-1), n >= 2). The denominators look like 1, seq(A120777(n-1), n >= 1). - Wolfdieter Lang, Aug 04 2014
The series of a(n)/A046161(n+1) is absolutely convergent to 1. - Ralf Steiner, Feb 09 2017

Examples

			1/(1 + sqrt(1-x)) = 1/2 + 1/8*x + 1/16*x^2 + 5/128*x^3 + 7/256*x^4 + ...
		

Crossrefs

Cf. Equals A000265(A000108(n)).
Essentially the absolute values of A002596. Cf. A000108, A001795.

Programs

  • Magma
    [Numerator(Catalan(n)/2^(2*n+1)):n in [0..30]]; // Vincenzo Librandi, Jan 14 2016
  • Maple
    a:= n-> abs(numer(binomial(1/2, n+1))): seq(a(n), n=0..50); # Alois P. Heinz, Apr 10 2009
  • Mathematica
    Table[Numerator[CatalanNumber[n]/2^(2n+1)],{n,0,30}] (* Harvey P. Dale, Jul 27 2011 *)
    A098597[n_] := With[{c = CatalanNumber[n]}, c / 2^IntegerExponent[c, 2]];
    Table[A098597[n], {n, 0, 29}]  (* Peter Luschny, Apr 16 2024 *)
  • PARI
    {a(n) = if( n < 0, 0, numerator(polcoeff(1 / (1 + sqrt(1 - x + x * O(x^n))), n)))};
    

Formula

Numerators of g.f.: 1/(1 + sqrt(1-x)).
a(n) = A000108(n) / 2^A048881(n).

Extensions

Edited by Ralf Stephan, Dec 28 2004

A161200 Numerators in expansion of (1-x)^(3/2).

Original entry on oeis.org

1, -3, 3, 1, 3, 3, 7, 9, 99, 143, 429, 663, 4199, 6783, 22287, 37145, 1002915, 1710855, 5892945, 10235115, 71645805, 126233085, 447553665, 797813055, 11435320455, 20583576819, 74417546961, 135054066707, 983965343151, 1798281489207, 6593698793759, 12123897782073
Offset: 0

Views

Author

Johannes W. Meijer, Jun 08 2009

Keywords

References

  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 6, equation 6:14:5 at page 50.

Crossrefs

Cf. A002596 ((1-x)^(1/2)) and A161202 ((1-x)^(5/2)).
Cf. A161199 (numerators in expansion of (1-x)^(-5/2)).
Cf. A161198 (triangle related to the series expansions of (1-x)^((-1-2*n)/2) for all values of n).
Cf. A046161 (denominators of the series expansions of (1-x)^(3/2)).

Programs

  • Mathematica
    Numerator[CoefficientList[Series[(1-x)^(3/2),{x,0,30}],x]] (* Harvey P. Dale, Aug 26 2016 *)
    a[n_]:= Numerator[3/(3-8*n+4*n^2)*Binomial[2*n,n]/(4^n)]; Array[a,28,0] (* Stefano Spezia, Dec 29 2024 *)

Formula

a(n) = numerator((3/(3-8*n+4*n^2))*binomial(2*n,n)/4^n).

A161199 Numerators in expansion of (1-x)^(-5/2).

Original entry on oeis.org

1, 5, 35, 105, 1155, 3003, 15015, 36465, 692835, 1616615, 7436429, 16900975, 152108775, 339319575, 1502700975, 3305942145, 115707975075, 251835004575, 1091285019825, 2354878200675, 20251952525805, 43397041126725, 185423721177825, 395033145117975
Offset: 0

Views

Author

Johannes W. Meijer, Jun 08 2009

Keywords

Crossrefs

Cf. A161198 (triangle for (1-x)^((-1-2*n)/2) for all values of n).
Cf. A046161 (denominators for (1-x)^(-5/2)).
Numerators of [x^n]( (1-x)^(p/2) ): A161202 (p=5), A161200 (p=3), A002596 (p=1), A001790 (p=-1), A001803 (p=-3), this sequence (p=-5), A161201 (p=-7).

Programs

  • Magma
    A161199:= func< n | Numerator( Binomial(n+3,3)*Catalan(n+2)/2^(2*n+1) ) >;
    [A161199(n): n in [0..30]]; // G. C. Greubel, Sep 24 2024
    
  • Mathematica
    Numerator[CoefficientList[Series[(1-x)^(-5/2),{x,0,30}],x]] (* or *) Numerator[Table[(4n^2+8n+3)/3 Binomial[2n,n]/4^n,{n,0,30}]] (* Harvey P. Dale, Oct 15 2011 *)
  • SageMath
    def A161199(n): return numerator((-1)^n*binomial(-5/2,n))
    [A161199(n) for n in range(31)] # G. C. Greubel, Sep 24 2024

Formula

a(n) = numerator(((3 + 8*n + 4*n^2)/3)*binomial(2*n,n)/(4^n)).
a(n) = denominator((3/2)*Integral_{x=0..1} x^n*sqrt(1-x) dx), where the integral is sqrt(Pi)*n!/Gamma(n+5/2) = n!/( (n+3/2)*(n+1/2)*(n-1/2)*...*(1/2)). - Groux Roland, Feb 23 2011

A161202 Numerators in expansion of (1-x)^(5/2).

Original entry on oeis.org

1, -5, 15, -5, -5, -3, -5, -5, -45, -55, -143, -195, -1105, -1615, -4845, -7429, -185725, -294975, -950475, -1550775, -10235115, -17058525, -57378675, -97294275, -1329688425, -2287064091, -7916760315, -13781027215
Offset: 0

Views

Author

Johannes W. Meijer, Jun 08 2009

Keywords

Crossrefs

Cf. A046161 (denominators).
Cf. A161198 (triangle of coefficients of (1-x)^((-1-2*n)/2)).
Numerators of [x^n]( (1-x)^(p/2) ): this sequence (p=5), A161200 (p=3), A002596 (p=1), A001790 (p=-1), A001803 (p=-3), A161199 (p=-5), A161201 (p=-7).

Programs

  • Magma
    A161202:= func< n | -Numerator(15*(n+1)*Catalan(n)/(4^n*(2*n-1)*(2*n-3)*(2*n-5))) >;
    [A161202(n): n in [0..30]]; // G. C. Greubel, Sep 24 2024
    
  • Mathematica
    Numerator[CoefficientList[Series[(1-x)^(5/2),{x,0,30}],x]] (* Harvey P. Dale, Aug 22 2011 *)
    Table[(-1)^n*Numerator[Binomial[5/2, n]], {n,0,30}] (* G. C. Greubel, Sep 24 2024 *)
  • SageMath
    def A161202(n): return (-1)^n*numerator(binomial(5/2,n))
    [A161202(n) for n in range(31)] # G. C. Greubel, Sep 24 2024

Formula

a(n) = numerator( (15/(15-46*n+36*n^2-8*n^3))*binomial(2*n,n)/(4^n) ).
a(n) = (-1)^n*numerator( binomial(5/2, n) ). - G. C. Greubel, Sep 24 2024

A161201 Numerators in expansion of (1-x)^(-7/2).

Original entry on oeis.org

1, 7, 63, 231, 3003, 9009, 51051, 138567, 2909907, 7436429, 37182145, 91265265, 882230895, 2103781365, 9917826435, 23141595015, 856239015555, 1964313035685, 8948537162565, 20251952525805, 182267572732245
Offset: 0

Views

Author

Johannes W. Meijer, Jun 08 2009

Keywords

Crossrefs

Cf. A046161 (denominators).
Cf. A161198 (triangle of coefficients of (1-x)^((-1-2*n)/2)).
Numerators of [x^n]( (1-x)^(p/2) ): A161202 (p=5), A161200 (p=3), A002596 (p=1), A001790 (p=-1), A001803 (p=-3), A161199 (p=-5), this sequence (p=-7).

Programs

  • Magma
    A161201:= func< n | Numerator((n+1)*(2*n+1)*(2*n+3)*(2*n+5)*Catalan(n)/(15*4^n)) >;
    [A161201(n): n in [0..30]]; // G. C. Greubel, Sep 24 2024
    
  • Mathematica
    CoefficientList[Series[(1-x)^(-7/2),{x,0,20}],x]//Numerator (* Harvey P. Dale, Jan 14 2020 *)
    Table[(-1)^n*Numerator[Binomial[-7/2, n]], {n, 0, 30}] (* G. C. Greubel, Sep 24 2024 *)
  • SageMath
    def A161201(n): return (-1)^n*numerator(binomial(-7/2,n))
    [A161201(n) for n in range(31)] # G. C. Greubel, Sep 24 2024

Formula

a(n) = numerator(((15+46*n+36*n^2+8*n^3)/15)*binomial(2*n,n)/(4^n)).
a(n) = (-1)^n*numerator( binomial(-7/2, n) ). - G. C. Greubel, Sep 24 2024
Showing 1-10 of 20 results. Next