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 43 results. Next

A123021 Triangle of coefficients of (1 - x)^n*B(x/(1 - x),n), where B(x,n) is the Morgan-Voyce polynomial related to A078812.

Original entry on oeis.org

1, 2, -1, 3, -2, 4, -2, -2, 1, 5, 0, -9, 6, -1, 6, 5, -24, 18, -4, 7, 14, -49, 36, -4, -4, 1, 8, 28, -84, 50, 20, -30, 10, -1, 9, 48, -126, 36, 115, -120, 45, -6, 10, 75, -168, -48, 358, -335, 120, -6, -6, 1, 11, 110, -198, -264, 847, -714, 175, 84, -63, 14
Offset: 0

Views

Author

Roger L. Bagula and Gary W. Adamson, Sep 24 2006

Keywords

Comments

The n-th row consists of the coefficients in the expansion of Sum_{j=0..n} A078812(n,j)*x^j*(1 - x)^(n - j).

Examples

			Triangle begins:
    1;
    2,  -1;
    3,  -2;
    4,  -2,   -2,    1;
    5,   0,   -9,    6,  -1;
    6,   5,  -24,   18,  -4;
    7,  14,  -49,   36,  -4,   -4,   1;
    8,  28,  -84,   50,  20,  -30,  10, -1;
    9,  48, -126,   36, 115, -120,  45, -6;
   10,  75, -168,  -48, 358, -335, 120, -6,  -6,  1;
   11, 110, -198, -264, 847, -714, 175, 84, -63, 14, -1;
   ... - _Franck Maminirina Ramaharo_, Oct 09 2018
		

Crossrefs

Programs

  • Mathematica
    Table[CoefficientList[Sum[Binomial[n+k+1, n-k]*x^k*(1-x)^(n-k), {k, 0, n}], x], {n, 0, 10}]//Flatten
  • Maxima
    t(n, k) := binomial(n + k + 1, n - k)$
    P(x, n) := expand(sum(t(n, j)*x^j*(1 - x)^(n - j), j, 0, n))$
    T(n, k) := ratcoef(P(x, n), x, k)$
    tabf(nn) := for n:0 thru nn do print(makelist(T(n, k), k, 0, hipow(P(x, n), x)))$ /* Franck Maminirina Ramaharo, Oct 09 2018 */
    
  • Sage
    def p(n,x): return sum( binomial(n+j+1, n-j)*x^j*(1-x)^(n-j) for j in (0..n) )
    def T(n): return ( p(n,x) ).full_simplify().coefficients(sparse=False)
    flatten([T(n) for n in (0..12)]) # G. C. Greubel, Jul 15 2021

Formula

From Franck Maminirina Ramaharo, Oct 09 2018: (Start)
Row n = coefficients in the expansion of (1/sqrt((4 - 3*x)*x))*(((2 - x + sqrt((4 - 3*x)*x))/2)^(n + 1) - ((2 - x - sqrt((4 - 3*x)*x))/2)^(n + 1)).
G.f.: 1/(1 - (2 - x)*y + (1 - x)^2*y^2).
E.g.f.: (1/sqrt((4 - 3*x)*x))*((2 - x + sqrt((4 - 3*x)*x))*exp(y*(2 - x + sqrt((4 - 3*x)*x))/2)/2 - (2 - x - sqrt((4 - 3*x)*x))*exp(y*(2 - x - sqrt((4 - 3*x)*x))/2)/2).
T(n,1) = -A254749(n+1). (End)

Extensions

Edited, new name, and offset corrected by Franck Maminirina Ramaharo, Oct 09 2018

A125274 Eigensequence of triangle A078812: a(n) = Sum_{k=0..n-1} A078812(n-1,k)*a(k) for n > 0 with a(0)=1.

Original entry on oeis.org

1, 1, 3, 10, 42, 210, 1199, 7670, 54224, 418744, 3499781, 31425207, 301324035, 3069644790, 33078375153, 375634524357, 4480492554993, 55971845014528, 730438139266281, 9935106417137098, 140553930403702487
Offset: 0

Views

Author

Paul D. Hanna, Nov 26 2006

Keywords

Examples

			a(3) = 3*(1) + 4*(1) + 1*(3) = 10;
a(4) = 4*(1) + 10*(1) + 6*(3) + 1*(10) = 42;
a(5) = 5*(1) + 20*(1) + 21*(3) + 8*(10) + 1*(42) = 210.
Triangle A078812(n,k) = binomial(n+k+1, n-k) begins:
  1;
  2,  1;
  3,  4,  1;
  4, 10,  6,  1;
  5, 20, 21,  8,  1;
  6, 35, 56, 36, 10,  1; ...
where g.f. of column k = 1/(1-x)^(2*k+2).
		

Crossrefs

Programs

  • Mathematica
    a[0] = 1; a[n_] := a[n] = Sum[Binomial[n+k, n-k-1] * a[k], {k, 0, n-1}]; Array[a, 20, 0] (* Amiram Eldar, Nov 24 2018 *)
  • PARI
    a(n)=if(n==0,1,sum(k=0,n-1, a(k)*binomial(n+k, n-k-1)))

Formula

a(n) = Sum_{k=0..n-1} binomial(n+k, n-k-1)*a(k) for n > 0 with a(0)=1.
G.f. satisfies A(x) = 1 + x/(1-x)^2*A(x/(1-x)^2). [Vladimir Kruchinin, Nov 28 2011]

A085643 Refines sequence A078812 using permutations and products of exponents on least prime signatures.

Original entry on oeis.org

1, 2, 1, 3, 4, 1, 4, 6, 4, 6, 1, 5, 8, 12, 9, 12, 8, 1, 6, 10, 16, 9, 12, 36, 8, 12, 24, 10, 1, 7, 12, 20, 24, 15, 48, 36, 27, 16, 72, 32, 15, 40, 12, 1, 8, 14, 24, 30, 16, 18, 60, 48, 72, 54, 20, 96, 144, 54, 16, 20, 120, 80, 18, 60, 14, 1
Offset: 1

Views

Author

Alford Arnold, Aug 15 2003

Keywords

Comments

The shape sequence for this table is A000041. Row sums in the example are A000079, A006906 and A001906.

Examples

			When the signatures (A025487) are listed in A036035 order the number of permutations of the exponents begin
1, 1 1, 1 2 1, 1 2 1 3 1, 1 2 2 3 3 4 1, ... b(n) and the products begin
1, 2 1, 3 2 1, 4 3 4 2 1, 5 4 6 3 4 2 1, ... c(n) thus we can write
1, 2 1, 3 4 1, 4 6 4 6 1, 5 8 12 9 12 8 1, ... a(n) = b(n)*c(n)
		

Crossrefs

Cf. A078812.

A144254 Eigentriangle by rows, termwise products of A078812 and its eigensequence, A125274.

Original entry on oeis.org

1, 2, 1, 3, 4, 3, 4, 10, 18, 10, 5, 20, 63, 80, 42, 6, 35, 168, 360, 420, 210, 7, 56, 378, 1200, 2310, 2520, 1199, 8, 84, 756, 3300, 9240, 16380, 16786, 7670, 9, 120, 1386, 7920, 30030, 76440, 125895, 122720, 54224, 10, 165, 2376, 17160, 84084, 286650
Offset: 1

Views

Author

Gary W. Adamson, Sep 16 2008

Keywords

Comments

Right border A144253 = A125274, the eigensequence of A078812: (1, 1, 3, 10, 42, 210, 1199,...).
Row sums = A125274 shifted.
Sum of row n terms = rightmost term of next row.

Examples

			First few rows of the triangle =
1;
2, 1;
3, 4, 3;
4, 10, 18, 10;
5, 20, 63, 80, 42;
6, 35, 168, 360, 420, 210;
7, 56, 378, 1200, 2310, 2520, 1199;
...
Triangle A078812 begins:
1;
2, 1;
3, 4, 1;
4, 10, 6, 1;
5, 20, 21, 8, 1;
...
Its eigensequence = A125274: (1, 1, 3, 10, 42, 210, 1199,...).
Row 3 of triangle A144253 = termwise products of (4, 10, 6, 1) and (1, 1, 3, 10) = (4*1, 10*1, 6*3, 1*10).
		

Crossrefs

Formula

Eigensequence by rows, T(n,k) = A078812(n,k) * A125274(k).

A049310 Triangle of coefficients of Chebyshev's S(n,x) := U(n,x/2) polynomials (exponents in increasing order).

Original entry on oeis.org

1, 0, 1, -1, 0, 1, 0, -2, 0, 1, 1, 0, -3, 0, 1, 0, 3, 0, -4, 0, 1, -1, 0, 6, 0, -5, 0, 1, 0, -4, 0, 10, 0, -6, 0, 1, 1, 0, -10, 0, 15, 0, -7, 0, 1, 0, 5, 0, -20, 0, 21, 0, -8, 0, 1, -1, 0, 15, 0, -35, 0, 28, 0, -9, 0, 1, 0, -6, 0, 35, 0, -56, 0, 36, 0, -10, 0, 1, 1, 0, -21, 0, 70, 0, -84, 0
Offset: 0

Views

Author

Keywords

Comments

G.f. for row polynomials S(n,x) (signed triangle): 1/(1-x*z+z^2). Unsigned triangle |a(n,m)| has Fibonacci polynomials F(n+1,x) as row polynomials with g.f. 1/(1-x*z-z^2). |a(n,m)| triangle has rows of Pascal's triangle A007318 in the even-numbered diagonals (odd-numbered ones have only 0's).
Row sums (unsigned triangle) A000045(n+1) (Fibonacci). Row sums (signed triangle) S(n,1) sequence = periodic(1,1,0,-1,-1,0) = A010892.
Alternating row sums A049347(n) = S(n,-1) = periodic(1,-1,0). - Wolfdieter Lang, Nov 04 2011
S(n,x) is the characteristic polynomial of the adjacency matrix of the n-path. - Michael Somos, Jun 24 2002
S(n,x) is also the matching polynomial of the n-path. - Eric W. Weisstein, Apr 10 2017
|T(n,k)| = number of compositions of n+1 into k+1 odd parts. Example: |T(7,3)| = 10 because we have (1,1,3,3), (1,3,1,3), (1,3,3,1), (3,1,1,3), (3,1,3,1), (3,3,1,1), (1,1,1,5), (1,1,5,1), (1,5,1,1) and (5,1,1,1). - Emeric Deutsch, Apr 09 2005
S(n,x)= R(n,x) + S(n-2,x), n >= 2, S(-1,x)=0, S(0,x)=1, R(n,x):=2*T(n,x/2) = Sum_{m=0..n} A127672(n,m)*x^m (monic integer Chebyshev T-Polynomials). This is the rewritten so-called trace of the transfer matrix formula for the T-polynomials. - Wolfdieter Lang, Dec 02 2010
In a regular N-gon inscribed in a unit circle, the side length is d(N,1) = 2*sin(Pi/N). The length ratio R(N,k):=d(N,k)/d(N,1) for the (k-1)-th diagonal, with k from {2,3,...,floor(N/2)}, N >= 4, equals S(k-1,x) = sin(k*Pi/N)/sin(Pi/N) with x=rho(N):=R(N,2) = 2*cos(Pi/N). Example: N=7 (heptagon), rho=R(7,2), sigma:=R(N,3) = S(2,rho) = rho^2 - 1. Motivated by the quoted paper by P. Steinbach. - Wolfdieter Lang, Dec 02 2010
From Wolfdieter Lang, Jul 12 2011: (Start)
In q- or basic analysis, q-numbers are [n]_q := S(n-1,q+1/q) = (q^n-(1/q)^n)/(q-1/q), with the row polynomials S(n,x), n >= 0.
The zeros of the row polynomials S(n-1,x) are (from those of Chebyshev U-polynomials):
x(n-1;k) = +- t(k,rho(n)), k = 1..ceiling((n-1)/2), n >= 2, with t(n,x) the row polynomials of A127672 and rho(n):= 2*cos(Pi/n). The simple vanishing zero for even n appears here as +0 and -0.
Factorization of the row polynomials S(n-1,x), x >= 1, in terms of the minimal polynomials of cos(2 Pi/2), called Psi(n,x), with coefficients given by A181875/A181876:
S(n-1,x) = (2^(n-1))*Product_{n>=1}(Psi(d,x/2), 2 < d | 2n).
(From the rewritten eq. (3) of the Watkins and Zeitlin reference, given under A181872.) [See the W. Lang ArXiv link, Proposition 9, eq. (62). - Wolfdieter Lang, Apr 14 2018]
(End)
The discriminants of the S(n,x) polynomials are found in A127670. - Wolfdieter Lang, Aug 03 2011
This is an example for a subclass of Riordan convolution arrays (lower triangular matrices) called Bell arrays. See the L. W. Shapiro et al. reference under A007318. If a Riordan array is named (G(z),F(z)) with F(z)=z*Fhat(z), the o.g.f. for the row polynomials is G(z)/(1-x*z*Fhat(z)), and it becomes a Bell array if G(z)=Fhat(z). For the present Bell type triangle G(z)=1/(1+z^2) (see the o.g.f. comment above). This leads to the o.g.f. for the column no. k, k >= 0, x^k/(1+x^2)^(k+1) (see the formula section), the one for the row sums and for the alternating row sums (see comments above). The Riordan (Bell) A- and Z-sequences (defined in a W. Lang link under A006232, with references) have o.g.f.s 1-x*c(x^2) and -x*c(x^2), with the o.g.f. of the Catalan numbers A000108. Together they lead to a recurrence given in the formula section. - Wolfdieter Lang, Nov 04 2011
The determinant of the N x N matrix S(N,[x[1], ..., x[N]]) with elements S(m-1,x[n]), for n, m = 1, 2, ..., N, and for any x[n], is identical with the determinant of V(N,[x[1], ..., x[N]]) with elements x[n]^(m-1) (a Vandermondian, which equals Product_{1 <= i < j<= N} (x[j] - x[i])). This is a special instance of a theorem valid for any N >= 1 and any monic polynomial system p(m,x), m>=0, with p(0,x) = 1. For this theorem see the Vein-Dale reference, p. 59. Thanks to L. Edson Jeffery for an email asking for a proof of the non-singularity of the matrix S(N,[x[1], ...., x[N]]) if and only if the x[j], j = 1..N, are pairwise distinct. - Wolfdieter Lang, Aug 26 2013
These S polynomials also appear in the context of modular forms. The rescaled Hecke operator T*n = n^((1-k)/2)*T_n acting on modular forms of weight k satisfies T*(p^n) = S(n, T*p), for each prime p and positive integer n. See the Koecher-Krieg reference, p. 223. - _Wolfdieter Lang, Jan 22 2016
For a shifted o.g.f. (mod signs), its compositional inverse, and connections to Motzkin and Fibonacci polynomials, non-crossing partitions and other combinatorial structures, see A097610. - Tom Copeland, Jan 23 2016
From M. Sinan Kul, Jan 30 2016; edited by Wolfdieter Lang, Jan 31 2016 and Feb 01 2016: (Start)
Solutions of the Diophantine equation u^2 + v^2 - k*u*v = 1 for integer k given by (u(k,n), v(k,n)) = (S(n,k), S(n-1,k)) because of the Cassini-Simson identity: S(n,x)^2 - S(n+1,x)*S(n-1, x) = 1, after use of the S-recurrence. Note that S(-n, x) = -S(-n-2, x), n >= 1, and the periodicity of some S(n, k) sequences.
Hence another way to obtain the row polynomials would be to take powers of the matrix [x, -1; 1,0]: S(n, x) = (([x, -1; 1, 0])^n)[1,1], n >= 0.
See also a Feb 01 2016 comment on A115139 for a well-known S(n, x) sum formula.
Then we have with the present T triangle
A039834(n) = -i^(n+1)*T(n-1, k) where i is the imaginary unit and n >= 0.
A051286(n) = Sum_{i=0..n} T(n,i)^2 (see the Philippe Deléham, Nov 21 2005 formula),
A181545(n) = Sum_{i=0..n+1} abs(T(n,i)^3),
A181546(n) = Sum_{i=0..n+1} T(n,i)^4,
A181547(n) = Sum_{i=0..n+1} abs(T(n,i)^5).
S(n, 0) = A056594(n), and for k = 1..10 the sequences S(n-1, k) with offset n = 0 are A128834, A001477, A001906, A001353, A004254, A001109, A004187, A001090, A018913, A004189.
(End)
For more on the Diophantine equation presented by Kul, see the Ismail paper. - Tom Copeland, Jan 31 2016
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 U(n,x), A053117, so Sum_{k=0..n} L(k,x/2) L(n-k,x/2) = S(n,x). This gives S(n,x) = L(n/2,x/2)^2 + 2*Sum_{k=0..n/2-1} L(k,x/2) L(n-k,x/2) for n even and S(n,x) = 2*Sum_{k=0..(n-1)/2} L(k,x/2) L(n-k,x/2) for odd n. For a connection to elliptic curves and modular forms, see A053117. For the normalized Legendre polynomials, see A100258. For other properties and relations to other polynomials, see Allouche et al. - Tom Copeland, Feb 04 2016
LG(x,h1,h2) = -log(1 - h1*x + h2*x^2) = Sum_{n>0} F(n,-h1,h2,0,..,0) x^n/n is a log series generator of the bivariate row polynomials of A127672 with A127672(0,0) = 0 and where F(n,b1,b2,..,bn) are the Faber polynomials of A263916. Exp(LG(x,h1,h2)) = 1 / (1 - h1*x + h2*x^2 ) is the o.g.f. of the bivariate row polynomials of this entry. - Tom Copeland, Feb 15 2016 (Instances of the bivariate o.g.f. for this entry are on pp. 5 and 18 of Sunada. - Tom Copeland, Jan 18 2021)
For distinct odd primes p and q the Legendre symbol can be written as Legendre(q,p) = Product_{k=1..P} S(q-1, 2*cos(2*Pi*k/p)), with P = (p-1)/2. See the Lemmermeyer reference, eq. (8.1) on p. 236. Using the zeros of S(q-1, x) (see above) one has S(q-1, x) = Product_{l=1..Q} (x^2 - (2*cos(Pi*l/q))^2), with Q = (q-1)/2. Thus S(q-1, 2*cos(2*Pi*k/p)) = ((-4)^Q)*Product_{l=1..Q} (sin^2(2*Pi*k/p) - sin^2(Pi*l/q)) = ((-4)^Q)*Product_{m=1..Q} (sin^2(2*Pi*k/p) - sin^2(2*Pi*m/q)). For the proof of the last equality see a W. Lang comment on the triangle A057059 for n = Q and an obvious function f. This leads to Eisenstein's proof of the quadratic reciprocity law Legendre(q,p) = ((-1)^(P*Q)) * Legendre(p,q), See the Lemmermeyer reference, pp. 236-237. - Wolfdieter Lang, Aug 28 2016
For connections to generalized Fibonacci polynomials, compare their generating function on p. 5 of the Amdeberhan et al. link with the o.g.f. given above for the bivariate row polynomials of this entry. - Tom Copeland, Jan 08 2017
The formula for Ramanujan's tau function (see A000594) for prime powers is tau(p^k) = p^(11*k/2)*S(k, p^(-11/2)*tau(p)) for k >= 1, and p = A000040(n), n >= 1. See the Hardy reference, p. 164, eqs. (10.3.4) and (10.3.6) rewritten in terms of S. - Wolfdieter Lang, Jan 27 2017
From Wolfdieter Lang, May 08 2017: (Start)
The number of zeros Z(n) of the S(n, x) polynomials in the open interval (-1,+1) is 2*b(n) for even n >= 0 and 1 + 2*b(n) for odd n >= 1, where b(n) = floor(n/2) - floor((n+1)/3). This b(n) is the number of integers k in the interval (n+1)/3 < k <= floor(n/2). See a comment on the zeros of S(n, x) above, and b(n) = A008615(n-2), n >= 0. The numbers Z(n) have been proposed (with a conjecture related to A008611) by Michel Lagneau, as the number of zeros of Fibonacci polynomials on the imaginary axis (-I,+I), with I=sqrt(-1). They are Z(n) = A008611(n-1), n >= 0, with A008611(-1) = 0. Also Z(n) = A194960(n-4), n >= 0. Proof using the A008611 version. A194960 follows from this.
In general the number of zeros Z(a;n) of S(n, x) for n >= 0 in the open interval (-a,+a) for a from the interval (0,2) (x >= 2 never has zeros, and a=0 is trivial: Z(0;n) = 0) is with b(a;n) = floor(n//2) - floor((n+1)*arccos(a/2)/Pi), as above Z(a;n) = 2*b(a;n) for even n >= 0 and 1 + 2*b(a;n) for odd n >= 1. For the closed interval [-a,+a] Z(0;n) = 1 and for a from (0,1) one uses for Z(a;n) the values b(a;n) = floor(n/2) - ceiling((n+1)*arccos(a/2)/Pi) + 1. (End)
The Riordan row polynomials S(n, x) (Chebyshev S) belong to the Boas-Buck class (see a comment and references in A046521), hence they satisfy the Boas-Buck identity: (E_x - n*1)*S(n, x) = (E_x + 1)*Sum_{p=0..n-1} (1 - (-1)^p)*(-1)^((p+1)/2)*S(n-1-p, x), for n >= 0, where E_x = x*d/dx (Euler operator). For the triangle T(n, k) this entails a recurrence for the sequence of column k, given in the formula section. - Wolfdieter Lang, Aug 11 2017
The e.g.f. E(x,t) := Sum_{n>=0} (t^n/n!)*S(n,x) for the row polynomials is obtained via inverse Laplace transformation from the above given o.g.f. as E(x,t) = ((1/xm)*exp(t/xm) - (1/xp)*exp(t/xp) )/(xp - xm) with xp = (x + sqrt(x^2-4))/2 and xm = (x - sqrt(x^2-4))/2. - Wolfdieter Lang, Nov 08 2017
From Wolfdieter Lang, Apr 12 2018: (Start)
Factorization of row polynomials S(n, x), for n >= 1, in terms of C polynomials (not Chebyshev C) with coefficients given in A187360. This is obtained from the factorization into Psi polynomials (see the Jul 12 2011 comment above) but written in terms of minimal polynomials of 2*cos(2*Pi/n) with coefficients in A232624:
S(2*k, x) = Product_{2 <= d | (2*k+1)} C(d, x)*(-1)^deg(d)*C(d, -x), with deg(d) = A055034(d) the degree of C(d, x).
S(2*k+1, x) = Product_{2 <= d | 2*(k+1)} C(d, x) * Product_{3 <= 2*d + 1 | (k+1)} (-1)^(deg(2*d+1))*C(2*d+1, -x).
Note that (-1)^(deg(2*d+1))*C(2*d+1, -x)*C(2*d+1, x) pairs always appear.
The number of C factors of S(2*k, x), for k >= 0, is 2*(tau(2*k+1) - 1) = 2*(A099774(k+1) - 1) = 2*A095374(k), and for S(2*k+1, x), for k >= 0, it is tau(2*(k+1)) + tau_{odd}(k+1) - 2 = A302707(k), with tau(2*k+1) = A099774(k+1), tau(n) = A000005 and tau(2*(k+1)) = A099777(k+1).
For the reverse problem, the factorization of C polynomials into S polynomials, see A255237. (End)
The S polynomials with general initial conditions S(a,b;n,x) = x*S(a,b;n-1,x) - S(a,b;n-2,x), for n >= 1, with S(a,b;-1,x) = a and S(a,b;0,x) = b are S(a,b;n,x) = b*S(n, x) - a*S(n-1, x), for n >= -1. Recall that S(-2, x) = -1 and S(-1, x) = 0. The o.g.f. is G(a,b;z,x) = (b - a*z)/(1 - x*z + z^2). - Wolfdieter Lang, Oct 18 2019
Also the convolution triangle of A101455. - Peter Luschny, Oct 06 2022
From Wolfdieter Lang, Apr 26 2023: (Start)
Multi-section of S-polynomials: S(m*n+k, x) = S(m+k, x)*S(n-1, R(m, x)) - S(k, x)*S(n-2, R(m, x)), with R(n, x) = S(n, x) - S(n-2, x) (see A127672), S(-2, x) = -1, and S(-1, x) = 0, for n >= 0, m >= 1, and k = 0, 1, ..., m-1.
O.g.f. of {S(m*n+k, y)}_{n>=0}: G(m,k,y,x) = (S(k, y) - (S(k, y)*R(m, y) - S(m+k, y))*x)/(1 - R(m,y)*x + x^2).
See eqs. (40) and (49), with r = x or y and s =-1, of the G. Detlefs and W. Lang link at A034807. (End)
S(n, x) for complex n and complex x: S(n, x) = ((-i/2)/sqrt(1 - (x/2)^2))*(q(x/2)*exp(+n*log(q(x/2))) - (1/q(x/2))*exp(-n*log(q(x/2)))), with q(x) = x + sqrt(1 - x^2)*i. Here log(z) = |z| + Arg(z)*i, with Arg(z) from [-Pi,+Pi) (principal branch). This satisfies the recurrence relation for S because it is derived from the Binet - de Moivre formula for S. Examples: S(n/m, 0) = cos((n/m)*Pi/4), for n >= 0 and m >= 1. S(n*i, 0) = (1/2)*(1 + exp(n*Pi))*exp(-(n/2)*Pi), for n >= 0. S(1+i, 2+i) = 0.6397424847... + 1.0355669490...*i. Thanks to Roberto Alfano for asking a question leading to this formula. - Wolfdieter Lang, Jun 05 2023
Lim_{n->oo} S(n, x)/S(n-1, x) = r(x) = (x - sqrt(x^2 -4))/2, for |x| >= 2. For x = +-2, this limit is +-1. - Wolfdieter Lang, Nov 15 2023

Examples

			The triangle T(n, k) begins:
  n\k  0  1   2   3   4   5   6    7   8   9  10  11
  0:   1
  1:   0  1
  2:  -1  0   1
  3:   0 -2   0   1
  4:   1  0  -3   0   1
  5:   0  3   0  -4   0   1
  6:  -1  0   6   0  -5   0   1
  7:   0 -4   0  10   0  -6   0    1
  8:   1  0 -10   0  15   0  -7    0   1
  9:   0  5   0 -20   0  21   0   -8   0   1
  10: -1  0  15   0 -35   0  28    0  -9   0   1
  11:  0 -6   0  35   0 -56   0   36   0 -10   0   1
  ... Reformatted and extended by _Wolfdieter Lang_, Oct 24 2012
For more rows see the link.
E.g., fourth row {0,-2,0,1} corresponds to polynomial S(3,x)= -2*x + x^3.
From _Wolfdieter Lang_, Jul 12 2011: (Start)
Zeros of S(3,x) with rho(4)= 2*cos(Pi/4) = sqrt(2):
  +- t(1,sqrt(2)) = +- sqrt(2) and
  +- t(2,sqrt(2)) = +- 0.
Factorization of S(3,x) in terms of Psi polynomials:
S(3,x) = (2^3)*Psi(4,x/2)*Psi(8,x/2) = x*(x^2-2).
(End)
From _Wolfdieter Lang_, Nov 04 2011: (Start)
A- and Z- sequence recurrence:
T(4,0) = - (C(0)*T(3,1) + C(1)*T(3,3)) = -(-2 + 1) = +1,
T(5,3) = -3 - 1*1 = -4.
(End)
Boas-Buck recurrence for column k = 2, n = 6: S(6, 2) = (3/4)*(0 - 2* S(4 ,2) + 0 + 2*S(2, 2)) = (3/4)*(-2*(-3) + 2) = 6. - _Wolfdieter Lang_, Aug 11 2017
From _Wolfdieter Lang_, Apr 12 2018: (Start)
Factorization into C polynomials (see the Apr 12 2018 comment):
S(4, x) = 1 - 3*x^2 + x^4 = (-1 + x + x^2)*(-1 - x + x^2) = (-C(5, -x)) * C(5, x); the number of factors is 2 = 2*A095374(2).
S(5, x) = 3*x - 4*x^3 + x^5 = x*(-1 + x)*(1 + x)*(-3 + x^2) = C(2, x)*C(3, x)*(-C(3, -x))*C(6, x); the number of factors is 4 = A302707(2). (End)
		

References

  • G. H. Hardy, Ramanujan: twelve lectures on subjects suggested by his life and work, AMS Chelsea Publishing, Providence, Rhode Island, 2002, p. 164.
  • Max Koecher and Aloys Krieg, Elliptische Funktionen und Modulformen, 2. Auflage, Springer, 2007, p. 223.
  • Franz Lemmermeyer, Reciprocity Laws. From Euler to Eisenstein, Springer, 2000.
  • D. S. Mitrinovic, Analytic Inequalities, Springer-Verlag, 1970; p. 232, Sect. 3.3.38.
  • Theodore J. Rivlin, Chebyshev polynomials: from approximation theory to algebra and number theory, 2. ed., Wiley, New York, 1990, pp. 60 - 61.
  • R. Vein and P. Dale, Determinants and Their Applications in Mathematical Physics, Springer, 1999.

Crossrefs

Cf. A000005, A000217, A000292, A000332, A000389, A001227, A007318, A008611, A008615, A101455, A010892, A011973, A053112 (without zeros), A053117, A053119 (reflection), A053121 (inverse triangle), A055034, A097610, A099774, A099777, A100258, A112552 (first column clipped), A127672, A168561 (absolute values), A187360. A194960, A232624, A255237.
Triangles of coefficients of Chebyshev's S(n,x+k) for k = 5, 4, 3, 2, 1, 0, -1, -2, -3, -4, -5: A207824, A207823, A125662, A078812, A101950, A049310, A104562, A053122, A207815, A159764, A123967.

Programs

  • Magma
    A049310:= func< n,k | ((n+k) mod 2) eq 0 select (-1)^(Floor((n+k)/2)+k)*Binomial(Floor((n+k)/2), k) else 0 >;
    [A049310(n,k): k in [0..n], n in [0..15]]; // G. C. Greubel, Jul 25 2022
  • Maple
    A049310 := proc(n,k): binomial((n+k)/2,(n-k)/2)*cos(Pi*(n-k)/2)*(1+(-1)^(n-k))/2 end: seq(seq(A049310(n,k), k=0..n),n=0..11); # Johannes W. Meijer, Aug 08 2011
    # Uses function PMatrix from A357368. Adds a row above and a column to the left.
    PMatrix(10, n -> ifelse(irem(n, 2) = 0, 0, (-1)^iquo(n-1, 2))); # Peter Luschny, Oct 06 2022
  • Mathematica
    t[n_, k_] /; EvenQ[n+k] = ((-1)^((n+k)/2+k))*Binomial[(n+k)/2, k]; t[n_, k_] /; OddQ[n+k] = 0; Flatten[Table[t[n, k], {n, 0, 12}, {k, 0, n}]][[;; 86]] (* Jean-François Alcover, Jul 05 2011 *)
    Table[Coefficient[(-I)^n Fibonacci[n + 1, - I x], x, k], {n, 0, 10}, {k, 0, n}] //Flatten (* Clark Kimberling, Aug 02 2011; corrected by Eric W. Weisstein, Apr 06 2017 *)
    CoefficientList[ChebyshevU[Range[0, 10], -x/2], x] // Flatten (* Eric W. Weisstein, Apr 06 2017 *)
    CoefficientList[Table[(-I)^n Fibonacci[n + 1, -I x], {n, 0, 10}], x] // Flatten (* Eric W. Weisstein, Apr 06 2017 *)
  • PARI
    {T(n, k) = if( k<0 || k>n || (n + k)%2, 0, (-1)^((n + k)/2 + k) * binomial((n + k)/2, k))} /* Michael Somos, Jun 24 2002 */
    
  • SageMath
    @CachedFunction
    def A049310(n,k):
        if n< 0: return 0
        if n==0: return 1 if k == 0 else 0
        return A049310(n-1,k-1) - A049310(n-2,k)
    for n in (0..9): [A049310(n,k) for k in (0..n)] # Peter Luschny, Nov 20 2012
    

Formula

T(n,k) := 0 if n < k or n+k odd, otherwise ((-1)^((n+k)/2+k))*binomial((n+k)/2, k); T(n, k) = -T(n-2, k)+T(n-1, k-1), T(n, -1) := 0 =: T(-1, k), T(0, 0)=1, T(n, k)= 0 if n < k or n+k odd; g.f. k-th column: (1 / (1 + x^2)^(k + 1)) * x^k. - Michael Somos, Jun 24 2002
T(n,k) = binomial((n+k)/2, (n-k)/2)*cos(Pi*(n-k)/2)*(1+(-1)^(n-k))/2. - Paul Barry, Aug 28 2005
Sum_{k=0..n} T(n,k)^2 = A051286(n). - Philippe Deléham, Nov 21 2005
Recurrence for the (unsigned) Fibonacci polynomials: F(1)=1, F(2)=x; for n > 2, F(n) = x*F(n-1) + F(n-2).
From Wolfdieter Lang, Nov 04 2011: (Start)
The Riordan A- and Z-sequences, given in a comment above, lead together to the recurrence:
T(n,k) = 0 if n < k, if k=0 then T(0,0)=1 and
T(n,0)= -Sum_{i=0..floor((n-1)/2)} C(i)*T(n-1,2*i+1), otherwise T(n,k) = T(n-1,k-1) - Sum_{i=1..floor((n-k)/2)} C(i)*T(n-1,k-1+2*i), with the Catalan numbers C(n)=A000108(n).
(End)
The row polynomials satisfy also S(n,x) = 2*(T(n+2, x/2) - T(n, x/2))/(x^2-4) with the Chebyshev T-polynomials. Proof: Use the trace formula 2*T(n, x/2) = S(n, x) - S(n-2, x) (see the Dec 02 2010 comment above) and the S-recurrence several times. This is a formula which expresses the S- in terms of the T-polynomials. - Wolfdieter Lang, Aug 07 2014
From Tom Copeland, Dec 06 2015: (Start)
The non-vanishing, unsigned subdiagonals Diag_(2n) contain the elements D(n,k) = Sum_{j=0..k} D(n-1,j) = (k+1) (k+2) ... (k+n) / n! = binomial(n+k,n), so the o.g.f. for the subdiagonal is (1-x)^(-(n+1)). E.g., Diag_4 contains D(2,3) = D(1,0) + D(1,1) + D(1,2) + D(1,3) = 1 + 2 + 3 + 4 = 10 = binomial(5,2). Diag_4 is shifted A000217; Diag_6, shifted A000292: Diag_8, shifted A000332; and Diag_10, A000389.
The non-vanishing antidiagonals are signed rows of the Pascal triangle A007318.
For a reversed, unsigned version with the zeros removed, see A011973. (End)
The Boas-Buck recurrence (see a comment above) for the sequence of column k is: S(n, k) = ((k+1)/(n-k))*Sum_{p=0..n-1-k} (1 - (-1)^p)*(-1)^((p+1)/2) * S(n-1-p, k), for n > k >= 0 and input S(k, k) = 1. - Wolfdieter Lang, Aug 11 2017
The m-th row consecutive nonzero entries in order are (-1)^c*(c+b)!/c!b! with c = m/2, m/2-1, ..., 0 and b = m-2c if m is even and with c = (m-1)/2, (m-1)/2-1, ..., 0 with b = m-2c if m is odd. For the 8th row starting at a(36) the 5 consecutive nonzero entries in order are 1,-10,15,-7,1 given by c = 4,3,2,1,0 and b = 0,2,4,6,8. - Richard Turk, Aug 20 2017
O.g.f.: exp( Sum_{n >= 0} 2*T(n,x/2)*t^n/n ) = 1 + x*t + (-1 + x^2)*t^2 + (-2*x + x^3)*t^3 + (1 - 3*x^2 + x^4)*t^4 + ..., where T(n,x) denotes the n-th Chebyshev polynomial of the first kind. - Peter Bala, Aug 15 2022

A001906 F(2n) = bisection of Fibonacci sequence: a(n) = 3*a(n-1) - a(n-2).

Original entry on oeis.org

0, 1, 3, 8, 21, 55, 144, 377, 987, 2584, 6765, 17711, 46368, 121393, 317811, 832040, 2178309, 5702887, 14930352, 39088169, 102334155, 267914296, 701408733, 1836311903, 4807526976, 12586269025, 32951280099, 86267571272, 225851433717, 591286729879, 1548008755920
Offset: 0

Views

Author

Keywords

Comments

Apart from initial term, same as A088305.
Second column of array A102310 and of A028412.
Numbers k such that 5*k^2 + 4 is a square. - Gregory V. Richardson, Oct 13 2002
Apart from initial terms, also Pisot sequences E(3,8), P(3,8), T(3,8). See A008776 for definitions of Pisot sequences.
Binomial transform of A000045. - Paul Barry, Apr 11 2003
Number of walks of length 2n+1 in the path graph P_4 from one end to the other one. Example: a(2)=3 because in the path ABCD we have ABABCD, ABCBCD and ABCDCD. - Emeric Deutsch, Apr 02 2004
Simplest example of a second-order recurrence with the sixth term a square.
Number of (s(0), s(1), ..., s(2n)) such that 0 < s(i) < 5 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n, s(0) = 1, s(2n) = 3. - Lekraj Beedassy, Jun 11 2004
a(n) (for n > 0) is the smallest positive integer that cannot be created by summing at most n values chosen among the previous terms (with repeats allowed). - Andrew Weimholt, Jul 20 2004
All nonnegative integer solutions of Pell equation b(n)^2 - 5*a(n)^2 = +4 together with b(n) = A005248(n), n >= 0. - Wolfdieter Lang, Aug 31 2004
a(n+1) is a Chebyshev transform of 3^n (A000244), where the sequence with g.f. G(x) is sent to the sequence with g.f. (1/(1+x^2))G(x/(1+x^2)). - Paul Barry, Oct 25 2004
a(n) is the number of distinct products of matrices A, B, C, in (A+B+C)^n where commutator [A,B] = 0 but C does not commute with A or B. - Paul D. Hanna and Max Alekseyev, Feb 01 2006
Number of binary words with exactly k-1 strictly increasing runs. Example: a(3)=F(6)=8 because we have 0|0,1|0,1|1,0|01,01|0,1|01,01|1 and 01|01. Column sums of A119900. - Emeric Deutsch, Jul 23 2006
See Table 1 on page 411 of Lukovits and Janezic paper. - Parthasarathy Nambi, Aug 22 2006
Inverse: With phi = (sqrt(5) + 1)/2, log_phi((sqrt(5) a(n) + sqrt(5 a(n)^2 + 4))/2) = n. - David W. Cantrell (DWCantrell(AT)sigmaxi.net), Feb 19 2007
[1,3,8,21,55,144,...] is the Hankel transform of [1,1,4,17,75,339,1558,...](see A026378). - Philippe Deléham, Apr 13 2007
The Diophantine equation a(n) = m has a solution (for m >= 1) if and only if floor(arcsinh(sqrt(5)*m/2)/log(phi)) <> floor(arccosh(sqrt(5)*m/2)/log(phi)) where phi is the golden ratio. An equivalent condition is A130259(m) = A130260(m). - Hieronymus Fischer, May 25 2007
a(n+1) = AB^(n)(1), n >= 0, with compositions of Wythoff's complementary A(n):=A000201(n) and B(n)=A001950(n) sequences. See the W. Lang link under A135817 for the Wythoff representation of numbers (with A as 1 and B as 0 and the argument 1 omitted). E.g., 1=`1`, 3=`10`, 8=`100`, 21=`1000`, ..., in Wythoff code.
Equals row sums of triangles A140069, A140736 and A140737. - Gary W. Adamson, May 25 2008
a(n) is also the number of idempotent order-preserving partial transformations (of an n-element chain) of width n (width(alpha) = max(Im(alpha))). Equivalently, it is the number of idempotent order-preserving full transformations (of an n-element chain). - Abdullahi Umar, Sep 08 2008
a(n) is the number of ways that a string of 0,1 and 2 of size (n-1) can be arranged with no 12-pairs. - Udita Katugampola, Sep 24 2008
Starting with offset 1 = row sums of triangle A175011. - Gary W. Adamson, Apr 03 2010
As a fraction: 1/71 = 0.01408450... or 1/9701 = 0.0001030821.... - Mark Dols, May 18 2010
Sum of the products of the elements in the compositions of n (example for n=3: the compositions are 1+1+1, 1+2, 2+1, and 3; a(3) = 1*1*1 + 1*2 + 2*1 + 3 = 8). - Dylon Hamilton, Jun 20 2010, Geoffrey Critzer, Joerg Arndt, Dec 06 2010
a(n) relates to regular polygons with even numbers of edges such that Product_{k=1..(n-2)/2} (1 + 4*cos^2 k*Pi/n) = even-indexed Fibonacci numbers with a(n) relating to the 2*n-gons. The constants as products = roots to even-indexed rows of triangle A152063. For example: a(5) = 55 satisfies the product formula relating to the 10-gon. - Gary W. Adamson, Aug 15 2010
Alternatively, product of roots to x^4 - 12x^3 + 51x^2 - 90x + 55, (10th row of triangle A152063) = (4.618...)*(3.618...)*(2.381...)*(1.381...) = 55. - Gary W. Adamson, Aug 15 2010
a(n) is the number of generalized compositions of n when there are i different types of i, (i=1,2,...). - Milan Janjic, Aug 26 2010
Starting with "1" = row sums of triangle A180339, and eigensequence of triangle A137710. - Gary W. Adamson, Aug 28 2010
a(2) = 3 is the only prime.
Number of nonisomorphic graded posets with 0 and uniform hasse graph of rank n > 0, with exactly 2 elements of each rank level above 0. (Uniform used in the sense of Retakh, Serconek, and Wilson. Graded used in Stanley's sense that every maximal chain has the same length n.) - David Nacin, Feb 13 2012
Pisano period lengths: 1, 3, 4, 3, 10, 12, 8, 6, 12, 30, 5, 12, 14, 24, 20, 12, 18, 12, 9, 30, ... - R. J. Mathar, Aug 10 2012
Solutions (x, y) = (a(n), a(n+1)) satisfying x^2 + y^2 = 3xy + 1. - Michel Lagneau, Feb 01 2014
For n >= 1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,2}. - Milan Janjic, Jan 25 2015
With a(0) = 0, for n > 1, a(n) is the smallest number not already in the sequence such that a(n)^2 - a(n-1)^2 is a Fibonacci number. - Derek Orr, Jun 08 2015
Let T be the tree generated by these rules: 0 is in T, and if p is in T, then p + 1 is in T and x*p is in T and y*p is in T. The n-th generation of T consists of A001906(n) polynomials, for n >= 0. - Clark Kimberling, Nov 24 2015
For n > 0, a(n) = exactly the maximum area of a quadrilateral with sides in order of lengths F(n), F(n), L(n), and L(n) with L(n)=A000032(n). - J. M. Bergot, Jan 20 2016
a(n) = twice the area of a triangle with vertices at (L(n+1), L(n+2)), (F(n+1), F(n+1)), and (L(n+2), L(n+1)), with L(n)=A000032(n). - J. M. Bergot, Apr 20 2016
Except for the initial 0, this is the p-INVERT of (1,1,1,1,1,...) for p(S) = 1 - S - S^2; see A291000. - Clark Kimberling, Aug 24 2017
a(n+1) is the number of spanning trees of the graph T_n, where T_n is a sequence of n triangles, where adjacent triangles share an edge. - Kevin Long, May 07 2018
a(n) is the number of ways to partition [n] such that each block is a run of consecutive numbers, and each block has a fixed point, e.g., for n=3, 12|3 with 1 and 3 as fixed points is valid, but 13|2 is not valid as 1 and 3 do not form a run. Consequently, a(n) also counts the spanning trees of the graph given by taking a path with n vertices and adding another vertex adjacent to all of them. - Kevin Long, May 11 2018
From Wolfdieter Lang, May 31 2018: (Start)
The preceding comment can be paraphrased as follows. a(n) is the row sum of the array A305309 for n >= 1. The array A305309(n, k) gives the sum of the products of the block lengths of the set partition of [n] := {1, 2, ..., n} with A048996(n, k) blocks of consecutive numbers, corresponding to the compositions obtained from the k-th partition of n in Abramowitz-Stegun order. See the comments and examples at A305309.
{a(n)} also gives the infinite sequence of nonnegative numbers k for which k * ||k*phi|| < 1/sqrt(5), where the irrational number phi = A001622 (golden section), and ||x|| is the absolute value of the difference between x and the nearest integer. See, e.g., the Havil reference, pp. 171-172. (End)
a(n) is the number of tilings of two n X 1 rectangles joined orthogonally at a common end-square (so to have 2n-1 squares in a right-angle V shape) with only 1 X 1 and 2 X 1 tiles. This is a consequence of F(2n) = F(n+1)*F(n) + F(n)*F(n-1). - Nathaniel Gregg, Oct 10 2021
These are the denominators of the upper convergents to the golden ratio, tau; they are also the numerators of the lower convergents (viz. 1/1 < 3/2 < 8/5 < 21/13 < ... < tau < ... 13/8 < 5/3 < 2/1). - Clark Kimberling, Jan 02 2022
For n > 1, a(n) is the smallest Fibonacci number of unit equilateral triangle tiles needed to make an isosceles trapezoid of height F(n) triangles. - Kiran Ananthpur Bacche, Sep 01 2024

Examples

			G.f. = x + 3*x^2 + 8*x^3 + 21*x^4 + 55*x^5 + 144*x^6 + 377*x^7 + 987*x^8 + ...
a(3) = 8 because there are exactly 8 idempotent order-preserving full transformations on a 3-element chain, namely: (1,2,3)->(1,1,1),(1,2,3)->(2,2,2),(1,2,3)->(3,3,3),(1,2,3)->(1,1,3),(1,2,3)->(2,2,3),(1,2,3)->(1,2,2),(1,2,3)->(1,3,3),(1,2,3)->(1,2,3)-mappings are coordinate-wise. - _Abdullahi Umar_, Sep 08 2008
		

References

  • Mohammad K. Azarian, The Generating Function for the Fibonacci Sequence, Missouri Journal of Mathematical Sciences, Vol. 2, No. 2, Spring 1990, pp. 78-79. Zentralblatt MATH, Zbl 1097.11516.
  • Mohammad K. Azarian, A Generalization of the Climbing Stairs Problem II, Missouri Journal of Mathematical Sciences, Vol. 16, No. 1, Winter 2004, pp. 12-17.
  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 2,5,6,14,33,55.
  • R. J. Douglas, Tournaments that admit exactly one Hamiltonian cycle, Proc. London Math. Soc., 21 (1970), 716-730.
  • G. Everest, A. van der Poorten, I. Shparlinski and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; see esp. p. 255.
  • A. Gerardin, Reply to Query 4389, L'Intermédiaire des Mathématiciens, 22 (1915), 23.
  • Julian Havil, The Irrationals, Princeton University Press, Princeton and Oxford, 2012, pp. 171-172.
  • Howie, J. M. Combinatorial and probabilistic results in transformation semigroups. Words, languages and combinatorics, II (Kyoto, 1992), 200--206, World Sci. Publ., River Edge, NJ, (1994).
  • Laradji, A. and Umar, A. Combinatorial results for semigroups of order-preserving full transformations. Semigroup Forum 72 (2006), 51-62.
  • I. Lukovits, A. Graovac, E. Kalman, G. Kaptay, P. Nagy, S. Nikolic, J. Sytchev and N. Trinajstich, "Nanotubes: Number of Kekulé Structures and Aromaticity", J. Chem. Inf. Comput. Sci, vol. 43 (2003), pp. 609-614. See Equation 6 on page 611.
  • T. Mansour, M. Shattuck, A statistic on n-color compositions and related sequences, Proc. Indian Acad. Sci. (Math. Sci.) Vol. 124, No. 2, May 2014, pp. 127-140.
  • H. Mathieu, Query 3932, L'Intermédiaire des Mathématiciens, 18 (1911), 222. - N. J. A. Sloane, Mar 08 2022
  • I. Niven and H. S. Zuckerman, An Introduction to the Theory of Numbers. 2nd ed., Wiley, NY, 1966, p. 101.
  • Paulo Ribenboim, Primes in Lucas sequences (Chap 4), in 'My Numbers, My Friends', Springer-Verlag 2000 NY, page 27.
  • 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).
  • R. Stanley, Enumerative combinatorics, Vol. 1, Cambridge University Press, Cambridge, 1997, pp. 96-100.

Crossrefs

Fibonacci A000045 = union of this sequence and A001519.
Inverse sequences A130259 and A130260.

Programs

  • Haskell
    a001906 n = a001906_list !! n
    a001906_list =
       0 : 1 : zipWith (-) (map (* 3) $ tail a001906_list) a001906_list
    -- Reinhard Zumkeller, Oct 03 2011
    
  • Magma
    [Fibonacci(2*n): n in [0..30]]; // Vincenzo Librandi, Sep 10 2014
  • Maple
    with(combstruct): SeqSeqSeqL := [T, {T=Sequence(S, card > 0), S=Sequence(U, card > 1), U=Sequence(Z, card >0)}, unlabeled]: seq(count(SeqSeqSeqL, size=n+1), n=0..28); # Zerinvary Lajos, Apr 04 2009
    H := (n, a, b) -> hypergeom([a - n/2, b - n/2], [1 - n], -4):
    a := n -> `if`(n = 0, 0, H(2*n, 1, 1/2)):
    seq(simplify(a(n)), n=0..30); # Peter Luschny, Sep 03 2019
    A001906 := proc(n)
        combinat[fibonacci](2*n) ;
    end proc:
    seq(A001906(n),n=0..20) ; # R. J. Mathar, Jan 11 2024
  • Mathematica
    f[n_] := Fibonacci[2n]; Array[f, 28, 0] (* or *)
    LinearRecurrence[{3, -1}, {0, 1}, 28] (* Robert G. Wilson v, Jul 13 2011 *)
    Take[Fibonacci[Range[0,60]],{1,-1,2}] (* Harvey P. Dale, May 23 2012 *)
    Table[ ChebyshevU[n-1, 3/2], {n, 0, 30}] (* Jean-François Alcover, Jan 25 2013, after Michael Somos *)
    CoefficientList[Series[(x)/(1 - 3x + x^2), {x, 0, 30}], x] (* Vincenzo Librandi, Sep 10 2014 *)
  • Maxima
    makelist(fib(2*n),n,0,30); /* Martin Ettl, Oct 21 2012 */
    
  • MuPAD
    numlib::fibonacci(2*n) $ n = 0..35; // Zerinvary Lajos, May 09 2008
    
  • PARI
    {a(n) = fibonacci(2*n)}; /* Michael Somos, Dec 06 2002 */
    
  • PARI
    {a(n) = subst( poltchebi(n+1)*4 - poltchebi(n)*6, x, 3/2)/5}; /* Michael Somos, Dec 06 2002 */
    
  • PARI
    {a(n) = polchebyshev( n-1, 2, 3/2)}; /* Michael Somos Jun 18 2011 */
    
  • PARI
    Vec(x/(1-3*x+x^2)+O(x^99)) \\ Charles R Greathouse IV, Oct 24 2012
    
  • Python
    def a(n, adict={0:0, 1:1}):
        if n in adict:
            return adict[n]
        adict[n]=3*a(n-1) - a(n-2)
        return adict[n] # David Nacin, Mar 04 2012
    
  • Sage
    [lucas_number1(n,3,1) for n in range(27)] # Zerinvary Lajos, Jun 25 2008
    
  • Sage
    [fibonacci(2*n) for n in range(0, 28)] # Zerinvary Lajos, May 15 2009
    

Formula

G.f.: x / (1 - 3*x + x^2). - Simon Plouffe in his 1992 dissertation
a(n) = 3*a(n-1) - a(n-2) = A000045(2*n).
a(n) = -a(-n).
a(n) = A060921(n-1, 0), n >= 1.
a(n) = sqrt((A005248(n)^2 - 4)/5).
a(n) = A007598(n) - A007598(n-2), n > 1.
a(n) = (ap^n - am^n)/(ap-am), with ap := (3+sqrt(5))/2, am := (3-sqrt(5))/2.
Invert transform of natural numbers: a(n) = Sum_{k=1..n} k*a(n-k), a(0) = 1. - Vladeta Jovovic, Apr 27 2001
a(n) = S(n-1, 3) with S(n, x) = U(n, x/2) Chebyshev's polynomials of the 2nd kind, see A049310.
a(n) = Sum_{k=0..n} binomial(n, k)*F(k). - Benoit Cloitre, Sep 03 2002
Limit_{n->infinity} a(n)/a(n-1) = 1 + phi = (3 + sqrt(5))/2. This sequence includes all of the elements of A033888 combined with A033890.
a(0)=0, a(1)=1, a(2)=3, a(n)*a(n-2) + 1 = a(n-1)^2. - Benoit Cloitre, Dec 06 2002
a(n) = n + Sum_{k=0..n-1} Sum_{i=0..k} a(i) = n + A054452(n). - Benoit Cloitre, Jan 26 2003
a(n) = Sum_{k=1..n} binomial(n+k-1, n-k). - Vladeta Jovovic, Mar 23 2003
E.g.f.: (2/sqrt(5))*exp(3*x/2)*sinh(sqrt(5)*x/2). - Paul Barry, Apr 11 2003
Second diagonal of array defined by T(i, 1) = T(1, j) = 1, T(i, j) = Max(T(i-1, j) + T(i-1, j-1); T(i-1, j-1) + T(i, j-1)). - Benoit Cloitre, Aug 05 2003
a(n) = F(n)*L(n) = A000045(n)*A000032(n). - Lekraj Beedassy, Nov 17 2003
F(2n+2) = 1, 3, 8, ... is the binomial transform of F(n+2). - Paul Barry, Apr 24 2004
Partial sums of A001519(n). - Lekraj Beedassy, Jun 11 2004
a(n) = Sum_{i=0..n-1} binomial(2*n-1-i, i)*5^(n-i-1)*(-1)^i. - Mario Catalani (mario.catalani(AT)unito.it), Jul 23 2004
a(n) = Sum_{k=0..n} binomial(n+k, n-k-1) = Sum_{k=0..n} binomial(n+k, 2k+1).
a(n+1) = Sum_{k=0..floor(n/2)} binomial(n-k, k)*(-1)^k*3^(n-2*k). - Paul Barry, Oct 25 2004
a(n) = (n*L(n) - F(n))/5 = Sum_{k=0..n-1} (-1)^n*L(2*n-2*k-1).
The i-th term of the sequence is the entry (1, 2) in the i-th power of the 2 X 2 matrix M = ((1, 1), (1, 2)). - Simone Severini, Oct 15 2005
Computation suggests that this sequence is the Hankel transform of A005807. The Hankel transform of {a(n)} is Det[{{a(1), ..., a(n)}, {a(2), ..., a(n+1)}, ..., {a(n), ..., a(2n-1)}}]. - John W. Layman, Jul 21 2000
a(n+1) = (A005248(n+1) - A001519(n))/2. - Creighton Dement, Aug 15 2004
a(n+1) = Sum_{i=0..n} Sum_{j=0..n} binomial(n-i, j)*binomial(n-j, i). - N. J. A. Sloane, Feb 20 2005
a(n) = (2/sqrt(5))*sinh(2*n*psi), where psi:=log(phi) and phi=(1+sqrt(5))/2. - Hieronymus Fischer, Apr 24 2007
a(n) = ((phi+1)^n - A001519(n))/phi with phi=(1+sqrt(5))/2. - Reinhard Zumkeller, Nov 22 2007
Row sums of triangle A135871. - Gary W. Adamson, Dec 02 2007
a(n)^2 = Sum_{k=1..n} a(2*k-1). This is a property of any sequence S(n) such that S(n) = B*S(n-1) - S(n-2) with S(0) = 0 and S(1) = 1 including {0,1,2,3,...} where B = 2. - Kenneth J Ramsey, Mar 23 2008
a(n) = 1/sqrt(5)*(phi^(2*n+2) - phi^(-2*n-2)), where phi = (1+sqrt(5))/2, the golden ratio. - Udita Katugampola (SIU), Sep 24 2008
If p[i] = i and if A is Hessenberg matrix of order n defined by: A[i,j] = p[j-i+1], (i<=j), A[i,j] = -1, (i = j+1), and A[i,j] = 0 otherwise. Then, for n >= 1, a(n) = det(A). - Milan Janjic, May 02 2010
If p[i] = Stirling2(i,2) and if A is the Hessenberg matrix of order n defined by: A[i,j] = p[j-i+1], (i<=j), A[i,j] = -1, (i = j+1), and A[i,j] = 0 otherwise. Then, for n >= 1, a(n-1) = det(A). - Milan Janjic, May 08 2010
a(n) = F(2*n+10) mod F(2*n+5).
a(n) = 1 + a(n-1) + Sum_{i=1..n-1} a(i), with a(0)=0. - Gary W. Adamson, Feb 19 2011
a(n) is equal to the permanent of the (n-1) X (n-1) Hessenberg matrix with 3's along the main diagonal, i's along the superdiagonal and the subdiagonal (i is the imaginary unit), and 0's everywhere else. - John M. Campbell, Jun 09 2011
a(n), n > 1 is equal to the determinant of an (n-x) X (n-1) tridiagonal matrix with 3's in the main diagonal, 1's in the super and subdiagonals, and the rest 0's. - Gary W. Adamson, Jun 27 2011
a(n) = b such that Integral_{x=0..Pi/2} sin(n*x)/(3/2-cos(x)) dx = c + b*log(3). - Francesco Daddi, Aug 01 2011
a(n+1) = Sum_{k=0..n} A101950(n,k)*2^k. - Philippe Deléham, Feb 10 2012
G.f.: A(x) = x/(1-3*x+x^2) = G(0)/sqrt(5); where G(k)= 1 -(a^k)/(1 - b*x/(b*x - 2*(a^k)/G(k+1))), a = (7-3*sqrt(5))/2, b = 3+sqrt(5), if |x|<(3-sqrt(5))/2 = 0.3819660...; (continued fraction 3 kind, 3-step ). - Sergei N. Gladkovskii, Jun 25 2012
a(n) = 2^n*b(n;1/2) = -b(n;-1), where b(n;d), n=0,1,...,d, denote the delta-Fibonacci numbers defined in comments to A000045 (see also Witula's et al. papers). - Roman Witula, Jul 12 2012
Product_{n>=1} (1 + 1/a(n)) = 1 + sqrt(5). - Peter Bala, Dec 23 2012
Product_{n>=2} (1 - 1/a(n)) = (1/6)*(1 + sqrt(5)). - Peter Bala, Dec 23 2012
G.f.: x/(1-2*x) + x^2/(1-2*x)/(Q(0)-x) where Q(k) = 1 - x/(x*k+1)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Feb 23 2013
G.f.: G(0)/2 - 1, where G(k) = 1 + 1/( 1 - x/(x + (1-x)^2/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 16 2013
G.f.: x*G(0)/(2-3*x), where G(k) = 1 + 1/( 1 - x*(5*k-9)/(x*(5*k-4) - 6/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 17 2013
Sum_{n>=1} 1/(a(n) + 1/a(n)) = 1. Compare with A001519, A049660 and A049670. - Peter Bala, Nov 29 2013
a(n) = U(n-1,3/2) where U(n-1,x) is Chebyshev polynomial of the second kind. - Milan Janjic, Jan 25 2015
The o.g.f. A(x) satisfies A(x) + A(-x) + 6*A(x)*A(-x) = 0. The o.g.f. for A004187 equals -A(sqrt(x))*A(-sqrt(x)). - Peter Bala, Apr 02 2015
For n > 1, a(n) = (3*F(n+1)^2 + 2*F(n-2)*F(n+1) - F(n-2)^2)/4. - J. M. Bergot, Feb 16 2016
For n > 3, a(n) = floor(MA) - 4 for n even and floor(MA) + 5 for n odd. MA is the maximum area of a quadrilateral with lengths of sides in order L(n), L(n), F(n-3), F(n+3), with L(n)=A000032(n). The ratio of the longer diagonal to the shorter approaches 5/3. - J. M. Bergot, Feb 16 2016
a(n+1) = Sum_{j=0..n} Sum_{k=0..j} binomial(n-j,k)*binomial(j,k)*2^(j-k). - Tony Foster III, Sep 18 2017
a(n) = Sum_{k=0..n-1} Sum_{i=0..n-1} C(k+i,k-i). - Wesley Ivan Hurt, Sep 21 2017
a(n) = Sum_{k=1..A000041(n)} A305309(n, k), n >= 1. Also row sums of triangle A078812.- Wolfdieter Lang, May 31 2018
a(n) = H(2*n, 1, 1/2) for n > 0 where H(n, a, b) -> hypergeom([a - n/2, b - n/2], [1 - n], -4). - Peter Luschny, Sep 03 2019
Sum_{n>=1} 1/a(n) = A153386. - Amiram Eldar, Oct 04 2020
a(n) = A249450(n) + 2. - Leo Tavares, Oct 10 2021
a(n) = -2/(sqrt(5)*tan(2*arctan(phi^(2*n)))), where phi = A001622 is the golden ratio. - Diego Rattaggi, Nov 21 2021
a(n) = sinh(2*n*arcsinh(1/2))/sqrt(5/4). - Peter Luschny, May 21 2022
From Amiram Eldar, Dec 02 2024: (Start)
Product_{n>=1} (1 - (-1)^n/a(n)) = 1 + 1/sqrt(5) (A344212).
Product_{n>=2} (1 + (-1)^n/a(n)) = (5/6) * (1 + 1/sqrt(5)). (End)
a(n) = Sum_{k>=0} Fibonacci(2*n*k)/(Lucas(2*n)^(k+1)). - Diego Rattaggi, Jan 12 2025
Sum_{n>=0} a(n)/3^n = 3. - Diego Rattaggi, Jan 20 2025

A208510 Triangle of coefficients of polynomials u(n,x) jointly generated with A029653; see the Formula section.

Original entry on oeis.org

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

Views

Author

Clark Kimberling, Feb 28 2012

Keywords

Comments

Row sums: A083329
Alternating row sums: 1,0,-1,-1,-1,-1,-1,-1,-1,-1,...
Antidiagonal sums: A000071 (-1+Fibonacci numbers)
col 1: A000012
col 2: A005408
col 3: A000290
col 4: A000330
col 5: A002415
col 6: A005585
col 7: A040977
col 8: A050486
col 9: A053347
col 10: A054333
col 11: A054334
col 12: A057788
col 2n-1 of A208510 is column n of A208508
col 2n of A208510 is column n of A208509.
...
GENERAL DISCUSSION:
A208510 typifies arrays generated by paired recurrence equations of the following form:
u(n,x)=a(n,x)*u(n-1,x)+b(n,x)*v(n-1,x)+c(n,x)
v(n,x)=d(n,x)*u(n-1,x)+e(n,x)*v(n-1,x)+f(n,x).
...
These first-order recurrences imply separate second-order recurrences. In order to show them, the six functions a(n,x),...,f(n,x) are abbreviated as a,b,c,d,e,f.
Then, starting with initial values u(1,x)=1 and u(2,x)=a+b+c: u(n,x) = (a+e)u(n-1,x) + (bd-ae)u(n-2,x) + bf-ce+c.
With initial values v(1,x)=1 and v(2,x)=d+e+f: v(n,x) = (a+e)v(n-1,x) + (bd-ae)v(n-2,x) + cd-af+f.
...
In the guide below, the last column codes certain sequences that occur in one of these ways: row, column, edge, row sum, alternating row sum. Coding:
A: 1,-1,1,-1,1,-1,1.... A033999
B: 1,2,4,8,16,32,64,... powers of 2
C: 1,1,1,1,1,1,1,1,.... A000012
D: 2,2,2,2,2,2,2,2,.... A007395
E: 2,4,6,8,10,12,14,... even numbers
F: 1,1,2,3,5,8,13,21,.. Fibonacci numbers
N: 1,2,3,4,5,6,7,8,.... A000027
O: 1,3,5,7,9,11,13,.... odd numbers
P: 1,3,9,27,81,243,.... powers of 3
S: 1,4,9,16,25,36,49,.. squares
T: 1,3,6,10,15,21,38,.. triangular numbers
Z: 1,0,0,0,0,0,0,0,0,.. A000007
*: (eventually) periodic alternating row sums
^: has a limiting row; i.e., the polynomials "approach" a power series
This coding includes indirect and repeated occurrences; e.g. F occurs thrice at A094441: in column 1 directly as Fibonacci numbers, in row sums as odd-indexed Fibonacci numbers, and in alternating row sums as signed Fibonacci numbers.
......... a....b....c....d....e....f....code
A034839 u 1....1....0....1....x....0....CCOT
A034867 v 1....1....0....1....x....0....CEN
A210221 u 1....1....0....1....2x...0....BBFF
A210596 v 1....1....0....1....2x...0....BBFF
A105070 v 1....2x...0....1....1....0....BN
A207605 u 1....1....0....1....x+1..0....BCFFN
A106195 v 1....1....0....1....x+1..0....BCFFN
A207606 u 1....1....0....x....x+1..0....DNT
A207607 v 1....1....0....x....x+1..0....DNT
A207608 u 1....1....0....2x...x+1..0....N
A207609 v 1....1....0....2x...x+1..0....C
A207610 u 1....1....0....1....x....1....CF
A207611 v 1....1....0....1....x....1....BCF
A207612 u 1....1....0....1....2x...1....BF
A207613 v 1....1....0....1....2x...1....BF
A207614 u 1....1....0....1....x+1..1....CN
A207615 v 1....1....0....1....x+1..1....CFN
A207616 u 1....1....0....x....1....1....CE
A207617 v 1....1....0....x....1....1....CNO
A029638 u 1....1....0....x....x....1....CDNO
A029635 v 1....1....0....x....x....1....CDNOZ
A207618 u 1....1....0....x....2x...1....N
A207619 v 1....1....0....x....2x...1....CFN
A207620 u 1....1....0....x....x+1..1....DET
A207621 v 1....1....0....x....x+1..1....DNO
A207622 u 1....1....0....2x...1....1....BT
A207623 v 1....1....0....2x...1....1....BN
A207624 u 1....1....0....2x...x....1....N
A102662 v 1....1....0....2x...x....1....CO
A207625 u 1....1....0....2x...x+1..1....T
A207626 v 1....1....0....2x...x+1..1....N
A207627 u 1....1....0....2x...2x...1....BN
A207628 v 1....1....0....2x...2x...1....BCE
A207629 u 1....1....0....x+1..1....1....CET
A207630 v 1....1....0....x+1..1....1....CO
A207631 u 1....1....0....x+1..x....1....DF
A207632 v 1....1....0....x+1..x....1....DEF
A207633 u 1....1....0....x+1..2x...1....F
A207634 v 1....1....0....x+1..2x...1....F
A207635 u 1....1....0....x+1..x+1..1....DN
A207636 v 1....1....0....x+1..x+1..1....CD
A160232 u 1....x....0....1....2x...0....BCFN
A208341 v 1....x....0....1....2x...0....BCFFN
A085478 u 1....x....0....1....x+1..0....CCOFT*
A078812 v 1....x....0....1....x+1..0....CEFN*
A208342 u 1....x....0....x....x....0....CCFNO
A208343 v 1....x....0....x....x....0....BBCDFZ
A208344 u 1....x....0....x....2x...0....CCFN
A208345 v 1....x....0....x....2x...0....CFZ
A094436 u 1....x....0....x....x+1..0....CFFN
A094437 v 1....x....0....x....x+1..0....CEFF
A117919 u 1....x....0....2x...1....0....BCNT
A135837 v 1....x....0....2x...1....0....BCET
A208328 u 1....x....0....2x...x....0....CCOP
A208329 v 1....x....0....2x...x....0....DPZ
A208330 u 1....x....0....2x...x+1..0....CNPT
A208331 v 1....x....0....2x...x+1..0....CN
A208332 u 1....x....0....2x...2x...0....CCE
A208333 v 1....x....0....2x...2x...0....DZ
A208334 u 1....x....0....x+1..1....0....CCNT
A208335 v 1....x....0....x+1..1....0....CCN*
A208336 u 1....x....0....x+1..x....0....CFNT*
A208337 v 1....x....0....x+1..x....0....ACFN*
A208338 u 1....x....0....x+1..2x...0....CNP
A208339 v 1....x....0....x+1..2x...0....BCNP
A202390 u 1....x....0....x+1..x+1..0....CFPTZ*
A208340 v 1....x....0....x+1..x+1..0....FNPZ*
A208508 u 1....x....0....1....1....1....CCES
A208509 v 1....x....0....1....1....1....BCO
A208510 u 1....x....0....1....x....1....CCCNOS*
A029653 v 1....x....0....1....x....1....BCDOSZ*
A208511 u 1....x....0....1....2x...1....BCFO
A208512 v 1....x....0....1....2x...1....BDFO
A208513 u 1....x....0....1....x+1..1....CCES*
A111125 v 1....x....0....1....x+1..1....COO*
A133567 u 1....x....0....x....1....1....CCOTT
A133084 v 1....x....0....x....1....1....BBCEN
A208514 u 1....x....0....x....x....1....CEFN
A208515 v 1....x....0....x....x....1....BCDFN
A208516 u 1....x....0....x....2x...1....CNN
A208517 v 1....x....0....x....2x...1....CCN
A208518 u 1....x....0....x....x+1..1....CFNT
A208519 v 1....x....0....x....x+1..1....NFFT
A208520 u 1....x....0....2x...1....1....BCTT
A208521 v 1....x....0....2x...1....1....BEN
A208522 u 1....x....0....2x...x....1....CCN
A208523 v 1....x....0....2x...x....1....CCO
A208524 u 1....x....0....2x...x+1..1....CT*
A208525 v 1....x....0....2x...x+1..1....ACNP*
A208526 u 1....x....0....2x...2x...1....CEN
A208527 v 1....x....0....2x...2x...1....CCE
A208606 u 1....x....0....x+1..1....1....CCS
A208607 v 1....x....0....x+1..1....1....CNO
A208608 u 1....x....0....x+1..x....1....CFOT
A208609 v 1....x....0....x+1..x....1....DEN*
A208610 u 1....x....0....x+1..2x...1....CO
A208611 v 1....x....0....x+1..2x...1....DE
A208612 u 1....x....0....x+1..x+1..1....CFNS
A208613 v 1....x....0....x+1..x+1..1....CFN*
A105070 u 1....2x...0....1....1....0....BN
A207536 u 1....2x...0....1....1....0....BCT
A208751 u 1....2x...0....1....x+1..0....CDPT
A208752 v 1....2x...0....1....x+1..0....CNP
A135837 u 1....2x...0....x....1....0....BCNT
A117919 v 1....2x...0....x....1....0....BCNT
A208755 u 1....2x...0....x....x....0....BCDEP
A208756 v 1....2x...0....x....x....0....BCCOZ
A208757 u 1....2x...0....x....2x...0....CDEP
A208758 v 1....2x...0....x....2x...0....CCEPZ
A208763 u 1....2x...0....2x...x....0....CDOP
A208764 v 1....2x...0....2x...x....0....CCCP
A208765 u 1....2x...0....2x...x+1..0....CE
A208766 v 1....2x...0....2x...x+1..0....CC
A208747 u 1....2x...0....2x...2x...0....CDE
A208748 v 1....2x...0....2x...2x...0....CCZ
A208749 u 1....2x...0....x+1..1....0....BCOPT
A208750 v 1....2x...0....x+1..1....0....BCNP*
A208759 u 1....2x...0....x+1..2x....0...CE
A208760 v 1....2x...0....x+1..2x....0...BCO
A208761 u 1....2x...0....x+1..x+1...0...BCCT*
A208762 v 1....2x...0....x+1..x+1...0...BNZ*
A208753 u 1....2x...0....1....1.....1...BCS
A208754 v 1....2x...0....1....1.....1...BO
A105045 u 1....2x...0....1....2x....1...BCCOS*
A208659 v 1....2x...0....1....2x....1...BDOSZ*
A208660 u 1....2x...0....1....x+1...1...CDS
A208904 v 1....2x...0....1....x+1...1...CNO
A208905 u 1....2x...0....x....1.....1...BCT
A208906 v 1....2x...0....x....1.....1...BNN
A208907 u 1....2x...0....x....x.....1...BCN
A208756 v 1....2x...0....x....x.....1...BCCE
A208755 u 1....2x...0....x....2x....1...CEN
A208910 v 1....2x...0....x....2x....1...CCE
A208911 u 1....2x...0....x....x+1...1...BCT
A208912 v 1....2x...0....x....x+1...1...BNT
A208913 u 1....2x...0....2x...1.....1...BCT
A208914 v 1....2x...0....2x...1.....1...BEN
A208915 u 1....2x...0....2x...x.....1...CE
A208916 v 1....2x...0....2x...x.....1...CCO
A208919 u 1....2x...0....2x...x+1...1...CT
A208920 v 1....2x...0....2x...x+1...1...N
A208917 u 1....2x...0....2x...2x....1...CEN
A208918 v 1....2x...0....2x...2x....1...CCNP
A208921 u 1....2x...0....x+1..1.....1...BC
A208922 v 1....2x...0....x+1..1.....1...BON
A208923 u 1....2x...0....x+1..x.....1...BCNO
A208908 v 1....2x...0....x+1..x.....1...BDN*
A208909 u 1....2x...0....x+1..2x....1...BN
A208930 v 1....2x...0....x+1..2x....1...DN
A208931 u 1....2x...0....x+1..x+1...1...BCOS
A208932 v 1....2x...0....x+1..x+1...1...BCO*
A207537 u 1....x+1..0....1....1.....0...BCO
A207538 v 1....x+1..0....1....1.....0...BCE
A122075 u 1....x+1..0....1....x.....0...CCFN*
A037027 v 1....x+1..0....1....x.....0...CCFN*
A209125 u 1....x+1..0....1....2x....0...BCFN*
A164975 v 1....x+1..0....1....2x....0...BF
A209126 u 1....x+1..0....x....x.....0...CDFO*
A209127 v 1....x+1..0....x....x.....0...DFOZ*
A209128 u 1....x+1..0....x....2x....0...CDE*
A209129 v 1....x+1..0....x....2x....0...DEZ
A102756 u 1....x+1..0....x....x+1...0...CFNP*
A209130 v 1....x+1..0....x....x+1...0...CCFNP*
A209131 u 1....x+1..0....2x...x.....0...CDEP*
A209132 v 1....x+1..0....2x...x.....0...CNPZ*
A209133 u 1....x+1..0....2x...2x....0...CDN
A209134 v 1....x+1..0....2x...2x....0...CCN*
A209135 u 1....x+1..0....2x...x+1...0...CN*
A209136 v 1....x+1..0....2x...x+1...0...CCS*
A209137 u 1....x+1..0....x+1..x.....0...CFFP*
A209138 v 1....x+1..0....x+1..x.....0...AFFP*
A209139 u 1....x+1..0....x+1..2x....0...CF*
A209140 v 1....x+1..0....x+1..2x....0...BF
A209141 u 1....x+1..0....x+1..x+1...0...BCF*
A209142 v 1....x+1..0....x+1..x+1...0...BFZ*
A209143 u 1....x+1..0....1....1.....1...CCE*
A209144 v 1....x+1..0....1....1.....1...COO*
A209145 u 1....x+1..0....1....x.....1...CCFN*
A122075 v 1....x+1..0....1....x.....1...CCFN*
A209146 u 1....x+1..0....1....2x....1...BCF*
A209147 v 1....x+1..0....1....2x....1...BF
A209148 u 1....x+1..0....1....x+1...1...CCO*
A209149 v 1....x+1..0....1....x+1...1...CDO*
A209150 u 1....x+1..0....x....1.....1...CCNT*
A208335 v 1....x+1..0....x....1.....1...CDNN*
A209151 u 1....x+1..0....x....x.....1...CFN*
A208337 v 1....x+1..0....x....x.....1...ACFN*
A209152 u 1....x+1..0....x....2x....1...CN*
A208339 v 1....x+1..0....x....x.....1...BCN
A209153 u 1....x+1..0....x....x+1...1...CFT*
A208340 v 1....x+1..0....x....x.....1...FNZ*
A209154 u 1....x+1..0....2x...1.....1...BCT*
A209157 v 1....x+1..0....2x...1.....1...BNN
A209158 u 1....x+1..0....2x...x.....1...CN*
A209159 v 1....x+1..0....2x...x.....1...CO*
A209160 u 1....x+1..0....2x...2x....1...CN*
A209161 v 1....x+1..0....2x...2x....1...CE
A209162 u 1....x+1..0....2x...x+1...1...CT*
A209163 v 1....x+1..0....2x...x+1...1...CO*
A209164 u 1....x+1..0....x+1..1.....1...CC*
A209165 v 1....x+1..0....x+1..1.....1...CCN
A209166 u 1....x+1..0....x+1..x.....1...CFF*
A209167 v 1....x+1..0....x+1..x.....1...FF*
A209168 u 1....x+1..0....x+1..2x....1...CF*
A209169 v 1....x+1..0....x+1..2x....1...CF
A209170 u 1....x+1..0....x+1..x+1...1...CF*
A209171 v 1....x+1..0....x+1..x+1...1...CF*
A053538 u x....1....0....1....1.....0...BBCCFN
A076791 v x....1....0....1....1.....0...BBCDF
A209172 u x....1....0....1....2x....0...BCCFF
A209413 v x....1....0....1....2x....0...BCCFF
A094441 u x....1....0....1....x+1...0...CFFFN
A094442 v x....1....0....1....x+1...0...CEFFF
A054142 u x....1....0....x....x+1...0...CCFOT*
A172431 v x....1....0....x....x+1...0...CEFN*
A008288 u x....1....0....2x...1.....0...CCOO*
A035607 v x....1....0....2x...1.....0...ACDE*
A209414 u x....1....0....2x...x+1...0...CCS
A112351 v x....1....0....2x...x+1...0...CON
A209415 u x....1....0....x+1..x.....0...CCTN
A209416 v x....1....0....x+1..x.....0...ACN*
A209417 u x....1....0....x+1..2x....0...CC
A209418 v x....1....0....x+1..2x....0...BBC
A209419 u x....1....0....x+1..x+1...0...CFTZ*
A209420 v x....1....0....x+1..x+1...0...FNZ*
A209421 u x....1....0....1....1.....1...CCN
A209422 v x....1....0....1....1.....1...CD
A209555 u x....1....0....1....x.....1...CNN
A209556 v x....1....0....1....x.....1...CNN
A209557 u x....1....0....1....2x....1...BCN
A209558 v x....1....0....1....2x....1...BN
A209559 u x....1....0....1....x+1...1...CN
A209560 v x....1....0....1....x+1...1...CN
A209561 u x....1....0....x....1.....1...CCNNT*
A209562 v x....1....0....x....1.....1...CDNNT*
A209563 u x....1....0....x....x.....1...CCFT^
A209564 v x....1....0....x....x.....1...CFN^
A209565 u x....1....0....x....2x....1...CC^
A209566 v x....1....0....x....2x....1...BC^
A209567 u x....1....0....x....x+1...1...CNT*
A209568 v x....1....0....x....x+1...1...NNS*
A209569 u x....1....0....2x...1.....1...CNO*
A209570 v x....1....0....2x...1.....1...DNN*
A209571 u x....1....0....2x...x.....1...CCS^
A209572 v x....1....0....2x...x.....1...CN^
A209573 u x....1....0....2x...x+1...1...CNS
A209574 v x....1....0....2x...x+1...1...NO
A209575 u x....1....0....2x...2x....1...CC
A209576 v x....1....0....2x...2x....1...C
A209577 u x....1....0....x+1..1.....1...CNNT
A209578 v x....1....0....x+1..1.....1...CNN
A209579 u x....1....0....x+1..x.....1...CNNT
A209580 v x....1....0....x+1..x.....1...NN*
A209581 u x....1....0....x+1..2x....1...CN
A209582 v x....1....0....x+1..2x....1...BN
A209583 u x....1....0....x+1..x+1...1...CT*
A209584 v x....1....0....x+1..x+1...1...CN*
A121462 u x....x....0....x....x+1...0...BCFFNZ
A208341 v x....x....0....x....x+1...0...BCFFN
A209687 u x....x....0....2x...x+1...0...BCNZ
A208339 v x....x....0....2x...x+1...0...BCN
A115241 u x....x....0....1....1.....1...CDNZ*
A209688 v x....x....0....1....1.....1...DDN*
A209689 u x....x....0....1....x.....1...FNZ^
A209690 v x....x....0....1....x.....1...FN^
A209691 u x....x....0....1....2x....1...BCZ^
A209692 v x....x....0....1....2x....1...BCC^
A209693 u x....x....0....1....x+1...1...NNZ*
A209694 v x....x....0....1....x+1...1...CN*
A209697 u x....x....0....x....x+1...1...BNZ
A209698 v x....x....0....x....x+1...1...BNT
A209699 u x....x....0....2x...1.....1...BNNZ
A209700 v x....x....0....2x...1.....1...BDN
A209701 u x....x....0....2x...x+1...1...NZ
A209702 v x....x....0....2x...x+1...1...N
A209703 u x....x....0....x+1..1.....1...FNTZ
A209704 v x....x....0....x+1..1.....1...FNNT
A209705 u x....x....0....x+1..x+1...1...BNZ*
A209706 v x....x....0....x+1..x+1...1...BCN*
A209695 u x....x+1..0....2x...x+1...0...ACN*
A209696 v x....x+1..0....2x...x+1...0...CDN*
A209830 u x....x+1..0....x+1..2x....0...ACF
A209831 v x....x+1..0....x+1..2x....0...BCF*
A209745 u x....x+1..0....x+1..x+1...0...ABF*
A209746 v x....x+1..0....x+1..x+1...0...BFZ*
A209747 u x....x+1..0....1....1.....1...ADE*
A209748 v x....x+1..0....1....1.....1...DEO
A209749 u x....x+1..0....1....x.....1...ANN*
A209750 v x....x+1..0....1....x.....1...CNO
A209751 u x....x+1..0....1....2x....1...ABN*
A209752 v x....x+1..0....1....2x....1...BN
A209753 u x....x+1..0....1....x+1...1...AN*
A209754 v x....x+1..0....1....x+1...1...NT*
A209755 u x....x+1..0....x....1.....1...AFN
A209756 v x....x+1..0....x....1.....1...FNO*
A209759 u x....x+1..0....x....2x....1...ACF^
A209760 v x....x+1..0....x....2x....1...CF^*
A209761 u x....x+1..0....x.....x+1..1...ABNS*
A209762 v x....x+1..0....x.....x+1..1...BNS*
A209763 u x....x+1..0....2x....1....1...ABN*
A209764 v x....x+1..0....2x....1....1...BNN
A209765 u x....x+1..0....2x....x....1...ACF^*
A209766 v x....x+1..0....2x....x....1...CF^
A209767 u x....x+1..0....2x....x+1..1...AN*
A209768 v x....x+1..0....2x....x+1..1...N*
A209769 u x....x+1..0....x+1...1....1...AF*
A209770 v x....x+1..0....x+1...1....1...FN
A209771 u x....x+1..0....x+1...x....1...ABN*
A209772 v x....x+1..0....x+1...x....1...BN*
A209773 u x....x+1..0....x+1...2x...1...AF
A209774 v x....x+1..0....x+1...2x...1...FN*
A209775 u x....x+1..0....x+1...x+1..1...AB*
A209776 v x....x+1..0....x+1...x+1..1...BC*
A210033 u 1....1....1....1.....x....1...BCN
A210034 v 1....1....1....1.....x....1...BCDFN
A210035 u 1....1....1....1.....2x...1...BBF
A210036 v 1....1....1....1.....2x...1...BBFF
A210037 u 1....1....1....1.....x+1..1...BCFFN
A210038 v 1....1....1....1.....x+1..1...BCFFN
A210039 u 1....1....1....x.....1....1...BCOT
A210040 v 1....1....1....x.....1....1...BCEN
A210042 u 1....1....1....x.....x....1...BCDEOT*
A124927 v 1....1....1....x.....x....1...BCDET*
A210041 u 1....1....1....x.....2x...1...BFO
A209758 v 1....1....1....x.....2x...1...BCFO
A210187 u 1....1....1....x.....x+1..1...DTF*
A210188 v 1....1....1....x.....x+1..1...DNF*
A210189 u 1....1....1....2x....1....1...BT
A210190 v 1....1....1....2x....1....1...BN
A210191 u 1....1....1....2x....x....1...CO*
A210192 v 1....1....1....2x....x....1...CCO*
A210193 u 1....1....1....2x....x+1..1...CPT
A210194 v 1....1....1....2x....x+1..1...CN
A210195 u 1....1....1....2x....2x...1...BOPT*
A210196 v 1....1....1....2x....2x...1...BCC*
A210197 u 1....1....1....x+1...1....1...BCOT
A210198 v 1....1....1....x+1...1....1...BCEN
A210199 u 1....1....1....x+1...x....1...DFT
A210200 v 1....1....1....x+1...x....1...DFO*
A210201 u 1....1....1....x+1...2x...1...BFP
A210202 v 1....1....1....x+1...2x...1...BF
A210203 u 1....1....1....x+1...x+1..1...BDOP
A210204 v 1....1....1....x+1...x+1..1...BCDN*
A210211 u x....1....1....1.....2x...1...BCFN
A210212 v x....1....1....1.....2x...1...BFN
A210213 u x....1....1....1.....x+1..1...CFFN
A210214 v x....1....1....1.....x+1..1...CFFO
A210215 u x....1....1....x.....x....1...BCDFT^
A210216 v x....1....1....x.....x....1...BCFO^
A210217 u x....1....1....x.....2x...1...CDF^
A210218 v x....1....1....x.....2x...1...BCF^
A210219 u x....1....1....x.....x+1..1...CNSTF*
A210220 v x....1....1....x.....x+1..1...FNNT*
A104698 u x....1....1....2x......1..1...CENS*
A210220 v x....1....1....2x....x+1..1...DNNT*
A210223 u x....1....1....2x....x....1...CD^
A210224 v x....1....1....2x....x....1...CO^
A210225 u x....1....1....2x....x+1..1...CNP
A210226 v x....1....1....2x....x+1..1...NOT
A210227 u x....1....1....2x....2x...1...CDP^
A210228 v x....1....1....2x....2x...1...C^
A210229 u x....1....1....x+1...1....1...CFNN
A210230 v x....1....1....x+1...1....1...CCN
A210231 u x....1....1....x+1...x....1...CNT
A210232 v x....1....1....x+1...x....1...NN*
A210233 u x....1....1....x+1...2x...1...CNP
A210234 v x....1....1....x+1...2x...1...BN
A210235 u x....1....1....x+1...x+1..1...CCFPT*
A210236 v x....1....1....x+1...x+1..1...CFN*
A124927 u x....x....1....1.....1....1...BCDEET*
A210042 v x....1....1....x+1...x+1..1...BDEOT*
A210216 u x....x....1....1.....x....1...BCFO^
A210215 v x....x....1....1.....x....1...BCDFT^
A210549 u x....x....1....1.....2x...1...BCF^
A210550 v x....x....1....1.....2x...1...BDF^
A172431 u x....x....1....1.....x+1..1...CEFN*
A210551 v x....x....1....1.....x+1..1...CFOT*
A210552 u x....x....1....x.....1....1...BBCFNO
A210553 v x....x....1....x.....1....1...BNNFB
A208341 u x....x....1....x.....x+1..1...BCFFN
A210554 v x....x....1....x.....x+1..1...BNFFT
A210555 u x....x....1....2x....1....1...BCNN
A210556 v x....x....1....2x....1....1...BENP
A210557 u x....x....1....2x....x+1..1...CNP
A210558 v x....x....1....2x....x+1..1...N
A210559 u x....x....1....x+1...1....1...CEF
A210560 v x....x....1....x+1...1....1...OFNS
A210561 u x....x....1....x+1...x....1...BCNP^
A210562 v x....x....1....x+1...x....1...BDP*^
A210563 u x....x....1....x+1...2x...1...CFP^
A210564 v x....x....1....x+1...2x...1...DF^
A013609 u x....x....1....x+1...x+1..1...BCEPT*
A209757 v x....x....1....x+1...x+1..1...BCOS*
A209819 u x....2x...1....x+1...x....1...CFN^
A209820 v x....2x...1....x+1...x....1...DF^
A209996 u x....2x...1....x+1...2x...1...CP^
A209998 v x....2x...1....x+1...2x...1...DP^
A209999 u x....x+1..1....1.....x+1..1...FN*
A210287 v x....x+1..1....1.....x+1..1...CFT*
A210565 u x....x+1..1....x.....1....1...FNT*
A210595 v x....x+1..1....x.....1....1...FNNT
A210598 u x....x+1..1....x+1...2x...1...FN*
A210599 v x....x+1..1....x+1...2x...1...FN
A210600 u x....x+1..1....x+1...x+1..1...BF*
A210601 v x....x+1..1....x+1...x+1..1...BF*
A210597 u 2x...1....1....x+1...1....1...BF
A210601 v 2x...1....1....x+1...1....1...BFN*
A210603 u 2x...1....1....x+1...x+1..1...BF
A210738 v 2x...1....1....x+1...x+1..1...CBF*
A210739 u 2x...x....1....x+1...x....1...CF^
A210740 v 2x...x....1....x+1...x....1...DF*^
A210741 u 2x...x....1....x+1...x+1..1...BCFO
A210742 v 2x...x....1....x+1...x+1..1...CFO*
A210743 u 2x...x+1..1....x+1...1....1...F
A210744 v 2x...x+1..1....x+1...1....1...FN
A210747 u 2x...x+1..1....x+1...x+1..1...FF
A210748 v 2x...x+1..1....x+1...x+1..1...CFF*
A210749 u x+1..1....1....x+1...2x...1...BCF
A210750 v x+1..1....1....x+1...2x...1...BF
A210751 u x+1..x....1....x+1...2x...1...FNT
A210752 v x+1..x....1....x+1...2x...1...FN
A210753 u x+1..x....1....x+1...x+1..1...BNZ*
A210754 v x+1..x....1....x+1...x+1..1...BCT*
A210755 u x+1..2x...1....x+1...x+1..1...N*
A210756 v x+1..2x...1....x+1...x+1..1...CT*
A210789 u 1....x....0....x+2...x-1..0...CFFN
A210790 v 1....x....0....x+2...x-1..0...CEFF
A210791 u 1....x....0....x-1...x+2..0...CFNP
A210792 v 1....x....0....x-1...x+2..0...CF
A210793 u 1....x+1..0....x+2...x-1..0...CFNP
A210794 v 1....x+1..0....x+2...x-1..0...FPP
A210795 u 1....x....1....x+2...x-1..0...FN
A210796 v 1....x....1....x+2...x-1..0...FO
A210797 u 1....x....0....x+2...x-1..1...CF
A210798 v 1....x....0....x+2...x-1..1...F
A210799 u 1....x+1..1....x+2...x-1..0...FN
A210800 v 1....x+1..1....x+2...x-1..0...F
A210801 u 1....x+1..1....x+2...x-1..1...FN
A210802 v 1....x+1..1....x+2...x-1..1...F
A210803 u 1....x....0....x-1...x+3..0...F*
A210804 v 1....x....0....x-1...x+3..0...F*
A210805 u 1....x....0....x+2...x-1.-1...CFFN
A210806 v 1....x....0....x+2...x-1.-1...FF
A210858 u 1....x....0....x+n...x....0...CFT*
A210859 v 1....x....0....x+n...x....0...FN*
A210860 u 1....x+1..0....x+n...x....0...F
A210861 v 1....x+1..0....x+n...x....0...F*
A210862 u 1....x....1....x+n-1.x....0...FN
A210863 v 1....x....1....x+n-1.x....0...FS
A210864 u 1....x....1....x+n...x....0...FN
A210865 v 1....x....1....x+n...x....0...FT
A210866 u 1....x....0....x+n...x...-x...CFT
A210867 v 1....x....0....x+n...x...-x...FN
A210868 u 1....x....0....x+1...x-1..0...BCFN
A210869 v 1....x....0....x+1...x-1..0...BBCFNZ
A210870 u 1....x....0....x+1...x-1..1...CFFN
A210871 v 1....x....0....x+1...x-1..1...CFF
A210872 u x....1...-1....x.....x....1...BDFZ^
A210873 v x....1...-1....x.....x....1...BCFN^
A210876 u x....1....1....x.....x....x...BCCF^
A210877 v x....1....1....x.....x....x...BDFNZ^
A210878 u x....2x...0....x+1...x....1...DFZ^
A210879 v x....2x...0....x+1...x....1...FC*^
Some of these triangles have irregular row lengths, making it difficult to retrieve individual rows/columns/diagonals without actually computing the recurrence. - Georg Fischer, Sep 04 2021

Examples

			First five rows:
1
1...1
1...3...1
1...5...4...1
1...7...9...5...1
First five polynomials u(n,x):
1
1 + x
1 + 3x + x^2
1 + 5x + 4x^2 + x^3
1 + 7x + 9x^2 + 5x^3 + x^4
		

Crossrefs

Programs

  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := u[n - 1, x] + x*v[n - 1, x];
    v[n_, x_] := u[n - 1, x] + x*v[n - 1, x] + 1;
    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[%]   (* A208510 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]   (* A029653 *)
  • Python
    from sympy import Poly
    from sympy.abc import x
    def u(n, x): return 1 if n==1 else u(n - 1, x) + x*v(n - 1, x)
    def v(n, x): return 1 if n==1 else u(n - 1, x) + x*v(n - 1, x) + 1
    def a(n): return Poly(u(n, x), x).all_coeffs()[::-1]
    for n in range(1, 13): print(a(n)) # Indranil Ghosh, May 27 2017

Formula

u(n,x)=u(n-1,x)+x*v(n-1,x),
v(n,x)=u(n-1,x)+x*v(n-1,x)+1,
where u(1,x)=1, v(1,x)=1.
Also, u(n,x)=(x+1)*u(n-1,x)+x for n>2, with u(n,2)=x+1.

Extensions

Corrected by Philippe Deléham, Apr 10 2012
Corrections and additions by Clark Kimberling, May 09 2012
Corrections in the overview by Georg Fischer, Sep 04 2021

A004736 Triangle read by rows: row n lists the first n positive integers in decreasing order.

Original entry on oeis.org

1, 2, 1, 3, 2, 1, 4, 3, 2, 1, 5, 4, 3, 2, 1, 6, 5, 4, 3, 2, 1, 7, 6, 5, 4, 3, 2, 1, 8, 7, 6, 5, 4, 3, 2, 1, 9, 8, 7, 6, 5, 4, 3, 2, 1, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1
Offset: 1

Views

Author

R. Muller

Keywords

Comments

Old name: Triangle T(n,k) = n-k, n >= 1, 0 <= k < n. Fractal sequence formed by repeatedly appending strings m, m-1, ..., 2, 1.
The PARI functions t1 (this sequence), t2 (A002260) can be used to read a square array T(n,k) (n >= 1, k >= 1) by antidiagonals upwards: n -> T(t1(n), t2(n)). - Michael Somos, Aug 23 2002, edited by M. F. Hasler, Mar 31 2020
A004736 is the mirror of the self-fission of the polynomial sequence (q(n,x)) given by q(n,x) = x^n+ x^(n-1) + ... + x + 1. See A193842 for the definition of fission. - Clark Kimberling, Aug 07 2011
Seen as flattened list: a(A000217(n)) = 1; a(A000124(n)) = n and a(m) <> n for m < A000124(n). - Reinhard Zumkeller, Jul 22 2012
Sequence B is called a reverse reluctant sequence of sequence A, if B is triangle array read by rows: row number k lists first k elements of the sequence A in reverse order. Sequence A004736 is the reverse reluctant sequence of sequence 1,2,3,... (A000027). - Boris Putievskiy, Dec 13 2012
The row sums equal A000217(n). The alternating row sums equal A004526(n+1). The antidiagonal sums equal A002620(n+1) respectively A008805(n-1). - Johannes W. Meijer, Sep 28 2013
From Peter Bala, Jul 29 2014: (Start)
Riordan array (1/(1-x)^2,x). Call this array M and for k = 0,1,2,... define M(k) to be the lower unit triangular block array
/I_k 0\
\ 0 M/
having the k X k identity matrix I_k as the upper left block; in particular, M(0) = M. Then the infinite matrix product M(0)*M(1)*M(2)*... is equal to A078812. (End)
T(n, k) gives the number of subsets of [n] := {1, 2, ..., n} with k consecutive numbers (consecutive k-subsets of [n]). - Wolfdieter Lang, May 30 2018
a(n) gives the distance from (n-1) to the smallest triangular number > (n-1). - Ctibor O. Zizka, Apr 09 2020
To construct the sequence, start from 1,2,,3,,,4,,,,5,,,,,6... where there are n commas after each "n". Then fill the empty places by the sequence itself. - Benoit Cloitre, Aug 17 2021
T(n,k) is the number of cycles of length 2*(k+1) in the (n+1)-ladder graph. There are no cycles of odd length. - Mohammed Yaseen, Jan 14 2023
The first 77 entries are also the signature sequence of log(3)=A002391. Then the two sequences start to differ. - R. J. Mathar, May 27 2024

Examples

			The triangle T(n, k) starts:
   n\k  1  2  3  4  5  6  7  8  9 10 11 12 ...
   1:   1
   2:   2  1
   3:   3  2  1
   4:   4  3  2  1
   5:   5  4  3  2  1
   6:   6  5  4  3  2  1
   7:   7  6  5  4  3  2  1
   8:   8  7  6  5  4  3  2  1
   9:   9  8  7  6  5  4  3  2  1
  10:  10  9  8  7  6  5  4  3  2  1
  11:  11 10  9  8  7  6  5  4  3  2  1
  12:  12 11 10  9  8  7  6  5  4  3  2  1
  ... Reformatted. - _Wolfdieter Lang_, Feb 04 2015
T(6, 3) = 4 because the four consecutive 3-subsets of [6] = {1, 2, ..., 6} are {1, 2, 3}, {2, 3, 4}, {3, 4, 5} and {4, 5, 6}. - _Wolfdieter Lang_, May 30 2018
		

References

  • H. S. M. Coxeter, Regular Polytopes, 3rd ed., Dover, NY, 1973, pp 159-162.

Crossrefs

Ordinal transform of A002260. See also A078812.
Cf. A141419 (partial sums per row).
Cf. A134546 (T * A051731, matrix product).
See A001511 for definition of ordinal transform.
Cf. A128174 (parity).

Programs

  • Excel
    =if(row()>=column();row()-column()+1;"") [Mats Granvik, Jan 19 2009]
    
  • Haskell
    a004736 n k = n - k + 1
    a004736_row n = a004736_tabl !! (n-1)
    a004736_tabl = map reverse a002260_tabl
    -- Reinhard Zumkeller, Aug 04 2014, Jul 22 2012
    
  • Maple
    A004736 := proc(n,m) n-m+1 ; end:
    T := (n, k) -> n-k+1: seq(seq(T(n,k), k=1..n), n=1..13); # Johannes W. Meijer, Sep 28 2013
  • Mathematica
    Flatten[ Table[ Reverse[ Range[n]], {n, 12}]] (* Robert G. Wilson v, Apr 27 2004 *)
    Table[Range[n,1,-1],{n,20}]//Flatten (* Harvey P. Dale, May 27 2020 *)
  • PARI
    {a(n) = 1 + binomial(1 + floor(1/2 + sqrt(2*n)), 2) - n}
    
  • PARI
    {t1(n) = binomial( floor(3/2 + sqrt(2*n)), 2) - n + 1} /* A004736 */
    
  • PARI
    {t2(n) = n - binomial( floor(1/2 + sqrt(2*n)), 2)} /* A002260 */
    
  • PARI
    apply( A004736(n)=1-n+(n=sqrtint(8*n)\/2)*(n+1)\2, [1..99]) \\ M. F. Hasler, Mar 31 2020
    
  • Python
    def agen(rows):
        for n in range(1, rows+1): yield from range(n, 0, -1)
    print([an for an in agen(13)]) # Michael S. Branicky, Aug 17 2021
    
  • Python
    from math import comb, isqrt
    def A004736(n): return comb((m:=isqrt(k:=n<<1))+(k>m*(m+1))+1,2)+1-n # Chai Wah Wu, Nov 08 2024

Formula

a(n+1) = 1 + A025581(n).
a(n) = (2 - 2*n + round(sqrt(2*n)) + round(sqrt(2*n))^2)/2. - Brian Tenneson, Oct 11 2003
G.f.: 1 / ((1-x)^2 * (1-x*y)). - Ralf Stephan, Jan 23 2005
Recursion: e(n,k) = (e(n - 1, k)*e(n, k - 1) + 1)/e(n - 1, k - 1). - Roger L. Bagula, Mar 25 2009
a(n) = (t*t+3*t+4)/2-n, where t = floor((-1+sqrt(8*n-7))/2). - Boris Putievskiy, Dec 13 2012
From Johannes W. Meijer, Sep 28 2013: (Start)
T(n, k) = n - k + 1, n >= 1 and 1 <= k <= n.
T(n, k) = A002260(n+k-1, n-k+1). (End)
a(n) = A000217(A002024(n)) - n + 1. - Enrique Pérez Herrero, Aug 29 2016

Extensions

New name from Omar E. Pol, Jul 15 2012

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

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 6, 5, 1, 1, 10, 15, 7, 1, 1, 15, 35, 28, 9, 1, 1, 21, 70, 84, 45, 11, 1, 1, 28, 126, 210, 165, 66, 13, 1, 1, 36, 210, 462, 495, 286, 91, 15, 1, 1, 45, 330, 924, 1287, 1001, 455, 120, 17, 1, 1, 55, 495, 1716, 3003, 3003, 1820, 680, 153, 19, 1
Offset: 0

Views

Author

Philippe Deléham, Aug 14 2003

Keywords

Comments

Coefficient array for Morgan-Voyce polynomial b(n,x). A053122 (unsigned) is the coefficient array for B(n,x). Reversal of A054142. - Paul Barry, Jan 19 2004
This triangle is formed from even-numbered rows of triangle A011973 read in reverse order. - Philippe Deléham, Feb 16 2004
T(n,k) is the number of nondecreasing Dyck paths of semilength n+1, having k+1 peaks. T(n,k) is the number of nondecreasing Dyck paths of semilength n+1, having k peaks at height >= 2. T(n,k) is the number of directed column-convex polyominoes of area n+1, having k+1 columns. - Emeric Deutsch, May 31 2004
Riordan array (1/(1-x), x/(1-x)^2). - Paul Barry, May 09 2005
The triangular matrix a(n,k) = (-1)^(n+k)*T(n,k) is the matrix inverse of A039599. - Philippe Deléham, May 26 2005
The n-th row gives absolute values of coefficients of reciprocal of g.f. of bottom-line of n-wave sequence. - Floor van Lamoen (fvlamoen(AT)planet.nl), Sep 24 2006
Unsigned version of A129818. - Philippe Deléham, Oct 25 2007
T(n, k) is also the number of idempotent order-preserving full transformations (of an n-chain) of height k >=1 (height(alpha) = |Im(alpha)|) and of waist n (waist(alpha) = max(Im(alpha))). - Abdullahi Umar, Oct 02 2008
A085478 is jointly generated with A078812 as a triangular array of coefficients of polynomials u(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
Per Kimberling's recursion relations, see A102426. - Tom Copeland, Jan 19 2016
Subtriangle of the triangle given by (0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Mar 26 2012
T(n,k) is also the number of compositions (ordered partitions) of 2*n+1 into 2*k+1 parts which are all odd. Proof: The o.g.f. of column k, x^k/(1-x)^(2*k+1) for k >= 0, is the o.g.f. of the odd-indexed members of the sequence with o.g.f. (x/(1-x^2))^(2*k+1) (bisection, odd part). Thus T(n,k) is obtained from the sum of the multinomial numbers A048996 for the partitions of 2*n+1 into 2*k+1 parts, all of which are odd. E.g., T(3,1) = 3 + 3 from the numbers for the partitions [1,1,5] and [1,3,3], namely 3!/(2!*1!) and 3!/(1!*2!), respectively. The number triangle with the number of these partitions as entries is A152157. - Wolfdieter Lang, Jul 09 2012
The matrix elements of the inverse are T^(-1)(n,k) = (-1)^(n+k)*A039599(n,k). - R. J. Mathar, Mar 12 2013
T(n,k) = A258993(n+1,k) for k = 0..n-1. - Reinhard Zumkeller, Jun 22 2015
The n-th row polynomial in descending powers of x is the n-th Taylor polynomial of the algebraic function F(x)*G(x)^n about 0, where F(x) = (1 + sqrt(1 + 4*x))/(2*sqrt(1 + 4*x)) and G(x) = ((1 + sqrt(1 + 4*x))/2)^2. For example, for n = 4, (1 + sqrt(1 + 4*x))/(2*sqrt(1 + 4*x)) * ((1 + sqrt(1 + 4*x))/2)^8 = (x^4 + 10*x^3 + 15*x^2 + 7*x + 1) + O(x^5). - Peter Bala, Feb 23 2018
Row n also gives the coefficients of the characteristc polynomial of the tridiagonal n X n matrix M_n given in A332602: Phi(n, x) := Det(M_n - x*1_n) = Sum_{k=0..n} T(n, k)*(-x)^k, for n >= 0, with Phi(0, x) := 1. - Wolfdieter Lang, Mar 25 2020
It appears that the largest root of the n-th degree polynomial is equal to the sum of the distinct diagonals of a (2*n+1)-gon including the edge, 1. The largest root of x^3 - 6*x^2 + 5*x - 1 is 5.048917... = the sum of (1 + 1.80193... + 2.24697...). Alternatively, the largest root of the n-th degree polynomial is equal to the square of sigma(2*n+1). Check: 5.048917... is the square of sigma(7), 2.24697.... Given N = 2*n+1, sigma(N) (N odd) can be defined as 1/(2*sin(Pi/(2*N))). Relating to the 9-gon, the largest root of x^4 - 10*x^3 + 15*x^2 - 7*x + 1 is 8.290859..., = the sum of (1 + 1.879385... + 2.532088... + 2.879385...), and is the square of sigma(9), 2.879385... Refer to A231187 for a further clarification of sigma(7). - Gary W. Adamson, Jun 28 2022
For n >=1, the n-th row is given by the coefficients of the minimal polynomial of -4*sin(Pi/(4*n + 2))^2. - Eric W. Weisstein, Jul 12 2023
Denoting this lower triangular array by L, then L * diag(binomial(2*k,k)^2) * transpose(L) is the LDU factorization of A143007, the square array of crystal ball sequences for the A_n X A_n lattices. - Peter Bala, Feb 06 2024
T(n, k) is the number of occurrences of the periodic substring (01)^k in the periodic string (01)^n (see Proposition 4.7 at page 7 in Fang). - Stefano Spezia, Jun 09 2024

Examples

			Triangle begins as:
  1;
  1    1;
  1    3    1;
  1    6    5    1;
  1   10   15    7    1;
  1   15   35   28    9    1;
  1   21   70   84   45   11    1;
  1   28  126  210  165   66   13    1;
  1   36  210  462  495  286   91   15    1;
  1   45  330  924 1287 1001  455  120   17    1;
  1   55  495 1716 3003 3003 1820  680  153   19    1;
...
From _Philippe Deléham_, Mar 26 2012: (Start)
(0, 1, 0, 1, 0, 0, 0, ...) DELTA (1, 0, 1, -1, 0, 0, 0, ...) begins:
  1
  0, 1
  0, 1,  1
  0, 1,  3,   1
  0, 1,  6,   5,   1
  0, 1, 10,  15,   7,   1
  0, 1, 15,  35,  28,   9,  1
  0, 1, 21,  70,  84,  45, 11,  1
  0, 1, 28, 126, 210, 165, 66, 13, 1. (End)
		

Crossrefs

Programs

  • GAP
    Flat(List([0..12], n-> List([0..n], k-> Binomial(n+k, 2*k) ))); # G. C. Greubel, Aug 01 2019
  • Haskell
    a085478 n k = a085478_tabl !! n !! k
    a085478_row n = a085478_tabl !! n
    a085478_tabl = zipWith (zipWith a007318) a051162_tabl a025581_tabl
    -- Reinhard Zumkeller, Jun 22 2015
    
  • Magma
    [Binomial(n+k, 2*k): k in [0..n], n in [0..12]]; // G. C. Greubel, Aug 01 2019
    
  • Maple
    T := (n,k) -> binomial(n+k,2*k): seq(seq(T(n,k), k=0..n), n=0..11);
  • 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, 2 k], {n, 0, 12}, {k, 0, n}] // Flatten (* G. C. Greubel, Aug 01 2019 *)
    CoefficientList[Table[Fibonacci[2 n + 1, Sqrt[x]], {n, 0, 10}], x] // Flatten (* Eric W. Weisstein, Jul 03 2023 *)
    Join[{{1}}, CoefficientList[Table[MinimalPolynomial[-4 Sin[Pi/(4 n + 2)]^2, x], {n, 20}], x]] (* Eric W. Weisstein, Jul 12 2023 *)
  • PARI
    T(n,k) = binomial(n+k,n-k)
    
  • Sage
    [[binomial(n+k,2*k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Aug 01 2019
    

Formula

T(n, k) = (n+k)!/((n-k)!*(2*k)!).
G.f.: (1-z)/((1-z)^2-tz). - Emeric Deutsch, May 31 2004
Row sums are A001519 (Fibonacci(2n+1)). Diagonal sums are A011782. Binomial transform of A026729 (product of lower triangular matrices). - Paul Barry, Jun 21 2004
T(n, 0) = 1, T(n, k) = 0 if n=0} T(n-1-j, k-1)*(j+1). T(0, 0) = 1, T(0, k) = 0 if k>0; T(n, k) = T(n-1, k-1) + T(n-1, k) + Sum_{j>=0} (-1)^j*T(n-1, k+j)*A000108(j). For the column k, g.f.: Sum_{n>=0} T(n, k)*x^n = (x^k) / (1-x)^(2*k+1). - Philippe Deléham, Feb 15 2004
Sum_{k=0..n} T(n,k)*x^(2*k) = A000012(n), A001519(n+1), A001653(n), A078922(n+1), A007805(n), A097835(n), A097315(n), A097838(n), A078988(n), A097841(n), A097727(n), A097843(n), A097730(n), A098244(n), A097733(n), A098247(n), A097736(n), A098250(n), A097739(n), A098253(n), A097742(n), A098256(n), A097767(n), A098259(n), A097770(n), A098262(n), A097773(n), A098292(n), A097776(n) for x=0,1,2,...,27,28 respectively. - Philippe Deléham, Dec 31 2007
T(2*n,n) = A005809(n). - Philippe Deléham, Sep 17 2009
A183160(n) = Sum_{k=0..n} T(n,k)*T(n,n-k). - Paul D. Hanna, Dec 27 2010
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k). - Philippe Deléham, Feb 06 2012
O.g.f. for column k: x^k/(1-x)^(2*k+1), k >= 0. [See the o.g.f. of the triangle above, and a comment on compositions. - Wolfdieter Lang, Jul 09 2012]
E.g.f.: (2/sqrt(x + 4))*sinh((1/2)*t*sqrt(x + 4))*cosh((1/2)*t*sqrt(x)) = t + (1 + x)*t^3/3! + (1 + 3*x + x^2)*t^5/5! + (1 + 6*x + 5*x^2 + x^3)*t^7/7! + .... Cf. A091042. - Peter Bala, Jul 29 2013
T(n, k) = A065941(n+3*k, 4*k) = A108299(n+3*k, 4*k) = A194005(n+3*k, 4*k). - Johannes W. Meijer, Sep 05 2013
Sum_{k=0..n} (-1)^k*T(n,k)*A000108(k) = A000007(n) for n >= 0. - Werner Schulte, Jul 12 2017
Sum_{k=0..floor(n/2)} T(n-k,k)*A000108(k) = A001006(n) for n >= 0. - Werner Schulte, Jul 12 2017
From Peter Bala, Jun 26 2025: (Start)
The n-th row polynomial b(n, x) = (-1)^n * U(2*n, (i/2)*sqrt(x)), where U(n,x) is the n-th Chebyshev polynomial of the second kind.
b(n, x) = (-1)^n * Dir(n, -1 - x/2), where Dir(n, x) is the n-th row polynomial of the triangle A244419.
b(n, -1 - x) is the n-th row polynomial of A098493. (End)

A029653 Numbers in (2,1)-Pascal triangle (by row).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Reverse of A029635. Row sums are A003945. Diagonal sums are Fibonacci(n+2) = Sum_{k=0..floor(n/2)} (2n-3k)*C(n-k,n-2k)/(n-k). - Paul Barry, Jan 30 2005
Riordan array ((1+x)/(1-x), x/(1-x)). The signed triangle (-1)^(n-k)T(n,k) or ((1-x)/(1+x), x/(1+x)) is the inverse of A055248. Row sums are A003945. Diagonal sums are F(n+2). - Paul Barry, Feb 03 2005
Row sums = A003945: (1, 3, 6, 12, 24, 48, 96, ...) = (1, 3, 7, 15, 31, 63, 127, ...) - (0, 0, 1, 3, 7, 15, 31, ...); where (1, 3, 7, 15, ...) = A000225. - Gary W. Adamson, Apr 22 2007
Triangle T(n,k), read by rows, given by (2,-1,0,0,0,0,0,0,0,...) DELTA (1,0,0,0,0,0,0,0,0,...) where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 17 2011
A029653 is jointly generated with A208510 as an 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*v(n-1,x)+1. See the Mathematica section. - Clark Kimberling, Feb 28 2012
For a closed-form formula for arbitrary left and right borders of Pascal like triangle, see A228196. - Boris Putievskiy, Aug 18 2013
For a closed-form formula for generalized Pascal's triangle, see A228576. - Boris Putievskiy, Sep 04 2013
The n-th row polynomial is (2 + x)*(1 + x)^(n-1) for n >= 1. More generally, the n-th row polynomial of the Riordan array ( (1-a*x)/(1-b*x), x/(1-b*x) ) is (b - a + x)*(b + x)^(n-1) for n >= 1. - Peter Bala, Feb 25 2018

Examples

			The triangle T(n,k) begins:
n\k 0  1  2   3   4   5   6   7  8  9 10 ...
0:  1
1:  2  1
2:  2  3  1
3:  2  5  4   1
4:  2  7  9   5   1
5:  2  9 16  14   6   1
6:  2 11 25  30  20   7   1
7:  2 13 36  55  50  27   8   1
8:  2 15 49  91 105  77  35   9  1
9:  2 17 64 140 196 182 112  44 10  1
10: 2 19 81 204 336 378 294 156 54 11  1
... Reformatted. - _Wolfdieter Lang_, Jan 09 2015
With the array M(k) as defined in the Formula section, the infinite product M(0)*M(1)*M(2)*... begins
/1        \/1         \/1        \      /1        \
|2 1      ||0 1       ||0 1      |      |2 1      |
|2 1 1    ||0 2 1     ||0 0 1    |... = |2 3 1    |
|2 1 1 1  ||0 2 1 1   ||0 0 2 1  |      |2 5 4 1  |
|2 1 1 1 1||0 2 1 1 1 ||0 0 2 1 1|      |2 7 9 5 1|
|...      ||...       ||...      |      |...      |
- _Peter Bala_, Dec 27 2014
		

References

  • Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8.

Crossrefs

(d, 1) Pascal triangles: A007318(d=1), A093560(3), A093561(4), A093562(5), A093563(6), A093564(7), A093565(8), A093644(9), A093645(10).

Programs

  • Haskell
    a029653 n k = a029653_tabl !! n !! k
    a029653_row n = a029653_tabl !! n
    a029653_tabl = [1] : iterate
                   (\xs -> zipWith (+) ([0] ++ xs) (xs ++ [0])) [2, 1]
    -- Reinhard Zumkeller, Dec 16 2013
    
  • Maple
    A029653 :=  proc(n,k)
    if n = 0 then
      1;
    else
      binomial(n-1, k)+binomial(n, k)
    fi
    end proc: # R. J. Mathar, Jun 30 2013
  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := u[n - 1, x] + x*v[n - 1, x];
    v[n_, x_] := u[n - 1, x] + x*v[n - 1, x] + 1;
    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[%]  (* A208510 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]  (* A029653 *)
    (* Clark Kimberling, Feb 28 2012 *)
  • Python
    from sympy import Poly
    from sympy.abc import x
    def u(n, x): return 1 if n==1 else u(n - 1, x) + x*v(n - 1, x)
    def v(n, x): return 1 if n==1 else u(n - 1, x) + x*v(n - 1, x) + 1
    def a(n): return Poly(v(n, x), x).all_coeffs()[::-1]
    for n in range(1, 13): print(a(n)) # Indranil Ghosh, May 27 2017
    
  • Python
    from math import comb, isqrt
    def A029653(n): return comb(r:=(m:=isqrt(k:=n+1<<1))-(k<=m*(m+1)),a:=n-comb(r+1,2))*((r<<1)-a)//r if n else 1 # Chai Wah Wu, Nov 12 2024

Formula

T(n, k) = C(n-2, k-1) + C(n-2, k) + C(n-1, k-1) + C(n-1, k) except for n=0.
G.f.: (1 + x + y + xy)/(1 - y - xy). - Ralf Stephan, May 17 2004
T(n, k) = (2n-k)*binomial(n, n-k)/n, n, k > 0. - Paul Barry, Jan 30 2005
Sum_{k=0..n} T(n, k)*x^k gives A003945-A003954 for x = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10. - Philippe Deléham, Jul 10 2005
T(n, k) = C(n-1, k) + C(n, k). - Philippe Deléham, Jul 10 2005
Equals A097806 * A007318, i.e., the pairwise operator * Pascal's Triangle as infinite lower triangular matrices. - Gary W. Adamson, Apr 22 2007
From Peter Bala, Dec 27 2014: (Start)
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(2 + 5*x + 4*x^2/2! + x^3/3!) = 2 + 7*x + 16*x^2/2! + 30*x^3/3! + 50*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ).
Let M denote the lower unit triangular array with 1's on the main diagonal and 1's everywhere else below the main diagonal except for the first column which consists of the sequence [1,2,2,2,...]. For k = 0,1,2,... define M(k) to be the lower unit triangular block array
/I_k 0\
\ 0 M/ having the k X k identity matrix I_k as the upper left block; in particular, M(0) = M. Then the present triangle equals the infinite product M(0)*M(1)*M(2)*... (which is clearly well-defined). See the Example section. (End)

Extensions

More terms from James Sellers
Showing 1-10 of 43 results. Next