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.

Previous Showing 11-20 of 52 results. Next

A002422 Expansion of (1-4*x)^(5/2).

Original entry on oeis.org

1, -10, 30, -20, -10, -12, -20, -40, -90, -220, -572, -1560, -4420, -12920, -38760, -118864, -371450, -1179900, -3801900, -12406200, -40940460, -136468200, -459029400, -1556708400, -5318753700, -18296512728, -63334082520
Offset: 0

Views

Author

Keywords

References

  • A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 55.
  • 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).
  • T. N. Thiele, Interpolationsrechnung. Teubner, Leipzig, 1909, p. 164.

Crossrefs

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( (1-4*x)^(5/2) )); // G. C. Greubel, Jul 03 2019
    
  • Maple
    A002422 := n -> -(15/8)*4^n*GAMMA(n-5/2)/(sqrt(Pi)*GAMMA(1+n)):
    seq(A002422(n), n=0..26); # Peter Luschny, Dec 14 2015
  • Mathematica
    CoefficientList[Series[(1-4x)^{5/2},{x,0,30}],x] (* Vincenzo Librandi, Jun 11 2012 *)
  • PARI
    vector(30, n, n--; (-4)^n*binomial(5/2, n)) \\ G. C. Greubel, Jul 03 2019
    
  • Sage
    [(-4)^n*binomial(5/2, n) for n in (0..30)] # G. C. Greubel, Jul 03 2019

Formula

a(n+3) = -2 * A007272(n).
a(n) = Sum_{m=0..n} binomial(n, m) * K_m(6), where K_m(x) = K_m(n, 2, x) is a Krawtchouk polynomial. - Alexander Barg (abarg(AT)research.bell-labs.com).
a(n) ~ -15/8*Pi^(-1/2)*n^(-7/2)*2^(2*n)*{1 + 35/8*n^-1 + ...}. - Joe Keane (jgk(AT)jgk.org), Nov 22 2001
a(n) = -(15/8)*4^n*Gamma(n-5/2)/(sqrt(Pi)*Gamma(1+n)). - Peter Luschny, Dec 14 2015
a(n) = (-4)^n*binomial(5/2, n). - Peter Luschny, Oct 22 2018
D-finite with recurrence: n*a(n) +2*(-2*n+7)*a(n-1)=0. - R. J. Mathar, Jan 16 2020
From Amiram Eldar, Mar 24 2022: (Start)
Sum_{n>=0} 1/a(n) = 32/45 - 14*Pi/(3^5*sqrt(3)).
Sum_{n>=0} (-1)^n/a(n) = 2144/1875 - 28*log(phi)/(5^4*sqrt(5)), where phi is the golden ratio (A001622). (End)

A002423 Expansion of (1-4*x)^(7/2).

Original entry on oeis.org

1, -14, 70, -140, 70, 28, 28, 40, 70, 140, 308, 728, 1820, 4760, 12920, 36176, 104006, 305900, 917700, 2801400, 8684340, 27293640, 86843400, 279409200, 908079900, 2978502072, 9851968392, 32839894640
Offset: 0

Views

Author

Keywords

References

  • A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 55.
  • 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).
  • T. N. Thiele, Interpolationsrechnung. Teubner, Leipzig, 1909, p. 164.

Crossrefs

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( (1-4*x)^(7/2) )); // G. C. Greubel, Jul 03 2019
    
  • Maple
    A002423 := n -> (105/16)*4^n*GAMMA(-7/2+n)/(sqrt(Pi)*GAMMA(1+n)):
    seq(A002423(n), n=0..27); # Peter Luschny, Dec 14 2015
  • Mathematica
    CoefficientList[Series[(1-4*x)^(7/2),{x,0,30}],x] (* Jean-François Alcover, Mar 21 2011 *)
    Table[(4^(-1+x) Pochhammer[-(7/2),-1+x])/Pochhammer[1,-1+x],{x,30}] (* Harvey P. Dale, Jul 13 2011 *)
  • PARI
    vector(30, n, n--; (-4)^n*binomial(7/2, n)) \\ G. C. Greubel, Jul 03 2019
    
  • Sage
    [(-4)^n*binomial(7/2, n) for n in (0..30)] # G. C. Greubel, Jul 03 2019

Formula

a(n) = Sum_{m=0..n} binomial(n, m) * K_m(8), where K_m(x) = K_m(n, 2, x) is a Krawtchouk polynomial. - Alexander Barg (abarg(AT)research.bell-labs.com)
a(n) ~ 105*4^(n-2)/(sqrt(Pi)*n^(9/2)). - Vaclav Kotesovec, Jul 28 2013
a(n) = (105/16)*4^n*Gamma(-7/2+n)/(sqrt(Pi)*Gamma(1+n)). - Peter Luschny, Dec 14 2015
a(n) = (-4)^n * binomial(7/2, n). - G. C. Greubel, Jul 03 2019
D-finite with recurrence: n*a(n) +2*(-2*n+9)*a(n-1)=0. - R. J. Mathar, Jan 16 2020
From Amiram Eldar, Mar 24 2022: (Start)
Sum_{n>=0} 1/a(n) = 36/35 + 2*Pi/(3^4*sqrt(3)).
Sum_{n>=0} (-1)^n/a(n) = 23932/21875 - 36*log(phi)/(5^5*sqrt(5)), where phi is the golden ratio (A001622). (End)

A115140 O.g.f. inverse of Catalan A000108 o.g.f.

Original entry on oeis.org

1, -1, -1, -2, -5, -14, -42, -132, -429, -1430, -4862, -16796, -58786, -208012, -742900, -2674440, -9694845, -35357670, -129644790, -477638700, -1767263190, -6564120420, -24466267020, -91482563640, -343059613650, -1289904147324, -4861946401452, -18367353072152
Offset: 0

Views

Author

Wolfdieter Lang, Jan 13 2006

Keywords

Crossrefs

See A034807 and A115149 for comments.
For convolutions of this sequence see A115141-A115149.

Formula

O.g.f.: 1/c(x) = 1-x*c(x) with the o.g.f. c(x):=(1-sqrt(1-4*x))/(2*x) of A000108 (Catalan numbers).
a(0) = 1, a(n) = -C(n-1), n>=1, with C(n):=A000108(n) (Catalan).
G.f.: (1 + sqrt(1-4*x))/2=U(0) where U(k)=1 - x/U(k+1) ; (continued fraction, 1-step). - Sergei N. Gladkovskii, Oct 29 2012
G.f.: 1/G(0) where G(k) = 1 - x/(x - 1/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Dec 12 2012
G.f.: G(0), where G(k)= 2*x*(2*k+1) + k + 1 - 2*x*(k+1)*(2*k+3)/G(k+1) ; (continued fraction). - Sergei N. Gladkovskii, Jul 14 2013
D-finite with recurrence n*a(n) +2*(-2*n+3)*a(n-1)=0. a(n) = A002420(n)/2, n>0. - R. J. Mathar, Aug 09 2015
a(n) ~ -2^(2*n-2) / (sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, May 06 2021

A068555 Triangle read by rows in which row n contains (2i)!*(2j)!/(i!*j!*(i+j)!) for i + j = n, i = 0..n.

Original entry on oeis.org

1, 2, 2, 6, 2, 6, 20, 4, 4, 20, 70, 10, 6, 10, 70, 252, 28, 12, 12, 28, 252, 924, 84, 28, 20, 28, 84, 924, 3432, 264, 72, 40, 40, 72, 264, 3432, 12870, 858, 198, 90, 70, 90, 198, 858, 12870, 48620, 2860, 572, 220, 140, 140, 220, 572, 2860, 48620, 184756, 9724
Offset: 0

Views

Author

N. J. A. Sloane, Mar 23 2002

Keywords

Comments

One of three infinite families of integral factorial ratio sequences of height 1 (see Bober, Theorem 1.2). The other two are A007318 and A046521. A related table is A182073. - Peter Bala, Apr 10 2012

Examples

			From _Bruno Berselli_, Apr 27 2012: (Start)
Triangle begins:
       1;
       2,    2;
       6,    2,    6;
      20,    4,    4,  20;
      70,   10,    6,  10,  70;
     252,   28,   12,  12,  28, 252;
     924,   84,   28,  20,  28,  84, 924;
    3432,  264,   72,  40,  40,  72, 264, 3432;
   12870,  858,  198,  90,  70,  90, 198,  858, 12870;
   48620, 2860,  572, 220, 140, 140, 220,  572,  2860, 48620;
  184756, 9724, 1716, 572, 308, 252, 308,  572,  1716,  9724, 184756; ...
(End)
T(4,0) = A000984(4) = 70, T(4,1) = 4*20 - 70 = 10, T(4,2) = 4*4 - 10 = 6, T(4,3) = 4*4 - 6 = 10, T(4,4) = 4*20 - 10 = 70. - _Philippe Deléham_, Mar 10 2014
		

References

  • R. K. Guy and Cal Long, Email to N. J. A. Sloane, Feb 22, 2002.
  • Peter J. Larcombe and David R. French, On the integrality of the Catalan-Larcombe-French sequence 1,8,80,896,10816,.... Proceedings of the Thirty-second Southeastern International Conference on Combinatorics, Graph Theory and Computing (Baton Rouge, LA, 2001). Congr. Numer. 148 (2001), 65-91. MR1887375
  • Umberto Scarpis, Sui numeri primi e sui problemi dell'analisi indeterminata in Questioni riguardanti le matematiche elementari, Nicola Zanichelli Editore (1924-1927, third edition), page 11.

Crossrefs

Apart perhaps from signs, diagonals give A000984, A002420, A078718.
Cf. A182411, A082590 (row sums).

Programs

  • Magma
    [Factorial(2*i)*Factorial(2*(n-i))/(Factorial(i)*Factorial(n)*Factorial(n-i)): i in [0..n], n in [0..10]]; // Bruno Berselli, Apr 27 2012
  • Maple
    A068555 := proc(n,i)
        j := n-i ;
        (2*i)!*(2*j)!/(i!*j!*(i+j)!) ;
    end proc: # R. J. Mathar, May 31 2016
  • Mathematica
    Flatten[ Table[ Table[ (2i)!*(2(n - i))!/(i!*(n - i)!*n!), {i, 0, n}], {n, 0, 9}]]
  • PARI
    a(n,k)=if(n<0 || k<0,0,(2*n)!*(2*k)!/n!/k!/(n+k)!);
    

Formula

The square array defined by f := (a, b)->add(binomial(2*a, k)*binomial(2*b, a+b-k)*(-1)^(a+b-k), k=0..2*a); and read by antidiagonals gives a signed version. See Sprugnoli, 3.38.
Let f(x) = 1/sqrt(1 - 4*x) denote the o.g.f for A000984. The o.g.f. for this table is (f(x) + f(y))*f(x)*f(y)*(1/(1 + f(x)*f(y))) = (1 + 2*x + 6*x^2 + 20*x^3 + ...) + (2 + 2*x + 4*x^2 + 10*x^3 + ...)*y + (6 + 4*x + 6*x^2 + 12*x^3 + ...)*y^2 + .... - Peter Bala, Apr 10 2012
T(n,0) = A000984(n), T(n,k) = 4*T(n-1,k-1) - T(n,k-1) for k = 1..n. - Philippe Deléham, Mar 10 2014

A241477 Triangle read by rows, number of orbitals classified with respect to the first zero crossing, n>=1, 1<=k<=n.

Original entry on oeis.org

1, 0, 2, 2, 2, 2, 0, 4, 0, 2, 6, 12, 4, 2, 6, 0, 12, 0, 4, 0, 4, 20, 60, 12, 12, 12, 4, 20, 0, 40, 0, 12, 0, 8, 0, 10, 70, 280, 40, 60, 36, 24, 40, 10, 70, 0, 140, 0, 40, 0, 24, 0, 20, 0, 28, 252, 1260, 140, 280, 120, 120, 120, 60, 140, 28, 252, 0, 504, 0
Offset: 1

Views

Author

Peter Luschny, Apr 23 2014

Keywords

Comments

For the combinatorial definitions see A232500. An orbital w over n sectors has its first zero crossing at k if k is the smallest j such that the partial sum(1<=i<=j, w(i))) = 0, where w(i) are the jumps of the orbital represented by -1, 0, 1.

Examples

			[1], [ 1]
[2], [ 0,  2]
[3], [ 2,  2,  2]
[4], [ 0,  4,  0,  2]
[5], [ 6, 12,  4,  2,  6]
[6], [ 0, 12,  0,  4,  0, 4]
[7], [20, 60, 12, 12, 12, 4, 20]
		

Crossrefs

Row sums: A056040.
Cf. A232500.

Programs

  • Maple
    A241477 := proc(n, k)
      if n = 0 then 1
    elif k = 0 then 0
    elif irem(n, 2) = 0 and irem(k, 2) = 1 then 0
    elif k = 1 then (n-1)!/iquo(n-1,2)!^2
    else 2*(n-k)!*(k-2)!/iquo(k,2)/(iquo(k-2,2)!*iquo(n-k,2)!)^2
      fi end:
    for n from 1 to 9 do seq(A241477(n, k), k=1..n) od;
  • Mathematica
    T[n_, k_] := Which[n == 0, 1, k == 0, 0, Mod[n, 2] == 0 && Mod[k, 2] == 1,  0, k == 1, (n-1)!/Quotient[n-1, 2]!^2, True, 2*(n-k)!*(k-2)!/Quotient[k, 2]/(Quotient[k-2, 2]!*Quotient[n-k, 2]!)^2];
    Table[T[n, k], {n, 1, 12}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 20 2018, from Maple *)
  • Sage
    def A241477_row(n):
        if n == 0: return [1]
        Z = [0]*n; T = [0] if is_odd(n) else []
        for i in (1..n//2): T.append(-1); T.append(1)
        for p in Permutations(T):
            i = 0; s = p[0]
            while s != 0: i += 1; s += p[i];
            Z[i] += 1
        return Z
    for n in (1..9): A241477_row(n)

Formula

If n is even and k is odd then T(n, k) = 0 else if k = 1 then T(n, 1) = A056040(n-1) else T(n, k) = 2*A057977(k-2)*A056040(n-k).
T(n, n) = A241543(n).
T(n+1, 1) = A126869(n).
T(2*n, 2*n) = |A002420(n)|.
T(2*n+1, 1) = A000984(n).
T(2*n+1, n+1) = A241530(n).
T(2*n+2, 2) = A028329(n).
T(4*n, 2*n) = |A010370(n)|.
T(4*n, 4*n) = |A024491(n)|.
T(4*n+1, 1) = A001448(n).
T(4*n+1, 2*n+1) = A002894(n).

A002424 Expansion of (1-4*x)^(9/2).

Original entry on oeis.org

1, -18, 126, -420, 630, -252, -84, -72, -90, -140, -252, -504, -1092, -2520, -6120, -15504, -40698, -110124, -305900, -869400, -2521260, -7443720, -22331160, -67964400, -209556900, -653817528, -2062039896, -6567978928, -21111360840
Offset: 0

Views

Author

Keywords

References

  • A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 55.
  • 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).
  • T. N. Thiele, Interpolationsrechnung. Teubner, Leipzig, 1909, p. 164.

Crossrefs

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( (1-4*x)^(9/2) )); // G. C. Greubel, Jul 03 2019
    
  • Maple
    A002424 := n -> -(945/32)*4^n*GAMMA(-9/2+n)/(sqrt(Pi)*GAMMA(1+n)):
    seq(A002424(n),n=0..28); # Peter Luschny, Dec 14 2015
  • Mathematica
    CoefficientList[Series[(1-4x)^(9/2),{x,0,30}],x] (* Harvey P. Dale, Dec 27 2011 *)
  • PARI
    my(x='x+O('x^30)); Vec((1-4*x)^(9/2)) \\ Altug Alkan, Dec 14 2015
    
  • PARI
    vector(30, n, n--; (-4)^n*binomial(9/2, n)) \\ G. C. Greubel, Jul 03 2019
    
  • Sage
    [(-4)^n*binomial(9/2, n) for n in (0..30)] # G. C. Greubel, Jul 03 2019

Formula

a(n) = Sum_{m=0..n} binomial(n, m) * K_m(10), where K_m(x) = K_m(n, 2, x) is a Krawtchouk polynomial. - Alexander Barg, abarg(AT)research.bell-labs.com.
a(n) = -(945/32)*4^n*Gamma(-9/2+n)/(sqrt(Pi)*Gamma(1+n)). - Peter Luschny, Dec 14 2015
a(n) = (-4)^n*binomial(9/2, n). - G. C. Greubel, Jul 03 2019
D-finite with recurrence: n*a(n) +2*(-2*n+11)*a(n-1)=0. - R. J. Mathar, Jan 16 2020
From Amiram Eldar, Mar 25 2022: (Start)
Sum_{n>=0} 1/a(n) = 32/35 - 22*Pi/(3^7*sqrt(3)).
Sum_{n>=0} (-1)^n/a(n) = 1050752/984375 - 44*log(phi)/(5^6*sqrt(5)), where phi is the golden ratio (A001622). (End)

A106191 Expansion of sqrt(1-4x)/(1-x).

Original entry on oeis.org

1, -1, -3, -7, -17, -45, -129, -393, -1251, -4111, -13835, -47427, -164999, -581023, -2066823, -7415703, -26805393, -97520733, -356810313, -1312087713, -4846614093, -17974854933, -66907388973, -249872516253, -935991743553, -3515800038201, -13239692841105
Offset: 0

Views

Author

Paul Barry, Apr 24 2005

Keywords

Comments

Row sums of number triangle A106190. Partial sums of A002420.
For n >= 1, the absolute values also give the iterates of A122237, starting from 0. (A122237(0), A122237(A122237(0)), A122237(A122237(A122237(0))), ...), this stems from the fact that the sequence gives the positions of terms with binary expansion 1(10){n-1}0 in A014486 (see A080675).

Crossrefs

|a(n)| = A080300(A080675(n)) = A075161(A001348(n)) (for n >= 1) = A075163(A000244(A008578(n-2))) = A014137(n-1)+A014138(n-2) = 2*A014137(n-1)-1, for n >= 2 (because binomial(2n+2, n+1)/(2n+1) = 2*A000108(n)).

Formula

a(n) = Sum_{k=0..n} binomial(2k, k)/(1-2k).
G.f.: (2/(1-x))/G(0), where G(k) = 1 + 1/(1 - 2*x*(2*k+1)/(2*x*(2*k+1) + (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 24 2013
D-finite with recurrence: a(0)=1, a(1)=-1; for n>1, a(n) = (1/n)*((5*n-6)*a(n-1) - (4*n-6)*a(n-2)). - Tani Akinari, Aug 25 2013

Extensions

Barry's formula made more succinct, as well as comments regarding interpretation as absolute values added by Antti Karttunen, Sep 14 2006

A010370 a(n) = binomial(2*n, n)^2 / (1-2*n).

Original entry on oeis.org

1, -4, -12, -80, -700, -7056, -77616, -906048, -11042460, -139053200, -1796567344, -23696871744, -317933029232, -4326899214400, -59605244280000, -829705000377600, -11654762427179100, -165021757273414800, -2353088020380174000, -33764531705178120000
Offset: 0

Views

Author

Joe Keane (jgk(AT)jgk.org)

Keywords

Comments

Expansion of hypergeometric function F(-1/2, 1/2; 1; 16*x).
Expansion of E(m)/(Pi/2) in powers of m/16 = (k/4)^2, where E(m) is the complete elliptic integral of the second kind evaluated at m. - Michael Somos, Mar 04 2003

Examples

			G.f. = 1 - 4*x - 12*x^2 - 80*x^3 - 700*x^4 - 7056*x^5 - 77616*x^6 - ...
		

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. 591.
  • J. M. Borwein and P. B. Borwein, Pi and the AGM, Wiley, 1987, p. 8.

Crossrefs

Programs

  • Maple
    seq(binomial(2*n,n)^2/(1-2*n), n=0..30); # Robert Israel, Jul 10 2017
  • Mathematica
    CoefficientList[Series[EllipticE[16x]2/Pi, {x, 0, 20}], x]
    Table[Binomial[2n,n]^2/(1-2n),{n,0,30}] (* Harvey P. Dale, Mar 07 2013 *)
  • PARI
    {a(n) = binomial(2*n, n)^2 / (1 - 2*n)}; /* Michael Somos, Dec 13 2002 */

Formula

a(n) ~ -1/2*Pi^-1*n^-2*2^(4*n). [corrected by Vaclav Kotesovec, Oct 04 2019]
a(n) = -4 * A000891(n-1), n>0. - Michael Somos, Dec 13 2002
G.f.: F(-1/2, 1/2; 1; 16x) = E(16*x) / (Pi/2). a(n) = binomial(2*n, n)^2 / (1 - 2*n). - Michael Somos, Mar 04 2003
E.g.f.: Sum_{n>=0} a(n) * (x/2)^(2n)/(2n)! = I0^2*(1-2*x^2) +2*x*I0*I1 +2*x^2*I1^2 where I0=BesselI(0, x), I1=BesselI(1, x). - Michael Somos, Jun 22 2005
n^2*a(n) -4*(2*n-1)*(2*n-3)*a(n-1)=0. - R. J. Mathar, Feb 15 2013
0 = a(n)*(+1048576*a(n+2) + 2695168*a(n+3) - 989568*a(n+4) + 65340*a(n+5)) + a(n+1)*(-8192*a(n+2) - 99840*a(n+3) + 52652*a(n+4) - 4236*a(n+5)) + a(n+2)*(-128*a(n+2) + 280*a(n+3) - 484*a(n+4) + 57*a(n+5)) for all n in Z. - Michael Somos, Jan 21 2017
a(n) = A002894(n) - 8 * A000894(n-1). - Michael Somos, Jul 10 2017

Extensions

Additional comments from Michael Somos, Dec 13 2002

A020923 Expansion of (1-4*x)^(11/2).

Original entry on oeis.org

1, -22, 198, -924, 2310, -2772, 924, 264, 198, 220, 308, 504, 924, 1848, 3960, 8976, 21318, 52668, 134596, 354200, 956340, 2641320, 7443720, 21360240, 62300700, 184410072, 553230216, 1680180656, 5160554872, 16015515120, 50181947376, 158639704608, 505664058438
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

  • Maple
    A002423 := n -> (10395/64)*4^n*GAMMA(-11/2+n)/(sqrt(Pi)*GAMMA(1+n)):
    seq(A002423(n),n=0..28); # Peter Luschny, Dec 14 2015
  • Mathematica
    CoefficientList[Series[(1 - 4*x)^(11/2), {x,0,50}], x] (* G. C. Greubel, Feb 15 2017 *)
  • PARI
    my(x='x+O('x^50)); Vec((1-4*x)^(11/2)) \\ G. C. Greubel, Feb 15 2017

Formula

a(n) = (10395/64)*4^n*Gamma(-11/2+n)/(sqrt(Pi)*Gamma(1+n)). - Peter Luschny, Dec 14 2015
D-finite with recurrence: n*a(n) +2*(-2*n+13)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
From Amiram Eldar, Mar 25 2022: (Start)
a(n) = (-4)^n*binomial(11/2, n).
Sum_{n>=0} 1/a(n) = 1124/1155 + 26*Pi/(3^8*sqrt(3)).
Sum_{n>=0} (-1)^n/a(n) = 56972276/54140625 - 52*log(phi)/(5^7*sqrt(5)), where phi is the golden ratio (A001622). (End)

A361397 Number A(n,k) of k-dimensional cubic lattice walks with 2n steps from origin to origin and avoiding early returns to the origin; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 0, 1, 2, 0, 1, 4, 2, 0, 1, 6, 20, 4, 0, 1, 8, 54, 176, 10, 0, 1, 10, 104, 996, 1876, 28, 0, 1, 12, 170, 2944, 22734, 22064, 84, 0, 1, 14, 252, 6500, 108136, 577692, 275568, 264, 0, 1, 16, 350, 12144, 332050, 4525888, 15680628, 3584064, 858, 0
Offset: 0

Views

Author

Alois P. Heinz, Mar 10 2023

Keywords

Comments

Column k is INVERTi transform of k-th row of A287318.

Examples

			Square array A(n,k) begins:
  1,  1,     1,      1,       1,        1,        1, ...
  0,  2,     4,      6,       8,       10,       12, ...
  0,  2,    20,     54,     104,      170,      252, ...
  0,  4,   176,    996,    2944,     6500,    12144, ...
  0, 10,  1876,  22734,  108136,   332050,   796860, ...
  0, 28, 22064, 577692, 4525888, 19784060, 62039088, ...
		

Crossrefs

Columns k=0-5 give: A000007, |A002420|, A054474, A049037, A359801, A361364.
Rows n=0-2 give: A000012, A005843, A139271.
Main diagonal gives A361297.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0 or i=1, 1,
          add(b(n-j, i-1)*binomial(n, j)^2, j=0..n))
        end:
    g:= proc(n, k) option remember; `if` (n<1, -1,
          -add(g(n-i, k)*(2*i)!*b(i, k)/i!^2, i=1..n))
        end:
    A:= (n,k)-> `if`(n=0, 1, `if`(k=0, 0, g(n, k))):
    seq(seq(A(n, d-n), n=0..d), d=0..10);
  • Mathematica
    b[n_, 0] = 0; b[n_, 1] = 1; b[0, k_] = 1;
    b[n_, k_] := b[n, k] = Sum[Binomial[n, i]^2*b[i, k - 1], {i, 0, n}]; (* A287316 *)
    g[n_, k_] := g[n, k] = b[n, k]*Binomial[2 n, n]; (* A287318 *)
    a[n_, k_] := a[n, k] = g[n, k] - Sum[a[i, k]*g[n - i, k], {i, 1, n - 1}];
    TableForm[Table[a[n, k], {k, 0, 10}, {n, 0, 10}]] (* Shel Kaphan, Mar 14 2023 *)

Formula

A(n,1)/2 = A000108(n-1) for n >= 1.
G.f. of column k: 2 - 1/Integral_{t=0..oo} exp(-t)*BesselI(0,2*t*sqrt(x))^k dt. - Shel Kaphan, Mar 19 2023
Previous Showing 11-20 of 52 results. Next