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.

A045622 Convolution of A000108 (Catalan numbers) with A045543.

Original entry on oeis.org

1, 25, 362, 3973, 36646, 299530, 2238676, 15613741, 103054094, 650194974, 3950996556, 23257207714, 133217073276, 745218012084, 4083224828328, 21966983072637, 116268166691358, 606474982072982, 3122157367765788
Offset: 1

Views

Author

Keywords

Comments

Also convolution of A045530 with A000984 (central binomial coefficients); also convolution of A045505 with A000302 (powers of 4).

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 40); Coefficients(R!( (1-Sqrt(1-4*x))/(2*(1-4*x)^6) )); // G. C. Greubel, Jan 13 2020
    
  • Maple
    seq(coeff(series((1-sqrt(1-4*x))/(2*(1-4*x)^6), x, n+1), x, n), n = 0..40); # G. C. Greubel, Jan 13 2020
  • Mathematica
    CoefficientList[Series[(1-Sqrt[1-4*x])/(2*x*(1-4*x)^6), {n,0,40}], x] (* G. C. Greubel, Jan 13 2020 *)
  • PARI
    my(x='x+O('x^40)); Vec((1-sqrt(1-4*x))/(2*(1-4*x)^6)) \\ G. C. Greubel, Jan 13 2020
    
  • Sage
    def A045622_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( (1-sqrt(1-4*x))/(2*(1-4*x)^6) ).list()
    A045622_list(40) # G. C. Greubel, Jan 13 2020

Formula

a(n) = binomial(n+6, 5)*(4^(n+1) - A000984(n+6)/A000984(5))/2, A000984(n) = binomial(2*n, n).
G.f.: x*c(x)/(1-4*x)^6, where c(x) = g.f. for Catalan numbers.

A054335 A convolution triangle of numbers based on A000984 (central binomial coefficients of even order).

Original entry on oeis.org

1, 2, 1, 6, 4, 1, 20, 16, 6, 1, 70, 64, 30, 8, 1, 252, 256, 140, 48, 10, 1, 924, 1024, 630, 256, 70, 12, 1, 3432, 4096, 2772, 1280, 420, 96, 14, 1, 12870, 16384, 12012, 6144, 2310, 640, 126, 16, 1, 48620, 65536, 51480, 28672, 12012, 3840, 924, 160, 18, 1
Offset: 0

Views

Author

Wolfdieter Lang, Mar 13 2000

Keywords

Comments

In the language of the Shapiro et al. reference (given in A053121) such a lower triangular (ordinary) convolution array, considered as a matrix, belongs to the Bell-subgroup of the Riordan-group. The g.f. for the row polynomials p(n,x) (increasing powers of x) is 1/(sqrt(1-4*z)-x*z).
Riordan array (1/sqrt(1-4*x),x/sqrt(1-4*x)). - Paul Barry, May 06 2009
The matrix inverse is apparently given by deleting the leftmost column from A206022. - R. J. Mathar, Mar 12 2013

Examples

			Triangle begins:
    1;
    2,    1;
    6,    4,   1;
   20,   16,   6,   1;
   70,   64,  30,   8,  1;
  252,  256, 140,  48, 10,  1;
  924, 1024, 630, 256, 70, 12, 1; ...
Fourth row polynomial (n=3): p(3,x) = 20 + 16*x + 6*x^2 + x^3.
From _Paul Barry_, May 06 2009: (Start)
Production matrix begins
    2,   1;
    2,   2,  1;
    0,   2,  2,  1;
   -2,   0,  2,  2,  1;
    0,  -2,  0,  2,  2,  1;
    4,   0, -2,  0,  2,  2, 1;
    0,   4,  0, -2,  0,  2, 2, 1;
  -10,   0,  4,  0, -2,  0, 2, 2, 1;
    0, -10,  0,  4,  0, -2, 0, 2, 2, 1; (End)
		

Crossrefs

Row sums: A026671.

Programs

  • GAP
    T:= function(n, k)
        if k mod 2=0 then return Binomial(2*n-k, n-Int(k/2))*Binomial(n-Int(k/2),Int(k/2))/Binomial(k,Int(k/2));
        else return 4^(n-k)*Binomial(n-Int((k-1)/2)-1, Int((k-1)/2));
        fi;
      end;
    Flat(List([0..10], n-> List([0..n], k-> T(n, k) ))); # G. C. Greubel, Jul 20 2019
  • Magma
    T:= func< n, k | (k mod 2) eq 0 select Binomial(2*n-k, n-Floor(k/2))* Binomial(n-Floor(k/2),Floor(k/2))/Binomial(k,Floor(k/2)) else 4^(n-k)*Binomial(n-Floor((k-1)/2)-1, Floor((k-1)/2)) >;
    [[T(n,k): k in [0..n]]: n in [0..10]]; // G. C. Greubel, Jul 20 2019
    
  • Maple
    A054335 := proc(n,k)
        if k <0 or k > n then
            0 ;
        elif type(k,odd) then
            kprime := floor(k/2) ;
            binomial(n-kprime-1,kprime)*4^(n-k) ;
        else
            kprime := k/2 ;
            binomial(2*n-k,n-kprime)*binomial(n-kprime,kprime)/binomial(k,kprime) ;
        end if;
    end proc: # R. J. Mathar, Mar 12 2013
    # Uses function PMatrix from A357368. Adds column 1,0,0,0,... to the left.
    PMatrix(10, n -> binomial(2*(n-1), n-1)); # Peter Luschny, Oct 19 2022
  • Mathematica
    Flatten[ CoefficientList[#1, x] & /@ CoefficientList[ Series[1/(Sqrt[1 - 4*z] - x*z), {z, 0, 9}], z]] (* or *)
    a[n_, k_?OddQ] := 4^(n-k)*Binomial[(2*n-k-1)/2, (k-1)/2]; a[n_, k_?EvenQ] := (Binomial[n-k/2, k/2]*Binomial[2*n-k, n-k/2])/Binomial[k, k/2]; Table[a[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Sep 08 2011, updated Jan 16 2014 *)
  • PARI
    T(n, k) = if(k%2==0, binomial(2*n-k, n-k/2)*binomial(n-k/2,k/2)/binomial(k,k/2), 4^(n-k)*binomial(n-(k-1)/2-1, (k-1)/2));
    for(n=0,10, for(k=0,n, print1(T(n,k), ", "))) \\ G. C. Greubel, Jul 20 2019
    
  • Sage
    def T(n, k):
        if (mod(k,2)==0): return binomial(2*n-k, n-k/2)*binomial(n-k/2,k/2)/binomial(k,k/2)
        else: return 4^(n-k)*binomial(n-(k-1)/2-1, (k-1)/2)
    [[T(n,k) for k in (0..n)] for n in (0..10)] # G. C. Greubel, Jul 20 2019
    

Formula

a(n, 2*k+1) = binomial(n-k-1, k)*4^(n-2*k-1), a(n, 2*k) = binomial(2*(n-k), n-k)*binomial(n-k, k)/binomial(2*k, k), k >= 0, n >= m >= 0; a(n, m) := 0 if n
Column recursion: a(n, m)=2*(2*n-m-1)*a(n-1, m)/(n-m), n>m >= 0, a(m, m) := 1.
G.f. for column m: cbie(x)*(x*cbie(x))^m, with cbie(x) := 1/sqrt(1-4*x).
G.f.: 1/(1-x*y-2*x/(1-x/(1-x/(1-x/(1-x/(1-... (continued fraction). - Paul Barry, May 06 2009
Sum_{k>=0} T(n,2*k)*(-1)^k*A000108(k) = A000108(n+1). - Philippe Deléham, Jan 30 2012
Sum_{k=0..floor(n/2)} T(n-k,n-2*k) = A098615(n). - Philippe Deléham, Feb 01 2012
T(n,k) = 4*T(n-1,k) + T(n-2,k-2) for k>=1. - Philippe Deléham, Feb 02 2012
Vertical recurrence: T(n,k) = 1*T(n-1,k-1) + 2*T(n-2,k-1) + 6*T(n-3,k-1) + 20*T(n-4,k-1) + ... for k >= 1 (the coefficients 1, 2, 6, 20, ... are the central binomial coefficients A000984). - Peter Bala, Oct 17 2015

A050982 5-idempotent numbers.

Original entry on oeis.org

1, 30, 525, 7000, 78750, 787500, 7218750, 61875000, 502734375, 3910156250, 29326171875, 213281250000, 1510742187500, 10458984375000, 70971679687500, 473144531250000, 3105010986328125, 20091247558593750, 128360748291015625, 810699462890625000
Offset: 5

Keywords

Comments

Number of n-permutations of 6 objects: t,u,v,z,x, y with repetition allowed, containing exactly five u's. Example: a(6)=30 because we have uuuuut, uuuutu, uuutuu, uutuuu, utuuuu, tuuuuu, uuuuuv, uuuuvu, uuuvuu, uuvuuu, uvuuuu, vuuuuu, uuuuuz, uuuuzu, uuuzuu, uuzuuu, uzuuuu, zuuuuu, uuuuux, uuuuxu, uuuxuu, uuxuuu, uxuuuu, xuuuuu, uuuuuy, uuuuyu, uuuyuu, uuyuuu, uyuuuu, yuuuuu. - Zerinvary Lajos, Jun 16 2008

References

  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 91, #43.

Programs

Formula

a(n) = C(n, 5)*5^(n-5).
G.f.: x^5/(1-5*x)^6. - Zerinvary Lajos, Aug 06 2008
From Amiram Eldar, Apr 17 2022: (Start)
Sum_{n>=5} 1/a(n) = 6400*log(5/4) - 17125/12.
Sum_{n>=5} (-1)^(n+1)/a(n) = 32400*log(6/5) - 23625/4. (End)

A140404 a(n) = binomial(n+5, 5)*7^n.

Original entry on oeis.org

1, 42, 1029, 19208, 302526, 4235364, 54353838, 652246056, 7419298887, 80787921214, 848273172747, 8636963213424, 85649885199788, 830145041167176, 7886377891088172, 73606193650156272, 676256904160810749, 6126091955339109138, 54794489156088698401, 484498640959100070072
Offset: 0

Author

Zerinvary Lajos, Jun 16 2008

Keywords

Comments

With a different offset, number of n-permutations of 8 objects:r,s,t,u,v,z,x,y with repetition allowed, containing exactly five (5) u's. Example: a(1)=42 because we have
uuuuur, uuuuru, uuuruu, uuruuu, uruuuu, ruuuuu
uuuuus, uuuusu, uuusuu, uusuuu, usuuuu, suuuuu,
uuuuut, uuuutu, uuutuu, uutuuu, utuuuu, tuuuuu,
uuuuuv, uuuuvu, uuuvuu, uuvuuu, uvuuuu, vuuuuu,
uuuuuz, uuuuzu, uuuzuu, uuzuuu, uzuuuu, zuuuuu,
uuuuux, uuuuxu, uuuxuu, uuxuuu, uxuuuu, xuuuuu,
uuuuuy, uuuuyu, uuuyuu, uuyuuu, uyuuuu, yuuuuu.

Programs

  • Magma
    [7^n* Binomial(n+5, 5): n in [0..20]]; // Vincenzo Librandi, Oct 12 2011
    
  • Maple
    seq(binomial(n+5,5)*7^n,n=0..17);
  • Mathematica
    Table[Binomial[n+5,5]7^n,{n,0,20}] (* or *) LinearRecurrence[ {42,-735,6860,-36015,100842,-117649},{1,42,1029,19208,302526,4235364},21] (* Harvey P. Dale, Sep 08 2011 *)
  • PARI
    a(n)=binomial(n+5,5)*7^n \\ Charles R Greathouse IV, Oct 07 2015

Formula

G.f.: 1/(1-7*x)^6. - Zerinvary Lajos, Aug 06 2008
a(n) = 42*a(n-1) - 735*a(n-2) + 6860*a(n-3) - 36015*a(n-4) + 100842*a(n-5) - 117649*a(n-6). - Harvey P. Dale, Sep 08 2011
From Amiram Eldar, Aug 28 2022: (Start)
Sum_{n>=0} 1/a(n) = 45360*log(7/6) - 27965/4.
Sum_{n>=0} (-1)^n/a(n) = 143360*log(8/7) - 229705/12. (End)

A172978 a(n) = binomial(n+10, 10)*4^n.

Original entry on oeis.org

1, 44, 1056, 18304, 256256, 3075072, 32800768, 318636032, 2867724288, 24216338432, 193730707456, 1479398129664, 10848919617536, 76776969601024, 526470648692736, 3509804324618240, 22813728110018560, 144934272698941440, 901813252348968960, 5505807224867389440
Offset: 0

Author

Zerinvary Lajos, Feb 06 2010

Keywords

Programs

  • Magma
    [Binomial(n+10, 10)*4^n: n in [0..30]]; // Vincenzo Librandi, Jun 06 2011
  • Mathematica
    Table[Binomial[n + 10, 10]*4^n, {n, 0, 20}]

Formula

From Amiram Eldar, Mar 27 2022: (Start)
G.f.: 1/(1 - 4*x)^11.
Sum_{n>=0} 1/a(n) = 14269429/63 - 787320*log(4/3).
Sum_{n>=0} (-1)^n/a(n) = 78125000*log(5/4) - 1098284605/63. (End)
Showing 1-5 of 5 results.