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

A011782 Coefficients of expansion of (1-x)/(1-2*x) in powers of x.

Original entry on oeis.org

1, 1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096, 8192, 16384, 32768, 65536, 131072, 262144, 524288, 1048576, 2097152, 4194304, 8388608, 16777216, 33554432, 67108864, 134217728, 268435456, 536870912, 1073741824, 2147483648, 4294967296, 8589934592
Offset: 0

Views

Author

Lee D. Killough (killough(AT)wagner.convex.com)

Keywords

Comments

Apart from initial term, same as A000079 (powers of 2).
Number of compositions (ordered partitions) of n. - Toby Bartels, Aug 27 2003
Number of ways of putting n unlabeled items into (any number of) labeled boxes where every box contains at least one item. Also "unimodal permutations of n items", i.e., those which rise then fall. (E.g., for three items: ABC, ACB, BCA and CBA are unimodal.) - Henry Bottomley, Jan 17 2001
Number of permutations in S_n avoiding the patterns 213 and 312. - Tuwani Albert Tshifhumulo, Apr 20 2001. More generally (see Simion and Schmidt), the number of permutations in S_n avoiding (i) the 123 and 132 patterns; (ii) the 123 and 213 patterns; (iii) the 132 and 213 patterns; (iv) the 132 and 231 patterns; (v) the 132 and 312 patterns; (vi) the 213 and 231 patterns; (vii) the 213 and 312 patterns; (viii) the 231 and 312 patterns; (ix) the 231 and 321 patterns; (x) the 312 and 321 patterns.
a(n+2) is the number of distinct Boolean functions of n variables under action of symmetric group.
Number of unlabeled (1+2)-free posets. - Detlef Pauly, May 25 2003
Image of the central binomial coefficients A000984 under the Riordan array ((1-x), x*(1-x)). - Paul Barry, Mar 18 2005
Binomial transform of (1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...); inverse binomial transform of A007051. - Philippe Deléham, Jul 04 2005
Also, number of rationals in [0, 1) whose binary expansions terminate after n bits. - Brad Chalfan, May 29 2006
Equals row sums of triangle A144157. - Gary W. Adamson, Sep 12 2008
Prepend A089067 with a 1, getting (1, 1, 3, 5, 13, 23, 51, ...) as polcoeff A(x); then (1, 1, 2, 4, 8, 16, ...) = A(x)/A(x^2). - Gary W. Adamson, Feb 18 2010
An elephant sequence, see A175655. For the central square four A[5] vectors, with decimal values 2, 8, 32 and 128, lead to this sequence. For the corner squares these vectors lead to the companion sequence A094373. - Johannes W. Meijer, Aug 15 2010
From Paul Curtz, Jul 20 2011: (Start)
Array T(m,n) = 2*T(m,n-1) + T(m-1,n):
1, 1, 2, 4, 8, 16, ... = a(n)
1, 3, 8, 20, 48, 112, ... = A001792,
1, 5, 18, 56, 160, 432, ... = A001793,
1, 7, 32, 120, 400, 1232, ... = A001794,
1, 9, 50, 220, 840, 2912, ... = A006974, followed with A006975, A006976, gives nonzero coefficients of Chebyshev polynomials of first kind A039991 =
1,
1, 0,
2, 0, -1,
4, 0, -3, 0,
8, 0, -8, 0, 1.
T(m,n) third vertical: 2*n^2, n positive (A001105).
Fourth vertical appears in Janet table even rows, last vertical (A168342 array, A138509, rank 3, 13, = A166911)). (End)
A131577(n) and differences are:
0, 1, 2, 4, 8, 16,
1, 1, 2, 4, 8, 16, = a(n),
0, 1, 2, 4, 8, 16,
1, 1, 2, 4, 8, 16.
Number of 2-color necklaces of length 2n equal to their complemented reversal. For length 2n+1, the number is 0. - David W. Wilson, Jan 01 2012
Edges and also central terms of triangle A198069: a(0) = A198069(0,0) and for n > 0: a(n) = A198069(n,0) = A198069(n,2^n) = A198069(n,2^(n-1)). - Reinhard Zumkeller, May 26 2013
These could be called the composition numbers (see the second comment) since the equivalent sequence for partitions is A000041, the partition numbers. - Omar E. Pol, Aug 28 2013
Number of self conjugate integer partitions with exactly n parts for n>=1. - David Christopher, Aug 18 2014
The sequence is the INVERT transform of (1, 0, 1, 1, 2, 3, 5, 8, 13, 21, 34, ...). - Gary W. Adamson, Jul 16 2015
Number of threshold graphs on n nodes [Hougardy]. - Falk Hüffner, Dec 03 2015
Number of ternary words of length n in which binary subwords appear in the form 10...0. - Milan Janjic, Jan 25 2017
a(n) is the number of words of length n over an alphabet of two letters, of which one letter appears an even number of times (the empty word of length 0 is included). See the analogous odd number case in A131577, and the Balakrishnan reference in A006516 (the 4-letter odd case), pp. 68-69, problems 2.66, 2.67 and 2.68. - Wolfdieter Lang, Jul 17 2017
Number of D-equivalence classes of Łukasiewicz paths. Łukasiewicz paths are D-equivalent iff the positions of pattern D are identical in these paths. - Sergey Kirgizov, Apr 08 2018
Number of color patterns (set partitions) for an oriented row of length n using two or fewer colors (subsets). Two color patterns are equivalent if we permute the colors. For a(4)=8, the 4 achiral patterns are AAAA, AABB, ABAB, and ABBA; the 4 chiral patterns are the 2 pairs AAAB-ABBB and AABA-ABAA. - Robert A. Russell, Oct 30 2018
The determinant of the symmetric n X n matrix M defined by M(i,j) = (-1)^max(i,j) for 1 <= i,j <= n is equal to a(n) * (-1)^(n*(n+1)/2). - Bernard Schott, Dec 29 2018
For n>=1, a(n) is the number of permutations of length n whose cyclic representations can be written in such a way that when the cycle parentheses are removed what remains is 1 through n in natural order. For example, a(4)=8 since there are exactly 8 permutations of this form, namely, (1 2 3 4), (1)(2 3 4), (1 2)(3 4), (1 2 3)(4), (1)(2)(3 4), (1)(2 3)(4), (1 2)(3)(4), and (1)(2)(3)(4). Our result follows readily by conditioning on k, the number of parentheses pairs of the form ")(" in the cyclic representation. Since there are C(n-1,k) ways to insert these in the cyclic representation and since k runs from 0 to n-1, we obtain a(n) = Sum_{k=0..n-1} C(n-1,k) = 2^(n-1). - Dennis P. Walsh, May 23 2020
Maximum number of preimages that a permutation of length n + 1 can have under the consecutive-231-avoiding stack-sorting map. - Colin Defant, Aug 28 2020
a(n) is the number of occurrences of the empty set {} in the von Neumann ordinals from 0 up to n. Each ordinal k is defined as the set of all smaller ordinals: 0 = {}, 1 = {0}, 2 = {0,1}, etc. Since {} is the foundational element of all ordinals, the total number of times it appears grows as powers of 2. - Kyle Wyonch, Mar 30 2025

Examples

			G.f. = 1 + x + 2*x^2 + 4*x^3 + 8*x^4 + 16*x^5 + 32*x^6 + 64*x^7 + 128*x^8 + ...
    ( -1   1  -1)
det (  1   1   1)  = 4
    ( -1  -1  -1)
		

References

  • Mohammad K. Azarian, A Generalization of the Climbing Stairs Problem, Mathematics and Computer Education Journal, Vol. 31, No. 1, pp. 24-28, Winter 1997.
  • S. Kitaev, Patterns in Permutations and Words, Springer-Verlag, 2011. see p. 399 Table A.7
  • Xavier Merlin, Methodix Algèbre, Ellipses, 1995, p. 153.

Crossrefs

Sequences with g.f.'s of the form ((1-x)/(1-2*x))^k: this sequence (k=1), A045623 (k=2), A058396 (k=3), A062109 (k=4), A169792 (k=5), A169793 (k=6), A169794 (k=7), A169795 (k=8), A169796 (k=9), A169797 (k=10).
Cf. A005418 (unoriented), A122746(n-3) (chiral), A016116 (achiral).
Row sums of triangle A100257.
A row of A160232.
Row 2 of A278984.

Programs

  • Haskell
    a011782 n = a011782_list !! n
    a011782_list = 1 : scanl1 (+) a011782_list
    -- Reinhard Zumkeller, Jul 21 2013
    
  • Magma
    [Floor((1+2^n)/2): n in [0..35]]; // Vincenzo Librandi, Aug 21 2011
    
  • Maple
    A011782:= n-> ceil(2^(n-1)): seq(A011782(n), n=0..50); # Wesley Ivan Hurt, Feb 21 2015
    with(PolynomialTools):  A011782:=seq(coeftayl((1-x)/(1-2*x), x = 0, k),k=0..10^2); # Muniru A Asiru, Sep 26 2017
  • Mathematica
    f[s_] := Append[s, Ceiling[Plus @@ s]]; Nest[f, {1}, 32] (* Robert G. Wilson v, Jul 07 2006 *)
    CoefficientList[ Series[(1-x)/(1-2x), {x, 0, 32}], x] (* Robert G. Wilson v, Jul 07 2006 *)
    Table[Sum[StirlingS2[n, k], {k,0,2}], {n, 0, 30}] (* Robert A. Russell, Apr 25 2018 *)
    Join[{1},NestList[2#&,1,40]] (* Harvey P. Dale, Dec 06 2018 *)
  • PARI
    {a(n) = if( n<1, n==0, 2^(n-1))};
    
  • PARI
    Vec((1-x)/(1-2*x) + O(x^30)) \\ Altug Alkan, Oct 31 2015
    
  • Python
    def A011782(n): return 1 if n == 0 else 2**(n-1) # Chai Wah Wu, May 11 2022
  • Sage
    [sum(stirling_number2(n,j) for j in (0..2)) for n in (0..35)] # G. C. Greubel, Jun 02 2020
    

Formula

a(0) = 1, a(n) = 2^(n-1).
G.f.: (1 - x) / (1 - 2*x) = 1 / (1 - x / (1 - x)). - Michael Somos, Apr 18 2012
E.g.f.: cosh(z)*exp(z) = (exp(2*z) + 1)/2.
a(0) = 1 and for n>0, a(n) = sum of all previous terms.
a(n) = Sum_{k=0..n} binomial(n, 2*k). - Paul Barry, Feb 25 2003
a(n) = Sum_{k=0..n} binomial(n,k)*(1+(-1)^k)/2. - Paul Barry, May 27 2003
a(n) = floor((1+2^n)/2). - Toby Bartels (toby+sloane(AT)math.ucr.edu), Aug 27 2003
G.f.: Sum_{i>=0} x^i/(1-x)^i. - Jon Perry, Jul 10 2004
a(n) = Sum_{k=0..n} (-1)^(n-k)*binomial(k+1, n-k)*binomial(2*k, k). - Paul Barry, Mar 18 2005
a(n) = Sum_{k=0..floor(n/2)} A055830(n-k, k). - Philippe Deléham, Oct 22 2006
a(n) = Sum_{k=0..n} A098158(n,k). - Philippe Deléham, Dec 04 2006
G.f.: 1/(1 - (x + x^2 + x^3 + ...)). - Geoffrey Critzer, Aug 30 2008
a(n) = A000079(n) - A131577(n).
a(n) = A173921(A000079(n)). - Reinhard Zumkeller, Mar 04 2010
a(n) = Sum_{k=2^n..2^(n+1)-1} A093873(k)/A093875(k), sums of rows of the full tree of Kepler's harmonic fractions. - Reinhard Zumkeller, Oct 17 2010
E.g.f.: (exp(2*x)+1)/2 = (G(0) + 1)/2; G(k) = 1 + 2*x/(2*k+1 - x*(2*k+1)/(x + (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Dec 03 2011
A051049(n) = p(n+1) where p(x) is the unique degree-n polynomial such that p(k) = a(k) for k = 0, 1, ..., n. - Michael Somos, Apr 18 2012
A008619(n) = p(-1) where p(x) is the unique degree-n polynomial such that p(k) = a(k) for k = 0, 1, ..., n. - Michael Somos, Apr 18 2012
INVERT transform is A122367. MOBIUS transform is A123707. EULER transform of A059966. PSUM transform is A000079. PSUMSIGN transform is A078008. BINOMIAL transform is A007051. REVERT transform is A105523. A002866(n) = a(n)*n!. - Michael Somos, Apr 18 2012
G.f.: U(0), where U(k) = 1 + x*(k+3) - x*(k+2)/U(k+1); (continued fraction, 1-step). - Sergei N. Gladkovskii, Oct 10 2012
a(n) = A000041(n) + A056823(n). - Omar E. Pol, Aug 31 2013
E.g.f.: E(0), where E(k) = 1 + x/( 2*k+1 - x/E(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Dec 25 2013
G.f.: 1 + x/(1 + x)*( 1 + 3*x/(1 + 3*x)*( 1 + 5*x/(1 + 5*x)*( 1 + 7*x/(1 + 7*x)*( 1 + ... )))). - Peter Bala, May 27 2017
a(n) = Sum_{k=0..2} stirling2(n, k).
G.f.: Sum_{j=0..k} A248925(k,j)*x^j / Product_{j=1..k} 1-j*x with k=2. - Robert A. Russell, Apr 25 2018
a(n) = A053120(n, n), n >= 0, (main diagonal of triangle of Chebyshev's T polynomials). - Wolfdieter Lang, Nov 26 2019

Extensions

Additional comments from Emeric Deutsch, May 14 2001
Typo corrected by Philippe Deléham, Oct 25 2008

A001792 a(n) = (n+2)*2^(n-1).

Original entry on oeis.org

1, 3, 8, 20, 48, 112, 256, 576, 1280, 2816, 6144, 13312, 28672, 61440, 131072, 278528, 589824, 1245184, 2621440, 5505024, 11534336, 24117248, 50331648, 104857600, 218103808, 452984832, 939524096, 1946157056, 4026531840, 8321499136, 17179869184, 35433480192
Offset: 0

Views

Author

Keywords

Comments

Number of parts in all compositions (ordered partitions) of n + 1. For example, a(2) = 8 because in 3 = 2 + 1 = 1 + 2 = 1 + 1 + 1 we have 8 parts. Also number of compositions (ordered partitions) of 2n + 1 with exactly 1 odd part. For example, a(2) = 8 because the only compositions of 5 with exactly 1 odd part are 5 = 1 + 4 = 2 + 3 = 3 + 2 = 4 + 1 = 1 + 2 + 2 = 2 + 1 + 2 = 2 + 2 + 1. - Emeric Deutsch, May 10 2001
Binomial transform of natural numbers [1, 2, 3, 4, ...].
For n >= 1 a(n) is also the determinant of the n X n matrix with 3's on the diagonal and 1's elsewhere. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 06 2001
The arithmetic mean of first n terms of the sequence is 2^(n-1). - Amarnath Murthy, Dec 25 2001, corrected by M. F. Hasler, Dec 17 2016
Also the number of "winning paths" of length n across an n X n Hex board. Satisfies the recursion a(n) = 2a(n-1) + 2^(n-2). - David Molnar (molnar(AT)stolaf.edu), Apr 10 2002
Diagonal in A053218. - Benoit Cloitre, May 08 2002
Let M_n be the n X n matrix m_(i, j) = 1 + abs(i-j) then det(M_n) = (-1)^(n-1)*a(n-1). - Benoit Cloitre, May 28 2002
Absolute value of determinant of n X n matrix of form: [1 2 3 4 5 / 2 1 2 3 4 / 3 2 1 2 3 / 4 3 2 1 2 / 5 4 3 2 1]. - Benoit Cloitre, Jun 20 2002
Number of ones in all (n+1)-bit integers (cf. A000120). - Ralf Stephan, Aug 02 2003
This sequence also emerges as a floretion force transform of powers of 2 (see program code). Define a(-1) = 0 (as the sequence is returned by FAMP). Then a(n-1) + A098156(n+1) = 2*a(n) (conjecture). - Creighton Dement, Mar 14 2005
This sequence gives the absolute value of the determinant of the Toeplitz matrix with first row containing the first n integers. - Paul Max Payton, May 23 2006
Equals sums of rows right of left edge of A102363 divided by three, + 2^K. - David G. Williams (davidwilliams(AT)paxway.com), Oct 08 2007
If X_1, X_2, ..., X_n are 2-blocks of a (2n+1)-set X then, for n >= 1, a(n) is the number of (n+1)-subsets of X intersecting each X_i, (i = 1, 2, ..., n). - Milan Janjic, Nov 18 2007
Also, a(n-1) is the determinant of the n X n matrix with A[i, j] = n - |i-j|. - M. F. Hasler, Dec 17 2008
1/2 the number of permutations of 1 .. (n+2) arranged in a circle with exactly one local maximum. - R. H. Hardin, Apr 19 2009
The first corrector line for transforming 2^n offset 0 with a leading 1 into the Fibonacci sequence. - Al Hakanson (hawkuu(AT)gmail.com), Jun 01 2009
a(n) is the number of runs of consecutive 1's in all binary sequences of length (n+1). - Geoffrey Critzer, Jul 02 2009
Let X_j (0 < j <= 2^n) all the subsets of N_n; m(i, j) := if {i} in X_j then 1 else 0. Let A = transpose(M).M; Then a(i, j) = (number of elements)|X_i intersect X_j|. Determinant(X*I-A) = (X-(n+1)*2^(n-2))*(X-2^(n-2))^(n-1)*X^(2^n-n).
Eigenvector for (n+1)*2^(n-2) is V_i=|X_i|.
Sum_{k=1..2^n} |X_i intersect X_k|*|X_k| = (n+1)*2^(n-2)*|X_i|.
Eigenvectors for 2^(n-2) are {line(M)[i] - line(M)[j], 1 <= i, j <= n}. - CLARISSE Philippe (clarissephilippe(AT)yahoo.fr), Mar 24 2010
The sequence b(n) = 2*A001792(n), for n >= 1 with b(0) = 1, is an elephant sequence, see A175655. For the central square four A[5] vectors, with decimal values 187, 190, 250 and 442, lead to the b(n) sequence. For the corner squares these vectors lead to the companion sequence A134401. - Johannes W. Meijer, Aug 15 2010
Equals partial sums of A045623: (1, 2, 5, 12, 28, ...); where A045623 = the convolution square of (1, 1, 2, 4, 8, 16, 32, ...). - Gary W. Adamson, Oct 26 2010
Equals (1, 2, 4, 8, 16, ...) convolved with (1, 1, 2, 4, 8, 16, ...); e.g., a(3) = 20 = (1, 1, 2, 4) dot (8, 4, 2, 1) = (8 + 4 + 4 + 4). - Gary W. Adamson, Oct 26 2010
This sequence seems to give the first x+1 nonzero terms in the sequence derived by subtracting the m-th term in the x_binacci sequence (where the first term is one and the y-th term is the sum of x terms immediately preceding it) from 2^(m-2). - Dylan Hamilton, Nov 03 2010
Recursive formulas for a(n) are in many cases derivable from its property wherein delta^k(a(n)) - a(n) = k*2^n where delta^k(a(n)) represents the k-th forward difference of a(n). Provable with a difference table and a little induction. - Ethan Beihl, May 02 2011
Let f(n,k) be the sum of numbers in the subsets of size k of {1, 2, ..., n}. Then a(n-1) is the average of the numbers f(n, 0), ... f(n, n). Example: (f(3, 1), f(3, 2), f(3, 3)) = (6, 12, 6), with average (6+12+6)/3. - Clark Kimberling, Feb 24 2012
a(n) is the number of length-2n binary sequences that contain a subsequence of ones with length n or more. To derive this result, note that there are 2^n sequences where the initial one of the subsequence occurs at entry one. If the initial one of the subsequence occurs at entry 2, 3, ..., or n + 1, there are 2^(n-1) sequences since a zero must precede the initial one. Hence a(n) = 2^n + n*2^(n-1)=(n+2)2^(n-1). An example is given in the example section below. - Dennis P. Walsh, Oct 25 2012
As the total number of parts in all compositions of n+1 (see the first line in Comments) the equivalent sequence for partitions is A006128. On the other hand, as the first differences of A001787 (see the first line in Crossrefs) the equivalent sequence for partitions is A138879. - Omar E. Pol, Aug 28 2013
a(n) is the number of spanning trees of the complete tripartite graph K_{n,1,1}. - James Mahoney, Oct 24 2013
a(n-1) = denominator of the mean (2n/(n+1), after reduction), of the compositions of n; numerator is given by A022998(n). - Clark Kimberling, Mar 11 2014
From Tom Copeland, Nov 09 2014: (Start)
The shifted array belongs to an interpolated family of arrays associated to the Catalan A000108 (t=1), and Riordan, or Motzkin sums A005043 (t=0), with the interpolating o.g.f. (1-sqrt(1-4x/(1+(1-t)x)))/2 and inverse x(1-x)/(1+(t-1)x(1-x)). See A091867 for more info on this family. Here the interpolation is t=-3 (mod signs in the results).
Let C(x) = (1 - sqrt(1-4x))/2, an o.g.f. for the Catalan numbers A000108, with inverse Cinv(x) = x*(1-x) and P(x,t) = x/(1+t*x) with inverse P(x,-t).
Shifted o.g.f: G(x) = x*(1-x)/(1 - 4x*(1-x)) = P[Cinv(x),-4].
Inverse o.g.f: Ginv(x) = [1 - sqrt(1 - 4*x/(1+4x))]/2 = C[P(x, 4)] (signed shifted A001700). Cf. A030528. (End)
For n > 0, element a(n) of the sequence is equal to the gradients of the (n-1)-th row of Pascal triangle multiplied with the square of the integers from n+1,...,1. I.e., row 3 of Pascal's triangle 1,3,3,1 has gradients 1,2,0,-2,-1, so a(4) = 1*(5^2) + 2*(4^2) + 0*(3^2) - 2*(2^2) - 1*(1^2) = 48. - Jens Martin Carlsson, May 18 2017
Number of self-avoiding paths connecting all the vertices of a convex (n+2)-gon. - Ivaylo Kortezov, Jan 19 2020
a(n-1) is the total number of elements of subsets of {1,2,..,n} that contain n. For example, for n = 3, a(2) = 8, and the subsets of {1,2,3} that contain 3 are {3}, {1,3}, {2,3}, {1,2,3}, with a total of 8 elements. - Enrique Navarrete, Aug 01 2020

Examples

			a(0) = 1, a(1) = 2*1 + 1 = 3, a(2) = 2*3 + 2 = 8, a(3) = 2*8 + 4 = 20, a(4) = 2*20 + 8 = 48, a(5) = 2*48 + 16 = 112, a(6) = 2*112 + 32 = 256, ... - _Philippe Deléham_, Apr 19 2009
a(2) = 8 since there are 8 length-4 binary sequences with a subsequence of ones of length 2 or more, namely, 1111, 1110, 1101, 1011, 0111, 1100, 0110, and 0011. - _Dennis P. Walsh_, Oct 25 2012
G.f. = 1 + 3*x + 8*x^2 + 20*x^3 + 48*x^4 + 112*x^5 + 256*x^6 + 576*x^7 + ...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 795.
  • 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).
  • A. M. Stepin and A. T. Tagi-Zade, Words with restrictions, pp. 67-74 of Kvant Selecta: Combinatorics I, Amer. Math. Soc., 2001 (G_n on p. 70).

Crossrefs

First differences of A001787.
a(n) = A049600(n, 1), a(n) = A030523(n + 1, 1).
Cf. A053113.
Row sums of triangles A008949 and A055248.
a(n) = -A039991(n+2, 2).
If the exponent E in a(n) = Sum_{m=0..n} (Sum_{k=0..m} C(n,k))^E is 1, 2, 3, 4, 5 we get A001792, A003583, A007403, A294435, A294436 respectively.

Programs

  • GAP
    List([0..35],n->(n+2)*2^(n-1)); # Muniru A Asiru, Sep 25 2018
    
  • Haskell
    a001792 n = a001792_list !! n
    a001792_list = scanl1 (+) a045623_list
    -- Reinhard Zumkeller, Jul 21 2013
    
  • Magma
    [(n+2)*2^(n-1): n in [0..40]]; // Vincenzo Librandi, Nov 10 2014
    
  • Maple
    A001792 := n-> (n+2)*2^(n-1);
    spec := [S, {B=Set(Z, 0 <= card), S=Prod(Z, B, B)}, labeled]: seq(combstruct[count](spec, size=n)/4, n=2..30); # Zerinvary Lajos, Oct 09 2006
    A001792:=-(-3+4*z)/(2*z-1)^2; # Simon Plouffe in his 1992 dissertation, which gives the sequence without the initial 1
    G(x):=1/exp(2*x)*(1-x): f[0]:=G(x): for n from 1 to 54 do f[n]:=diff(f[n-1],x) od: x:=0: seq(abs(f[n]),n=0..28 ); # Zerinvary Lajos, Apr 17 2009
    a := n -> hypergeom([-n, 2], [1], -1);
    seq(round(evalf(a(n),32)), n=0..31); # Peter Luschny, Aug 02 2014
  • Mathematica
    matrix[n_Integer /; n >= 1] := Table[Abs[p - q] + 1, {q, n}, {p, n}]; a[n_Integer /; n >= 1] := Abs[Det[matrix[n]]] (* Josh Locker (joshlocker(AT)macfora.com), Apr 29 2004 *)
    g[n_,m_,r_] := Binomial[n - 1, r - 1] Binomial[m + 1, r] r; Table[1 + Sum[g[n, k - n, r], {r, 1, k}, {n, 1, k - 1}], {k, 1, 29}] (* Geoffrey Critzer, Jul 02 2009 *)
    a[n_] := (n + 2)*2^(n - 1); a[Range[0, 40]] (* Vladimir Joseph Stephan Orlovsky, Feb 09 2011 *)
    LinearRecurrence[{4, -4}, {1, 3}, 40] (* Harvey P. Dale, Aug 29 2011 *)
    CoefficientList[Series[(1 - x) / (1 - 2 x)^2, {x, 0, 40}], x] (* Vincenzo Librandi, Nov 10 2014 *)
    b[i_]:=i; a[n_]:=Abs[Det[ToeplitzMatrix[Array[b, n], Array[b, n]]]]; Array[a, 40] (* Stefano Spezia, Sep 25 2018 *)
    a[n_]:=Hypergeometric2F1[2,-n+1,1,-1];Array[a,32] (* Giorgos Kalogeropoulos, Jan 04 2022 *)
  • PARI
    A001792(n)=(n+2)<<(n-1) \\ M. F. Hasler, Dec 17 2008
    
  • Python
    for n in range(0,40): print(int((n+2)*2**(n-1)), end=' ') # Stefano Spezia, Oct 16 2018

Formula

a(n) = (n+2)*2^(n-1).
G.f.: (1 - x)/(1 - 2*x)^2 = 2F1(1,3;2;2x).
a(n) = 4*a(n-1) - 4*a(n-2).
G.f. (-1 + (1-2*x)^(-2))/(x*2^2). - Wolfdieter Lang
a(n) = A018804(2^n). - Matthew Vandermast, Mar 01 2003
a(n) = Sum_{k=0..n+2} binomial(n+2, 2k)*k. - Paul Barry, Mar 06 2003
a(n) = (1/4)*A001787(n+2). - Emeric Deutsch, May 24 2003
With a leading 0, this is ((n+1)2^n - 0^n)/4 = Sum_{m=0..n} binomial(n - 1, m - 1)*m, the binomial transform of A004526(n+1). - Paul Barry, Jun 05 2003
a(n) = Sum_{k=0..n} binomial(n, k)*(k + 1). - Lekraj Beedassy, Jun 24 2004
a(n) = A000244(n) - A066810(n). - Ross La Haye, Apr 29 2006
Row sums of triangle A130585. - Gary W. Adamson, Jun 06 2007
Equals A125092 * [1/1, 1/2, 1/3, ...]. - Gary W. Adamson, Nov 16 2007
a(n) = (n+1)*2^n - n*2^(n-1). Equals A128064 * A000079. - Gary W. Adamson, Dec 28 2007
G.f.: F(3, 1; 2; 2x). - Paul Barry, Sep 03 2008
a(n) = 1 + Sum_{k=1..n} (n - k + 4)2^(n - k - 1). This follows from the result that the number of parts equal to k in all compositions of n is (n - k + 3)2^(n - k - 2) for 0 < k < n. - Geoffrey Critzer, Sep 21 2008
a(n) = 2^(n-1) + 2 a(n-1) ; a(n-1) = det(n - |i - j|){i, j = 1..n}. - _M. F. Hasler, Dec 17 2008
a(n) = 2*a(n-1) + 2^(n-1). - Philippe Deléham, Apr 19 2009
a(n) = A164910(2^n). - Gary W. Adamson, Aug 30 2009
a(n) = Sum_{i=1..2^n} gcd(i, 2^n) = A018804(2^n). So we have: 2^0 * phi(2^n) + ... + 2^n * phi(2^0) = (n + 2)*2^(n-1), where phi is the Euler totient function. - Jeffrey R. Goodwin, Nov 11 2011
a(n) = Sum_{j=0..n} Sum_{i=0..n} binomial(n, i + j). - Yalcin Aktar, Jan 17 2012
Eigensequence of an infinite lower triangular matrix with 2^n as the left border and the rest 1's. - Gary W. Adamson, Jan 30 2012
G.f.: 1 + 2*x*U(0) where U(k) = 1 + (k + 1)/(2 - 8*x/(4*x + (k + 1)/U(k + 1))); (continued fraction, 3 - step). - Sergei N. Gladkovskii, Oct 19 2012
a(n) = Sum_{k=0..n} Sum_{j=0..k} binomial(n,j). - Peter Luschny, Dec 03 2013
a(n) = Hyper2F1([-n, 2], [1], -1). - Peter Luschny, Aug 02 2014
G.f.: 1 / (1 - 3*x / (1 + x / (3 - 4*x))). - Michael Somos, Aug 26 2015
a(n) = -A053120(2+n, n), n >= 0, the negative of the third (sub)diagonal of the triangle of Chebyshev's T polynomials. - Wolfdieter Lang, Nov 26 2019
From Amiram Eldar, Jan 12 2021: (Start)
Sum_{n>=0} 1/a(n) = 8*log(2) - 4.
Sum_{n>=0} (-1)^n/a(n) = 4 - 8*log(3/2). (End)
E.g.f.: exp(2*x)*(1 + x). - Stefano Spezia, Jun 11 2021

A053120 Triangle of coefficients of Chebyshev's T(n,x) polynomials (powers of x in increasing order).

Original entry on oeis.org

1, 0, 1, -1, 0, 2, 0, -3, 0, 4, 1, 0, -8, 0, 8, 0, 5, 0, -20, 0, 16, -1, 0, 18, 0, -48, 0, 32, 0, -7, 0, 56, 0, -112, 0, 64, 1, 0, -32, 0, 160, 0, -256, 0, 128, 0, 9, 0, -120, 0, 432, 0, -576, 0, 256, -1, 0, 50, 0, -400, 0, 1120, 0, -1280, 0, 512, 0, -11, 0, 220, 0, -1232, 0, 2816, 0, -2816, 0, 1024
Offset: 0

Views

Author

Keywords

Comments

Row sums (signed triangle): A000012 (powers of 1). Row sums (unsigned triangle): A001333(n).
From Wolfdieter Lang, Oct 21 2013: (Start)
The row polynomials T(n,x) equal (S(n,2*x) - S(n-2,2*x))/2, n >= 0, with the row polynomials S from A049310, with S(-1,x) = 0, and S(-2,x) = -1.
The zeros of T(n,x) are x(n,k) = cos((2*k+1)*Pi/(2*n)), k = 0, 1, ..., n-1, n >= 1. (End)
From Wolfdieter Lang, Jan 03 2020 and Paul Weisenhorn: (Start)
The (sub)diagonal sequences {D_{2*k}(m)}{m >= 0}, for k >= 0, have o.g.f. GD{2*k}(x) = (-1)^k*(1-x)/(1-2*x)^(k+1), for k >= 0, and GD_{2*k+1}(x) = 0, for k >= 0. This follows from their o.g.f. GGD(z, x) := Sum_{k>=0} GD_k(x)*z^n which is obtained from the o.g.f. of the T-triangle GT(z, x) = (1-x*z)/(1 - 2*x + z^2) (see the formula section) by GGD(z, x) = GT(z, x/z).
The explicit form is then D_{2*k}(m) = (-1)^k, for m = 0, and
(-1)^k*(2*k+m)*2^(m-1)*risefac(k+1, m-1)/m!, for m >= 1, with the rising factorial risefac(x, n). (End)

Examples

			The triangle a(n,m) begins:
n\m  0  1   2    3     4    5     6     7      8    9   10...
0:   1
1:   0  1
2:  -1  0   2
3:   0 -3   0    4
4:   1  0  -8    0     8
5:   0  5   0  -20     0   16
6:  -1  0  18    0   -48    0    32
7:   0 -7   0   56     0 -112     0    64
8:   1  0 -32    0   160    0  -256     0    128
9:   0  9   0 -120     0  432     0  -576      0  256
10: -1  0  50    0  -400    0  1120     0  -1280    0  512
... Reformatted and extended - _Wolfdieter Lang_, Oct 21 2013
E.g., the fourth row (n=3) corresponds to the polynomial T(3,x) = -3*x + 4*x^3.
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964. Tenth printing, Wiley, 2002 (also electronically available), p. 795.
  • F. Hirzebruch et al., Manifolds and Modular Forms, Vieweg 1994 pp. 77, 105.
  • 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.
  • TableCurve 2D, Automated curve fitting and equation discovery, Version 5.01 for Windows, User's Manual, Chebyshev Series Polynomials and Rationals, pages 12-21 - 12-24, SYSTAT Software, Inc., Richmond, WA, 2002.

Crossrefs

The first nonzero (sub)diagonal sequences are A011782, -A001792, A001793(n+1), -A001794, A006974, -A006975, A006976, -A209404.

Programs

  • Julia
    using Nemo
    function A053120Row(n)
        R, x = PolynomialRing(ZZ, "x")
        p = chebyshev_t(n, x)
        [coeff(p, j) for j in 0:n] end
    for n in 0:6 A053120Row(n) |> println end # Peter Luschny, Mar 13 2018
    
  • Magma
    &cat[ Coefficients(ChebyshevT(n)): n in [0..11] ]; // Klaus Brockhaus, Mar 08 2008
    
  • Maple
    with(orthopoly) ;
    A053120 := proc(n,k)
        T(n,x) ;
        coeftayl(%,x=0,k) ;
    end proc: # R. J. Mathar, Jun 30 2013
    T := (n, x) -> `if`(n = 0, 1, add((-1)^(n - k) * (n/(2*k))*binomial(k, n - k) *(2*x)^(2*k - n), k = 1 ..n)):
    seq(seq(coeff(T(n, x), x, k), k = 0..n), n = 0..11); # Peter Luschny, Sep 20 2022
  • Mathematica
    t[n_, k_] := Coefficient[ ChebyshevT[n, x], x, k]; Flatten[ Table[ t[n, k], {n, 0, 11}, {k, 0, n}]] (* Jean-François Alcover, Jan 16 2012 *)
  • PARI
    for(n=0,5,P=polchebyshev(n);for(k=0,n,print1(polcoeff(P,k)", "))) \\ Charles R Greathouse IV, Jan 16 2012
    
  • SageMath
    def f(n,k): # f = A039991
        if (n<2 and k==0): return 1
        elif (k<0 or k>n): return 0
        else: return 2*f(n-1, k) - f(n-2, k-2)
    def A053120(n,k): return f(n, n-k)
    flatten([[A053120(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Aug 10 2022

Formula

T(n, m) = A039991(n, n-m).
G.f. for row polynomials T(n,x) (signed triangle): (1-x*z)/(1-2*x*z+z^2). If unsigned: (1-x*z)/(1-2*x*z-z^2).
T(n, m) := 0 if n < m or n+m odd; T(n, m) = (-1)^(n/2) if m=0 (n even); otherwise T(n, m) = ((-1)^((n+m)/2 + m))*(2^(m-1))*n*binomial((n+m)/2-1, m-1)/m.
Recursion for n >= 2: T(n, m) = T*a(n-1, m-1) - T(n-2, m), T(n, m)=0 if n < m, T(n, -1) := 0, T(0, 0) = T(1, 1) = 1.
G.f. for m-th column (signed triangle): 1/(1+x^2) if m=0, otherwise (2^(m-1))*(x^m)*(1-x^2)/(1+x^2)^(m+1).
From G. C. Greubel, Aug 10 2022: (Start)
Sum_{k=0..floor(n/2)} T(n-k, k) = A000007(n).
T(2*n, n) = i^n * A036909(n/2) * (1+(-1)^n)/2 + [n=0]/3. (End)
T(n, k) = [x^k] T(n, x) for n >= 1, where T(n, x) = Sum_{k=1..n}(-1)^(n - k)*(n/ (2*k))*binomial(k, n - k)*(2*x)^(2*k - n). - Peter Luschny, Sep 20 2022

A053763 a(n) = 2^(n^2 - n).

Original entry on oeis.org

1, 1, 4, 64, 4096, 1048576, 1073741824, 4398046511104, 72057594037927936, 4722366482869645213696, 1237940039285380274899124224, 1298074214633706907132624082305024, 5444517870735015415413993718908291383296, 91343852333181432387730302044767688728495783936
Offset: 0

Views

Author

Stephen G Penrice, Mar 29 2000

Keywords

Comments

Nilpotent n X n matrices over GF(2). Also number of simple digraphs (without self-loops) on n labeled nodes (see also A002416).
For n >= 1 a(n) is the size of the Sylow 2-subgroup of the Chevalley group A_n(4) (sequence A053291). - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 30 2001
(-1)^ceiling(n/2) * resultant of the Chebyshev polynomial of first kind of degree n and Chebyshev polynomial of first kind of degree (n+1) (cf. A039991). - Benoit Cloitre, Jan 26 2003
The number of reflexive binary relations on an n-element set. - Justin Witt (justinmwitt(AT)gmail.com), Jul 12 2005
From Rick L. Shepherd, Dec 24 2008: (Start)
Number of gift exchange scenarios where, for each person k of n people,
i) k gives gifts to g(k) of the others, where 0 <= g(k) <= n-1,
ii) k gives no more than one gift to any specific person,
iii) k gives no single gift to two or more people and
iv) there is no other person j such that j and k jointly give a single gift.
(In other words -- but less precisely -- each person k either gives no gifts or gives exactly one gift per person to 1 <= g(k) <= n-1 others.) (End)
In general, sequences of the form m^((n^2 - n)/2) enumerate the graphs with n labeled nodes with m types of edge. a(n) therefore is the number of labeled graphs with n nodes with 4 types of edge. To clarify the comment from Benoit Cloitre, dated Jan 26 2003, in this context: simple digraphs (without self-loops) have four types of edge. These types of edges are as follows: the absent edge, the directed edge from A -> B, the directed edge from B -> A and the bidirectional edge, A <-> B. - Mark Stander, Apr 11 2019

Examples

			a(2)=4 because there are four 2 x 2 nilpotent matrices over GF(2):{{0,0},{0,0}},{{0,1},{0,0}},{{0,0},{1,0}},{{1,1,},{1,1}} where 1+1=0. - _Geoffrey Critzer_, Oct 05 2012
		

References

  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 521.
  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 5, Eq. (1.1.5).

Crossrefs

Programs

Formula

Sequence given by the Hankel transform (see A001906 for definition) of A059231 = {1, 1, 5, 29, 185, 1257, 8925, 65445, 491825, ...}; example: det([1, 1, 5, 29; 1, 5, 29, 185; 5, 29, 185, 1257; 29, 185, 1257, 8925]) = 4^6 = 4096. - Philippe Deléham, Aug 20 2005
a(n) = 4^binomial(n, n-2). - Zerinvary Lajos, Jun 16 2007
a(n) = Sum_{i=0..n^2-n} binomial(n^2-n, i). - Rick L. Shepherd, Dec 24 2008
G.f. A(x) satisfies: A(x) = 1 + x * A(4*x). - Ilya Gutkovskiy, Jun 04 2020
Sum_{n>=1} 1/a(n) = A319016. - Amiram Eldar, Oct 27 2020
Sum_{n>=0} a(n)*u^n/A002884(n) = Product_{r>=1} 1/(1-u/q^r). - Geoffrey Critzer, Oct 28 2021

A001793 a(n) = n*(n+3)*2^(n-3).

Original entry on oeis.org

1, 5, 18, 56, 160, 432, 1120, 2816, 6912, 16640, 39424, 92160, 212992, 487424, 1105920, 2490368, 5570560, 12386304, 27394048, 60293120, 132120576, 288358400, 627048448, 1358954496, 2936012800, 6325010432, 13589544960, 29125246976
Offset: 1

Views

Author

Keywords

Comments

Coefficients of Chebyshev T polynomials: the subdiagonal A053120(n+3, n-1), for n > = 1. [rewritten by Wolfdieter Lang, Nov 25 2019]
Number of 132-avoiding permutations of [n+3] containing exactly two 123 patterns. - Emeric Deutsch, Jul 13 2001
Number of Dyck paths of semilength n+2 having pyramid weight n+1 (for pyramid weight see Denise and Simion). Example: a(2)=5 because the Dyck paths of semilength 4 having pyramid weight 3 are: (ud)u(ud)(ud)d, u(ud)(ud)d(ud), u(ud)(ud)(ud)d, u(ud)(uudd)d and u(uudd)(ud)d [here u=(1,1), d=(1,-1) and the maximal pyramids, of total length 3, are shown between parentheses]. - Emeric Deutsch, Mar 10 2004
a(n) is the number of dissections of a regular (n+3)-gon using n-1 noncrossing diagonals such that every piece of the dissection contains at least one non-base side of the (n+3)-gon. (One side of the (n+3)-gon is designated the base.) - David Callan, Mar 23 2004
If X_1,X_2,...,X_n are 2-blocks of a (2n+1)-set X then a(n) is the number of (n+2)-subsets of X intersecting each X_i, (i=1..n). - Milan Janjic, Nov 18 2007
The second corrector line for transforming 2^n offset 0 with a leading 1 into the Fibonacci sequence. - Al Hakanson (hawkuu(AT)gmail.com), Jun 01 2009
Sum of all nodes of all integer compositions of n, see example. - Olivier Gérard, Oct 22 2011
Number of compositions of 2n with exactly two odd summands (see example). - Mamuka Jibladze, Sep 04 2013
4*a(n) is the number of North-East paths from (0,0) to (n+2,n+2) with exactly two east steps below y = x-1 or above y = x+1. It is related to paired pattern P_1 and P_6 in Pan and Remmel's link. - Ran Pan, Feb 04 2016
From Paul Weisenhorn, Oct 18 2019: (Start)
The polynomials S(n,x)= Sum_(k>=1) b(n,k)*x^k has the recurrence relation S(n+2,x)=2*S(n+1,x))-x*S(n) with S(1,x)=1, S(2,x)=2-x and are generated by the coefficients b(n,k). b(n,k) is defined by b(n,k)=Sum_(j=1..k) binomials(k+1,j)*b(n-j,k) or by b(n,k)=((n-2+k)!*(n-1+2k)*2^n)/(4*(n-1)!*k!). b(n,1)=A001792, b(n,2)=A001793, b(n,3)=A001794, b(n,4)=A006974, b(n,5)=A006975, b(n,6)=A006976, b(n,7)=A209404.
The general formula for the sequences with k>=1: a(n)=((n-2+k)!*(n-1+2k)*2^n)/(4*(n-1)!*k!) with n >= 1. (End) [See a comment in A053120 on subdiagonal sequences. - Wolfdieter Lang, Jan 03 2020]

Examples

			a(2)=5 since 32415, 32451, 34125, 42135 and 52134 are the only 132-avoiding permutations of 12345 containing exactly two increasing subsequences of length 3.
a(4)=56: the compositions of 4 are 4, 3+1, 1+3, 2+2, 2+1+1, 1+2+1, 1+1+2, 1+1+1+1, the corresponding nodes (partial sums) are {0, 4}, {0, 3, 4}, {0, 1, 4}, {0, 2, 4}, {0, 2, 3, 4}, {0, 1, 3, 4}, {0, 1, 2, 4}, {0, 1, 2, 3, 4}, with individual sums {4, 7, 5, 6, 9, 8, 7, 10} and total 56. - _Olivier Gérard_, Oct 22 2011
The a(3)=18 compositions of 2*3=6 with two odd summands are 5+1, 1+5, 3+3, 4+1+1, 1+4+1, 1+1+4, 3+2+1, 3+1+2, 2+3+1, 2+1+3, 1+3+2, 1+2+3, 2+2+1+1, 2+1+2+1, 2+1+1+2, 1+2+2+1, 1+2+1+2, 1+1+2+2. - _Mamuka Jibladze_, Sep 04 2013
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 795.
  • 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).

Crossrefs

a(n) = A039991(n+3, 4) = A055252(n, 1).
Cf. A053120.

Programs

Formula

G.f.: x*(1-x)/(1-2*x)^3. Binomial transform of squares [1, 4, 9, ...].
a(n) = Sum_{k=0..floor((n+4)/2)} C(n+4, 2k)*C(k, 2). - Paul Barry, May 15 2003
With two leading zeros, binomial transform of quarter-squares A002620. - Paul Barry, May 27 2003
a(n) = Sum_{k=0..n+2} C(n+2, k) * floor(k^2/4). - Paul Barry, May 27 2003
a(n) = Sum_{i=0..j} binomial(i+1, 2)*binomial(j, i). - Jon Perry, Feb 26 2004
With one leading zero, binomial transform of triangular numbers A000217. - Philippe Deléham, Aug 02 2005
a(n) = Sum_{k=0..n+1} (-1)^(n-k+1)*C(k, n-k+1)*k*C(2k, k)/2. - Paul Barry, Oct 07 2005
Left-shifted sequence is binomial transform of left-shifted squares (A000290). - Franklin T. Adams-Watters, Nov 29 2006
Binomial transform of a(n) = n^2 offset 1. a(3)=18. - Al Hakanson (hawkuu(AT)gmail.com), Jun 01 2009
a(n) = (1/n) * Sum_{k=0..n} binomial(n,k)*k^3. - Gary Detlefs, Nov 26 2011
For n > 1, a(n) = Sum_{k=0..n-1} Sum_{i=0..n} (k+2) * C(n-2,i). - Wesley Ivan Hurt, Sep 20 2017
a(n) = a(-3-n)*2^(2*n+3), a(n)*(n+3) = -A058645(-3-n)*2^(2*n+3) for all n in Z. - Michael Somos, Apr 19 2019
E.g.f.: (1/2)*exp(2*x)*x*(2 + x). - Stefano Spezia, Aug 17 2019
From Amiram Eldar, Jan 05 2022: (Start)
Sum_{n>=1} 1/a(n) = 128/9 - 56*log(2)/3.
Sum_{n>=1} (-1)^(n+1)/a(n) = 24*log(3/2) - 80/9. (End)

A008310 Triangle of coefficients of Chebyshev polynomials T_n(x).

Original entry on oeis.org

1, 1, -1, 2, -3, 4, 1, -8, 8, 5, -20, 16, -1, 18, -48, 32, -7, 56, -112, 64, 1, -32, 160, -256, 128, 9, -120, 432, -576, 256, -1, 50, -400, 1120, -1280, 512, -11, 220, -1232, 2816, -2816, 1024, 1, -72, 840, -3584, 6912, -6144, 2048, 13, -364, 2912, -9984, 16640, -13312, 4096
Offset: 0

Views

Author

Keywords

Comments

The row length sequence of this irregular array is A008619(n), n >= 0. Even or odd powers appear in increasing order starting with 1 or x for even or odd row numbers n, respectively. This is the standard triangle A053120 with 0 deleted. - Wolfdieter Lang, Aug 02 2014
Let T* denote the triangle obtained by replacing each number in this triangle by its absolute value. Then T* gives the coefficients for cos(nx) as a polynomial in cos x. - Clark Kimberling, Aug 04 2024

Examples

			Rows are: (1), (1), (-1,2), (-3,4), (1,-8,8), (5,-20,16) etc., since if c = cos(x): cos(0x) = 1, cos(1x) = 1c; cos(2x) = -1+2c^2; cos(3x) = -3c+4c^3, cos(4x) = 1-8c^2+8c^4, cos(5x) = 5c-20c^3+16c^5, etc.
From _Wolfdieter Lang_, Aug 02 2014: (Start)
This irregular triangle a(n,k) begins:
  n\k   0    1     2      3      4      5      6      7 ...
  0:    1
  1:    1
  2:   -1    2
  3:   -3    4
  4:    1   -8     8
  5:    5  -20    16
  6:   -1   18   -48     32
  7:   -7   56  -112     64
  8:    1  -32   160   -256    128
  9:    9 -120   432   -576    256
 10:   -1   50  -400   1120  -1280    512
 11:  -11  220 -1232   2816  -2816   1024
 12:    1  -72   840  -3584   6912  -6144   2048
 13:   13 -364  2912  -9984  16640 -13312   4096
 14:   -1   98 -1568   9408 -26880  39424 -28672   8192
 15:  -15  560 -6048  28800 -70400  92160 -61440  16384
  ...
T(4,x) = 1 - 8*x^2 + 8*x^4, T(5,x) = 5*x - 20*x^3 +16*x^5.
(End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 795.
  • E. A. Guilleman, Synthesis of Passive Networks, Wiley, 1957, p. 593.
  • Yaroslav Zolotaryuk, J. Chris Eilbeck, "Analytical approach to the Davydov-Scott theory with on-site potential", Physical Review B 63, p543402, Jan. 2001. The authors write, "Since the algebra of these is 'hyperbolic', contrary to the usual Chebyshev polynomials defined on the interval 0 <= x <= 1, we call the set of functions (21) the hyperbolic Chebyshev polynomials." (This refers to the triangle T* described in Comments.)

Crossrefs

A039991 is a row reversed version, but has zeros which enable the triangle to be seen. Columns/diagonals are A011782, A001792, A001793, A001794, A006974, A006975, A006976 etc.
Reflection of A028297. Cf. A008312, A053112.
Row sums are one. Polynomial evaluations include A001075 (x=2), A001541 (x=3), A001091, A001079, A023038, A011943, A001081, A023039, A001085, A077422, A077424, A097308, A097310, A068203.
Cf. A053120.

Programs

  • Maple
    A008310 := proc(n,m) local x ; coeftayl(simplify(ChebyshevT(n,x),'ChebyshevT'),x=0,m) ; end: i := 0 : for n from 0 to 100 do for m from n mod 2 to n by 2 do printf("%d %d ",i,A008310(n,m)) ; i := i+1 ; od ; od ; # R. J. Mathar, Apr 20 2007
    # second Maple program:
    b:= proc(n) b(n):= `if`(n<2, 1, expand(2*b(n-1)-x*b(n-2))) end:
    T:= n-> (p-> (d-> seq(coeff(p, x, d-i), i=0..d))(degree(p)))(b(n)):
    seq(T(n), n=0..15);  # Alois P. Heinz, Sep 04 2019
  • Mathematica
    Flatten[{1, Table[CoefficientList[ChebyshevT[n, x], x], {n, 1, 13}]}]//DeleteCases[#, 0, Infinity]& (* or *) Flatten[{1, Table[Table[((-1)^k*2^(n-2 k-1)*n*Binomial[n-k, k])/(n-k), {k, Floor[n/2], 0, -1}], {n, 1, 13}]}] (* Eugeniy Sokol, Sep 04 2019 *)

Formula

a(n,m) = 2^(m-1) * n * (-1)^((n-m)/2) * ((n+m)/2-1)! / (((n-m)/2)! * m!) if n>0. - R. J. Mathar, Apr 20 2007
From Paul Weisenhorn, Oct 02 2019: (Start)
T_n(x) = 2*x*T_(n-1)(x) - T_(n-2)(x), T_0(x) = 1, T_1(x) = x.
T_n(x) = ((x+sqrt(x^2-1))^n + (x-sqrt(x^2-1))^n)/2. (End)
From Peter Bala, Aug 15 2022: (Start)
T(n,x) = [z^n] ( z*x + sqrt(1 + z^2*(x^2 - 1)) )^n.
Sum_{k = 0..2*n} binomial(2*n,k)*T(k,x) = (2^n)*(1 + x)^n*T(n,x).
exp( Sum_{n >= 1} T(n,x)*t^n/n ) = Sum_{n >= 0} P(n,x)*t^n, where P(n,x) denotes the n-th Legendre polynomial. (End)

Extensions

Additional comments and more terms from Henry Bottomley, Dec 13 2000
Edited: Corrected Cf. A039991 statement. Cf. A053120 added. - Wolfdieter Lang, Aug 06 2014

A028297 Coefficients of Chebyshev polynomials of the first kind: triangle of coefficients in expansion of cos(n*x) in descending powers of cos(x).

Original entry on oeis.org

1, 1, 2, -1, 4, -3, 8, -8, 1, 16, -20, 5, 32, -48, 18, -1, 64, -112, 56, -7, 128, -256, 160, -32, 1, 256, -576, 432, -120, 9, 512, -1280, 1120, -400, 50, -1, 1024, -2816, 2816, -1232, 220, -11, 2048, -6144, 6912, -3584, 840, -72, 1, 4096, -13312, 16640, -9984
Offset: 0

Views

Author

Keywords

Comments

Rows are of lengths 1, 1, 2, 2, 3, 3, ... (A008619).
This triangle is generated from A118800 by shifting down columns to allow for (1, 1, 2, 2, 3, 3, ...) terms in each row. - Gary W. Adamson, Dec 16 2007
Unsigned triangle = A034839 * A007318. - Gary W. Adamson, Nov 28 2008
Triangle, with zeros omitted, given by (1, 1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, -1, 1, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 16 2011
From Wolfdieter Lang, Aug 02 2014: (Start)
This irregular triangle is the row reversed version of the Chebyshev T-triangle A053120 given by A039991 with vanishing odd-indexed columns removed.
If zeros are appended in each row n >= 1, in order to obtain a regular triangle (see the Philippe Deléham comment, g.f. and example) this becomes the Riordan triangle (1-x)/(1-2*x), -x^2/(1-2*x). See also the unsigned version A201701 of this regular triangle.
(End)
Apparently, unsigned diagonals of this array are rows of A200139. - Tom Copeland, Oct 11 2014
It appears that the coefficients are generated by the following: Let SM_k = Sum( d_(t_1, t_2)* eM_1^t_1 * eM_2^t_2) summed over all length 2 integer partitions of k, i.e., 1*t_1 + 2*t_2 = k, where SM_k are the averaged k-th power sum symmetric polynomials in 2 data (i.e., SM_k = S_k/2 where S_k are the k-th power sum symmetric polynomials, and where eM_k are the averaged k-th elementary symmetric polynomials, eM_k = e_k/binomial(2,k) with e_k being the k-th elementary symmetric polynomials. The data d_(t_1, t_2) form an irregular triangle, with one row for each k value starting with k=1. Thus this procedure and associated OEIS sequences A287768, A288199, A288207, A288211, A288245, A288188 are generalizations of Chebyshev polynomials of the first kind. - Gregory Gerard Wojnar, Jul 01 2017

Examples

			Letting c = cos x, we have: cos 0x = 1, cos 1x = 1c; cos 2x = 2c^2-1; cos 3x = 4c^3-3c, cos 4x = 8c^4-8c^2+1, etc.
T4 = 8x^4 - 8x^2 + 1 = 8, -8, +1 = 2^(3) - (4)(2) + [2^(-1)](4)/2.
From _Wolfdieter Lang_, Aug 02 2014: (Start)
The irregular triangle T(n,k) begins:
n\k     1      2     3      4     5     6   7   8 ....
0:      1
1:      1
2:      2     -1
3:      4     -3
4:      8     -8     1
5:     16    -20     5
6:     32    -48    18     -1
7:     64   -112    56     -7
8:    128   -256   160    -32     1
9:    256   -576   432   -120     9
10:   512  -1280  1120   -400    50    -1
11:  1024  -2816  2816  -1232   220   -11
12:  2048  -6144  6912  -3584   840   -72   1
13:  4096 -13312 16640  -9984  2912  -364  13
14:  8192 -28672 39424 -26880  9408 -1568  98  -1
15: 16384 -61440 92160 -70400 28800 -6048 560 -15
...
T(4,x) = 8*x^4 -8*x^2 + 1*x^0, T(5,x) = 16*x^5 - 20*x^3 + 5*x^1, with Chebyshev's T-polynomials (A053120). (End)
From _Philippe Deléham_, Dec 16 2011: (Start)
The triangle (1,1,0,0,0,0,...) DELTA (0,-1,1,0,0,0,0,...) includes zeros and begins:
   1;
   1,   0;
   2,  -1,  0;
   4,  -3,  0,  0;
   8,  -8,  1,  0, 0;
  16, -20,  5,  0, 0, 0;
  32, -48, 18, -1, 0, 0, 0; (End)
		

References

  • I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series and Products, 5th ed., Section 1.335, p. 35.
  • S. Selby, editor, CRC Basic Mathematical Tables, CRC Press, 1970, p. 106. [From Rick L. Shepherd, Jul 06 2010]

Crossrefs

Cf. A028298.
Reflection of A008310, the main entry. With zeros: A039991.
Cf. A053120 (row reversed table including zeros).
Cf. A001333 (row sums 1), A001333 (alternating row sums). - Wolfdieter Lang, Aug 02 2014

Programs

  • Maple
    b:= proc(n) b(n):= `if`(n<2, 1, expand(2*b(n-1)-x*b(n-2))) end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(n)):
    seq(T(n), n=0..15);  # Alois P. Heinz, Sep 04 2019
  • Mathematica
    t[n_] := (Cos[n x] // TrigExpand) /. Sin[x]^m_ /; EvenQ[m] -> (1 - Cos[x]^2)^(m/2) // Expand; Flatten[Table[ r = Reverse @ CoefficientList[t[n], Cos[x]]; If[OddQ[Length[r]], AppendTo[r,0]]; Partition[r,2][[All, 1]],{n, 0, 13}] ][[1 ;; 53]] (* Jean-François Alcover, May 06 2011 *)
    Tpoly[n_] := HypergeometricPFQ[{(1 - n)/2, -n/2}, {1/2}, 1 - x];
    Table[CoefficientList[Tpoly[n], x], {n, 0, 12}] // Flatten (* Peter Luschny, Feb 03 2021 *)

Formula

cos(n*x) = 2 * cos((n-1)*x) * cos(x) - cos((n-2)*x) (from CRC's Multiple-angle relations). - Rick L. Shepherd, Jul 06 2010
G.f.: (1-x) / (1-2x+y*x^2). - Philippe Deléham, Dec 16 2011
Sum_{k=0..n} T(n,k)*x^k = A011782(n), A000012(n), A146559(n), A087455(n), A138230(n), A006495(n), A138229(n) for x = 0, 1, 2, 3, 4, 5, 6, respectively. - Philippe Deléham, Dec 16 2011
T(n,k) = [x^k] hypergeom([1/2 - n/2, -n/2], [1/2], 1 - x). - Peter Luschny, Feb 03 2021
T(n,k) = (-1)^k * 2^(n-1-2*k) * A034807(n,k). - Hoang Xuan Thanh, Jun 21 2025

Extensions

More terms from David W. Wilson
Row length sequence and link to Abramowitz-Stegun added by Wolfdieter Lang, Aug 02 2014

A001794 Negated coefficients of Chebyshev T polynomials: [x^n](-T(n+6, x)), n >= 0.

Original entry on oeis.org

1, 7, 32, 120, 400, 1232, 3584, 9984, 26880, 70400, 180224, 452608, 1118208, 2723840, 6553600, 15597568, 36765696, 85917696, 199229440, 458752000, 1049624576, 2387607552, 5402263552, 12163481600, 27262976000, 60850962432
Offset: 0

Views

Author

Keywords

Comments

A negated subdiagonal of A053120.
If X_1,X_2,...,X_n are 2-blocks of a (2n+1)-set X then a(n-2) is the number of (n+3)-subsets of X intersecting each X_i, (i=1,2,...,n). - Milan Janjic, Nov 18 2007
The third corrector line for transforming 2^n offset 0 with a leading 1 into the Fibonacci sequence. - Al Hakanson (hawkuu(AT)gmail.com), Jun 01 2009
a(n-2) is the number of strings of length n defined on {0, 1, 2, 3} that have exactly two 2's and no 3's or exactly three 3's and no 2's. For example, for n=4, a(2)=32 since the strings are the 6 permutations of 2200, the 12 permutations of 2201, the 6 permutations of 2211, the 4 permutations of 3330, and the 4 permutations of 3331. - Enrique Navarrete, Jun 03 2025

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 795.
  • 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).

Crossrefs

Cf. A039991 (negative of column 6), A028297, A008310, A053120.
With alternating signs, the o.g.f. (with offset 1) is the inverse of the o.g.f. of A065097.
Cf. A001789 (partial sums), A081279 (binomial transform), A005900 (0 followed by inverse binomial transform).

Programs

  • GAP
    List([0..25],n->2^(n-2)*(n+1)*(n+2)*(n+6)/3); # Muniru A Asiru, Mar 20 2018
  • Magma
    [2^(n-1)/3*Binomial(n+2,2)*(n+6) : n in [0..25]]; // Brad Clardy, Mar 08 2012
    
  • Maple
    [seq(coeftayl((1-x)/(1-2*x)^4, x = 0, k), k=0..25)]; # Muniru A Asiru, Mar 20 2018
  • Mathematica
    a[n_] := 2^(n-2)*(n+1)*(n+2)*(n+6)/3; a /@ Range[0, 20] (* Giovanni Resta, Mar 25 2017 *)
    LinearRecurrence[{8,-24,32,-16},{1,7,32,120},30] (* Harvey P. Dale, Oct 08 2024 *)
  • PARI
    a(n) = sum(i=0, n+1, sum(k=0, i, k^2*binomial(n+1, i))); \\ Michel Marcus, Mar 25 2017
    
  • PARI
    a(n) = - polcoeff(polchebyshev(n+6), n); \\ Michel Marcus, Mar 20 2018
    

Formula

a(n) = 2^(n-2)*(n+1)*(n+2)*(n+6)/3. [See a comment in A053120 on subdiagonals. - Wolfdieter Lang, Jan 03 2020]
G.f.: (1-x)/(1-2*x)^4. - Simon Plouffe in his 1992 dissertation
a(n) = Sum_{k=0..floor((n+6)/2)} C(n+6, 2*k)*C(k, 3). - Paul Barry, May 15 2003
With a leading zero, the binomial transform of A000330. - Paul Barry, Jul 19 2003
a(n) = Sum_{i=0..n+1} (Sum{k=0..i} (k^2*binomial(n+1, i))). - Jon Perry, Feb 26 2004 [corrected by Michel Marcus, Mar 25 2017]
Binomial transform of a(n) = (2*n^3 + 6*n^2 + 7*n + 3)/3 offset 0. a(3)=120. - Al Hakanson (hawkuu(AT)gmail.com), Jun 01 2009
a(n) = (2^(n-1)/3)*binomial(n+2,2)*(n+6). - Brad Clardy, Mar 08 2012
E.g.f.: (1/3)*exp(2*x)*(3 + 15*x + 12*x^2 + 2*x^3). - Stefano Spezia, Jan 03 2020
From Amiram Eldar, Jan 05 2022: (Start)
Sum_{n>=0} 1/a(n) = 156*log(2)/5 - 511/25.
Sum_{n>=0} (-1)^n/a(n) = 241/25 - 108*log(3/2)/5. (End)
E.g.f.: exp(2*x)*(x^2/2 + x^3/6) (with two leading zeros). - Enrique Navarrete, Jun 03 2025

Extensions

Name clarified by Wolfdieter Lang, Nov 26 2019

A006974 Coefficients of Chebyshev T polynomials: a(n) = A053120(n+8, n), n >= 0.

Original entry on oeis.org

1, 9, 50, 220, 840, 2912, 9408, 28800, 84480, 239360, 658944, 1770496, 4659200, 12042240, 30638080, 76873728, 190513152, 466944000, 1133117440, 2724986880, 6499598336, 15386804224, 36175872000, 84515225600, 196293427200, 453437816832
Offset: 0

Views

Author

Keywords

Comments

If X_1,X_2,...,X_n are 2-blocks of a (2n+1)-set X then, for n>=3, a(n-3) is the number of (n+4)-subsets of X intersecting each X_i, (i=1,2,...,n). - Milan Janjic, Nov 18 2007
The fourth corrector line for transforming 2^n offset 0 with a leading 1 into the Fibonacci sequence. [Al Hakanson (hawkuu(AT)gmail.com), Jun 01 2009]

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 795.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A039991 (see column 8), A003472 (partial sums), A053120.

Programs

  • Magma
    [2^(n-1)/4*Binomial(n+3,3)*(n+8) : n in [0..25]]; // Brad Clardy, Mar 08 2012
  • Maple
    a := n->n*(n+1)*(n+2)*(n+7)*2^(n-5)/3;

Formula

G.f.: (1-x)/(1-2*x)^5.
a(n) = Sum_{k=0..floor((n+8)/2)} C(n+8, 2k)*C(k, 4). - Paul Barry, May 15 2003
Binomial transform of a(n)=(24*n^4-134*n^3+261*n^2-130*n+3)/3 offset 0. a(3)=220. [Al Hakanson (hawkuu(AT)gmail.com), Jun 01 2009]
a(n) = 2^(n-3)*binomial(n+3, 3)*(n+8). - Brad Clardy, Mar 08 2012 [See a comment in A053120 on subdiagonals. - Wolfdieter Lang, Jan 03 2020]
E.g.f.: (1/3)*exp(2*x)*(3 + 21*x + 27*x^2 + 10*x^3 + x^4). - Stefano Spezia, Aug 17 2019

Extensions

Name clarified by Wolfdieter Lang, Nov 26 2019

A006975 Negated coefficients of Chebyshev T polynomials: a(n) = -A053120(n+10, n), n >= 0.

Original entry on oeis.org

1, 11, 72, 364, 1568, 6048, 21504, 71808, 228096, 695552, 2050048, 5870592, 16400384, 44843008, 120324096, 317521920, 825556992, 2118057984, 5369233408, 13463453696, 33426505728, 82239815680, 200655503360, 485826232320
Offset: 0

Views

Author

Keywords

Comments

Binomial transform of A069038. - Paul Barry, Feb 19 2003
If X_1, X_2, ..., X_n are 2-blocks of a (2n+1)-set X then, for n>=4, a(n-4) is the number of (n+5)-subsets of X intersecting each X_i, (i=1,2,...,n). - Milan Janjic, Nov 18 2007
The 5th corrector line for transforming 2^n offset 0 with a leading 1 into the Fibonacci sequence. - Al Hakanson (hawkuu(AT)gmail.com), Jun 01 2009

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 795.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

First differences of A054849.

Programs

  • Magma
    [2^(n-1)/5*Binomial(n+4,4)*(n+10): n in [0..25]]; // Brad Clardy, Mar 10 2012
  • Mathematica
    Table[2^(n-1)/5*Binomial[n+4, 4]*(n+10), {n, 0, 30}] (* Paolo Xausa, Jun 26 2024 *)

Formula

G.f.: (1-x)/(1-2*x)^6. a(n) = 2^(n-1)*binomial(n+4, 4)*(n+10)/5, for n >= 0. [a(n) from Mar 06 2000 rewritten. See the Brad Clardy formula below, and a comment in A053120 on subdiagonals. - Wolfdieter Lang, Jan 03 2020]
a(n) = 2^(n-4)*(n+1)(n+2)(n+3)(n+4)(n+10)/15. - Paul Barry, Feb 19 2003
a(n) = Sum_{k=0..floor((n+10)/2)} C(n+10, 2k)*C(k, 5). - Paul Barry, May 15 2003
a(n) = -A039991(n+10, 10). - N. J. A. Sloane, May 16 2003
a(n) = binomial transform of b(n)= (2*n^5 + 10*n^4 + 30*n^3 + 50*n^2 + 43*n + 15) / 15 offset 0. a(3) = 364. - Al Hakanson (hawkuu(AT)gmail.com), Jun 01 2009
a(n) = 2^(n-1)/5*binomial(n+4,4)*(n+10). - Brad Clardy, Mar 10 2012
E.g.f.: (1/15)*exp(2*x)*(15+135*x+240*x^2+140*x^3+30*x^4+2*x^5). - Stefano Spezia, Jan 03 2020

Extensions

Name clarified by Wolfdieter Lang, Nov 26 2019
Showing 1-10 of 17 results. Next