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 21-30 of 56 results. Next

A003947 Expansion of (1+x)/(1-4*x).

Original entry on oeis.org

1, 5, 20, 80, 320, 1280, 5120, 20480, 81920, 327680, 1310720, 5242880, 20971520, 83886080, 335544320, 1342177280, 5368709120, 21474836480, 85899345920, 343597383680, 1374389534720, 5497558138880, 21990232555520, 87960930222080, 351843720888320
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 5.
For n>=1, a(n+1) is equal to the number of functions f:{1,2,...,n+1}->{1,2,3,4,5} such that for fixed, different x_1, x_2,...,x_n in {1,2,...,n+1} and fixed y_1, y_2,...,y_n in {1,2,3,4,5} we have f(x_i)<>y_i, (i=1,2,...,n). - Milan Janjic, May 10 2007
Number of length-n strings of 5 letters with no two adjacent letters identical. The general case (strings of r letters) is the sequence with g.f. (1+x)/(1-(r-1)*x). - Joerg Arndt, Oct 11 2012
Create a rectangular prism with edges of lengths 2^(n-2), 2^(n-1), and 2^(n) starting at n=2; then the surface area = a(n). - J. M. Bergot, Aug 08 2013

Crossrefs

Cf. A003948, A003949. Column 5 in A265583.

Programs

  • GAP
    Concatenation([1], List([1..30], n-> 5*4^(n-1) )); # G. C. Greubel, Aug 10 2019
  • Magma
    [1] cat [5*4^(n-1): n in [1..30]]; // G. C. Greubel, Aug 10 2019
    
  • Maple
    k := 5; if n = 0 then 1 else k*(k-1)^(n-1); fi;
  • Mathematica
    q = 5; Join[{a = 1}, Table[If[n != 0, a = q*a - a, a = q*a], {n, 0, 25}]] (* and *) Join[{1}, 5*4^Range[0, 25]] (* Vladimir Joseph Stephan Orlovsky, Jul 11 2011 *)
    LinearRecurrence[{4},{1,5},30] (* Harvey P. Dale, Apr 19 2015 *)
  • PARI
    a(n)=5*4^n\4 \\ Charles R Greathouse IV, Sep 08 2011
    
  • Sage
    [1]+[5*4^(n-1) for n in (1..30)] # G. C. Greubel, Aug 10 2019
    

Formula

Binomial transform of A060925. Its binomial transform is A003463 (without leading zero). - Paul Barry, May 19 2003
From Paul Barry, May 19 2003: (Start)
a(n) = (5*4^n - 0^n)/4.
G.f.: (1+x)/(1-4*x).
E.g.f.: (5*exp(4*x) - exp(0))/4. (End)
a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 3. - Philippe Deléham, Jul 10 2005
a(n) = A146523(n)*A011782(n). - R. J. Mathar, Jul 08 2009
a(n) = 5*A000302(n-1), n>0.
a(n) = 4*a(n-1), n>1. - Vincenzo Librandi, Dec 31 2010
G.f.: 2+x- 2/G(0), where G(k)= 1 + 1/(1 - x*(5*k-4)/(x*(5*k+1) - 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 04 2013

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A029635 The (1,2)-Pascal triangle (or Lucas triangle) read by rows.

Original entry on oeis.org

2, 1, 2, 1, 3, 2, 1, 4, 5, 2, 1, 5, 9, 7, 2, 1, 6, 14, 16, 9, 2, 1, 7, 20, 30, 25, 11, 2, 1, 8, 27, 50, 55, 36, 13, 2, 1, 9, 35, 77, 105, 91, 49, 15, 2, 1, 10, 44, 112, 182, 196, 140, 64, 17, 2, 1, 11, 54, 156, 294, 378, 336, 204, 81, 19, 2, 1, 12, 65, 210, 450, 672, 714, 540, 285, 100
Offset: 0

Views

Author

Keywords

Comments

This is also called Vieta's array. - N. J. A. Sloane, Nov 22 2017
Dropping the first term and changing the boundary conditions to T(n,1)=n, T(n,n-1)=2 (n>=2), T(n,n)=1 yields the number of nonterminal symbols (which generate strings of length k) in a certain context-free grammar in Chomsky normal form that generates all permutations of n symbols. Summation over k (1<=k<=n) results in A003945. For the number of productions of this grammar: see A090327. Example: 1; 2, 1; 3, 2, 1; 4, 5, 2, 1; 5, 9, 7, 2, 1; 6, 14, 16, 9, 2, 1; In addition to the example of A090327 we have T(3,3)=#{S}=1, T(3,2)=#{D,E}=2 and T(3,1)=#{A,B,C}=3. - Peter R. J. Asveld, Jan 29 2004
Much as the original Pascal triangle gives the Fibonacci numbers as sums of its diagonals, this triangle gives the Lucas numbers (A000032) as sums of its diagonals; see Posamentier & Lehmann (2007). - Alonso del Arte, Apr 09 2012
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 04 2013
It appears that for the infinite set of (1,N) Pascal's triangles, the binomial transform of the n-th row (n>0), followed by zeros, is equal to the n-th partial sum of (1, N, N, N, ...). Example: for the (1,2) Pascal's triangle, the binomial transform of the second row followed by zeros, i.e., of (1, 3, 2, 0, 0, 0, ...), is equal to the second partial sum of (1, 2, 2, 2, ...) = (1, 4, 9, 16, ...). - Gary W. Adamson, Aug 11 2015
Given any (1,N) Pascal triangle, let the binomial transform of the n-th row (n>1) followed by zeros be Q(x). It appears that the binomial transform of the (n-1)-th row prefaced by a zero is Q(n-1). Example: In the (1,2) Pascal triangle the binomial transform of row 3: (1, 4, 5, 2, 0, 0, 0, ...) is A000330 starting with 1: (1, 5, 14, 30, 55, 91, ...). The binomial transform of row 2 prefaced by a zero and followed by zeros, i.e., of (0, 1, 3, 2, 0, 0, 0, ...) is (0, 1, 5, 14, 30, 55, ...). - Gary W. Adamson, Sep 28 2015
It appears that in the array accompanying each (1,N) Pascal triangle (diagonals of the triangle), the binomial transform of (..., 1, N, 0, 0, 0, ...) preceded by (n-1) zeros generates the n-th row of the array (n>0). Then delete the zeros in the result. Example: in the (1,2) Pascal triangle, row 3 (1, 5, 14, 30, ...) is the binomial transform of (0, 0, 1, 2, 0, 0, 0, ...) with the resulting zeros deleted. - Gary W. Adamson, Oct 11 2015
Read as a square array (similar to the Example section Sq(m,j), but with Sq(0,0)=0 and Sq(m,j)=P(m+1,j) otherwise), P(n,k) are the multiplicities of the eigenvalues, lambda_n = n(n+k-1), of the Laplacians on the unit k-hypersphere, given by Teo (and Choi) as P(n,k) = (2n-k+1)(n+k-2)!/(n!(k-1)!). P(n,k) is also the numerator of a Dirichlet series for the Minakashisundarum-Pleijel zeta function for the sphere. Also P(n,k) is the dimension of the space of homogeneous, harmonic polynomials of degree k in n variables (Shubin, p. 169). For relations to Chebyshev polynomials and simple Lie algebras, see A034807. - Tom Copeland, Jan 10 2016
For a relation to a formulation for a universal Lie Weyl algebra for su(1,1), see page 16 of Durov et al. - Tom Copeland, Jan 15 2016

Examples

			Triangle begins:
  [0] [2]
  [1] [1, 2]
  [2] [1, 3,  2]
  [3] [1, 4,  5,  2]
  [4] [1, 5,  9,  7,   2]
  [5] [1, 6, 14, 16,   9,  2]
  [6] [1, 7, 20, 30,  25, 11,  2]
  [7] [1, 8, 27, 50,  55, 36, 13,  2]
  [8] [1, 9, 35, 77, 105, 91, 49, 15, 2]
.
Read as a square, the array begins:
  n\k| 0  1   2    3    4    5
  --------------------------------------
  0 |  2  2   2    2    2    2   A040000
  1 |  1  3   5    7    9   11   A005408
  2 |  1  4   9   16   25   36   A000290
  3 |  1  5  14   30   55   91   A000330
  4 |  1  6  20   50  105  196   A002415
  5 |  1  7  27   77  182  378   A005585
  6 |  1  8  35  112  294  672   A040977
		

References

  • Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8.
  • Alfred S. Posamentier & Ingmar Lehmann, The (Fabulous) Fibonacci Numbers. New York: Prometheus Books (2007): 97 - 105.
  • M. Shubin and S. Andersson, Pseudodifferential Operators and Spectral Theory, Springer Series in Soviet Mathematics, 1987.

Crossrefs

Cf. A003945 (row sums), A007318, A034807, A061896, A029653 (row-reversed), A157000.
Sums along ascending antidiagonals give Lucas numbers, n>0.

Programs

  • Haskell
    a029635 n k = a029635_tabl !! n !! k
    a029635_row n = a029635_tabl !! n
    a029635_tabl = [2] : iterate
       (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [1,2]
    -- Reinhard Zumkeller, Mar 12 2012, Feb 23 2012
    
  • Maple
    T := proc(n, k) option remember;
    if n = k then 2 elif k = 0 then 1 else T(n-1, k-1) + T(n-1, k) fi end:
    for n from 0 to 8 do seq(T(n, k), k = 0..n) od;  # Peter Luschny, Dec 22 2024
  • Mathematica
    t[0, 0] = 2; t[n_, k_] := If[k < 0 || k > n, 0, Binomial[n, k] + Binomial[n-1, k-1]]; Flatten[Table[t[n, k], {n, 0, 11}, {k, 0, n}]] (* Jean-François Alcover, May 03 2011 *)
    (* The next program cogenerates A029635 and A029638. *)
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := u[n - 1, x] + v[n - 1, x]
    v[n_, x_] := x*u[n - 1, x] + x*v[n - 1, x] + 1
    Table[Factor[u[n, x]], {n, 1, z}]
    Table[Factor[v[n, x]], {n, 1, z}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%]   (* A029638  *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]   (* A029635 *)
    (* Clark Kimberling, Feb 20 2012 *)
    Table[Binomial[n,k]+Binomial[n-1,k-1],{n,0,20},{k,0,n}]//Flatten (* Harvey P. Dale, Feb 08 2024 *)
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, (n==0) + binomial(n, k) + binomial(n-1, k-1))}; /* Michael Somos, Jul 15 2003 */
    
  • Sage
    # uses[riordan_array from A256893]
    riordan_array((2-x)/(1-x), x/(1-x), 8) # Peter Luschny, Nov 09 2019

Formula

From Henry Bottomley, Apr 26 2002; (Start)
T(n, k) = T(n-1, k-1) + T(n-1, k).
T(n, k) = C(n, k) + C(n-1, k-1).
T(n, k) = C(n, k)*(n + k)/n.
T(n, k) = A007318(n, k) + A007318(n-1, k-1).
T(n, k) = A061896(n + k, k) but with T(0, 0) = 1 and T(1, 1) = 2.
Row sum is floor(3^2(n-1)) i.e., A003945. (End)
G.f.: 1 + (1 + x*y) / (1 - x - x*y). - Michael Somos, Jul 15 2003
G.f. for n-th row: (x+2*y)*(x+y)^(n-1).
O.g.f. for row n: (1+x)/(1-x)^(n+1). The entries in row n are the nonzero entries in column n of A053120 divided by 2^(n-1). - Peter Bala, Aug 14 2008
T(2n, n) - T(2n, n+1)= Catalan(n)= A000108(n). - Philippe Deléham, Mar 19 2009
With T(0, 0) = 1 : Triangle T(n, k), read by rows, given by [1,0,0,0,0,0,...] DELTA [2,-1,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Oct 10 2011
With T(0, 0) = 1, as in the Example section below, this is known as Vieta's array. The LU factorization of the square array is given by Yang and Leida, equation 20. - Peter Bala, Feb 11 2012
For n > 0: T(n, k) = A097207(n-1, k), 0 <= k < n. - Reinhard Zumkeller, Mar 12 2012
For n > 0: T(n, k) = A029600(n, k) - A007318(n, k), 0 <= k <= n. - Reinhard Zumkeller, Apr 16 2012
Riordan array ((2-x)/(1-x), x/(1-x)). - Philippe Deléham, Mar 15 2013
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(1 + 4*x + 5*x^2/2! + 2*x^3/3!) = 1 + 5*x + 14*x^2/2! + 30*x^3/3! + 55*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ). - Peter Bala, Dec 22 2014
For n>=1: T(n, 0) + T(n, 1) + T(n, 2) = A000217(n+1). T(n, n-2) = (n-1)^2. - Bob Selcoe, Mar 29 2016:

Extensions

More terms from David W. Wilson
a(0) changed to 2 (was 1) by Daniel Forgues, Jul 06 2010

A003948 Expansion of (1+x)/(1-5*x).

Original entry on oeis.org

1, 6, 30, 150, 750, 3750, 18750, 93750, 468750, 2343750, 11718750, 58593750, 292968750, 1464843750, 7324218750, 36621093750, 183105468750, 915527343750, 4577636718750, 22888183593750, 114440917968750, 572204589843750, 2861022949218750, 14305114746093750
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 6.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954, m is 2, 3, 4, 5, 6. - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001
Hamiltonian path in S_4 X P_2n.
For n>=1, a(n+1) is equal to the number of functions f:{1,2,...,n+1}->{1,2,3,4,5,6} such that for fixed, different x_1, x_2,...,x_n in {1,2,...,n+1} and fixed y_1, y_2,...,y_n in {1,2,3,4,5,6} we have f(x_i)<>y_i, (i=1..n). - Milan Janjic, May 10 2007
For n>=1, a(n) equals the numbers of words of length n over the alphabet {0..5} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015 [Corrected by David Nacin, May 30 2017]
a(n) equals the numbers of sequences of length n on {0,...,5} where no two adjacent terms differ by three. - David Nacin, May 30 2017
It appears that these are the only n>1 for which alpha(n)=2n, where alpha(n) is the entry point of n in the Fibonacci sequence, see A001177. - Philippe Schnoebelen, Apr 11 2024

Crossrefs

Programs

  • GAP
    Concatenation([1], List([1..30], n-> 6*5^(n-1) )); # G. C. Greubel, Sep 24 2019
  • Magma
    [1] cat [6*5^(n-1): n in [1..30]]; // G. C. Greubel, Sep 24 2019
    
  • Maple
    k := 6; if n = 0 then 1 else k*(k-1)^(n-1); fi;
  • Mathematica
    q = 6; Join[{a = 1}, Table[If[n != 0, a = q*a - a, a = q*a], {n, 0, 25}]] (* and *) Join[{1}, 6*5^Range[0, 25]] (* Vladimir Joseph Stephan Orlovsky, Jul 11 2011 *)
    Join[{1},NestList[5#&,6,30]] (* Harvey P. Dale, Dec 31 2013 *)
    CoefficientList[Series[(1+x)/(1-5x), {x,0,30}], x] (* Michael De Vlieger, Dec 10 2016 *)
  • PARI
    Vec((1+x)/(1-5*x)+O(x^30)) \\ Charles R Greathouse IV, Nov 20 2012
    
  • Sage
    [1]+[6*5^(n-1) for n in (1..30)] # G. C. Greubel, Sep 24 2019
    

Formula

G.f.: (1+x)/(1-5*x).
a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 4. - Philippe Deléham, Jul 10 2005
The Hankel transform of this sequence is [1,-6,0,0,0,0,0,0,0,0,...]. - Philippe Deléham, Nov 21 2007
a(n) = 6*5^(n-1) for n>0, a(0)=1. - Vincenzo Librandi, Nov 18 2010
G.f.: 2/x - 5 - 8/(x*U(0)) where U(k)= 1 + 2/(3^k - 3^k/(2 + 1 - 12*x*3^k/(6*x*3^k + 1/U(k+1)))) ; (continued fraction, 4-step). - Sergei N. Gladkovskii, Oct 30 2012
E.g.f.: (6*exp(5*x) - 1)/5. - Ilya Gutkovskiy, Dec 10 2016
Sum_{n>=0} 1/a(n) = 29/24. - Bernard Schott, Oct 25 2021

Extensions

Definition corrected by Frans J. Faase, Feb 07 2009
Edited by N. J. A. Sloane, Dec 04 2009

A003953 Expansion of g.f.: (1+x)/(1-10*x).

Original entry on oeis.org

1, 11, 110, 1100, 11000, 110000, 1100000, 11000000, 110000000, 1100000000, 11000000000, 110000000000, 1100000000000, 11000000000000, 110000000000000, 1100000000000000, 11000000000000000
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 11.
a(n) is sequence A003945(n-1) written in base 2: a(0)=1, a(n) for n >= 1: 2 times 1, (n-1) times 0. a(n) is also A007283(n-1) and A042950(n) for n >= 1 written in base 2. a(n) is also A098011(n+3) and A101229(n+10) for n >= 1 written in base 2. a(n) is also abs(A110164(n+1)) for n >= 1 written in base 2. - Jaroslav Krizek, Aug 17 2009
a(n) equals the numbers of words of length n on alphabet {0,1,...,10} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015 [Corrected by David Nacin, Jun 02 2017]

Crossrefs

Programs

Formula

a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 9. - Philippe Deléham, Jul 10 2005
G.f.: (1+x)/(1-10*x). - Paul Barry, Mar 22 2006
a(0) = 1, a(n) = 10^n + 10^(n-1) = 11*10^(n-1) for n >= 1. - Jaroslav Krizek, Aug 17 2009
E.g.f.: (11*exp(10*x) - 1)/10. - G. C. Greubel, Sep 24 2019

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A003950 Expansion of g.f.: (1+x)/(1-7*x).

Original entry on oeis.org

1, 8, 56, 392, 2744, 19208, 134456, 941192, 6588344, 46118408, 322828856, 2259801992, 15818613944, 110730297608, 775112083256, 5425784582792, 37980492079544, 265863444556808, 1861044111897656, 13027308783283592, 91191161482985144, 638338130380896008
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 8.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954 m is 2, 3, 4, 5, 6 . - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001.
For n>=1, a(n) equals the number of words of length n on the alphabet {0,1,...,7} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015 [Corrected by David Nacin, May 31 2017]
a(n) is the number of octonary sequences of length n such that no two consecutive terms have distance 4. - David Nacin, May 31 2017

Crossrefs

Programs

Formula

a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 6. - Philippe Deléham, Jul 10 2005
From Philippe Deléham, Nov 21 2007: (Start)
a(n) = 8*7^(n-1) for n>=1, a(0)=1 .
G.f.: (1+x)/(1-7x).
The Hankel transform of this sequence is [1,-8,0,0,0,0,0,0,0,0,...]. (End)
a(0)=1, a(1)=8, a(n) = 7*a(n-1). - Vincenzo Librandi, Dec 10 2012
E.g.f.: (8*exp(7*x) - 1)/7. - G. C. Greubel, Sep 24 2019

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A003949 Expansion of g.f. (1+x)/(1-6*x).

Original entry on oeis.org

1, 7, 42, 252, 1512, 9072, 54432, 326592, 1959552, 11757312, 70543872, 423263232, 2539579392, 15237476352, 91424858112, 548549148672, 3291294892032, 19747769352192, 118486616113152, 710919696678912, 4265518180073472, 25593109080440832, 153558654482644992
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 7.
For n >= 1, a(n+1) is equal to the number of functions f:{1,2,...,n+1}->{1,2,3,4,5,6,7} such that for fixed, different x_1, x_2,...,x_n in {1,2,...,n+1} and fixed y_1, y_2,...,y_n in {1,2,3,4,5,6,7} we have f(x_i)<>y_i, (i=1,2,...,n). - Milan Janjic, May 10 2007
For n >= 1, a(n) equals the numbers of words of length n-1 on alphabet {0,1,2,3,5,6} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015

Crossrefs

Programs

  • GAP
    k:=7;; Concatenation([1], List([1..25], n-> k*(k-1)^(n-1) )); # G. C. Greubel, Sep 24 2019
  • Magma
    k:=7; [1] cat [k*(k-1)^(n-1): n in [1..25]]; // G. C. Greubel, Sep 24 2019
    
  • Magma
    R:=PowerSeriesRing(Rationals(), 25); Coefficients(R!( (1+x)/(1-6*x))); // Marius A. Burtea, Jan 20 2020
    
  • Maple
    k:=7; seq(`if`(n=0, 1, k*(k-1)^(n-1)), n = 0..25); # modified by G. C. Greubel, Sep 24 2019
  • Mathematica
    q = 7; Join[{a = 1}, Table[If[n != 0, a = q*a - a, a = q*a], {n, 0, 25}]] (* or *) Join[{1}, 7*6^Range[0, 25]] (* Vladimir Joseph Stephan Orlovsky, Jul 11 2011 *)
    CoefficientList[Series[(1+x)/(1-6*x), {x, 0, 30}], x] (* Vincenzo Librandi, Dec 10 2012 *)
    LinearRecurrence[{6},{1,7},30] (* or *) Join[{1},NestList[6#&,7,30]] (* Harvey P. Dale, May 03 2025 *)
  • PARI
    a(n)=if(n,7*6^(n-1),1) \\ Charles R Greathouse IV, Mar 22 2016
    
  • Sage
    k=7; [1]+[k*(k-1)^(n-1) for n in (1..25)] # G. C. Greubel, Sep 24 2019
    

Formula

G.f.: (1+x)/(1-6*x).
a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 5. - Philippe Deléham, Jul 10 2005
a(0)=1; for n > 0, a(n) = 7*6^(n-1). - Vincenzo Librandi, Nov 18 2010
a(0)=1, a(1)=7, a(n) = 6*a(n-1). - Vincenzo Librandi, Dec 10 2012
E.g.f.: (7*exp(6*x) - 1)/6. - G. C. Greubel, Sep 24 2019

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A003952 Expansion of g.f.: (1+x)/(1-9*x).

Original entry on oeis.org

1, 10, 90, 810, 7290, 65610, 590490, 5314410, 47829690, 430467210, 3874204890, 34867844010, 313810596090, 2824295364810, 25418658283290, 228767924549610, 2058911320946490, 18530201888518410, 166771816996665690, 1500946352969991210, 13508517176729920890
Offset: 0

Views

Author

N. J. A. Sloane, Mar 15 1996

Keywords

Comments

Coordination sequence for infinite tree with valency 10.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954 m is 2, 3, 4, 5, 6 . - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001
Except 1, all terms are in A033583. - Vincenzo Librandi, May 26 2014
For n>=1, a(n) equals the number of words of length n on alphabet {0,1,...,9} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015 [Corrected by David Nacin, May 31 2017]
a(n) is the number of sequences over the alphabet {0,1,...,9} of length n such that no two consecutive terms have distance 5. - David Nacin, May 31 2017

Crossrefs

Programs

Formula

a(n) = (10*9^n - 0^n)/9. Binomial transform is A000042. - Paul Barry, Jan 29 2004
G.f.: (1+x)/(1-9*x). - Philippe Deléham, Jan 31 2004
a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 8. - Philippe Deléham, Jul 10 2005
The Hankel transform of this sequence is: [1,-10,0,0,0,0,0,0,0,...]. - Philippe Deléham, Nov 21 2007
E.g.f.: (10*exp(9*x) - 1)/9. - G. C. Greubel, Sep 24 2019

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A003954 Expansion of g.f.: (1+x)/(1-11*x).

Original entry on oeis.org

1, 12, 132, 1452, 15972, 175692, 1932612, 21258732, 233846052, 2572306572, 28295372292, 311249095212, 3423740047332, 37661140520652, 414272545727172, 4556998002998892, 50126978032987812, 551396758362865932, 6065364341991525252, 66719007761906777772, 733909085380974555492
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 12.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954 m is 2, 3, 4, 5, 6 . - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001.
For n>=1, a(n) equals the numbers of words of length n-1 on alphabet {0,1,...,11} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015

Crossrefs

Programs

Formula

a(n) = Sum_{k=0..n} A029653(n,k)*x^k for x = 10. - Philippe Deléham, Jul 10 2005
G.f.: (1+x)/(1-11*x). The Hankel transform of this sequence is [1,-12,0,0,0,0,0,0,0,...]. - Philippe Deléham, Nov 21 2007
a(0) = 1; for n>0, a(n) = 12*11^(n-1). - Vincenzo Librandi, Nov 18 2010
a(0) = 1, a(1)=12, a(n) = 11*a(n-1). - Vincenzo Librandi, Dec 10 2012
E.g.f.: (12*exp(11*x) - 1)/11. - Elmo R. Oliveira, Mar 24 2025

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A003951 Expansion of g.f.: (1+x)/(1-8*x).

Original entry on oeis.org

1, 9, 72, 576, 4608, 36864, 294912, 2359296, 18874368, 150994944, 1207959552, 9663676416, 77309411328, 618475290624, 4947802324992, 39582418599936, 316659348799488, 2533274790395904, 20266198323167232, 162129586585337856, 1297036692682702848
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 9.
Binomial transform is {1, 10, 91, 820, 7381, ...}, see A002452. - Philippe Deléham, Jul 22 2005
a(n) equals the number of words of length n on alphabet {0,1,...,8} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015 [Corrected by David Nacin, May 31 2017]

Crossrefs

Cf. A003945.

Programs

Formula

a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 7. - Philippe Deléham, Jul 10 2005
a(0) = 1; for n>0, a(n) = 9*8^(n-1). - Vincenzo Librandi, Nov 18 2010
a(0) = 1, a(1) = 9, a(n) = 8*a(n-1). - Vincenzo Librandi, Dec 10 2012
E.g.f.: (9*exp(8*x) -1)/8. - G. C. Greubel, Sep 24 2019

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A078812 Triangle read by rows: T(n, k) = binomial(n+k-1, 2*k-1).

Original entry on oeis.org

1, 2, 1, 3, 4, 1, 4, 10, 6, 1, 5, 20, 21, 8, 1, 6, 35, 56, 36, 10, 1, 7, 56, 126, 120, 55, 12, 1, 8, 84, 252, 330, 220, 78, 14, 1, 9, 120, 462, 792, 715, 364, 105, 16, 1, 10, 165, 792, 1716, 2002, 1365, 560, 136, 18, 1, 11, 220, 1287, 3432, 5005, 4368, 2380, 816, 171, 20, 1
Offset: 0

Views

Author

Michael Somos, Dec 05 2002

Keywords

Comments

Warning: formulas and programs sometimes refer to offset 0 and sometimes to offset 1.
Apart from signs, identical to A053122.
Coefficient array for Morgan-Voyce polynomial B(n,x); see A085478 for references. - Philippe Deléham, Feb 16 2004
T(n,k) is the number of compositions of n having k parts when there are q kinds of part q (q=1,2,...). Example: T(4,2) = 10 because we have (1,3),(1,3'),(1,3"), (3,1),(3',1),(3",1),(2,2),(2,2'),(2',2) and (2',2'). - Emeric Deutsch, Apr 09 2005
T(n, k) is also the number of idempotent order-preserving full transformations (of an n-chain) of height k (height(alpha) = |Im(alpha)|). - Abdullahi Umar, Oct 02 2008
This sequence is jointly generated with A085478 as a triangular array of coefficients of polynomials v(n,x): initially, u(1,x) = v(1,x) = 1; for n > 1, u(n,x) = u(n-1,x) + x*v(n-1)x and v(n,x) = u(n-1,x) + (x+1)*v(n-1,x). See the Mathematica section. - Clark Kimberling, Feb 25 2012
Concerning Kimberling's recursion relations, see A102426. - Tom Copeland, Jan 19 2016
Subtriangle of the triangle T(n,k), 0 <= k <= n, read by rows, given by (0, 2, -1/2, 1/2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Mar 27 2012
From Wolfdieter Lang, Aug 30 2012: (Start)
With offset [0,0] the triangle with entries R(n,k) = T(n+1,k+1):= binomial(n+k+1, 2*k+1), n >= k >= 0, and zero otherwise, becomes the Riordan lower triangular convolution matrix R = (G(x)/x, G(x)) with G(x):=x/(1-x)^2 (o.g.f. of A000027). This means that the o.g.f. of column number k of R is (G(x)^(k+1))/x. This matrix R is the inverse of the signed Riordan lower triangular matrix A039598, called in a comment there S.
The Riordan matrix with entries R(n,k), just defined, provides the transition matrix between the sequence entry F(4*m*(n+1))/L(2*l), with m >= 0, for n=0,1,... and the sequence entries 5^k*F(2*m)^(2*k+1) for k = 0,1,...,n, with F=A000045 (Fibonacci) and L=A000032 (Lucas). Proof: from the inverse of the signed triangle Riordan matrix S used in a comment on A039598.
For the transition matrix R (T with offset [0,0]) defined above, row n=2: F(12*m) /L(2*m) = 3*5^0*F(2*m)^1 + 4*5^1*F(2*m)^3 + 1*5^2*F(2*m)^5, m >= 0. (End)
From R. Bagula's comment in A053122 (cf. Damianou link p. 10), this array gives the coefficients (mod sign) of the characteristic polynomials for the Cartan matrix of the root system A_n. - Tom Copeland, Oct 11 2014
For 1 <= k <= n, T(n,k) equals the number of (n-1)-length ternary words containing k-1 letters equal 2 and avoiding 01. - Milan Janjic, Dec 20 2016
The infinite sum (Sum_{i >= 0} (T(s+i,1+i) / 2^(s+2*i)) * zeta(s+1+2*i)) = 1 allows any zeta(s+1) to be expressed as a sum of rational multiples of zeta(s+1+2*i) having higher arguments. For example, zeta(3) can be expressed as a sum involving zeta(5), zeta(7), etc. The summation for each s >= 1 uses the s-th diagonal of the triangle. - Robert B Fowler, Feb 23 2022
The convolution triangle of the nonnegative integers. - Peter Luschny, Oct 07 2022

Examples

			Triangle begins, 1 <= k <= n:
                          1
                        2   1
                      3   4   1
                    4  10   6   1
                  5  20  21   8   1
                6  35  56  36  10   1
              7  56 126 120  55  12   1
            8  84 252 330 220  78  14   1
From _Peter Bala_, Feb 11 2025: (Start)
The array factorizes as an infinite product of lower triangular arrays:
  / 1               \    / 1              \ / 1              \ / 1             \
  | 2    1           |   | 2   1          | | 0  1           | | 0  1          |
  | 3    4   1       | = | 3   2   1      | | 0  2   1       | | 0  0  1       | ...
  | 4   10   6   1   |   | 4   3   2  1   | | 0  3   2  1    | | 0  0  2  1    |
  | 5   20  21   8  1|   | 5   4   3  2  1| | 0  4   3  2  1 | | 0  0  3  2  1 |
  |...               |   |...             | |...             | |...            |
Cf. A092276. (End)
		

Crossrefs

This triangle is formed from odd-numbered rows of triangle A011973 read in reverse order.
Row sums give A001906. With signs: A053122.
The column sequences are A000027, A000292, A000389, A000580, A000582, A001288 for k=1..6, resp. For k=7..24 they are A010966..(+2)..A011000 and for k=25..50 they are A017713..(+2)..A017763.

Programs

  • GAP
    Flat(List([0..12], n-> List([0..n], k-> Binomial(n+k+1, 2*k+1) ))); # G. C. Greubel, Aug 01 2019
  • Haskell
    a078812 n k = a078812_tabl !! n !! k
    a078812_row n = a078812_tabl !! n
    a078812_tabl = [1] : [2, 1] : f [1] [2, 1] where
       f us vs = ws : f vs ws where
         ws = zipWith (-) (zipWith (+) ([0] ++ vs) (map (* 2) vs ++ [0]))
                          (us ++ [0, 0])
    -- Reinhard Zumkeller, Dec 16 2013
    
  • Magma
    /* As triangle */ [[Binomial(n+k-1, 2*k-1): k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Jun 01 2018
    
  • Maple
    for n from 1 to 11 do seq(binomial(n+k-1,2*k-1),k=1..n) od; # yields sequence in triangular form; Emeric Deutsch, Apr 09 2005
    # Uses function PMatrix from A357368. Adds a row and column above and to the left.
    PMatrix(10, n -> n); # Peter Luschny, Oct 07 2022
  • Mathematica
    (* First program *)
    u[1, x_]:= 1; v[1, x_]:= 1; z = 13;
    u[n_, x_]:= u[n-1, x] + x*v[n-1, x];
    v[n_, x_]:= u[n-1, x] + (x+1)*v[n-1, x];
    Table[Expand[u[n, x]], {n, 1, z/2}]
    Table[Expand[v[n, x]], {n, 1, z/2}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%] (* A085478 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%] (* A078812 *) (* Clark Kimberling, Feb 25 2012 *)
    (* Second program *)
    Table[Binomial[n+k+1, 2*k+1], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Aug 01 2019 *)
  • Maxima
    T(n,m):=sum(binomial(2*k,n-m)*binomial(m+k,k)*(-1)^(n-m+k)*binomial(n+1,m+k+1),k,0,n-m); /* Vladimir Kruchinin, Apr 13 2016 */
    
  • PARI
    {T(n, k) = if( n<0, 0, binomial(n+k-1, 2*k-1))};
    
  • PARI
    {T(n, k) = polcoeff( polcoeff( x*y / (1 - (2 + y) * x + x^2) + x * O(x^n), n), k)};
    
  • Sage
    @cached_function
    def T(k,n):
        if k==n: return 1
        if k==0: return 0
        return sum(i*T(k-1,n-i) for i in (1..n-k+1))
    A078812 = lambda n,k: T(k,n)
    [[A078812(n,k) for k in (1..n)] for n in (1..8)] # Peter Luschny, Mar 12 2016
    
  • Sage
    [[binomial(n+k+1, 2*k+1) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Aug 01 2019
    

Formula

G.f.: x*y / (1 - (2 + y)*x + x^2). To get row n, expand this in powers of x then expand the coefficient of x^n in increasing powers of y.
From Philippe Deléham, Feb 16 2004: (Start)
If indexing begins at 0 we have
T(n,k) = (n+k+1)!/((n-k)!*(2k+1))!.
T(n,k) = Sum_{j>=0} T(n-1-j, k-1)*(j+1) with T(n, 0) = n+1, T(n, k) = 0 if n < k.
T(n,k) = T(n-1, k-1) + T(n-1, k) + Sum_{j>=0} (-1)^j*T(n-1, k+j)*A000108(j) with T(n,k) = 0 if k < 0, T(0, 0)=1 and T(0, k) = 0 for k > 0.
G.f. for the column k: Sum_{n>=0} T(n, k)*x^n = (x^k)/(1-x)^(2k+2).
Row sums: Sum_{k>=0} T(n, k) = A001906(n+1). (End)
Antidiagonal sums are A000079(n) = Sum_{k=0..floor(n/2)} binomial(n+k+1, n-k). - Paul Barry, Jun 21 2004
Riordan array (1/(1-x)^2, x/(1-x)^2). - Paul Barry, Oct 22 2006
T(0,0) = 1, T(n,k) = 0 if k < 0 or if k > n, T(n,k) = T(n-1,k-1) + 2*T(n-1,k) - T(n-2,k). - Philippe Deléham, Jan 26 2010
For another version see A128908. - Philippe Deléham, Mar 27 2012
T(n,m) = Sum_{k=0..n-m} (binomial(2*k,n-m)*binomial(m+k,k)*(-1)^(n-m+k)* binomial(n+1,m+k+1)). - Vladimir Kruchinin, Apr 13 2016
T(n, k) = T(n-1, k) + (T(n-1, k-1) + T(n-2, k-1) + T(n-3, k-1) + ...) for k >= 2 with T(n, 1) = n. - Peter Bala, Feb 11 2025
From Peter Bala, May 04 2025: (Start)
With the column offset starting at 0, the n-th row polynomial B(n, x) = 1/sqrt(x + 4) * Chebyshev_U(2*n+1, (1/2)*sqrt(x + 4)) = (-1)^n * Chebyshev_U(n, -(1/2)*(x + 2)).
B(n, x) / Product_{k = 1..2*n} (1 + 1/B(k, x)) = b(n, x), the n-th row polynomial of A085478. (End)

Extensions

Edited by N. J. A. Sloane, Apr 28 2008
Previous Showing 21-30 of 56 results. Next