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-9 of 9 results.

A038207 Triangle whose (i,j)-th entry is binomial(i,j)*2^(i-j).

Original entry on oeis.org

1, 2, 1, 4, 4, 1, 8, 12, 6, 1, 16, 32, 24, 8, 1, 32, 80, 80, 40, 10, 1, 64, 192, 240, 160, 60, 12, 1, 128, 448, 672, 560, 280, 84, 14, 1, 256, 1024, 1792, 1792, 1120, 448, 112, 16, 1, 512, 2304, 4608, 5376, 4032, 2016, 672, 144, 18, 1, 1024, 5120, 11520, 15360, 13440, 8064, 3360, 960, 180, 20, 1
Offset: 0

Views

Author

Keywords

Comments

This infinite matrix is the square of the Pascal matrix (A007318) whose rows are [ 1,0,... ], [ 1,1,0,... ], [ 1,2,1,0,... ], ...
As an upper right triangle, table rows give number of points, edges, faces, cubes,
4D hypercubes etc. in hypercubes of increasing dimension by column. - Henry Bottomley, Apr 14 2000. More precisely, the (i,j)-th entry is the number of j-dimensional subspaces of an i-dimensional hypercube (see the Coxeter reference). - Christof Weber, May 08 2009
Number of different partial sums of 1+[1,1,2]+[2,2,3]+[3,3,4]+[4,4,5]+... with entries that are zero removed. - Jon Perry, Jan 01 2004
Row sums are powers of 3 (A000244), antidiagonal sums are Pell numbers (A000129). - Gerald McGarvey, May 17 2005
Riordan array (1/(1-2x), x/(1-2x)). - Paul Barry, Jul 28 2005
T(n,k) is the number of elements of the Coxeter group B_n with descent set contained in {s_k}, 0<=k<=n-1. For T(n,n), we interpret this as the number of elements of B_n with empty descent set (since s_n does not exist). - Elizabeth Morris (epmorris(AT)math.washington.edu), Mar 01 2006
Let S be a binary relation on the power set P(A) of a set A having n = |A| elements such that for every element x, y of P(A), xSy if x is a subset of y. Then T(n,k) = the number of elements (x,y) of S for which y has exactly k more elements than x. - Ross La Haye, Oct 12 2007
T(n,k) is number of paths in the first quadrant going from (0,0) to (n,k) using only steps B=(1,0) colored blue, R=(1,0) colored red and U=(1,1). Example: T(3,2)=6 because we have BUU, RUU, UBU, URU, UUB and UUR. - Emeric Deutsch, Nov 04 2007
T(n,k) is the number of lattice paths from (0,0) to (n,k) using steps (0,1), and two kinds of step (1,0). - Joerg Arndt, Jul 01 2011
T(i,j) is the number of i-permutations of {1,2,3} containing j 1's. Example: T(2,1)=4 because we have 12, 13, 21 and 31; T(3,2)=6 because we have 112, 113, 121, 131, 211 and 311. - Zerinvary Lajos, Dec 21 2007
Triangle of coefficients in expansion of (2+x)^n. - N-E. Fahssi, Apr 13 2008
Sum of diagonals are Jacobsthal-numbers: A001045. - Mark Dols, Aug 31 2009
Triangle T(n,k), read by rows, given by [2,0,0,0,0,0,0,0,...] DELTA [1,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 15 2009
Eigensequence of the triangle = A004211: (1, 3, 11, 49, 257, 1539, ...). - Gary W. Adamson, Feb 07 2010
f-vectors ("face"-vectors) for n-dimensional cubes [see e.g., Hoare]. (This is a restatement of Bottomley's above.) - Tom Copeland, Oct 19 2012
With P = Pascal matrix, the sequence of matrices I, A007318, A038207, A027465, A038231, A038243, A038255, A027466 ... = P^0, P^1, P^2, ... are related by Copeland's formula below to the evolution at integral time steps n= 0, 1, 2, ... of an exponential distribution exp(-x*z) governed by the Fokker-Planck equation as given in the Dattoli et al. ref. below. - Tom Copeland, Oct 26 2012
The matrix elements of the inverse are T^(-1)(n,k) = (-1)^(n+k)*T(n,k). - R. J. Mathar, Mar 12 2013
Unsigned diagonals of A133156 are rows of this array. - Tom Copeland, Oct 11 2014
Omitting the first row, this is the production matrix for A039683, where an equivalent differential operator can be found. - Tom Copeland, Oct 11 2016
T(n,k) is the number of functions f:[n]->[3] with exactly k elements mapped to 3. Note that there are C(n,k) ways to choose the k elements mapped to 3, and there are 2^(n-k) ways to map the other (n-k) elements to {1,2}. Hence, by summing T(n,k) as k runs from 0 to n, we obtain 3^n = Sum_{k=0..n} T(n,k). - Dennis P. Walsh, Sep 26 2017
Since this array is the square of the Pascal lower triangular matrix, the row polynomials of this array are obtained as the umbral composition of the row polynomials P_n(x) of the Pascal matrix with themselves. E.g., P_3(P.(x)) = 1 P_3(x) + 3 P_2(x) + 3 P_1(x) + 1 = (x^3 + 3 x^2 + 3 x + 1) + 3 (x^2 + 2 x + 1) + 3 (x + 1) + 1 = x^3 + 6 x^2 + 12 x + 8. - Tom Copeland, Nov 12 2018
T(n,k) is the number of 2-compositions of n+1 with some zeros allowed that have k zeros; see the Hopkins & Ouvry reference. - Brian Hopkins, Aug 16 2020
Also the convolution triangle of A000079. - Peter Luschny, Oct 09 2022

Examples

			Triangle begins with T(0,0):
   1;
   2,  1;
   4,  4,  1;
   8, 12,  6,  1;
  16, 32, 24,  8,  1;
  32, 80, 80, 40, 10,  1;
  ... -  corrected by _Clark Kimberling_, Aug 05 2011
Seen as an array read by descending antidiagonals:
[0] 1, 2,  4,   8,    16,    32,    64,     128,     256, ...     [A000079]
[1] 1, 4,  12,  32,   80,    192,   448,    1024,    2304, ...    [A001787]
[2] 1, 6,  24,  80,   240,   672,   1792,   4608,    11520, ...   [A001788]
[3] 1, 8,  40,  160,  560,   1792,  5376,   15360,   42240, ...   [A001789]
[4] 1, 10, 60,  280,  1120,  4032,  13440,  42240,   126720, ...  [A003472]
[5] 1, 12, 84,  448,  2016,  8064,  29568,  101376,  329472, ...  [A054849]
[6] 1, 14, 112, 672,  3360,  14784, 59136,  219648,  768768, ...  [A002409]
[7] 1, 16, 144, 960,  5280,  25344, 109824, 439296,  1647360, ... [A054851]
[8] 1, 18, 180, 1320, 7920,  41184, 192192, 823680,  3294720, ... [A140325]
[9] 1, 20, 220, 1760, 11440, 64064, 320320, 1464320, 6223360, ... [A140354]
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 155.
  • H. S. M. Coxeter, Regular Polytopes, Dover Publications, New York (1973), p. 122.

Crossrefs

Programs

  • GAP
    Flat(List([0..15], n->List([0..n], k->Binomial(n, k)*2^(n-k)))); # Stefano Spezia, Nov 21 2018
  • Haskell
    a038207 n = a038207_list !! n
    a038207_list = concat $ iterate ([2,1] *) [1]
    instance Num a => Num [a] where
       fromInteger k = [fromInteger k]
       (p:ps) + (q:qs) = p + q : ps + qs
       ps + qs         = ps ++ qs
       (p:ps) * qs'@(q:qs) = p * q : ps * qs' + [p] * qs
        *                = []
    -- Reinhard Zumkeller, Apr 02 2011
    
  • Haskell
    a038207' n k = a038207_tabl !! n !! k
    a038207_row n = a038207_tabl !! n
    a038207_tabl = iterate f [1] where
       f row = zipWith (+) ([0] ++ row) (map (* 2) row ++ [0])
    -- Reinhard Zumkeller, Feb 27 2013
    
  • Magma
    /* As triangle */ [[(&+[Binomial(n,i)*Binomial(i,k): i in [k..n]]): k in [0..n]]: n in [0..15]]; // Vincenzo Librandi, Nov 16 2018
    
  • Maple
    for i from 0 to 12 do seq(binomial(i, j)*2^(i-j), j = 0 .. i) end do; # yields sequence in triangular form - Emeric Deutsch, Nov 04 2007
    # Uses function PMatrix from A357368. Adds column 1, 0, 0, ... to the left.
    PMatrix(10, n -> 2^(n-1)); # Peter Luschny, Oct 09 2022
  • Mathematica
    Table[CoefficientList[Expand[(y + x + x^2)^n], y] /. x -> 1, {n, 0,10}] // TableForm (* Geoffrey Critzer, Nov 20 2011 *)
    Table[Binomial[n,k]2^(n-k),{n,0,10},{k,0,n}]//Flatten (* Harvey P. Dale, May 22 2020 *)
  • PARI
    {T(n, k) = polcoeff((x+2)^n, k)}; /* Michael Somos, Apr 27 2000 */
    
  • Sage
    def A038207_triangle(dim):
        M = matrix(ZZ,dim,dim)
        for n in range(dim): M[n,n] = 1
        for n in (1..dim-1):
            for k in (0..n-1):
                M[n,k] = M[n-1,k-1]+2*M[n-1,k]
        return M
    A038207_triangle(9)  # Peter Luschny, Sep 20 2012
    

Formula

T(n, k) = Sum_{i=0..n} binomial(n,i)*binomial(i,k).
T(n, k) = (-1)^k*A065109(n,k).
G.f.: 1/(1-2*z-t*z). - Emeric Deutsch, Nov 04 2007
Rows of the triangle are generated by taking successive iterates of (A135387)^n * [1, 0, 0, 0, ...]. - Gary W. Adamson, Dec 09 2007
From the formalism of A133314, the e.g.f. for the row polynomials of A038207 is exp(x*t)*exp(2x). The e.g.f. for the row polynomials of the inverse matrix is exp(x*t)*exp(-2x). p iterates of the matrix give the matrix with e.g.f. exp(x*t)*exp(p*2x). The results generalize for 2 replaced by any number. - Tom Copeland, Aug 18 2008
Sum_{k=0..n} T(n,k)*x^k = (2+x)^n. - Philippe Deléham, Dec 15 2009
n-th row is obtained by taking pairwise sums of triangle A112857 terms starting from the right. - Gary W. Adamson, Feb 06 2012
T(n,n) = 1 and T(n,k) = T(n-1,k-1) + 2*T(n-1,k) for kJon Perry, Oct 11 2012
The e.g.f. for the n-th row is given by umbral composition of the normalized Laguerre polynomials A021009 as p(n,x) = L(n, -L(.,-x))/n! = 2^n L(n, -x/2)/n!. E.g., L(2,x) = 2 -4*x +x^2, so p(2,x)= (1/2)*L(2, -L(.,-x)) = (1/2)*(2*L(0,-x) + 4*L(1,-x) + L(2,-x)) = (1/2)*(2 + 4*(1+x) + (2+4*x+x^2)) = 4 + 4*x + x^2/2. - Tom Copeland, Oct 20 2012
From Tom Copeland, Oct 26 2012: (Start)
From the formalism of A132440 and A218272:
Let P and P^T be the Pascal matrix and its transpose and H= P^2= A038207.
Then with D the derivative operator,
exp(x*z/(1-2*z))/(1-2*z)= exp(2*z D_z z) e^(x*z)= exp(2*D_x (x D_x)) e^(z*x)
= (1 z z^2 z^3 ...) H (1 x x^2/2! x^3/3! ...)^T
= (1 x x^2/2! x^3/3! ...) H^T (1 z z^2 z^3 ...)^T
= Sum_{n>=0} z^n * 2^n Lag_n(-x/2)= exp[z*EF(.,x)], an o.g.f. for the f-vectors (rows) of A038207 where EF(n,x) is an e.g.f. for the n-th f-vector. (Lag_n(x) are the un-normalized Laguerre polynomials.)
Conversely,
exp(z*(2+x))= exp(2D_x) exp(x*z)= exp(2x) exp(x*z)
= (1 x x^2 x^3 ...) H^T (1 z z^2/2! z^3/3! ...)^T
= (1 z z^2/2! z^3/3! ...) H (1 x x^2 x^3 ...)^T
= exp(z*OF(.,x)), an e.g.f for the f-vectors of A038207 where
OF(n,x)= (2+x)^n is an o.g.f. for the n-th f-vector.
(End)
G.f.: R(0)/2, where R(k) = 1 + 1/(1 - (2*k+1+ (1+y))*x/((2*k+2+ (1+y))*x + 1/R(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Nov 09 2013
A038207 = exp[M*B(.,2)] where M = A238385-I and (B(.,x))^n = B(n,x) are the Bell polynomials (cf. A008277). B(n,2) = A001861(n). - Tom Copeland, Apr 17 2014
T = (A007318)^2 = A112857*|A167374| = |A118801|*|A167374| = |A118801*A167374| = |P*A167374*P^(-1)*A167374| = |P*NpdP*A167374|. Cf. A118801. - Tom Copeland, Nov 17 2016
E.g.f. for the n-th subdiagonal, n = 0,1,2,..., equals exp(x)*P(n,x), where P(n,x) is the polynomial 2^n*Sum_{k = 0..n} binomial(n,k)*x^k/k!. For example, the e.g.f. for the third subdiagonal is exp(x)*(8 + 24*x + 12*x^2 + 4*x^3/3) = 8 + 32*x + 80*x^2/2! + 160*x^3/3! + .... - Peter Bala, Mar 05 2017
T(3*k+2,k) = T(3*k+2,k+1), T(2*k+1,k) = 2*T(2*k+1,k+1). - Yuchun Ji, May 26 2020
From Robert A. Russell, Aug 05 2020: (Start)
G.f. for column k: x^k / (1-2*x)^(k+1).
E.g.f. for column k: exp(2*x) * x^k / k!. (End)
Also the array A(n, k) read by descending antidiagonals, where A(n, k) = (-1)^n*Sum_{j= 0..n+k} binomial(n + k, j)*hypergeom([-n, j+1], [1], 1). - Peter Luschny, Nov 09 2021

A090017 a(n) = 4*a(n-1) + 2*a(n-2) for n>1, a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 4, 18, 80, 356, 1584, 7048, 31360, 139536, 620864, 2762528, 12291840, 54692416, 243353344, 1082798208, 4817899520, 21437194496, 95384577024, 424412697088, 1888419942400, 8402505163776, 37386860539904, 166352452487168
Offset: 0

Views

Author

Paul Barry, Nov 19 2003

Keywords

Comments

Starting with "1" = INVERT transform of A007482: (1, 3, 11, 39, 139, ...). - Gary W. Adamson, Aug 06 2010
This is the Lucas sequence U(4,-2). - Bruno Berselli, Jan 09 2013
Lower left term in matrix powers of [(1,5); (1,3)]. Convolved with (1, 2, 0, 0, 0, ...) the result is A164549: (1, 6, 26, 116, ...). - Gary W. Adamson, Aug 10 2016
For n>0, a(n) equals the number of words of length n-1 over {0,1,2,3,4,5} in which 0 and 1 avoid runs of odd lengths. - Milan Janjic, Jan 08 2017

Crossrefs

Programs

  • Magma
    I:=[0,1]; [n le 2 select I[n] else 4*Self(n-1)+2*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Oct 12 2011
    
  • Mathematica
    a[n_Integer] := (-I Sqrt[2])^(n - 1) ChebyshevU[ n - 1, I Sqrt[2] ]
    a[n_]:=(MatrixPower[{{1,5},{1,3}},n].{{1},{1}})[[2,1]]; Table[Abs[a[n]],{n,-1,40}] (* Vladimir Joseph Stephan Orlovsky, Feb 19 2010 *)
    t={0,1};Do[AppendTo[t,4*t[[-1]]+2*t[[-2]]],{n,2,23}];t (* or *) LinearRecurrence[{4,2},{0,1},24] (* Indranil Ghosh, Feb 21 2017 *)
  • PARI
    Vec(x/(1-4*x-2*x^2)+O(x^99)) \\ Charles R Greathouse IV, Oct 12 2011
  • Sage
    [lucas_number1(n, 4, -2) for n in range(0, 23)] # Zerinvary Lajos, Apr 23 2009
    

Formula

G.f.: x/(1 - 4*x - 2*x^2).
a(n) = (-i*sqrt(2))^(n-1) U(n-1, i*sqrt(2)) where U is the Chebyshev polynomial of the second kind and i^2 = -1.
a(n) = ((2+sqrt(6))^n - (2-sqrt(6))^n)/(2 sqrt(6)). - Al Hakanson (hawkuu(AT)gmail.com), Jan 05 2009, Jan 07 2009
a(n+1) = Sum_{k=0..n} A099089(n,k)*2^k. - Philippe Deléham, Nov 21 2011
From Ilya Gutkovskiy, Aug 22 2016: (Start)
E.g.f.: sinh(sqrt(6)*x)*exp(2*x)/sqrt(6).
Number of zeros in substitution system {0 -> 11, 1 -> 11011} at step n from initial string "1" (1 -> 11011 -> 1101111011111101111011 -> ...). (End)

Extensions

Edited by Stuart Clary, Oct 25 2009

A053117 Triangle read by rows of coefficients of Chebyshev's U(n,x) polynomials (exponents in increasing order).

Original entry on oeis.org

1, 0, 2, -1, 0, 4, 0, -4, 0, 8, 1, 0, -12, 0, 16, 0, 6, 0, -32, 0, 32, -1, 0, 24, 0, -80, 0, 64, 0, -8, 0, 80, 0, -192, 0, 128, 1, 0, -40, 0, 240, 0, -448, 0, 256, 0, 10, 0, -160, 0, 672, 0, -1024, 0, 512, -1, 0, 60, 0, -560, 0, 1792, 0, -2304, 0, 1024, 0, -12, 0, 280, 0, -1792, 0, 4608, 0, -5120, 0, 2048, 1, 0, -84, 0, 1120, 0, -5376, 0, 11520, 0, -11264, 0, 4096
Offset: 0

Views

Author

Keywords

Comments

G.f. for row polynomials U(n,x) (signed triangle): 1/(1-2*x*z+z^2). Unsigned triangle |a(n,m)| has Fibonacci polynomials F(n+1,2*x) as row polynomials with g.f. 1/(1-2*x*z-z^2).
Row sums (unsigned triangle) A000129(n+1) (Pell). Row sums (signed triangle) A000027(n+1) (natural numbers).
The o.g.f. for the Legendre polynomials L(n,x) is 1 / sqrt(1- 2x*z + z^2), and squaring it gives the o.g.f. of this entry, so Sum_{k=0..n} L(k,x) L(n-k,x) = U(n,x). This reduces to U(n,x) = L(n/2,x)^2 + 2*Sum_{k=0...n/2-1} L(k,x) L(n-k,x) for n even and U(n,x) = 2*Sum_{k=0..(n-1)/2} L(k,x) L(n-k.x) for odd n. (Cf. also Allouche et al.) For a connection through the Legendre polynomials to elliptic curves and modular forms, see the MathOverflow question below. For the normalized Legendre polynomials, see A100258. (Cf. A097610 with h1 = -2x and h2 = 1, A207538, A099089 and A133156.) - Tom Copeland, Feb 04 2016
The compositional inverse of the shifted o.g.f. x / (1 + 2xz + z^2) for differently signed row polynomials of this entry is the shifted o.g.f. of A121448. The unsigned, non-vanishing antidiagonals (top to bottom) of this triangle are the rows of A038207. - Tom Copeland, Feb 08 2016

Examples

			Triangle begins:
   1;
   0,  2;
  -1,  0,   4;
   0, -4,   0, 8;
   1,  0, -12, 0, 16;
  ...
E.g., fourth row (n=3) {0,-4,0,8} corresponds to polynomial U(3,x) = -4*x + 8*x^3.
		

References

  • Theodore J. Rivlin, Chebyshev polynomials: from approximation theory to algebra and number theory, 2. ed., Wiley, New York, 1990.
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 22, page 196.

Crossrefs

Programs

  • Julia
    using Nemo
    function A053117Row(n)
        R, x = PolynomialRing(ZZ, "x")
        p = chebyshev_u(n, x)
        [coeff(p, j) for j in 0:n] end
    for n in 0:6 A053117Row(n) |> println end # Peter Luschny, Mar 13 2018
  • Maple
    seq(seq(coeff(orthopoly[U](n,x),x,j),j=0..n),n=0..16); # Robert Israel, Feb 09 2016
  • Mathematica
    Flatten[ Table[ CoefficientList[ ChebyshevU[n, x], x], {n, 0, 12}]](* Jean-François Alcover, Nov 24 2011 *)
  • PARI
    T(n, k) = polcoeff(polchebyshev(n,2), k); \\ Michel Marcus, Feb 10 2016
    

Formula

a(n, m) = (2^m)*A049310(n,m).
a(n, m) := 0 if n
If n and k are of the same parity then a(n,k)=(-1)^((n-k)/2)*sum(binomial((n+k)/2,i)*binomial((n+k)/2-i,(n-k)/2),i=0..k) and a(n,k)=0 otherwise. - Milan Janjic, Apr 13 2008

A133156 Irregular triangle read by rows: coefficients of U(n,x), Chebyshev polynomials of the second kind with exponents in decreasing order.

Original entry on oeis.org

1, 2, 4, -1, 8, -4, 16, -12, 1, 32, -32, 6, 64, -80, 24, -1, 128, -192, 80, -8, 256, -448, 240, -40, 1, 512, -1024, 672, -160, 10, 1024, -2304, 1792, -560, 60, -1, 2048, -5120, 4608, -1792, 280, -12, 4096, -11264, 11520, -5376, 1120, -84, 1
Offset: 0

Author

Gary W. Adamson, Dec 16 2007

Keywords

Comments

The Chebyshev polynomials of the second kind are defined by the recurrence relation: U(0,x) = 1; U(1,x) = 2x; U(n+1,x) = 2x*U(n,x) - U(n-1,x).
From Gary W. Adamson, Nov 28 2008: (Start)
Triangle read by rows, unsigned = A000012 * A028297.
Row sums of absolute values give the Pell series, A000129.
(End)
The row sums are {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, ...}.
Triangle, with zeros omitted, given by (2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, -1/2, 1/2, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 27 2011
Coefficients in the expansion of sin((n+1)*x)/sin(x) in descending powers of cos(x). The length of the n-th row is A008619(n). - Jianing Song, Nov 02 2018

Examples

			The first few Chebyshev polynomials of the second kind are
    1;
    2x;
    4x^2 -    1;
    8x^3 -    4x;
   16x^4 -   12x^2 +   1;
   32x^5 -   32x^3 +   6x;
   64x^6 -   80x^4 +  24x^2 -   1;
  128x^7 -  192x^5 +  80x^3 -   8x;
  256x^8 -  448x^6 + 240x^4 -  40x^2 +  1;
  512x^9 - 1024x^7 + 672x^5 - 160x^3 + 10x;
  ...
From _Roger L. Bagula_ and _Gary W. Adamson_: (Start)
     1;
     2;
     4,    -1;
     8,    -4;
    16,   -12,    1;
    32,   -32,    6;
    64,   -80,   24,   -1;
   128,  -192,   80,   -8;
   256,  -448,  240,  -40,  1;
   512, -1024,  672, -160, 10;
  1024, -2304, 1792, -560, 60, -1; (End)
From  _Philippe Deléham_, Dec 27 2011: (Start)
Triangle (2, 0, 0, 0, 0, ...) DELTA (0, -1/2, 1/2, 0, 0, 0, 0, 0, ...) begins:
   1;
   2,   0;
   4,  -1,  0;
   8,  -4,  0,  0;
  16, -12,  1,  0,  0;
  32, -32,  6,  0,  0,  0;
  64, -80, 24, -1,  0,  0,  0; (End)
		

Programs

  • Mathematica
    t[n_, m_] = (-1)^m*Binomial[n - m, m]*2^(n - 2*m);
    Table[Table[t[n, m], {m, 0, Floor[n/2]}], {n, 0, 10}];
    Flatten[%] (* Roger L. Bagula, Dec 19 2008 *)

Formula

A generating function for U(n) is 1/(1 - 2tx + t^2). Given A038207, shift down columns to allow for (1, 1, 2, 2, 3, 3, ...) terms in each row, then insert alternate signs.
T(n,m) = (-1)^m*binomial(n - m, m)*2^(n - 2*m). - Roger L. Bagula and Gary W. Adamson, Dec 19 2008
From Tom Copeland, Feb 11 2016: (Start)
Shifted o.g.f.: G(x,t) = x/(1 - 2x + tx^2).
A053117 is a reflected, aerated version of this entry; A207538, an unsigned version; and A099089, a reflected, shifted version.
The compositional inverse of G(x,t) is Ginv(x,t) = ((1 + 2x) - sqrt((1 + 2x)^2 - 4tx^2))/(2tx) = x - 2x^2 + (4 + t)x^3 - (8 + 6t)x^4 + ..., a shifted o.g.f. for A091894 (mod signs with A091894(0,0) = 0.). Cf. A097610 with h_1 = -2 and h_2 = t. (End)

Extensions

More terms from Philippe Deléham, Sep 12 2009

A207538 Triangle of coefficients of polynomials v(n,x) jointly generated with A207537; see Formula section.

Original entry on oeis.org

1, 2, 4, 1, 8, 4, 16, 12, 1, 32, 32, 6, 64, 80, 24, 1, 128, 192, 80, 8, 256, 448, 240, 40, 1, 512, 1024, 672, 160, 10, 1024, 2304, 1792, 560, 60, 1, 2048, 5120, 4608, 1792, 280, 12, 4096, 11264, 11520, 5376, 1120, 84, 1, 8192, 24576, 28160, 15360
Offset: 1

Author

Clark Kimberling, Feb 18 2012

Keywords

Comments

As triangle T(n,k) with 0<=k<=n and with zeros omitted, it is the triangle given by (2, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, 1/2, -1/2, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Mar 04 2012
The numbers in rows of the triangle are along "first layer" skew diagonals pointing top-left in center-justified triangle given in A013609 ((1+2*x)^n) and along (first layer) skew diagonals pointing top-right in center-justified triangle given in A038207 ((2+x)^n), see links. - Zagros Lalo, Jul 31 2018
If s(n) is the row sum at n, then the ratio s(n)/s(n-1) is approximately 2.414213562373095... (A014176: Decimal expansion of the silver mean, 1+sqrt(2)), when n approaches infinity. - Zagros Lalo, Jul 31 2018

Examples

			First seven rows:
1
2
4...1
8...4
16..12..1
32..32..6
64..80..24..1
(2, 0, 0, 0, 0, ...) DELTA (0, 1/2, -1/2, 0, 0, 0, ...) begins:
    1
    2,   0
    4,   1,  0
    8,   4,  0, 0
   16,  12,  1, 0, 0
   32,  32,  6, 0, 0, 0
   64,  80, 24, 1, 0, 0, 0
  128, 192, 80, 8, 0, 0, 0, 0
		

References

  • Shara Lalo and Zagros Lalo, Polynomial Expansion Theorems and Number Triangles, Zana Publishing, 2018, ISBN: 978-1-9995914-0-3, pp. 80-83, 357-358.

Programs

  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := u[n - 1, x] + (x + 1)*v[n - 1, x]
    v[n_, x_] := u[n - 1, x] + v[n - 1, x]
    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[%]  (* A207537, |A028297| *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]  (* A207538, |A133156| *)
    t[0, 0] = 1; t[n_, k_] := t[n, k] = If[n < 0 || k < 0, 0, 2 t[n - 1, k] + t[n - 2, k - 1]]; Table[t[n, k], {n, 0, 15}, {k, 0, Floor[n/2]}] // Flatten (* Zagros Lalo, Jul 31 2018 *)
    t[n_, k_] := t[n, k] = 2^(n - 2 k) * (n -  k)!/((n - 2 k)! k!) ; Table[t[n, k], {n, 0, 15}, {k, 0, Floor[n/2]} ]  // Flatten (* Zagros Lalo, Jul 31 2018 *)

Formula

u(n,x) = u(n-1,x)+(x+1)*v(n-1,x), v(n,x) = u(n-1,x)+v(n-1,x), where u(1,x) = 1, v(1,x) = 1. Also, A207538 = |A133156|.
From Philippe Deléham, Mar 04 2012: (Start)
With 0<=k<=n:
Mirror image of triangle in A099089.
Skew version of A038207.
Riordan array (1/(1-2*x), x^2/(1-2*x)).
G.f.: 1/(1-2*x-y*x^2).
Sum_{k, 0<=k<=n} T(n,k)*x^k = A190958(n+1), A127357(n), A090591(n), A089181(n+1), A088139(n+1), A045873(n+1), A088138(n+1), A088137(n+1), A099087(n), A000027(n+1), A000079(n), A000129(n+1), A002605(n+1), A015518(n+1), A063727(n), A002532(n+1), A083099(n+1), A015519(n+1), A003683(n+1), A002534(n+1), A083102(n), A015520(n+1), A091914(n) for x = -10, -9, -8, -7, -6, -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 respectively.
T(n,k) = 2*T(n-1,k) + T(-2,k-1) with T(0,0) = 1, T(1,0) = 2, T(1,1) = 0 and T(n, k) = 0 if k<0 or if k>n. (End)
T(n,k) = A013609(n-k, n-2*k+1). - Johannes W. Meijer, Sep 05 2013
From Tom Copeland, Feb 11 2016: (Start)
A053117 is a reflected, aerated and signed version of this entry. This entry belongs to a family discussed in A097610 with parameters h1 = -2 and h2 = -y.
Shifted o.g.f.: G(x,t) = x / (1 - 2 x - t x^2).
The compositional inverse of G(x,t) is Ginv(x,t) = -[(1 + 2x) - sqrt[(1+2x)^2 + 4t x^2]] / (2tx) = x - 2 x^2 + (4-t) x^3 - (8-6t) x^4 + ..., a shifted o.g.f. for A091894 (mod signs with A091894(0,0) = 0).
(End)

A090018 a(n) = 6*a(n-1) + 3*a(n-2) for n > 2, a(0)=1, a(1)=6.

Original entry on oeis.org

1, 6, 39, 252, 1629, 10530, 68067, 439992, 2844153, 18384894, 118841823, 768205620, 4965759189, 32099171994, 207492309531, 1341251373168, 8669985167601, 56043665125110, 362271946253463, 2341762672896108, 15137391876137037, 97849639275510546, 632510011281474387
Offset: 0

Author

Paul Barry, Nov 19 2003

Keywords

Comments

From Johannes W. Meijer, Aug 09 2010: (Start)
a(n) represents the number of n-move routes of a fairy chess piece starting in a given corner or side square on a 3 X 3 chessboard. This fairy chess piece behaves like a white queen on the eight side and corner squares but on the central square the queen explodes with fury and turns into a red queen, see A180032. The central square leads to A180028. (End)

Crossrefs

Sequences with g.f. of the form 1/(1 - 6*x - k*x^2): A106392 (k=-10), A027471 (k=-9), A006516 (k=-8), A081179 (k=-7), A030192 (k=-6), A003463 (k=-5), A084326 (k=-4), A138395 (k=-3), A154244 (k=-2), A001109 (k=-1), A000400 (k=0), A005668 (k=1), A135030 (k=2), this sequence (k=3), A135032 (k=4), A015551 (k=5), A057089 (k=6), A015552 (k=7), A189800 (k=8), A189801 (k=9), A190005 (k=10), A015553 (k=11).

Programs

  • Magma
    [n le 2 select 6^(n-1) else 6*Self(n-1)+3*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 15 2011
    
  • Maple
    a:= n-> (<<0|1>, <3|6>>^n. <<1,6>>)[1,1]:
    seq(a(n), n=0..30);  # Alois P. Heinz, Jan 17 2011
  • Mathematica
    Join[{a=1,b=6},Table[c=6*b+3*a;a=b;b=c,{n,100}]] (* Vladimir Joseph Stephan Orlovsky, Jan 16 2011 *)
    LinearRecurrence[{6,3}, {1,6}, 41] (* G. C. Greubel, Oct 10 2022 *)
  • PARI
    my(x='x+O('x^30)); Vec(1/(1-6*x-3*x^2)) \\ G. C. Greubel, Jan 24 2018
  • Sage
    [lucas_number1(n,6,-3) for n in range(1, 31)] # Zerinvary Lajos, Apr 24 2009
    

Formula

a(n) = (3+2*sqrt(3))^n*(sqrt(3)/4+1/2) + (1/2-sqrt(3)/4)*(3-2*sqrt(3))^n.
a(n) = (-i*sqrt(3))^n * ChebyshevU(n, isqrt(3)), i^2=-1.
From Johannes W. Meijer, Aug 09 2010: (Start)
G.f.: 1/(1 - 6*x - 3*x^2).
Limit_{k->oo} a(n+k)/a(k) = A141041(n) + A090018(n-1)*sqrt(12) for n >= 1.
Limit_{n->oo} A141041(n)/A090018(n-1) = sqrt(12). (End)
a(n) = Sum_{k=0..n} A099089(n,k)*3^k. - Philippe Deléham, Nov 21 2011
E.g.f.: exp(3*x)*(2*cosh(2*sqrt(3)*x) + sqrt(3)*sinh(2*sqrt(3)*x))/2. - Stefano Spezia, Apr 23 2025

Extensions

Typo in Mathematica program corrected by Vincenzo Librandi, Nov 15 2011

A190510 a(n) = 8*a(n-1) + 4*a(n-2), with a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 8, 68, 576, 4880, 41344, 350272, 2967552, 25141504, 213002240, 1804583936, 15288680448, 129527779328, 1097376956416, 9297126768640, 78766521974784, 667320682872832, 5653631550881792, 47898335138545664, 405801207311892480, 3438002999049322496
Offset: 0

Author

Keywords

Programs

  • Magma
    I:=[0, 1]; [n le 2 select I[n] else 8*Self(n-1)+4*Self(n-2): n in [1..20]]; // Vincenzo Librandi, Nov 23 2011
    
  • Mathematica
    LinearRecurrence[{8,4}, {0,1}, 50]
  • PARI
    x='x+O('x^30); concat([0], Vec(x/(1-8*x-4*x^2))) \\ G. C. Greubel, Jan 24 2018

Formula

G.f.: x/(1 - 8*x - 4*x^2). - R. J. Mathar, Nov 21 2011
a(n+1) = Sum_{k, 0<=k<=n} A099089(n,k)*4^k. - Philippe Deléham, Nov 21 2011

A190955 a(n) = 10*a(n-1) + 5*a(n-2), with a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 10, 105, 1100, 11525, 120750, 1265125, 13255000, 138875625, 1455031250, 15244690625, 159722062500, 1673444078125, 17533051093750, 183697731328125, 1924642568750000, 20164914344140625, 211272356285156250, 2213548134572265625, 23191843127148437500
Offset: 0

Author

Keywords

Programs

  • Magma
    I:=[0,1]; [n le 2 select I[n] else 10*Self(n-1) + 5*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 25 2018
  • Mathematica
    LinearRecurrence[{10,5}, {0,1}, 50]
  • PARI
    x='x+O('x^30); concat([0], Vec(x/(1-10*x-5*x^2))) \\ G. C. Greubel, Jan 25 2018
    

Formula

G.f.: x/(1 - 10*x - 5*x^2). - R. J. Mathar, Nov 21 2011
a(n+1) = Sum_{k=0..n} A099089(n,k)*5^k. - Philippe Deléham, Nov 21 2011

A099096 Riordan array (1,2-x).

Original entry on oeis.org

1, 0, 2, 0, -1, 4, 0, 0, -4, 8, 0, 0, 1, -12, 16, 0, 0, 0, 6, -32, 32, 0, 0, 0, -1, 24, -80, 64, 0, 0, 0, 0, -8, 80, -192, 128, 0, 0, 0, 0, 1, -40, 240, -448, 256, 0, 0, 0, 0, 0, 10, -160, 672, -1024, 512, 0, 0, 0, 0, 0, -1, 60, -560, 1792, -2304, 1024, 0, 0, 0, 0, 0
Offset: 0

Author

Paul Barry, Sep 25 2004

Keywords

Comments

Row sums are n+1 = Sum_{k=0..n} binomial(k,n-k)*2^(2k-n)*(-1)^(n-k). Diagonal sums are (-1)^n*A008346(n). The Riordan array (1,s+tx) defines T(n,k) = binomial(k,n-k)*s^k*(t/s)^(n-k). The row sums satisfy a(n) = s*a(n-1) + t*a(n-2) and the diagonal sums satisfy a(n) = s*a(n-2) + t*a(n-3).
Triangle T(n,k), 0 <= k <= n, read by rows given by [0, -1/2, 1/2, 0, 0, 0, 0, ...] DELTA [2, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 10 2008

Examples

			Rows begin
  1;
  0,   2;
  0,  -1,   4;
  0,   0,  -4,   8;
  0,   0,   1, -12,  16;
  ...
		

Crossrefs

Cf. A099089.

Programs

  • PARI
    /* Matrix power T^m formula: [T^m](n,k) = */ {T(n,k,m=1)=polcoeff((1 - (1-x +x*O(x^n))^(2^m) )^k,n)} \\ Paul D. Hanna, Nov 15 2007

Formula

Number triangle T(n, k) = binomial(k, n-k)*2^k*(-1/2)^(n-k); columns have g.f. (2x-x^2)^k.
G.f. of column k of matrix power T^m = (1 - (1-x)^(2^m))^k for k >= 0, when including the leading zeros that appear above the diagonal. - Paul D. Hanna, Nov 15 2007
T(n,k) = 2*T(n-1,k-1) - T(n-2,k-1), with T(0,0)=1, T(n,k)=0 if k < 0 or if k > n. - Philippe Deléham, Nov 25 2013
G.f.: 1/(1-2*x*y+x^2*y). - R. J. Mathar, Aug 12 2015
Showing 1-9 of 9 results.