A007318 Pascal's triangle read by rows: C(n,k) = binomial(n,k) = n!/(k!*(n-k)!), 0 <= k <= n.
1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 6, 4, 1, 1, 5, 10, 10, 5, 1, 1, 6, 15, 20, 15, 6, 1, 1, 7, 21, 35, 35, 21, 7, 1, 1, 8, 28, 56, 70, 56, 28, 8, 1, 1, 9, 36, 84, 126, 126, 84, 36, 9, 1, 1, 10, 45, 120, 210, 252, 210, 120, 45, 10, 1, 1, 11, 55, 165, 330, 462, 462, 330, 165, 55, 11, 1
Offset: 0
Examples
Triangle T(n,k) begins: n\k 0 1 2 3 4 5 6 7 8 9 10 11 ... 0 1 1 1 1 2 1 2 1 3 1 3 3 1 4 1 4 6 4 1 5 1 5 10 10 5 1 6 1 6 15 20 15 6 1 7 1 7 21 35 35 21 7 1 8 1 8 28 56 70 56 28 8 1 9 1 9 36 84 126 126 84 36 9 1 10 1 10 45 120 210 252 210 120 45 10 1 11 1 11 55 165 330 462 462 330 165 55 11 1 ... There are C(4,2)=6 ways to distribute 5 balls BBBBB, among 3 different urns, < > ( ) [ ], so that each urn gets at least one ball, namely, <BBB>(B)[B], <B>(BBB)[B], <B>(B)[BBB], <BB>(BB)[B], <BB>(B)[BB], and <B>(BB)[BB]. There are C(4,2)=6 increasing functions from {1,2} to {1,2,3,4}, namely, {(1,1),(2,2)},{(1,1),(2,3)}, {(1,1),(2,4)}, {(1,2),(2,3)}, {(1,2),(2,4)}, and {(1,3),(2,4)}. - _Dennis P. Walsh_, Apr 07 2011 There are C(4,2)=6 subsets of {1,2,3,4,5} with median element 3, namely, {3}, {1,3,4}, {1,3,5}, {2,3,4}, {2,3,5}, and {1,2,3,4,5}. - _Dennis P. Walsh_, Dec 15 2011 The successive k-iterations of {A(0)} = E are E;E;E;...; the corresponding number of elements are 1,1,1,... The successive k-iterations of {A(1)} = {a} are (omitting brackets) a;a,E; a,E,E;...; the corresponding number of elements are 1,2,3,... The successive k-iterations of {A(2)} = {a,a} are aa; aa,a,E; aa, a, E and a,E and E;...; the corresponding number of elements are 1,3,6,... - _Gregory L. Simay_, Aug 06 2018 Boas-Buck type recurrence for column k = 4: T(8, 4) = (5/4)*(1 + 5 + 15 + 35) = 70. See the Boas-Buck comment above. - _Wolfdieter Lang_, Nov 12 2018
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. 828.
- Amulya Kumar Bag, Binomial theorem in ancient India, Indian Journal of History of Science, vol. 1 (1966), pp. 68-74.
- Arthur T. Benjamin and Jennifer Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, p. 63ff.
- Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8.
- Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 306.
- John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 68-74.
- Paul Curtz, Intégration numérique des systèmes différentiels à conditions initiales, Centre de Calcul Scientifique de l'Armement, Arcueil, 1969.
- A. W. F. Edwards, Pascal's Arithmetical Triangle, 2002.
- William Feller, An Introduction to Probability Theory and Its Application, Vol. 1, 2nd ed. New York: Wiley, p. 36, 1968.
- Ronald L. Graham, Donald E. Knuth, and Oren Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 2nd. ed., 1994, p. 155.
- Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §4.4 Powers and Roots, pp. 140-141.
- David Hök, Parvisa mönster i permutationer [Swedish], 2007.
- Donald E. Knuth, The Art of Computer Programming, Vol. 1, 2nd ed., p. 52.
- Sergei K. Lando, Lecture on Generating Functions, Amer. Math. Soc., Providence, R.I., 2003, pp. 60-61.
- Blaise Pascal, Traité du triangle arithmétique, avec quelques autres petits traitez sur la mesme matière, Desprez, Paris, 1665.
- Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 71.
- Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 271-275.
- A. P. Prudnikov, Yu. A. Brychkov, and O. I. Marichev, "Integrals and Series", Volume 1: "Elementary Functions", Chapter 4: "Finite Sums", New York, Gordon and Breach Science Publishers, 1986-1992.
- John Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 6.
- John Riordan, Combinatorial Identities, Wiley, 1968, p. 2.
- Robert Sedgewick and Philippe Flajolet, An Introduction to the Analysis of Algorithms, Addison-Wesley, Reading, MA, 1996, p. 143.
- N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
- Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 6, pages 43-52.
- James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 13, 30-33.
- David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, 1987, pp. 115-118.
- Douglas B. West, Combinatorial Mathematics, Cambridge, 2021, p. 25.
Links
- N. J. A. Sloane, First 141 rows of Pascal's triangle, formatted as a simple linear sequence: (n, a(n)), n=0..10152.
- M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
- Tewodros Amdeberhan, Moa Apagodu, and Doron Zeilberger, Wilf's "Snake Oil" Method Proves an Identity in The Motzkin Triangle, arXiv:1507.07660 [math.CO], 2015.
- Said Amrouche and Hacène Belbachir, Asymmetric extension of Pascal-Dellanoy triangles, arXiv:2001.11665 [math.CO], 2020.
- Shaun V. Ault and Charles Kicey, Counting paths in corridors using circular Pascal arrays, Discrete Mathematics, Vol. 332, No. 6 (2014), pp. 45-54.
- Mohammad K. Azarian, Fibonacci Identities as Binomial Sums, International Journal of Contemporary Mathematical Sciences, Vol. 7, No. 38 (2012), pp. 1871-1876.
- Mohammad K. Azarian, Fibonacci Identities as Binomial Sums II, International Journal of Contemporary Mathematical Sciences, Vol. 7, No. 42 (2012), pp. 2053-2059.
- Amulya Kumar Bag, Binomial theorem in ancient India, Indian Journal of History of Science, Vol. 1 (1966), pp. 68-74.
- Armen G. Bagdasaryan and Ovidiu Bagdasar, On some results concerning generalized arithmetic triangles, Electronic Notes in Discrete Mathematics, Vol. 67 (2018), pp. 71-77.
- Peter Bala, A combinatorial interpretation for the binomial coefficients, 2013.
- Cyril Banderier and Donatella Merlini, Lattice paths with an infinite set of jumps, Proceedings of the 14th International Conference on Formal Power Series and Algebraic Combinatorics, Melbourne, Australia. 2002.
- J. Fernando Barbero G., Jesús Salas, and Eduardo J. S. Villaseñor, Bivariate Generating Functions for a Class of Linear Recurrences. I. General Structure, arXiv:1307.2010 [math.CO], 2013.
- Paul Barry, On Integer-Sequence-Based Constructions of Generalized Pascal Triangles, Journal of Integer Sequences, Vol. 9 (2006), Article 06.2.4.
- Paul Barry, Symmetric Third-Order Recurring Sequences, Chebyshev Polynomials, and Riordan Arrays , JIS, Vol. 12 (2009) Article 09.8.6.
- Paul Barry, Eulerian polynomials as moments, via exponential Riordan arrays, arXiv:1105.3043 [math.CO], 2011.
- Paul Barry, Combinatorial polynomials as moments, Hankel transforms and exponential Riordan arrays, arXiv:1105.3044 [math.CO], 2011.
- Paul Barry, On the Central Coefficients of Bell Matrices, J. Int. Seq., Vol. 14 (2011) Article 11.4.3, example 2.
- Paul Barry, Riordan-Bernstein Polynomials, Hankel Transforms and Somos Sequences, Journal of Integer Sequences, Vol. 15 (2012), Article 12.8.2.
- Paul Barry, On the Central Coefficients of Riordan Matrices, Journal of Integer Sequences, Vol. 16 (2013), Article 13.5.1.
- Paul Barry, A Note on a Family of Generalized Pascal Matrices Defined by Riordan Arrays, Journal of Integer Sequences, Vol. 16 (2013), Article 13.5.4.
- Paul Barry, On the Inverses of a Family of Pascal-Like Matrices Defined by Riordan Arrays, Journal of Integer Sequences, Vol. 16 (2013), Article 13.5.6.
- Paul Barry, On the Connection Coefficients of the Chebyshev-Boubaker polynomials, The Scientific World Journal, Vol. 2013 (2013), Article ID 657806, 10 pages.
- Paul Barry, General Eulerian Polynomials as Moments Using Exponential Riordan Arrays, Journal of Integer Sequences, Vol. 16 (2013), Article 13.9.6.
- Paul Barry, Riordan arrays, generalized Narayana triangles, and series reversion, Linear Algebra and its Applications, Vol. 491 (2016), pp. 343-385.
- Paul Barry, The Gamma-Vectors of Pascal-like Triangles Defined by Riordan Arrays, arXiv:1804.05027 [math.CO], 2018.
- Paul Barry, On the f-Matrices of Pascal-like Triangles Defined by Riordan Arrays, arXiv:1805.02274 [math.CO], 2018.
- Paul Barry, The Central Coefficients of a Family of Pascal-like Triangles and Colored Lattice Paths, J. Int. Seq., Vol. 22 (2019), Article 19.1.3.
- Paul Barry, On the halves of a Riordan array and their antecedents, arXiv:1906.06373 [math.CO], 2019.
- Paul Barry, On the r-shifted central triangles of a Riordan array, arXiv:1906.01328 [math.CO], 2019.
- Paul Barry, Generalized Catalan Numbers Associated with a Family of Pascal-like Triangles, J. Int. Seq., Vol. 22 (2019), Article 19.5.8.
- Paul Barry, A Note on Riordan Arrays with Catalan Halves, arXiv:1912.01124 [math.CO], 2019.
- Paul Barry, Chebyshev moments and Riordan involutions, arXiv:1912.11845 [math.CO], 2019.
- Paul Barry, Characterizations of the Borel triangle and Borel polynomials, arXiv:2001.08799 [math.CO], 2020.
- Paul Barry, On a Central Transform of Integer Sequences, arXiv:2004.04577 [math.CO], 2020.
- Paul Barry, Extensions of Riordan Arrays and Their Applications, Mathematics (2025) Vol. 13, No. 2, 242. See p. 13.
- Paul Barry, Notes on Riordan arrays and lattice paths, arXiv:2504.09719 [math.CO], 2025. See p. 2.
- Paul Barry and Aoife Hennessy, Four-term Recurrences, Orthogonal Polynomials and Riordan Arrays, Journal of Integer Sequences, Vol. 15 (2012), Article 12.4.2.
- Jonathan W. Bober, Factorial ratios, hypergeometric series, and a family of step functions, arXiv:0709.1977v1 [math.NT], J. London Math. Soc. (2), Vol. 79 (2009), pp. 422-444.
- Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids, English translation published by Fibonacci Association, Santa Clara Univ., Santa Clara, CA, 1993; see p. 4.
- Michael Bukata, Ryan Kulwicki, Nicholas Lewandowski, Lara Pudwell, Jacob Roth and Teresa Wheeland, Distributions of Statistics over Pattern-Avoiding Permutations, arXiv preprint arXiv:1812.07112 [math.CO], 2018.
- Douglas Butler, Pascal's Triangle.
- Isabel Cação, Helmuth R. Malonek, Maria Irene Falcão, and Graça Tomaz, Intrinsic Properties of a Non-Symmetric Number Triangle, J. Int. Seq., Vol. 26 (2023), Article 23.4.8.
- Naiomi T. Cameron and Asamoah Nkwanta, On Some (Pseudo) Involutions in the Riordan Group, Journal of Integer Sequences, Vol. 8 (2005), Article 05.3.7.
- Dario T. de Castro, p-adic Order of Positive Integers via Binomial Coefficients, INTEGERS, Electronic J. of Combinatorial Number Theory, Vol. 22, Paper A61, 2022.
- Ji Young Choi, Digit Sums Generalizing Binomial Coefficients, J. Int. Seq., Vol. 22 (2019), Article 19.8.3.
- Cristian Cobeli and Alexandru Zaharescu, Promenade around Pascal Triangle - Number Motives, Bull. Math. Soc. Sci. Math. Roumanie, Tome 56(104) No. 1 (2013), pp. 73-98.
- CombOS - Combinatorial Object Server, Generate combinations.
- J. H. Conway and N. J. A. Sloane, Low-dimensional lattices. VII Coordination sequences, Proc. R. Soc. Lond. A, Vo. 453, No. 1966 (1997), pp. 2369-2389.
- Tom Copeland, Infinigens, the Pascal Triangle, and the Witt and Virasoro Algebras.
- Persi Diaconis, The distribution of leading digits and uniform distribution mod 1, Ann. Probability, Vol. 5 (1977), pp. 72-81.
- Karl Dilcher and Kenneth B. Stolarsky, A Pascal-Type Triangle Characterizing Twin Primes, The American Mathematical Monthly, Vol. 112, No. 8 (Oct 2005), pp. 673-681.
- Tomislav Došlic and Darko Veljan, Logarithmic behavior of some combinatorial sequences, Discrete Math., Vol. 308, No. 11 (2008), pp. 2182-2212. MR2404544 (2009j:05019).
- Steffen Eger, Some Elementary Congruences for the Number of Weighted Integer Compositions, J. Int. Seq., Vol. 18 (2015), Article 15.4.1.
- Leonhard Euler, On the expansion of the power of any polynomial (1+x+x^2+x^3+x^4+etc.)^n, arXiv:math/0505425 [math.HO], 2005. See also The Euler Archive, item E709.
- Jackson Evoniuk, Steven Klee, and Van Magnan, Enumerating Minimal Length Lattice Paths, J. Int. Seq., Vol. 21 (2018), Article 18.3.6.
- A. Farina, S. Giompapa, A. Graziano, A. Liburdi, M. Ravanelli, and F. Zirilli, Tartaglia-Pascal's triangle: a historical perspective with applications, Signal, Image and Video Processing, Vol. 7, No. 1 (January 2013), pp. 173-188.
- Steven Finch, Pascal Sebah, and Zai-Qiao Bai, Odd Entries in Pascal's Trinomial Triangle, arXiv:0802.2654 [math.NT], 2008.
- David Fowler, The binomial coefficient function, Amer. Math. Monthly, Vol. 103, No. 1 (1996), pp. 1-17.
- Shishuo Fu and Yaling Wang, Bijective recurrences concerning two Schröder triangles, arXiv:1908.03912 [math.CO], 2019.
- Tom Halverson and Theodore N. Jacobson, Set-partition tableaux and representations of diagram algebras, arXiv:1808.08118 [math.RT], 2018.
- T. Han and S. Kitaev, Joint distributions of statistics over permutations avoiding two patterns of length 3, arXiv:2311.02974 [math.CO], 2023
- Brady Haran and Casandra Monroe, Pascal's Triangle, Numberphile video (2017).
- Tian-Xiao He and Renzo Sprugnoli, Sequence characterization of Riordan arrays, Discrete Math., Vol. 309, No. 12 (2009), pp. 3962-3974.
- Nick Hobson, Python program for A007318.
- V. E. Hoggatt, Jr. and Marjorie Bicknell, Catalan and related sequences arising from inverses of Pascal's triangle matrices, Fib. Quart., Vol. 14, No. 5 (1976), pp. 395-405.
- Matthew Hubbard and Tom Roby, Pascal's Triangle From Top to Bottom. [archived page]
- Charles Jordan, Calculus of Finite Differences (p. 65).
- Subhash Kak, The golden mean and the physics of aesthetics, in: B. Yadav and M. Mohan (eds.), Ancient Indian Leaps into Mathematics, Birkhäuser, Boston, MA, 2009, pp. 111-119; arXiv preprint, arXiv:physics/0411195 [physics.hist-ph], 2004.
- Petro Kolosov, Polynomial identities involving Pascal's triangle rows, 2022.
- Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seq., Vol. 3 (2000), Article 00.2.4.
- Eitan Y. Levine, GCD formula proof.
- Meng Li and Ron Goldman, Limits of sums for binomial and Eulerian numbers and their associated distributions, Discrete mathematics, Vol. 343, No. 7 (2020), 111870.
- P. A. MacMahon, Memoir on the Theory of the Compositions of Numbers, Phil. Trans. Royal Soc. London A, Vol. 184 (1893), pp. 835-901.
- Mathforum, Pascal's Triangle
- Carl McTague, On the Greatest Common Divisor of C(q*n,n), C(q*n,2*n), ...C(q*n,q*n-q), arXiv:1510.06696 [math.CO], 2015.
- D. Merlini, R. Sprugnoli, and M. C. Verri, An algebra for proper generating trees, in: D. Gardy and A. Mokkadem (eds.), Mathematics and Computer Science, Trends in Mathematics, Birkhäuser, Basel, 2000, pp. 127-139; alternative link.
- Donatella Merlini, Francesca Uncini, and M. Cecilia Verri, A unified approach to the study of general and palindromic compositions, Integers, Vol. 4 (2004), A23, 26 pp.
- Ângela Mestre and José Agapito, A Family of Riordan Group Automorphisms, J. Int. Seq., Vol. 22 (2019), Article 19.8.5.
- Pierre Remond de Montmort, Essay d'analyse sur les jeux de hazard, Paris: Chez Jacque Quillau, 1708, p. 80.
- Yossi Moshe, The density of 0's in recurrence double sequences, J. Number Theory, Vol. 103 (2003), pp. 109-121.
- Lili Mu and Sai-nan Zheng, On the Total Positivity of Delannoy-Like Triangles, Journal of Integer Sequences, Vol. 20 (2017), Article 17.1.6.
- Abdelkader Necer, Séries formelles et produit de Hadamard, Journal de théorie des nombres de Bordeaux, Vol. 9, No. 2 (1997), pp. 319-335.
- Asamoah Nkwanta and Earl R. Barnes, Two Catalan-type Riordan Arrays and their Connections to the Chebyshev Polynomials of the First Kind, Journal of Integer Sequences, Vol. 15 (2012), Article 12.3.3.
- Asamoah Nkwanta and Akalu Tefera, Curious Relations and Identities Involving the Catalan Generating Function and Numbers, Journal of Integer Sequences, Vol. 16 (2013), Article 13.9.5.
- Mustafa A. A. Obaid, S. Khalid Nauman, Wafaa M. Fakieh, and Claus Michael Ringel, The numbers of support-tilting modules for a Dynkin algebra, 2014.
- OEIS Wiki, Binomial coefficients
- Richard L. Ollerton and Anthony G. Shannon, Some properties of generalized Pascal squares and triangles, Fib. Q., Vol. 36, No. 2 (1998), pp. 98-109.
- Ed Pegg, Jr., Sequence Pictures, Math Games column, Dec 08 2003.
- Ed Pegg, Jr., Sequence Pictures, Math Games column, Dec 08 2003. [Cached copy, with permission (pdf only)]
- Balak Ram, Common factors of n!/(m!(n-m)!), (m = 1, 2, ... n-1), Journal of the Indian Mathematical Club (Madras) 1 (1909), pp. 39-43.
- Franck Ramaharo, Statistics on some classes of knot shadows, arXiv:1802.07701 [math.CO], 2018.
- Franck Ramaharo, A generating polynomial for the pretzel knot, arXiv:1805.10680 [math.CO], 2018.
- Franck Ramaharo, A generating polynomial for the two-bridge knot with Conway's notation C(n,r), arXiv:1902.08989 [math.CO], 2019.
- Franck Ramaharo, A bracket polynomial for 2-tangle shadows, arXiv:2002.06672 [math.CO], 2020.
- Jack Ramsay, On Arithmetical Triangles, The Pulse of Long Island, June 1965 [Mentions application to design of antenna arrays. Annotated scan.]
- Thomas M. Richardson, The Reciprocal Pascal Matrix, arXiv preprint arXiv:1405.6315 [math.CO], 2014.
- Yuriy Shablya, Dmitry Kruchinin, and Vladimir Kruchinin, Method for Developing Combinatorial Generation Algorithms Based on AND/OR Trees and Its Application, Mathematics, Vol. 8, No. 6 (2020), 962.
- Louis W. Shapiro, Seyoum Getu, Wen-Jin Woan, and Leon C. Woodson, The Riordan group, Discrete Applied Math., Vol. 34 (1991), pp. 229-239.
- N. J. A. Sloane, My favorite integer sequences, in Sequences and their Applications (Proceedings of SETA '98).
- N. J. A. Sloane, Triangle showing silhouette of first 30 rows of Pascal's triangle (after Cobeli and Zaharescu)
- N. J. A. Sloane, The OEIS: A Fingerprint File for Mathematics, arXiv:2105.05111 [math.HO], 2021.
- N. J. A. Sloane, "A Handbook of Integer Sequences" Fifty Years Later, arXiv:2301.03149 [math.NT], 2023, p. 5.
- Hermann Stamm-Wilbrandt, Compute C(n+m,...) based on C(n,...) and C(m,...) values animation.
- Igor Victorovich Statsenko, On the ordinal numbers of triangles of generalized special numbers, Innovation science No 2-2, State Ufa, Aeterna Publishing House, 2024, pp. 15-19. In Russian.
- Christopher Stover and Eric W. Weisstein, Composition. From MathWorld - A Wolfram Web Resource.
- Gérard Villemin's Almanach of Numbers, Triangle de Pascal.
- Eric Weisstein's World of Mathematics, Pascal's Triangle.
- Wikipedia, Pascal's triangle.
- Herbert S. Wilf, Generatingfunctionology, 2nd edn., Academic Press, NY, 1994, pp. 12ff.
- Ken Williams, Mathforum, Interactive Pascal's Triangle.
- Doron Zeilberger, The Combinatorial Astrology of Rabbi Abraham Ibn Ezra, arXiv:math/9809136 [math.CO], 1998.
- Chris Zheng and Jeffrey Zheng, Triangular Numbers and Their Inherent Properties, Variant Construction from Theoretical Foundation to Applications, Springer, Singapore, 51-65.
- Index entries for triangles and arrays related to Pascal's triangle.
- Index entries for "core" sequences.
- Index entries for sequences related to Benford's law.
Crossrefs
Equals differences between consecutive terms of A102363. - David G. Williams (davidwilliams(AT)Paxway.com), Jan 23 2006
Row sums give A000079 (powers of 2).
Partial sums of rows give triangle A008949.
The triangle of the antidiagonals is A011973.
Another version: A108044.
Cf. A008277, A132311, A132312, A052216, A052217, A052218, A052219, A052220, A052221, A052222, A052223, A144225, A202750, A211226, A047999, A026729, A052553, A051920, A193242.
Triangle sums (see the comments): A000079 (Row1); A000007 (Row2); A000045 (Kn11 & Kn21); A000071 (Kn12 & Kn22); A001924 (Kn13 & Kn23); A014162 (Kn14 & Kn24); A014166 (Kn15 & Kn25); A053739 (Kn16 & Kn26); A053295 (Kn17 & Kn27); A053296 (Kn18 & Kn28); A053308 (Kn19 & Kn29); A053309 (Kn110 & Kn210); A001519 (Kn3 & Kn4); A011782 (Fi1 & Fi2); A000930 (Ca1 & Ca2); A052544 (Ca3 & Ca4); A003269 (Gi1 & Gi2); A055988 (Gi3 & Gi4); A034943 (Ze1 & Ze2); A005251 (Ze3 & Ze4). - Johannes W. Meijer, Sep 22 2010
Fibonacci-Pascal triangles: A027926, A036355, A037027, A074829, A105809, A109906, A111006, A114197, A162741, A228074, A228196, A228576.
Cf. A115940 (pandigital binomial coefficients C(m,k) with k>1).
Programs
-
Axiom
-- (start) )set expose add constructor OutputForm pascal(0,n) == 1 pascal(n,n) == 1 pascal(i,j | 0 < i and i < j) == pascal(i-1,j-1) + pascal(i,j-1) pascalRow(n) == [pascal(i,n) for i in 0..n] displayRow(n) == output center blankSeparate pascalRow(n) for i in 0..20 repeat displayRow i -- (end)
-
GAP
Flat(List([0..12],n->List([0..n],k->Binomial(n,k)))); # Stefano Spezia, Dec 22 2018
-
Haskell
a007318 n k = a007318_tabl !! n !! k a007318_row n = a007318_tabl !! n a007318_list = concat a007318_tabl a007318_tabl = iterate (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [1] -- Cf. http://www.haskell.org/haskellwiki/Blow_your_mind#Mathematical_sequences -- Reinhard Zumkeller, Nov 09 2011, Oct 22 2010
-
Magma
/* As triangle: */ [[Binomial(n, k): k in [0..n]]: n in [0.. 10]]; // Vincenzo Librandi, Jul 29 2015
-
Maple
A007318 := (n,k)->binomial(n,k);
-
Mathematica
Flatten[Table[Binomial[n, k], {n, 0, 11}, {k, 0, n}]] (* Robert G. Wilson v, Jan 19 2004 *) Flatten[CoefficientList[CoefficientList[Series[1/(1 - x - x*y), {x, 0, 12}], x], y]] (* Mats Granvik, Jul 08 2014 *)
-
Maxima
create_list(binomial(n,k),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */
-
PARI
C(n,k)=binomial(n,k) \\ Charles R Greathouse IV, Jun 08 2011
-
Python
# See Hobson link. Further programs: from math import prod,factorial def C(n,k): return prod(range(n,n-k,-1))//factorial(k) # M. F. Hasler, Dec 13 2019, updated Apr 29 2022, Feb 17 2023
-
Python
from math import comb, isqrt def A007318(n): return comb(r:=(m:=isqrt(k:=n+1<<1))-(k<=m*(m+1)),n-comb(r+1,2)) # Chai Wah Wu, Nov 11 2024
-
Sage
def C(n,k): return Subsets(range(n), k).cardinality() # Ralf Stephan, Jan 21 2014
Formula
a(n, k) = C(n,k) = binomial(n, k).
C(n, k) = C(n-1, k) + C(n-1, k-1).
The triangle is symmetric: C(n,k) = C(n,n-k).
a(n+1, m) = a(n, m) + a(n, m-1), a(n, -1) := 0, a(n, m) := 0, n
C(n, k) = n!/(k!(n-k)!) if 0<=k<=n, otherwise 0.
C(n, k) = ((n-k+1)/k) * C(n, k-1) with C(n, 0) = 1. - Michael B. Porter, Mar 23 2025
G.f.: 1/(1-y-x*y) = Sum_(C(n, k)*x^k*y^n, n, k>=0)
G.f.: 1/(1-x-y) = Sum_(C(n+k, k)*x^k*y^n, n, k>=0).
G.f. for row n: (1+x)^n = Sum_{k=0..n} C(n, k)*x^k.
G.f. for column k: x^k/(1-x)^(k+1); [corrected by Werner Schulte, Jun 15 2022].
E.g.f.: A(x, y) = exp(x+x*y).
E.g.f. for column n: x^n*exp(x)/n!.
In general the m-th power of A007318 is given by: T(0, 0) = 1, T(n, k) = T(n-1, k-1) + m*T(n-1, k), where n is the row-index and k is the column; also T(n, k) = m^(n-k)*C(n, k).
Triangle T(n, k) read by rows; given by A000007 DELTA A000007, where DELTA is Deléham's operator defined in A084938.
Let P(n+1) = the number of integer partitions of (n+1); let p(i) = the number of parts of the i-th partition of (n+1); let d(i) = the number of different parts of the i-th partition of (n+1); let m(i, j) = multiplicity of the j-th part of the i-th partition of (n+1). Define the operator Sum_{i=1..P(n+1), p(i)=k+1} as the sum running from i=1 to i=P(n+1) but taking only partitions with p(i)=(k+1) parts into account. Define the operator Product_{j=1..d(i)} = product running from j=1 to j=d(i). Then C(n, k) = Sum_{p(i)=(k+1), i=1..P(n+1)} p(i)! / [Product_{j=1..d(i)} m(i, j)!]. E.g., C(5, 3) = 10 because n=6 has the following partitions with m=3 parts: (114), (123), (222). For their multiplicities one has: (114): 3!/(2!*1!) = 3; (123): 3!/(1!*1!*1!) = 6; (222): 3!/3! = 1. The sum is 3 + 6 + 1 = 10 = C(5, 3). - Thomas Wieder, Jun 03 2005
C(n, k) = Sum_{j=0..k} (-1)^j*C(n+1+j, k-j)*A000108(j). - Philippe Deléham, Oct 10 2005
G.f.: 1 + x*(1 + x) + x^3*(1 + x)^2 + x^6*(1 + x)^3 + ... . - Michael Somos, Sep 16 2006
Sum_{k=0..floor(n/2)} x^(n-k)*T(n-k,k) = A000007(n), A000045(n+1), A002605(n), A030195(n+1), A057087(n), A057088(n), A057089(n), A057090(n), A057091(n), A057092(n), A057093(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, respectively. Sum_{k=0..floor(n/2)} (-1)^k*x^(n-k)*T(n-k,k) = A000007(n), A010892(n), A009545(n+1), A057083(n), A001787(n+1), A030191(n), A030192(n), A030240(n), A057084(n), A057085(n+1), A057086(n), A084329(n+1) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, respectively. - Philippe Deléham, Sep 16 2006
C(n,k) <= A062758(n) for n > 1. - Reinhard Zumkeller, Mar 04 2008
C(t+p-1, t) = Sum_{i=0..t} C(i+p-2, i) = Sum_{i=1..p} C(i+t-2, t-1). A binomial number is the sum of its left parent and all its right ancestors, which equals the sum of its right parent and all its left ancestors. - Lee Naish (lee(AT)cs.mu.oz.au), Mar 07 2008
From Paul D. Hanna, Mar 24 2011: (Start)
Let A(x) = Sum_{n>=0} x^(n*(n+1)/2)*(1+x)^n be the g.f. of the flattened triangle:
A(x) = 1 + (x + x^2) + (x^3 + 2*x^4 + x^5) + (x^6 + 3*x^7 + 3*x^8 + x^9) + ...
then A(x) equals the series Sum_{n>=0} (1+x)^n*x^n*Product_{k=1..n} (1-(1+x)*x^(2*k-1))/(1-(1+x)*x^(2*k));
also, A(x) equals the continued fraction 1/(1- x*(1+x)/(1+ x*(1-x)*(1+x)/(1- x^3*(1+x)/(1+ x^2*(1-x^2)*(1+x)/(1- x^5*(1+x)/(1+ x^3*(1-x^3)*(1+x)/(1- x^7*(1+x)/(1+ x^4*(1-x^4)*(1+x)/(1- ...))))))))).
These formulas are due to (1) a q-series identity and (2) a partial elliptic theta function expression. (End)
Row n of the triangle is the result of applying the ConvOffs transform to the first n terms of the natural numbers (1, 2, 3, ..., n). See A001263 or A214281 for a definition of this transformation. - Gary W. Adamson, Jul 12 2012
From L. Edson Jeffery, Aug 02 2012: (Start)
Row n (n >= 0) of the triangle is given by the n-th antidiagonal of the infinite matrix P^n, where P = (p_{i,j}), i,j >= 0, is the production matrix
0, 1,
1, 0, 1,
0, 1, 0, 1,
0, 0, 1, 0, 1,
0, 0, 0, 1, 0, 1,
0, 0, 0, 0, 1, 0, 1,
0, 0, 0, 0, 0, 1, 0, 1,
0, 0, 0, 0, 0, 0, 1, 0, 1,
... (End)
Row n of the triangle is also given by the n+1 coefficients of the polynomial P_n(x) defined by the recurrence P_0(x) = 1, P_1(x) = x + 1, P_n(x) = x*P_{n-1}(x) + P_{n-2}(x), n > 1. - L. Edson Jeffery, Aug 12 2013
For a closed-form formula for arbitrary left and right borders of Pascal-like triangles see A228196. - Boris Putievskiy, Aug 18 2013
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 04 2013
(1+x)^n = Sum_{k=0..n} (-1)^(n-k)*binomial(n,k)*Sum_{i=0..k} k^(n-i)*binomial(k,i)*x^(n-i)/(n-i)!. - Vladimir Kruchinin, Oct 21 2013
E.g.f.: A(x,y) = exp(x+x*y) = 1 + (x+y*x)/( E(0)-(x+y*x)), where E(k) = 1 + (x+y*x)/(1 + (k+1)/E(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Nov 08 2013
E.g.f.: E(0) -1, where E(k) = 2 + x*(1+y)/(2*k+1 - x*(1+y)/E(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Dec 24 2013
G.f.: 1 + x*(1+x)*(1+x^2*(1+x)/(W(0)-x^2-x^3)), where W(k) = 1 + (1+x)*x^(k+2) - (1+x)*x^(k+3)/W(k+1); (continued fraction). - Sergei N. Gladkovskii, Dec 24 2013
Sum_{n>=0} C(n,k)/n! = e/k!, where e = exp(1), while allowing n < k where C(n,k) = 0. Also Sum_{n>=0} C(n+k-1,k)/n! = e * A000262(k)/k!, and for k>=1 equals e * A067764(k)/A067653(k). - Richard R. Forberg, Jan 01 2014
Sum_{n>=k} 1/C(n,k) = k/(k-1) for k>=1. - Richard R. Forberg, Feb 10 2014
From Tom Copeland, Apr 26 2014: (Start)
Multiply each n-th diagonal of the Pascal lower triangular matrix by x^n and designate the result by A007318(x) = P(x). Then with :xD:^n = x^n*(d/dx)^n and B(n,x), the Bell polynomials (A008277),
A) P(x)= exp(x*dP) = exp[x*(e^M-I)] = exp[M*B(.,x)] = (I+dP)^B(.,x)
B) P(:xD:) = exp(dP:xD:) = exp[(e^M-I):xD:] = exp[M*B(.,:xD:)] = exp[M*xD] = (I+dP)^(xD) with action P(:xD:)g(x) = exp(dP:xD:)g(x) = g[(I+dP)*x] (cf. also A238363).
C) P(x)^y = P(y*x). P(2x) = A038207(x) = exp[M*B(.,2x)], the face vectors of the n-dim hypercubes.
D) P(x) = [St2]*exp(x*M)*[St1] = [St2]*(I+dP)^x*[St1]
E) = [St1]^(-1)*(I+dP)^x*[St1] = [St2]*(I+dP)^x*[St2]^(-1)
where [St1]=padded A008275 just as [St2]=A048993=padded A008277 and exp(x*M) = (I+dP)^x = Sum_{k>=0} C(x,k) dP^k. (End)
From Peter Bala, Dec 21 2014: (Start)
Recurrence equation: T(n,k) = T(n-1,k)*(n + k)/(n - k) - T(n-1,k-1) for n >= 2 and 1 <= k < n, with boundary conditions T(n,0) = T(n,n) = 1. Note, changing the minus sign in the recurrence to a plus sign gives a recurrence for the square of the binomial coefficients - see A008459.
There is a relation between the e.g.f.'s of the rows and the diagonals of the triangle, namely, exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(1 + 3*x + 3*x^2/2! + x^3/3!) = 1 + 4*x + 10*x^2/2! + 20*x^3/3! + 35*x^4/4! + .... This property holds more generally for the Riordan arrays of the form ( f(x), x/(1 - x) ), where f(x) is an o.g.f. of the form 1 + f_1*x + f_2*x^2 + .... See, for example, A055248 and A106516.
Let P denote the present triangle. For k = 0,1,2,... define P(k) to be the lower unit triangular block array
/I_k 0\
\ 0 P/ having the k X k identity matrix I_k as the upper left block; in particular, P(0) = P. The infinite product P(0)*P(1)*P(2)*..., which is clearly well-defined, is equal to the triangle of Stirling numbers of the second kind A008277. The infinite product in the reverse order, that is, ...*P(2)*P(1)*P(0), is equal to the triangle of Stirling cycle numbers A130534. (End)
C(a+b,c) = Sum_{k=0..a} C(a,k)*C(b,b-c+k). This is a generalization of equation 1 from section 4.2.5 of the Prudnikov et al. reference, for a=b=c=n: C(2*n,n) = Sum_{k=0..n} C(n,k)^2. See Links section for animation of new formula. - Hermann Stamm-Wilbrandt, Aug 26 2015
The row polynomials of the Pascal matrix P(n,x) = (1+x)^n are related to the Bernoulli polynomials Br(n,x) and their umbral compositional inverses Bv(n,x) by the umbral relation P(n,x) = (-Br(.,-Bv(.,x)))^n = (-1)^n Br(n,-Bv(.,x)), which translates into the matrix relation P = M * Br * M * Bv, where P is the Pascal matrix, M is the diagonal matrix diag(1,-1,1,-1,...), Br is the matrix for the coefficients of the Bernoulli polynomials, and Bv that for the umbral inverse polynomials defined umbrally by Br(n,Bv(.,x)) = x^n = Bv(n,Br(.,x)). Note M = M^(-1). - Tom Copeland, Sep 05 2015
1/(1-x)^k = (r(x) * r(x^2) * r(x^4) * ...) where r(x) = (1+x)^k. - Gary W. Adamson, Oct 17 2016
Boas-Buck type recurrence for column k for Riordan arrays (see the Aug 10 2017 remark in A046521, also for the reference) with the Boas-Buck sequence b(n) = {repeat(1)}. T(n, k) = ((k+1)/(n-k))*Sum_{j=k..n-1} T(j, k), for n >= 1, with T(n, n) = 1. This reduces, with T(n, k) = binomial(n, k), to a known binomial identity (e.g, Graham et al. p. 161). - Wolfdieter Lang, Nov 12 2018
C((p-1)/a, b) == (-1)^b * fact_a(a*b-a+1)/fact_a(a*b) (mod p), where fact_n denotes the n-th multifactorial, a divides p-1, and the denominator of the fraction on the right side of the equation represents the modular inverse. - Isaac Saffold, Jan 07 2019
C(n,k-1) = A325002(n,k) - [k==n+1] = (A325002(n,k) + A325003(n,k)) / 2 = [k==n+1] + A325003(n,k). - Robert A. Russell, Oct 20 2020
From Hermann Stamm-Wilbrandt, May 13 2021: (Start)
Binomial sums are Fibonacci numbers A000045:
Sum_{k=0..n} C(n + k, 2*k + 1) = F(2*n).
Sum_{k=0..n} C(n + k, 2*k) = F(2*n + 1). (End)
C(n,k) = Sum_{i=0..k} A000108(i) * C(n-2i-1, k-i), for 0 <= k <= floor(n/2)-1. - Tushar Bansal, May 17 2025
Extensions
Checked all links, deleted 8 that seemed lost forever and were probably not of great importance. - N. J. A. Sloane, May 08 2018
A027649 a(n) = 2*(3^n) - 2^n.
1, 4, 14, 46, 146, 454, 1394, 4246, 12866, 38854, 117074, 352246, 1058786, 3180454, 9549554, 28665046, 86027906, 258149254, 774578834, 2323998646, 6972520226, 20918609254, 62757924914, 188277969046, 564842295746, 1694543664454, 5083664547794, 15251060752246
Offset: 0
Comments
Poly-Bernoulli numbers B_n^(k) with k=-2.
Binomial transform of A007051, if both sequences start at 0. Binomial transform of A000225(n+1). - Paul Barry, Mar 24 2003
Euler expands (1-z)/(1-5z+6z^2) and finds the general term. Section 226 of the Introductio indicates that he could have written down the recursion relation: a(n) = 5 a(n-1)-6 a(n-2). - V. Frederick Rickey (fred-rickey(AT)usma.edu), Feb 10 2006
Let R be a binary relation on the power set P(A) of a set A having n = |A| elements such that for every element x, y of P(A), xRy if x is a subset of y or y is a subset of x. Then a(n) = |R|. - Ross La Haye, Dec 22 2006
With regard to the comment by Ross La Haye: For proper subsets see A056182. - For nonempty subsets see A091344. - For nonempty proper subsets see a(n+1) in A260217. - Manfred Boergens, Aug 02 2023
If x, y are two n-bit binary strings then a(n) gives the number of pairs (x,y) such that XOR(x, y) = ABS(x - y). - Ramasamy Chandramouli, Feb 15 2009
Equals row sums of the triangular version of A038573. - Gary W. Adamson, Jun 04 2009
Inverse binomial transform of A085350. - Paul Curtz, Nov 14 2009
Related to the number of even a's in a nontrivial cycle (should one exist) in the 3x+1 Problem, where a <= floor(log_2(2*(3^n) - 2^n)). The value n correlates to the number of odds in such a nontrivial cycle. See page 1288 of Crandall's paper. Also, this relation gives another proof that the number of odds divided by the number of evens in a nontrivial cycle is bounded by log 2 / log 3 (this observation does not resolve the finite cycles conjecture as the value could be arbitrarily close to this bound). However, the same argument gives that log 2 / log 3 is less than or equal to the number of odds divided by the number of evens in a divergent sequence (should one exist), as log 2 / log 3 is the limit value for a cycle of an arbitrarily large length, where the length is given by the value n. - Jeffrey R. Goodwin, Aug 04 2011
Row sums of Riordan triangle A106516. - Wolfdieter Lang, Jan 09 2015
Number of restricted barred preferential arrangements having 3 bars in which the sections are all restricted sections such that (for fixed sections i and j) section i or section j is empty. - Sithembele Nkonkobe, Oct 12 2015
This is also row 2 of A281891: for n >= 1, when consecutive positive integers are written as a product of primes in nondecreasing order, a factor of 2 or 3 occurs in n-th position a(n) times out of every 6^n. - Peter Munn, May 18 2017
Also row sums of A124929. - Omar E. Pol, Jun 15 2017
This is the sum of A318921(n) for n in the range 2^(k+1) to 2^(k+2)-1. See A318921 for proof. - N. J. A. Sloane, Sep 25 2018
a(n) is also the number of acyclic orientations of the complete bipartite graph K_{2,n}. - Vincent Pilaud, Sep 15 2020
a(n-1) is also the number of n-digit numbers whose largest decimal digit is 2. - Stefano Spezia, Nov 15 2023
References
- Leonhard Euler, Introductio in analysin infinitorum (1748), section 216.
Links
- T. D. Noe, Table of n, a(n) for n = 0..200
- Taylor Brysiewicz, Holger Eble, and Lukas Kühne, Enumerating chambers of hyperplane arrangements with symmetry, arXiv:2105.14542 [math.CO], 2021.
- R. E. Crandall, On the 3x+1 problem, Math. Comp., 32 (1978) 1281-1292.
- Shalosh B. Ekhad, N. J. A. Sloane, and Doron Zeilberger, A Meta-Algorithm for Creating Fast Algorithms for Counting ON Cells in Odd-Rule Cellular Automata, arXiv:1503.01796 [math.CO], 2015; see also the Accompanying Maple Package.
- Shalosh B. Ekhad, N. J. A. Sloane, and Doron Zeilberger, Odd-Rule Cellular Automata on the Square Grid, arXiv:1503.04249 [math.CO], 2015.
- John Elias, Illustration of initial terms: Reflected Sierpinski triangle
- K. Imatomi, M. Kaneko, and E. Takeda, Multi-Poly-Bernoulli Numbers and Finite Multiple Zeta Values, J. Int. Seq. 17 (2014) # 14.4.5
- Ken Kamano, Sums of Products of Poly-Bernoulli Numbers of Negative Index, Journal of Integer Sequences, Vol. 15 (2012), #12.1.3.
- K. Kamano, Sums of Products of Bernoulli Numbers, Including Poly-Bernoulli Numbers, J. Int. Seq. 13 (2010), 10.5.2.
- Masanobu Kaneko, Poly-Bernoulli numbers, Journal de théorie des nombres de Bordeaux, 9 no. 1 (1997), Pages 221-228.
- Takao Komatsu, Some recurrence relations of poly-Cauchy numbers, J. Nonlinear Sci. Appl., (2019) Vol. 12, Issue 12, 829-845.
- Ross La Haye, Binary Relations on the Power Set of an n-Element Set, Journal of Integer Sequences, Vol. 12 (2009), Article 09.2.6.
- S. Nkonkobe and V. Murali, On some properties and relations between restricted barred preferential arrangements, multi-poly-bernoulli numbers and related numbers, arXiv:1509.07352 [math.CO], 2015.
- Eric Weisstein's World of Mathematics, Sierpinski Sieve
- Wikipedia, Sierpinski triangle
- Index entries for sequences related to Bernoulli numbers.
- Index entries for linear recurrences with constant coefficients, signature (5,-6).
Crossrefs
Programs
-
Haskell
a027649 n = a027649_list !! n a027649_list = map fst $ iterate (\(u, v) -> (3 * u + v, 2 * v)) (1, 1) -- Reinhard Zumkeller, Jun 09 2013
-
Magma
[2*(3^n)-2^n: n in [0..30]]; // Vincenzo Librandi, Jul 17 2011
-
Maple
a(n, k):= (-1)^n*sum( (-1)^'m'*'m'!*Stirling2(n,'m')/('m'+1)^k,'m'=0..n); seq(a(n, -2), n=0..30);
-
Mathematica
Table[2(3^n)-2^n,{n,0,30}] (* or *) LinearRecurrence[ {5,-6},{1,4},31] (* Harvey P. Dale, Apr 22 2011 *)
-
PARI
a(n)=2*(3^n)-2^n \\ Charles R Greathouse IV, Jul 16 2011
-
PARI
Vec((1-x)/((1-2*x)*(1-3*x)) + O(x^50)) \\ Altug Alkan, Oct 12 2015
-
SageMath
[2*(3^n - 2^(n-1)) for n in (0..30)] # G. C. Greubel, Aug 01 2022
Formula
G.f.: (1-x)/((1-2*x)*(1-3*x)).
a(n) = 3*a(n-1) + 2^(n-1), with a(0) = 1.
a(n) = Sum_{k=0..n} binomial(n, k)*(2^(k+1) - 1). - Paul Barry, Mar 24 2003
Partial sums of A053581. - Paul Barry, Jun 26 2003
Main diagonal of array (A085870) defined by T(i, 1) = 2^i - 1, T(1, j) = 2^j - 1, T(i, j) = T(i-1, j) + T(i-1, j-1). - Benoit Cloitre, Aug 05 2003
a(n) = A090888(n, 3). - Ross La Haye, Sep 21 2004
a(n) = Sum_{k=0..n} binomial(n+2, k+1)*Sum_{j=0..floor(k/2)} A001045(k-2j). - Paul Barry, Apr 17 2005
a(n) = Sum_{k=0..n} Sum_{j=0..n} binomial(n,j)*binomial(j+1,k+1). - Paul Barry, Sep 18 2006
a(n) = A166060(n+1)/6. - Philippe Deléham, Oct 21 2009
a(n) = 5*a(n-1) - 6*a(n-2), a(0)=1, a(1)=4. - Harvey P. Dale, Apr 22 2011
a(n) = A217764(n,2). - Ross La Haye, Mar 27 2013
For n>0, a(n) = 3 * a(n-1) + 2^(n-1) = 2 * (a(n-1) + 3^(n-1)). - J. Conrad, Oct 29 2015
for n>0, a(n) = 2 * (1 + 2^(n-2) + Sum_{x=1..n-2} Sum_{k=0..x-1} (binomial(x-1,k)*(2^(k+1) + 2^(n-x+k)))). - J. Conrad, Dec 10 2015
E.g.f.: exp(2*x)*(2*exp(x) - 1). - Stefano Spezia, May 18 2024
Extensions
Better formulas from David W. Wilson and Michael Somos
Incorrect formula removed by Charles R Greathouse IV, Mar 18 2010
Duplications (due to corrections to A numbers) removed by Peter Munn, Jun 15 2017
A228196 A triangle formed like Pascal's triangle, but with n^2 on the left border and 2^n on the right border instead of 1.
0, 1, 2, 4, 3, 4, 9, 7, 7, 8, 16, 16, 14, 15, 16, 25, 32, 30, 29, 31, 32, 36, 57, 62, 59, 60, 63, 64, 49, 93, 119, 121, 119, 123, 127, 128, 64, 142, 212, 240, 240, 242, 250, 255, 256, 81, 206, 354, 452, 480, 482, 492, 505, 511, 512, 100, 287, 560, 806, 932, 962, 974, 997, 1016, 1023, 1024
Offset: 1
Comments
The third row is (n^4 - n^2 + 24*n + 24)/12.
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 04 2013
Examples
The start of the sequence as a triangular array read by rows: 0; 1, 2; 4, 3, 4; 9, 7, 7, 8; 16, 16, 14, 15, 16; 25, 32, 30, 29, 31, 32; 36, 57, 62, 59, 60, 63, 64;
Links
- Boris Putievskiy, Rows n = 1..140 of triangle, flattened
- Rely Pellicer and David Alvo, Modified Pascal Triangle and Pascal Surfaces p.4
- Index entries for triangles and arrays related to Pascal's triangle
Crossrefs
Cf. We denote Pascal-like triangle with L(n) on the left border and R(n) on the right border by (L(n),R(n)). A007318 (1,1), A008949 (1,2^n), A029600 (2,3), A029618 (3,2), A029635 (1,2), A029653 (2,1), A037027 (Fibonacci(n),1), A051601 (n,n) n>=0, A051597 (n,n) n>0, A051666 (n^2,n^2), A071919 (1,0), A074829 (Fibonacci(n), Fibonacci(n)), A074909 (1,n), A093560 (3,1), A093561 (4,1), A093562 (5,1), A093563 (6,1), A093564 (7,1), A093565 (8,1), A093644 (9,1), A093645 (10,1), A095660 (1,3), A095666 (1,4), A096940 (1,5), A096956 (1,6), A106516 (3^n,1), A108561(1,(-1)^n), A132200 (4,4), A134636 (2n+1,2n+1), A137688 (2^n,2^n), A160760 (3^(n-1),1), A164844(1,10^n), A164847 (100^n,1), A164855 (101*100^n,1), A164866 (101^n,1), A172171 (1,9), A172185 (9,11), A172283 (-9,11), A177954 (int(n/2),1), A193820 (1,2^n), A214292 (n,-n), A227074 (4^n,4^n), A227075 (3^n,3^n), A227076 (5^n,5^n), A227550 (n!,n!), A228053 ((-1)^n,(-1)^n), A228074 (Fibonacci(n), n).
Programs
-
GAP
T:= function(n,k) if k=0 then return n^2; elif k=n then return 2^n; else return T(n-1,k-1) + T(n-1,k); fi; end; Flat(List([0..12], n-> List([0..n], k-> T(n,k) ))); # G. C. Greubel, Nov 12 2019
-
Maple
T:= proc(n, k) option remember; if k=0 then n^2 elif k=n then 2^k else T(n-1, k-1) + T(n-1, k) fi end: seq(seq(T(n, k), k=0..n), n=0..10); # G. C. Greubel, Nov 12 2019
-
Mathematica
T[n_, k_]:= T[n, k] = If[k==0, n^2, If[k==n, 2^k, T[n-1, k-1] + T[n-1, k]]]; Table[T[n, k], {n,0,10}, {k,0,n}]//Flatten (* G. C. Greubel, Nov 12 2019 *) Flatten[Table[Sum[i^2 Binomial[n-1-i, n-k-i], {i,1,n-k}] + Sum[2^i Binomial[n-1-i, k-i], {i,1,k}], {n,0,10}, {k,0,n}]] (* Greg Dresden, Aug 06 2022 *)
-
PARI
T(n,k) = if(k==0, n^2, if(k==n, 2^k, T(n-1, k-1) + T(n-1, k) )); \\ G. C. Greubel, Nov 12 2019
-
Python
def funcL(n): q = n**2 return q def funcR(n): q = 2**n return q for n in range (1,9871): t=int((math.sqrt(8*n-7) - 1)/ 2) i=n-t*(t+1)/2-1 j=(t*t+3*t+4)/2-n-1 sum1=0 sum2=0 for m1 in range (1,i+1): sum1=sum1+funcR(m1)*binomial(i+j-m1-1,i-m1) for m2 in range (1,j+1): sum2=sum2+funcL(m2)*binomial(i+j-m2-1,j-m2) sum=sum1+sum2
-
Sage
@CachedFunction def T(n, k): if (k==0): return n^2 elif (k==n): return 2^n else: return T(n-1, k-1) + T(n-1, k) [[T(n, k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Nov 12 2019
Formula
T(n,0) = n^2, n>0; T(0,k) = 2^k; T(n, k) = T(n-1, k-1) + T(n-1, k) for n,k > 0. [corrected by G. C. Greubel, Nov 12 2019]
Closed-form formula for general case. Let L(m) and R(m) be the left border and the right border of Pascal like triangle, respectively. We denote binomial(n,k) by C(n,k).
As table read by antidiagonals T(n,k) = Sum_{m1=1..n} R(m1)*C(n+k-m1-1, n-m1) + Sum_{m2=1..k} L(m2)*C(n+k-m2-1, k-m2); n,k >=0.
As linear sequence a(n) = Sum_{m1=1..i} R(m1)*C(i+j-m1-1, i-m1) + Sum_{m2=1..j} L(m2)*C(i+j-m2-1, j-m2), where i=n-t*(t+1)/2-1, j=(t*t+3*t+4)/2-n-1, t=floor((-1+sqrt(8*n-7))/2); n>0.
Some special cases. If L(m)={b,b,b...} b*A000012, then the second sum takes form b*C(n+k-1,j). If L(m) is {0,b,2b,...} b*A001477, then the second sum takes form b*C(n+k,n-1). Similarly for R(m) and the first sum.
For this sequence L(m)=m^2 and R(m)=2^m.
As table read by antidiagonals T(n,k) = Sum_{m1=1..n} (2^m1)*C(n+k-m1-1, n-m1) + Sum_{m2=1..k} (m2^2)*C(n+k-m2-1, k-m2); n,k >=0.
As linear sequence a(n) = Sum_{m1=1..i} (2^m1)*C(i+j-m1-1, i-m1) + Sum_{m2=1..j} (m2^2)*C(i+j-m2-1, j-m2), where i=n-t*(t+1)/2-1, j=(t*t+3*t+4)/2-n-1, t=floor((-1+sqrt(8*n-7))/2).
As a triangular array read by rows, T(n,k) = Sum_{i=1..n-k} i^2*C(n-1-i, n-k-i) + Sum_{i=1..k} 2^i*C(n-1-i, k-i); n,k >=0. - Greg Dresden, Aug 06 2022
Extensions
Cross-references corrected and extended by Philippe Deléham, Dec 27 2013
A055248 Triangle of partial row sums of triangle A007318(n,m) (Pascal's triangle). Triangle A008949 read backwards. Riordan (1/(1-2x), x/(1-x)).
1, 2, 1, 4, 3, 1, 8, 7, 4, 1, 16, 15, 11, 5, 1, 32, 31, 26, 16, 6, 1, 64, 63, 57, 42, 22, 7, 1, 128, 127, 120, 99, 64, 29, 8, 1, 256, 255, 247, 219, 163, 93, 37, 9, 1, 512, 511, 502, 466, 382, 256, 130, 46, 10, 1, 1024, 1023, 1013, 968, 848, 638, 386, 176, 56, 11, 1
Offset: 0
Comments
In the language of the Shapiro et al. reference (also given in A053121) such a lower triangular (ordinary) convolution array, considered as matrix, belongs to the Riordan-group. The g.f. for the row polynomials p(n,x) (increasing powers of x) is 1/((1-2*z)*(1-x*z/(1-z))).
Binomial transform of the all 1's triangle: as a Riordan array, it factors to give (1/(1-x),x/(1-x))(1/(1-x),x). Viewed as a number square read by antidiagonals, it has T(n,k) = Sum_{j=0..n} binomial(n+k,n-j) and is then the binomial transform of the Whitney square A004070. - Paul Barry, Feb 03 2005
Riordan array (1/(1-2x), x/(1-x)). Antidiagonal sums are A027934(n+1), n >= 0. - Paul Barry, Jan 30 2005; edited by Wolfdieter Lang, Jan 09 2015
Eigensequence of the triangle = A005493: (1, 3, 10, 37, 151, 674, ...); row sums of triangles A011971 and A159573. - Gary W. Adamson, Apr 16 2009
Read as a square array, this is the generalized Riordan array ( 1/(1 - 2*x), 1/(1 - x) ) as defined in the Bala link (p. 5), which factorizes as ( 1/(1 - x), x/(1 - x) )*( 1/(1 - x), x )*( 1, 1 + x ) = P*U*transpose(P), where P denotes Pascal's triangle, A007318, and U is the lower unit triangular array with 1's on or below the main diagonal. - Peter Bala, Jan 13 2016
Examples
The triangle a(n,m) begins: n\m 0 1 2 3 4 5 6 7 8 9 10 ... 0: 1 1: 2 1 2: 4 3 1 3: 8 7 4 1 4: 16 15 11 5 1 5: 32 31 26 16 6 1 6: 64 63 57 42 22 7 1 7: 128 127 120 99 64 29 8 1 8: 256 255 247 219 163 93 37 9 1 9: 512 511 502 466 382 256 130 46 10 1 10: 1024 1023 1013 968 848 638 386 176 56 11 1 ... Reformatted. - _Wolfdieter Lang_, Jan 09 2015 Fourth row polynomial (n=3): p(3,x)= 8 + 7*x + 4*x^2 + x^3. The matrix inverse starts 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; ... which may be related to A029653. - _R. J. Mathar_, Mar 29 2013 From _Peter Bala_, Dec 23 2014: (Start) 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 | |4 3 1 ||0 2 1 ||0 0 1 |... = |4 5 1 | |8 7 4 1 ||0 4 3 1 ||0 0 2 1 | |8 19 9 1 | |... ||0 8 7 4 1 ||0 0 4 3 1| |... | |... ||... ||... | | | = A143494. (End) Matrix factorization of square array as P*U*transpose(P): /1 \ /1 \ /1 1 1 1 ...\ /1 1 1 1 ...\ |1 1 ||1 1 ||0 1 2 3 ... | |2 3 4 5 ... | |1 2 1 ||1 1 1 ||0 0 1 3 ... | = |4 7 11 16 ... | |1 3 3 1 ||1 1 1 1 ||0 0 0 1 ... | |8 15 26 42 ... | |... ||... ||... | |... | - _Peter Bala_, Jan 13 2016
Links
- Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
- Peter Bala, Notes on generalized Riordan arrays
- Peter Bala, A055248: Rapidly converging series for log(2) and Pi
- Jean-Luc Baril, Javier F. González, and José L. Ramírez, Last symbol distribution in pattern avoiding Catalan words, Univ. Bourgogne (France, 2022).
- Paul Barry, On a Central Transform of Integer Sequences, arXiv:2004.04577 [math.CO], 2020.
- Norman Lindquist and Gerard Sierksma, Extensions of set partitions, Journal of Combinatorial Theory, Series A 31.2 (1981): 190-198. See Table I.
- L. W. Shapiro, S. Getu, Wen-Jin Woan and L. C. Woodson, The Riordan Group, Discrete Appl. Maths. 34 (1991) 229-239.
Crossrefs
Programs
-
Haskell
a055248 n k = a055248_tabl !! n !! k a055248_row n = a055248_tabl !! n a055248_tabl = map reverse a008949_tabl -- Reinhard Zumkeller, Jun 20 2015
-
Maple
T := (n,k) -> 2^n - (1/2)*binomial(n, k-1)*hypergeom([1, n + 1], [n-k + 2], 1/2). seq(seq(simplify(T(n,k)), k=0..n),n=0..10); # Peter Luschny, Oct 10 2019
-
Mathematica
a[n_, m_] := Sum[ Binomial[n, m + j], {j, 0, n}]; Table[a[n, m], {n, 0, 10}, {m, 0, n}] // Flatten (* Jean-François Alcover, Jul 05 2013, after Paul Barry *) T[n_, k_] := Binomial[n, k] * Hypergeometric2F1[1, k - n, k + 1, -1]; Flatten[Table[T[n, k], {n, 0, 7}, {k, 0, n}]] (* Peter Luschny, Oct 06 2023 *)
Formula
a(n, m) = A008949(n, n-m), if n > m >= 0.
a(n, m) = Sum_{k=m..n} A007318(n, k) (partial row sums in columns m).
Column m recursion: a(n, m) = Sum_{j=m..n-1} a(j, m) + A007318(n, m) if n >= m >= 0, a(n, m) := 0 if n
G.f. for column m: (1/(1-2*x))*(x/(1-x))^m, m >= 0.
a(n, m) = Sum_{j=0..n} binomial(n, m+j). - Paul Barry, Feb 03 2005
Inverse binomial transform (by columns) of A112626. - Ross La Haye, Dec 31 2006
T(2n,n) = A032443(n). - Philippe Deléham, Sep 16 2009
From Peter Bala, Dec 23 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)*(8 + 7*x + 4*x^2/2! + x^3/3!) = 8 + 15*x + 26*x^2/2! + 42*x^3/3! + 64*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ).
Let M denote the present triangle. 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. The infinite product M(0)*M(1)*M(2)*..., which is clearly well-defined, is equal to A143494 (but with a different offset). See the Example section. Cf. A106516. (End)
a(n,m) = Sum_{p=m..n} 2^(n-p)*binomial(p-1,m-1), n >= m >= 0, else 0. - Wolfdieter Lang, Jan 09 2015
T(n, k) = 2^n - (1/2)*binomial(n, k-1)*hypergeom([1, n+1], [n-k+2], 1/2). - Peter Luschny, Oct 10 2019
T(n, k) = binomial(n, k)*hypergeom([1, k - n], [k + 1], -1). - Peter Luschny, Oct 06 2023
n-th row polynomial R(n, x) = (2^n - x*(1 + x)^n)/(1 - x). These polynomials can be used to find series acceleration formulas for the constants log(2) and Pi. - Peter Bala, Mar 03 2025
A000340 a(0)=1, a(n) = 3*a(n-1) + n + 1.
1, 5, 18, 58, 179, 543, 1636, 4916, 14757, 44281, 132854, 398574, 1195735, 3587219, 10761672, 32285032, 96855113, 290565357, 871696090, 2615088290, 7845264891, 23535794695, 70607384108, 211822152348, 635466457069
Offset: 0
Comments
From Johannes W. Meijer, Feb 20 2009: (Start)
Second right hand column (n-m=1) of the A156920 triangle.
The generating function of this sequence enabled the analysis of the polynomials A156921 and A156925.
(End)
Partial sums of A003462, and thus the second partial sums of A000244 (3^n). Also column k=2 of A106516. - John Keith, Jan 04 2022
Examples
G.f. = 1 + 5*x + 18*x^2 + 58*x^3 + 179*x^4 + 543*x^5 + 1636*x^6 + ...
References
- F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 260.
- 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).
Links
- Vincenzo Librandi, Table of n, a(n) for n = 0..1000
- Dillan Agrawal, Selena Ge, Jate Greene, Tanya Khovanova, Dohun Kim, Rajarshi Mandal, Tanish Parida, Anirudh Pulugurtha, Gordon Redwine, Soham Samanta, and Albert Xu, Chip-Firing on Infinite k-ary Trees, arXiv:2501.06675 [math.CO], 2025. See pp. 9, 18.
- Shaoshi Chen, Hanqian Fang, Sergey Kitaev, and Candice X.T. Zhang, Patterns in Multi-dimensional Permutations, arXiv:2411.02897 [math.CO], 2024. See p. 7.
- INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 389
- Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
- Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
- László Tóth, On Schizophrenic Patterns in b-ary Expansions of Some Irrational Numbers, arXiv:2002.06584 [math.NT], 2020. Mentions this sequence. See also Proc. Amer. Math. Soc. 148 (2020), 461-469.
- Index entries for linear recurrences with constant coefficients, signature (5,-7,3).
Crossrefs
From Johannes W. Meijer, Feb 20 2009: (Start)
Equals A156920 second right hand column.
Equals A142963 second right hand column divided by 2^n.
Equals A156919 second right hand column divided by 2.
(End)
Cf. A014915.
Equals column k=1 of A008971 (shifted). - Jeremy Dover, Jul 11 2021
Programs
-
Magma
[(3^(n+2)-2*n-5)/4: n in [0..30]]; // Vincenzo Librandi, Aug 15 2011
-
Maple
a[ -1]:=0:a[0]:=1:for n from 1 to 50 do a[n]:=4*a[n-1]-3*a[n-2]+1 od: seq(a[n],n=0..50); # Miklos Kristof, Mar 09 2005 A000340:=-1/(3*z-1)/(z-1)**2; # conjectured by Simon Plouffe in his 1992 dissertation
-
Mathematica
a[ n_] := MatrixPower[ {{1, 0, 0}, {1, 1, 0}, {1, 1, 3}}, n + 1][[3, 1]]; (* Michael Somos, May 28 2014 *) RecurrenceTable[{a[0]==1,a[n]==3a[n-1]+n+1},a,{n,30}] (* or *) LinearRecurrence[{5,-7,3},{1,5,18},30] (* Harvey P. Dale, Jan 31 2017 *)
Formula
G.f.: 1/((1-3*x)*(1-x)^2).
a(n) = (3^(n+2) - 2*n - 5)/4.
a(n) = Sum_{k=0..n+1} (n-k+1)*3^k = Sum_{k=0..n+1} k*3^(n-k+1). - Paul Barry, Jul 30 2004
a(n) = Sum_{k=0..n} binomial(n+2, k+2)*2^k. - Paul Barry, Jul 30 2004
a(-1)=0, a(0)=1, a(n) = 4*a(n-1) - 3*a(n-2) + 1. - Miklos Kristof, Mar 09 2005
a(n) = 5*a(n-1) - 7*a(n-2) + 3*a(n-3). - Johannes W. Meijer, Feb 20 2009
a(-2 - n) = 3^-n * A014915(n). - Michael Somos, May 28 2014
E.g.f.: exp(x)*(9*exp(2*x) - 2*x - 5)/4. - Stefano Spezia, Nov 09 2024
A143495 Triangle read by rows: 3-Stirling numbers of the second kind.
1, 3, 1, 9, 7, 1, 27, 37, 12, 1, 81, 175, 97, 18, 1, 243, 781, 660, 205, 25, 1, 729, 3367, 4081, 1890, 380, 33, 1, 2187, 14197, 23772, 15421, 4550, 644, 42, 1, 6561, 58975, 133057, 116298, 47271, 9702, 1022, 52, 1, 19683, 242461, 724260, 830845, 447195
Offset: 3
Comments
This is the case r = 3 of the r-Stirling numbers of the second kind. The 3-Stirling numbers of the second kind give the number of ways of partitioning the set {1,2,...,n} into k nonempty disjoint subsets with the restriction that the elements 1, 2 and 3 belong to distinct subsets. For remarks on the general case see A143494 (r = 2). The corresponding array of 3-Stirling numbers of the first kind is A143492. The theory of r-Stirling numbers of both kinds is developed in [Broder]. For 3-Lah numbers refer to A143498.
From Peter Bala, Sep 19 2008: (Start)
Let D be the derivative operator d/dx and E the Euler operator x*d/dx. Then x^(-3)*E^n*x^3 = Sum_{k = 0..n} T(n+3,k+3)*x^k*D^k.
The row generating polynomials R_n(x) := Sum_{k= 3..n} T(n,k)*x^k satisfy the recurrence R_(n+1)(x) = x*R_n(x) + x*d/dx(R_n(x)) with R_3(x) = x^3. It follows that the polynomials R_n(x) have only real zeros (apply Corollary 1.2. of [Liu and Wang]).
Relation with the 3-Eulerian numbers E_3(n,j) := A144697(n,j): T(n,k) = (3!/k!)*Sum_{j = n-k..n-3} E_3(n,j)*binomial(j,n-k) for n >= k >= 3.
(End)
T(n,k) = S(n,k,3), n>=k>=3, in Mikhailov's first paper, eq.(28) or (A3). E.g.f. column k from (A20) with k->3, r->k. Therefore, with offset [0,0], this triangle is the Sheffer triangle (exp(3*x),exp(x)-1) with e.g.f. of column no. m>=0: exp(3*x)*((exp(x)-1)^m)/m!. See one of the formulas given below. For Sheffer matrices see the W. Lang link under A006232 with the S. Roman reference, also found in A132393. - Wolfdieter Lang, Sep 29 2011
Examples
Triangle begins n\k|....3....4....5....6....7....8 ================================== 3..|....1 4..|....3....1 5..|....9....7....1 6..|...27...37...12....1 7..|...81..175...97...18....1 8..|..243..781..660..205...25....1 ... T(5,4) = 7. The set {1,2,3,4,5} can be partitioned into four subsets such that 1, 2 and 3 belong to different subsets in 7 ways: {{1}{2}{3}{4,5}}, {{1}{2}{5}{3,4}}, {{1}{2}{4}{3,5}}, {{1}{3}{4}{2,5}}, {{1}{3}{5}{2,4}}, {{2}{3}{4}{1,5}} and {{2}{3}{5}{1,4}}. From _Peter Bala_, Feb 23 2025: (Start) The array factorizes as / 1 \ /1 \ /1 \ /1 \ | 3 1 | | 3 1 ||0 1 ||0 1 | | 9 7 1 | = | 9 4 1 ||0 3 1 ||0 0 1 | ... |27 37 12 1 | |27 13 5 1 ||0 9 4 1 ||0 0 3 1 | |81 175 97 18 1| |81 40 18 6 1||0 27 13 5 1 ||0 0 9 4 1 | |... | |... ||... ||... | where, in the infinite product on the right-hand side, the first array is the Riordan array (1/(1 - 3*x), x/(1 - x)). See A106516. (End)
Links
- Peter Bala, Factorising (r,b)-Stirling arrays
- Andrei Z. Broder, The r-Stirling numbers, Report Number: CS-TR-82-949, Stanford University, Department of Computer Science; see also, Discrete Math. 49, 241-259 (1984).
- A. Dzhumadildaev and D. Yeliussizov, Path decompositions of digraphs and their applications to Weyl algebra, arXiv preprint arXiv:1408.6764v1 [math.CO], 2014. [Version 1 contained many references to the OEIS, which were removed in Version 2. - _N. J. A. Sloane_, Mar 28 2015]
- Askar Dzhumadil’daev and Damir Yeliussizov, Walks, partitions, and normal ordering, Electronic Journal of Combinatorics, 22(4) (2015), #P4.10.
- V. V. Mikhailov, Ordering of some boson operator functions, J. Phys A: Math. Gen. 16 (1983) 3817-3827.
- V. V. Mikhailov, Normal ordering and generalised Stirling numbers, J. Phys A: Math. Gen. 18 (1985) 231-235.
- Erich Neuwirth, Recursively defined combinatorial functions: Extending Galton's board, Discrete Math. 239 No. 1-3, 33-51 (2001).
- L. Liu and Y. Wang, A unified approach to polynomial sequences with only real zeros, arXiv:math/0509207 [math.CO], 2005-2006.
- Michael J. Schlosser and Meesue Yoo, Elliptic Rook and File Numbers, Electronic Journal of Combinatorics, 24(1) (2017), #P1.31.
Crossrefs
Programs
-
Maple
A143495 := (n, k) -> (1/(k-3)!)*add((-1)^(k-i-1)*binomial(k-3,i)*(i+3)^(n-3), i = 0..k-3): for n from 3 to 12 do seq(A143495(n, k), k = 3..n) end do;
-
Mathematica
nmax = 12; t[n_, k_] := 1/(k-3)!* Sum[ (-1)^(k-j-1)*Binomial[k-3, j]*(j+3)^(n-3), {j, 0, k-3}]; Flatten[ Table[ t[n, k], {n, 3, nmax}, {k, 3, n}]] (* Jean-François Alcover, Dec 07 2011, after Maple *)
-
Sage
@CachedFunction def stirling2r(n, k, r) : if n < r: return 0 if n == r: return 1 if k == r else 0 return stirling2r(n-1, k-1, r) + k*stirling2r(n-1, k, r) A143495 = lambda n, k: stirling2r(n, k, 3) for n in (3..8): [A143495(n, k) for k in (3..n)] # Peter Luschny, Nov 19 2012
Formula
T(n+3,k+3) = (1/k!)*Sum_{i = 0..k} (-1)^(k-i)*binomial(k,i)*(i+3)^n, n,k >= 0.
T(n,k) = Stirling2(n,k) - 3*Stirling2(n-1,k) + 2*Stirling2(n-2,k), n,k >= 3.
Recurrence relation: T(n,k) = T(n-1,k-1) + k*T(n-1,k) for n > 3, with boundary conditions: T(n,2) = T(2,n) = 0 for all n; T(3,3) = 1; T(3,k) = 0 for k > 3.
Special cases: T(n,3) = 3^(n-3); T(n,4) = 4^(n-3) - 3^(n-3).
E.g.f. (k+3) column (with offset 3): (1/k!)*exp(3x)*(exp(x)-1)^k.
O.g.f. k-th column: Sum_{n >= k} T(n,k)*x^n = x^k/((1-3*x)*(1-4*x)*...*(1-k*x)).
E.g.f.: exp(3*t + x*(exp(t)-1)) = Sum_{n >= 0} Sum_{k = 0..n} T(n+3,k+3)*x^k*t^n/n! = Sum_{n >= 0} B_n(3;x)*t^n/n! = 1 + (3+x)*t/1! + (9+7*x+x^2)*t^2/2! + ..., where the row polynomials, B_n(3;x) := Sum_{k = 0..n} T(n+3,k+3)*x^k, may be called the 3-Bell polynomials.
Dobinski-type identities: Row polynomial B_n(3;x) = exp(-x)*Sum_{i >= 0} (i+3)^n*x^i/i!; Sum_{k = 0..n} k!*T(n+3,k+3)*x^k = Sum_{i >= 0} (i+3)^n*x^i/(1+x)^(i+1).
The T(n,k) are the connection coefficients between the falling factorials and the shifted monomials (x+3)^(n-3). For example, 9 + 7*x + x*(x-1) = (x+3)^2 and 27 + 37*x + 12x*(x-1) + x*(x-1)*(x-2) = (x+3)^3.
A093560 (3,1) Pascal triangle.
1, 3, 1, 3, 4, 1, 3, 7, 5, 1, 3, 10, 12, 6, 1, 3, 13, 22, 18, 7, 1, 3, 16, 35, 40, 25, 8, 1, 3, 19, 51, 75, 65, 33, 9, 1, 3, 22, 70, 126, 140, 98, 42, 10, 1, 3, 25, 92, 196, 266, 238, 140, 52, 11, 1, 3, 28, 117, 288, 462, 504, 378, 192, 63, 12, 1, 3, 31, 145, 405, 750, 966, 882, 570, 255, 75, 13, 1
Offset: 0
Comments
The array F(3;n,m) gives in the columns m >= 1 the figurate numbers based on A016777, including the pentagonal numbers A000326 (see the W. Lang link).
This is the third member, d=3, in the family of triangles of figurate numbers, called (d,1) Pascal triangles: A007318 (Pascal (d=1), A029653 (d=2).
This is an example of a Riordan triangle (see A053121 for a comment and the 1991 Shapiro et al. reference on the Riordan group) with o.g.f. of column no. m of the type g(x)*(x*f(x))^m with f(0)=1. Therefore the o.g.f. for the row polynomials p(n,x):=Sum_{m=0..n} a(n,m)*x^m is G(z,x)=g(z)/(1-x*z*f(z)). Here: g(x)=(1+2*x)/(1-x), f(x)=1/(1-x), hence G(z,x)=(1+2*z)/(1-(1+x)*z).
The SW-NE diagonals give the Lucas numbers A000032: L(n) = Sum_{k=0..ceiling((n-1)/2)} a(n-1-k,k), n >= 1, with L(0)=2. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs.
Triangle T(n,k), read by rows, given by [3,-2,0,0,0,0,0,0,...] DELTA [1,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Sep 17 2009
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 09 2013
From Wolfdieter Lang, Jan 09 2015: (Start)
The signed lower triangular matrix (-1)^(n-1)*a(n,m) is the inverse of the Riordan matrix A106516; that is Riordan ((1-2*x)/(1+x),x/(1+x)).
See the Peter Bala comment from Dec 23 2014 in A106516 for general Riordan triangles of the type (g(x), x/(1-x)): exp(x)*r(n,x) = d(n,x) with the e.g.f. r(n,x) of row n and the e.g.f. of diagonal n.
Similarly, for general Riordan triangles of the type (g(x), x/(1+x)): exp(x)*r(n,-x) = d(n,x). (End)
The n-th row polynomial is (3 + 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, Mar 02 2018
Binomial(n-2,k)+2*Binomial(n-3,k) is also the number of permutations avoiding both 123 and 132 with k double descents, i.e., positions with w[i]>w[i+1]>w[i+2]. - Lara Pudwell, Dec 19 2018
Examples
Triangle begins 1, 3, 1, 3, 4, 1, 3, 7, 5, 1, 3, 10, 12, 6, 1, 3, 13, 22, 18, 7, 1, 3, 16, 35, 40, 25, 8, 1, 3, 19, 51, 75, 65, 33, 9, 1, 3, 22, 70, 126, 140, 98, 42, 10, 1, 3, 25, 92, 196, 266, 238, 140, 52, 11, 1,
References
- Kurt Hawlitschek, Johann Faulhaber 1580-1635, Veroeffentlichung der Stadtbibliothek Ulm, Band 18, Ulm, Germany, 1995, Ch. 2.1.4. Figurierte Zahlen.
- Ivo Schneider, Johannes Faulhaber 1580-1635, Birkhäuser, Basel, Boston, Berlin, 1993, ch.5, pp. 109-122.
Links
- Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
- Peter Bala, A note on the diagonals of a proper Riordan Array
- M. Bukata, R. Kulwicki, N. Lewandowski, L. Pudwell, J. Roth, and T. Wheeland, Distributions of Statistics over Pattern-Avoiding Permutations, arXiv preprint arXiv:1812.07112 [math.CO], 2018.
- W. Lang, First 10 rows and array of figurate numbers .
Crossrefs
Programs
-
GAP
Concatenation([1],Flat(List([1..11],n->List([0..n],k->Binomial(n,k)+2*Binomial(n-1,k))))); # Muniru A Asiru, Dec 20 2018
-
Haskell
a093560 n k = a093560_tabl !! n !! k a093560_row n = a093560_tabl !! n a093560_tabl = [1] : iterate (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [3, 1] -- Reinhard Zumkeller, Aug 31 2014
-
Python
from math import comb, isqrt def A093560(n): return comb(r:=(m:=isqrt(k:=n+1<<1))-(k<=m*(m+1)),a:=n-comb(r+1,2))*(r+(r-a<<1))//r if n else 1 # Chai Wah Wu, Nov 12 2024
Formula
a(n, m)=F(3;n-m, m) for 0<= m <= n, otherwise 0, with F(3;0, 0)=1, F(3;n, 0)=3 if n>=1 and F(3;n, m):=(3*n+m)*binomial(n+m-1, m-1)/m if m>=1.
G.f. column m (without leading zeros): (1+2*x)/(1-x)^(m+1), m>=0.
Recursion: a(n, m)=0 if m>n, a(0, 0)= 1; a(n, 0)=3 if n>=1; a(n, m)= a(n-1, m) + a(n-1, m-1).
T(n, k) = C(n, k) + 2*C(n-1, k). - Philippe Deléham, Aug 28 2005
Equals M * A007318, where M = an infinite triangular matrix with all 1's in the main diagonal and all 2's in the subdiagonal. - Gary W. Adamson, Dec 01 2007
Sum_{k=0..n} T(n,k) = A151821(n+1). - Philippe Deléham, Sep 17 2009
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(3 + 7*x + 5*x^2/2! + x^3/3!) = 3 + 10*x + 22*x^2/2! + 40*x^3/3! + 65*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ). - Peter Bala, Dec 22 2014
G.f.: (-1-2*x)/(-1+x+x*y). - R. J. Mathar, Aug 11 2015
Extensions
Incorrect connection with A046055 deleted by N. J. A. Sloane, Jul 08 2009
A106517 Convolution of Fibonacci(n-1) and 3^n.
1, 3, 10, 31, 95, 288, 869, 2615, 7858, 23595, 70819, 212512, 637625, 1913019, 5739290, 17218247, 51655351, 154967040, 464902717, 1394710735, 4184136386, 12552415923, 37657258715, 112971793856, 338915410225, 1016746277043
Offset: 0
Links
- G. C. Greubel, Table of n, a(n) for n = 0..1000
- Index entries for linear recurrences with constant coefficients, signature (4,-2,-3).
Crossrefs
Programs
-
Magma
I:=[1,3,10]; [n le 3 select I[n] else 4*Self(n-1) -2*Self(n-2) -3*Self(n-3): n in [1..41]]; // G. C. Greubel, Aug 05 2021
-
Mathematica
LinearRecurrence[{4,-2,-3},{1,3,10},30] (* Harvey P. Dale, Oct 08 2014 *)
-
PARI
a(n) = sum(k=0, n, fibonacci(n-k-1) * 3^k); \\ Michel Marcus, Aug 06 2021
-
Sage
[(2*3^(n+1) - lucas_number2(n+1, 1, -1))/5 for n in (0..40)] # G. C. Greubel, Aug 05 2021
Formula
G.f.: (1-x)/((1-x-x^2)*(1-3*x)).
a(n) = Sum_{k=0..n} Fibonacci(n-k-1) * 3^k.
a(n) = A101220(2, 3, n+1). - Ross La Haye, Jul 25 2005
a(n) = (1/5)*(6*3^n - Lucas(n+1)). - Ralf Stephan, Nov 16 2010
Sum_{k=0..n} a(k) = A094688(n+1). - G. C. Greubel, Aug 05 2021
Comments