A013609 -id:A013609 - cp's OEIS frontend

A001845 Centered octahedral numbers (crystal ball sequence for cubic lattice).

1, 7, 25, 63, 129, 231, 377, 575, 833, 1159, 1561, 2047, 2625, 3303, 4089, 4991, 6017, 7175, 8473, 9919, 11521, 13287, 15225, 17343, 19649, 22151, 24857, 27775, 30913, 34279, 37881, 41727, 45825, 50183, 54809, 59711, 64897, 70375, 76153, 82239

Offset: 0

N. J. A. Sloane

Number of points in simple cubic lattice at most n steps from origin.
If X is an n-set and Y_i (i=1,2,3) mutually disjoint 2-subsets of X then a(n-6) is equal to the number of 6-subsets of X intersecting each Y_i (i=1,2,3). - Milan Janjic, Aug 26 2007
Equals binomial transform of [1, 6, 12, 8, 0, 0, 0, ...] where (1, 6, 12, 8) = row 3 of the Chebyshev triangle A013609. - Gary W. Adamson, Jul 19 2008
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=2,(i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n >= 4, a(n-2) = -coeff(charpoly(A,x),x^(n-3)). - Milan Janjic, Jan 26 2010
a(n) = A005408(n) * A097080(n-1) / 3. - Reinhard Zumkeller, Dec 15 2013
a(n) = D(3,n) where D are the Delannoy numbers (A008288). As such, a(n) gives the number of grid paths from (0,0) to (3,n) using steps that move one unit north, east, or northeast. - David Eppstein, Sep 07 2014
The first comment above can be re-expressed and generalized as follows: a(n) is the number of points in Z^3 that are L1 (Manhattan) distance <= n from any given point. Equivalently, due to a symmetry that is easier to see in the Delannoy numbers array (A008288), as a special case of Dmitry Zaitsev's Dec 10 2015 comment on A008288, a(n) is the number of points in Z^n that are L1 (Manhattan) distance <= 3 from any given point. - Shel Kaphan, Jan 02 2023

L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81.
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).

T. D. Noe, Table of n, a(n) for n = 0..1000
Luciano Ancora, The Square Pyramidal Number and other figurate numbers, ch. 4.
Bela Bajnok, Additive Combinatorics: A Menu of Research Problems, arXiv:1705.07444 [math.NT], May 2017. See Section 2.3.
D. Bump, K. Choi, P. Kurlberg, and J. Vaaler, A local Riemann hypothesis, I pages 16 and 17
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
Milan Janjić, Two Enumerative Functions
Milan Janjić, Hessenberg Matrices and Integer Sequences, J. Int. Seq. 13 (2010) # 10.7.8.
Milan Janjić, On Restricted Ternary Words and Insets, arXiv:1905.04465 [math.CO], 2019.
G. Kreweras, Sur les hiérarchies de segments, Cahiers Bureau Universitaire Recherche Opérationnelle, Cahier 20, Inst. Statistiques, Univ. Paris, 1973.
G. Kreweras, Sur les hiérarchies de segments, Cahiers du Bureau Universitaire de Recherche Opérationnelle, Institut de Statistique, Université de Paris, #20 (1973). (Annotated scanned copy)
T. P. Martin, Shells of atoms, Phys. Reports, 273 (1996), 199-241, eq. (10).
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
R. G. Stanton and D. D. Cowan, Note on a "square" functional equation, SIAM Rev., 12 (1970), 277-279.
Eric Weisstein's World of Mathematics, Haüy Construction
Eric Weisstein's World of Mathematics, Octahedral Number
Index entries for crystal ball sequences
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Sums of 2 consecutive terms give A008412.
(1/12)*t*(2*n^3 - 3*n^2 + n) + 2*n - 1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A005917, A063493, A063494, A063495, A063496.
Partial sums of A005899.
Cf. A001846, A001847, A001848, etc., A014820, A013609.
Cf. also A240876, A002412, A016061, A005899.
The 28 uniform 3D tilings: cab: A299266, A299267; crs: A299268, A299269; fcu: A005901, A005902; fee: A299259, A299265; flu-e: A299272, A299273; fst: A299258, A299264; hal: A299274, A299275; hcp: A007899, A007202; hex: A005897, A005898; kag: A299256, A299262; lta: A008137, A299276; pcu: A005899, A001845; pcu-i: A299277, A299278; reo: A299279, A299280; reo-e:  A299281, A299282; rho: A008137, A299276; sod: A005893, A005894; sve: A299255, A299261; svh: A299283, A299284; svj: A299254, A299260; svk: A010001, A063489; tca: A299285, A299286; tcd: A299287, A299288; tfs: A005899, A001845; tsi: A299289, A299290; ttw: A299257, A299263; ubt: A299291, A299292; bnn: A007899, A007202. See the Proserpio link in A299266 for overview.
Row/column 3 of A008288.

a001845 n = (2 * n + 1) * (2 * n ^ 2 + 2 * n + 3) `div` 3
-- Reinhard Zumkeller, Dec 15 2013

Table[(4 n^3 - 6 n^2 + 8 n - 3)/3, {n, 100}] (* Vladimir Joseph Stephan Orlovsky, Jan 15 2011 *)
LinearRecurrence[{4, -6, 4, -1}, {1, 7, 25, 63}, 40] (* Harvey P. Dale, Jun 05 2013 *)
CoefficientList[Series[(1 + x)^3/(-1 + x)^4, {x, 0, 20}], x] (* Eric W. Weisstein, Sep 27 2017 *)

a(n)=(2*n+1)*(2*n^2+2*n+3)/3 \\ Charles R Greathouse IV, Dec 06 2011

G.f.: (1+x)^3 /(1-x)^4. [conjectured (correctly) by Simon Plouffe in his 1992 dissertation]
a(n) = (2*n+1)*(2*n^2 + 2*n + 3)/3.
First differences of A014820(n). - Alexander Adamchuk, May 23 2006
a(n) = a(n-1) + 4*n^2 + 2, a(0)=1. - Vincenzo Librandi, Mar 27 2011
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4), with a(0)=1, a(1)=7, a(2)=25, a(3)=63. - Harvey P. Dale, Jun 05 2013
a(n) = Sum_{k=0..min(3,n)} 2^k * binomial(3,k) * binomial(n,k). See Bump et al. - Tom Copeland, Sep 05 2014
From Luciano Ancora, Jan 08 2015: (Start)
a(n) = 2 * A000330(n) + A000330(n+1) + A000330(n-1).
a(n) = A005900(n) + A005900(n+1).
a(n) = A005900(n) + A000330(n) + A000330(n+1).
a(n) = A000330(n-1) + A000330(n) + A005900(n+1). (End)
a(n) = A002412(n+1) + A016061(n-1) for n > 0. - Bruce J. Nicholson, Nov 12 2017
E.g.f.: exp(x)*(3 + 18*x + 18*x^2 + 4*x^3)/3. - Stefano Spezia, Mar 14 2024
Sum_{n >= 1} (-1)^(n+1)/(n*a(n-1)*a(n)) = 5/6 - log(2) = (1 - 1/2 + 1/3) - log(2). - Peter Bala, Mar 21 2024

A024023 a(n) = 3^n - 1.

Original entry on oeis.org

0, 2, 8, 26, 80, 242, 728, 2186, 6560, 19682, 59048, 177146, 531440, 1594322, 4782968, 14348906, 43046720, 129140162, 387420488, 1162261466, 3486784400, 10460353202, 31381059608, 94143178826, 282429536480, 847288609442, 2541865828328, 7625597484986, 22876792454960

Offset: 0

N. J. A. Sloane

Number of different directions along lines and hyper-diagonals in an n-dimensional cubic lattice for the attacking queens problem (A036464 in n=2, A068940 in n=3 and A068941 in n=4). The n-dimensional direction vectors have the a(n)+1 Cartesian coordinates (i,j,k,l,...) where i,j,k,l,... = -1, 0, or +1, excluding the zero-vector i=j=k=l=...=0. The corresponding hyper-line count is A003462. - R. J. Mathar, May 01 2006
Total number of sequences of length m=1,...,n with nonzero integer elements satisfying the condition Sum_{k=1..m} |n_k| <= n. See the K. A. Meissner link p. 6 (with a typo: it should be 3^([2a]-1)-1). - Wolfdieter Lang, Jan 21 2008
Let P(A) be the power set of an n-element set A and R be a relation on P(A) such that for all x, y of P(A), xRy if x and y are disjoint and either 0) x is a proper subset of y or y is a proper subset of x, or 1) x is not a subset of y and y is not a subset of x. Then a(n) = |R|. - Ross La Haye, Mar 19 2009
Number of neighbors in Moore's neighborhood in n dimensions. - Dmitry Zaitsev, Nov 30 2015
Number of terms in conjunctive normal form of Boolean expression with n variables. E.g., a(2) = 8: [~x, ~y, x, y, ~x|~y, ~x|y, x|~y, x|y]. - Yuchun Ji, May 12 2023
Number of rays of the Coxeter arrangement of type B_n. Equivalently, number of facets of the n-dimensional type B permutahedron. - Jose Bastidas, Sep 12 2023

			From _Zerinvary Lajos_, Jan 14 2007: (Start)
Ternary......decimal:
0...............0
2...............2
22..............8
222............26
2222...........80
22222.........242
222222........728
2222222......2186
22222222.....6560
222222222...19682
2222222222..59048
etc...........etc.
(End)
Sequence combinatorics: n=3: With length m=1: [1],[2],[3] each with 2 signs, with m=2: [1,1], [1,2], [2,1], each 2^2 = 4 times from choosing signs; m=3: [1,1,1] coming in 2^3 signed versions: 3*2 + 3*4 + 1*8 = 26 = a(3). The order is important, hence the M_0 multinomials A048996 enter as factors.
A027902 gives the 384 divisors of a(24). - _Reinhard Zumkeller_, Mar 11 2010

Mordechai Ben-Ari, Mathematical Logic for Computer Science, Third edition, 173-203.

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Omran Ahmadi and Robert Granger, An efficient deterministic test for Kloosterman sum zeros, Mathematics of Computation, Vol. 83, No. 285 (2014), pp. 347-363; arXiv preprint, arXiv:1104.3882 [math.NT], 2011-2012. See 1st column of Table 2, p. 9.
Feryal Alayont and Evan Henning, Edge Covers of Caterpillars, Cycles with Pendants, and Spider Graphs, J. Int. Seq. (2023) Vol. 26, Art. 23.9.4.
Michael Baake, Franz Gähler, and Uwe Grimm, Examples of Substitution Systems and Their Factors, Journal of Integer Sequences, Vol. 16 (2013), #13.2.14.
R. Samuel Buss, Herbrand's Theorem, University of California, Logic and Computational Complexity pp. 195-209, Lecture Notes in Computer Science, vol 960. Springer.
Jan Draisma, Tyrrell B. McAllister and Benjamin Nill, Lattice width directions and Minkowski's 3^d-theorem, SIAM J. Discrete Math., Vol. 26, No. 3 (2012), pp. 1104-1107; arXiv preprint, arXiv:0901.1375 [math.CO], Jan 10 2009.
Alessandro Farinelli, Herbrand Universe and Herbrand Base.
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.
Krzysztof A. Meissner, Black hole entropy in Loop Quantum Gravity, Classical and Quantum Gravity, Vol. 21, No. 22 (2004), pp. 5245--5251; arXiv preprint, arXiv:gr-qc/0407052, 2004.
Amir Sapir, The Tower of Hanoi with Forbidden Moves, The Computer J. 47 (1) (2004) 20, case three-in-a row, sequence b(n).
Steven Schlicker, Roman Vasquez, and Rachel Wofford, Integer Sequences from Configurations in the Hausdorff Metric Geometry via Edge Covers of Bipartite Graphs, J. Int. Seq. (2023) Vol. 26, Art. 23.6.6.
Amelia Carolina Sparavigna, The groupoids of Mersenne, Fermat, Cullen, Woodall and other Numbers and their representations by means of integer sequences, Politecnico di Torino, Italy (2019), [math.NT].
Amelia Carolina Sparavigna, Some Groupoids and their Representations by Means of Integer Sequences, International Journal of Sciences (2019) Vol. 8, No. 10.
Wikipedia, Herbrand Structure.
Damiano Zanardini, Computational Logic, Slides, UPM European Master in Computational Logic (EMCL) School of Computer Science Technical University of Madrid, 2009-2010.
Index entries for linear recurrences with constant coefficients, signature (4,-3).

Cf. triangle A013609.
Cf. A003462, A007051, A034472, A214369.
Cf. second column of A145901.

a024023 = subtract 1 . a000244  -- Reinhard Zumkeller, Jun 30 2013

[3^n-1: n in [0..35]]; // Vincenzo Librandi, Apr 30 2011

3^Range[0,30]-1 (* Paolo Xausa, Jul 15 2023 *)

a(n)=3^n-1 \\ Charles R Greathouse IV, Sep 24 2015

vector(50, n, sum(k=0, n, 2^k*binomial(n-1, k))-1) \\ Altug Alkan, Oct 04 2015

my(x='x+O('x^100)); concat([0], Vec(2*x/(-1+x)/(-1+3*x))) \\ Altug Alkan, Oct 16 2015

a(n) = A000244(n) - 1.
a(n) = 2*A003462(n). - R. J. Mathar, May 01 2006
A128760(a(n)) > 0. - Reinhard Zumkeller, Mar 25 2007
G.f.: 2*x/((-1+x)*(-1+3*x)) = 1/(-1+x) - 1/(-1+3*x). - R. J. Mathar, Nov 19 2007
a(n) = Sum_{k=1..n} Sum_{m=1..k} binomial(k-1,m-1)*2^m, n >= 1. a(0)=0. From the sequence combinatorics mentioned above. Twice partial sums of powers of 3.
E.g.f.: e^(3*x) - e^x. - Mohammad K. Azarian, Jan 14 2009
a(n) = A024101(n)/A034472(n). - Reinhard Zumkeller, Feb 14 2009
a(n) = 3*a(n-1) + 2 (with a(0)=0). - Vincenzo Librandi, Nov 19 2010
E.g.f.: -E(0) where E(k) = 1 - 3^k/(1 - x/(x - 3^k*(k+1)/E(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Dec 06 2012
a(n) = A227048(n,A020914(n)). - Reinhard Zumkeller, Jun 30 2013
Sum_{n>=1} 1/a(n) = A214369. - Amiram Eldar, Nov 11 2020
a(n) = Sum_{k=1..n} 2^k*binomial(n,k). - Ridouane Oudra, Jun 15 2025
From Peter Bala, Jul 01 2025: (Start)
For n >= 1, a(2*n)/a(n) = A034472(n) and a(3*n)/a(n) = A034513(n).
Modulo differences in offsets, exp( Sum_{n >= 1} a(k*n)/a(n)*x^n/n ) is the o.g.f. of A003462 (k = 2), A006100 (k = 3), A006101 (k = 4), A006102 (k = 5), A022196 (k = 6), A022197 (k = 7), A022198 (k = 8), A022199 (k = 9), A022200 (k = 10), A022201 (k = 11), A022202 (k = 12) and A022203 (k = 13).
The following are all examples of telescoping series:
Sum_{n >= 1} 3^n/(a(n)*a(n+1)) = 1/2^2; Sum_{n >= 1} 3^n/(a(n)*a(n+1)*a(n+2)) = 1/(2*8^2).
In general, for k >= 1, Sum_{n >= 1} 3^n/(a(n)*a(n+1)*...*a(n+k)) = 1/(a(1)*a(2)*...*a(k)*a(k)).
Sum_{n >= 1} 3^n/(a(n)*a(n+2)) = 5/64; Sum_{n >= 1} (-3)^n/(a(n)*a(n+2)) = -3/64.
Sum_{n >= 1} 3^n/(a(n)*a(n+4)) = 703/83200; Sum_{n >= 1} (-3)^n/(a(n)*a(n+4)) = - 417/83200. (End)

A059304 a(n) = 2^n * (2*n)! / (n!)^2.

Original entry on oeis.org

1, 4, 24, 160, 1120, 8064, 59136, 439296, 3294720, 24893440, 189190144, 1444724736, 11076222976, 85201715200, 657270374400, 5082890895360, 39392404439040, 305870434467840, 2378992268083200, 18531097667174400

Offset: 0

Henry Bottomley, Jan 25 2001

Number of lattice paths from (0,0) to (n,n) using steps (0,1), and two kinds of steps (1,0). - Joerg Arndt, Jul 01 2011
The convolution square root of this sequence is A004981. - T. D. Noe, Jun 11 2002
Also main diagonal of array: T(i,1)=2^(i-1), T(1,j)=1, T(i,j) = T(i,j-1) + 2*T(i-1,j). - Benoit Cloitre, Feb 26 2003
The Hankel transform (see A001906 for definition) of this sequence with interpolated zeros(1, 0, 4, 0, 24, 0, 160, 0, 1120, ...) = is A036442: 1, 4, 32, 512, 16384, ... . - Philippe Deléham, Jul 03 2005
The Hankel transform of this sequence gives A103488. - Philippe Deléham, Dec 02 2007
Equals the central column of the triangle A038207. - Zerinvary Lajos, Dec 08 2007
Equals number of permutations whose reverse shares the same recording tableau in the Robinson-Schensted correspondence with n=(k-1)/2 for k odd. - Dang-Son Nguyen, Jul 02 2024
Number of ternary strings of length 2*n that have the same number of 0's as the combined number of 1's and 2's. For example, a(2)=24 since the strings of length 4 are the 6 permutations of 0011, the 12 permutations of 0012, and the 6 permutations of 0022. - Enrique Navarrete, Jul 30 2025

Harry J. Smith, Table of n, a(n) for n = 0..200
Paul Barry and Arnauld Mesinga Mwafise, Classical and Semi-Classical Orthogonal Polynomials Defined by Riordan Arrays, and Their Moment Sequences, Journal of Integer Sequences, Vol. 21 (2018), Article 18.1.5.
Hacène Belbachir and Abdelghani Mehdaoui, Recurrence relation associated with the sums of square binomial coefficients, Quaestiones Mathematicae (2021) Vol. 44, Issue 5, 615-624.
Hacène Belbachir, Abdelghani Mehdaoui, and László Szalay, Diagonal Sums in the Pascal Pyramid, II: Applications, J. Int. Seq., Vol. 22 (2019), Article 19.3.5.
H. J. Brothers, Pascal's Prism: Supplementary Material.
Tucker J. Ervin, Blake Jackson, Jay Lane, Kyungyong Lee, Son Dang Nguyen, Jack O'Donohue and Michael Vaughan, Permutations whose Reverse Shares the Same Recording Tableau in the Robinson-Schensted Correspondence, Séminaire Lotharingien de Combinatoire 86 (2022), Article B86a.
Tony D. Noe, On the Divisibility of Generalized Central Trinomial Coefficients, Journal of Integer Sequences, Vol. 9 (2006), Article 06.2.7.
Michael Z. Spivey and Laura L. Steil, The k-Binomial Transforms and the Hankel Transform, Journal of Integer Sequences, Vol. 9 (2006), Article 06.1.1.

Diagonal of A013609.
Column k=0 of A067001.
Cf. A000079, A000984, A001906, A036442, A038207, A103488.

[2^n*Binomial(2*n,n): n in [0..25]]; // Vincenzo Librandi, Oct 08 2015

seq(binomial(2*n,n)*2^n,n=0..19); # Zerinvary Lajos, Dec 08 2007

Table[2^n Binomial[2n,n],{n,0,30}] (* Harvey P. Dale, Dec 16 2014 *)

{a(n)=if(n<0, 0, 2^n*(2*n)!/n!^2)} /* Michael Somos, Jan 31 2007 */

{ for (n = 0, 200, write("b059304.txt", n, " ", 2^n * (2*n)! / n!^2); ) } \\ Harry J. Smith, Jun 25 2009

/* as lattice paths: same as in A092566 but use */
steps=[[1, 0], [1, 0], [0, 1]]; /* note the double [1, 0] */
/* Joerg Arndt, Jul 01 2011 */

def A059304(n): return pow(2,n)*binomial(2*n,n)
print([A059304(n) for n in range(41)]) # G. C. Greubel, Jan 18 2025

a(n) = 2^n * C(2*n,n).
D-finite with recurrence a(n) = 4*(2-1/n)*a(n-1).
a(n) = A000079(n)*A000984(n)
G.f.: 1/sqrt(1-8*x) - T. D. Noe, Jun 11 2002
E.g.f.: exp(4*x)*BesselI(0, 4*x). - Vladeta Jovovic, Aug 20 2003
a(n) = A038207(n,n). - Joerg Arndt, Jul 01 2011
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - 4*x*(2*k+1)/(4*x*(2*k+1) + (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 24 2013
E.g.f.: E(0)/2, where E(k) = 1 + 1/(1 - 4*x/(4*x + (k+1)^2/(2*k+1)/E(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 01 2013
G.f.: Q(0)/(1+2*sqrt(x)), where Q(k) = 1 + 2*sqrt(x)/(1 - 2*sqrt(x)*(2*k+1)/(2*sqrt(x)*(2*k+1) + (k+1)/Q(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 09 2013
O.g.f.: hypergeom([1/2], [], 8*x). - Peter Luschny, Oct 08 2015
a(n) = Sum_{k = 0..2*n} (-1)^(n+k)*binomial(2*n,k)*binomial(3*n-2*k,n)* binomial(n+k,n). - Peter Bala, Aug 04 2016
a(n) ~ 8^n/sqrt(Pi*n). - Ilya Gutkovskiy, Aug 04 2016
From Amiram Eldar, Jul 21 2020: (Start)
Sum_{n>=0} 1/a(n) = 8/7 + 8*sqrt(7)*arcsin(1/sqrt(8))/49.
Sum_{n>=0} (-1)^n/a(n) = (8/27)*(3 - arcsinh(1/sqrt(8))). (End)
a(n) = Sum_{k = n..2*n} binomial(2*n,k)*binomial(k,n). In general, for m >= 1, Sum_{k = n..m*n} binomial(m*n,k)*binomial(k,n) = 2^((m-1)*n)*binomial(m*n,n). - Peter Bala, Mar 25 2023
Conjecture: a(n) = Sum_{0 <= j, k <= n} binomial(n, j)*binomial(n, k)* binomial(k+j, n). - Peter Bala, Jul 16 2024

A046816 Pascal's tetrahedron: entries in 3-dimensional version of Pascal triangle, or trinomial coefficients.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 2, 1, 2, 1, 1, 3, 3, 3, 6, 3, 1, 3, 3, 1, 1, 4, 4, 6, 12, 6, 4, 12, 12, 4, 1, 4, 6, 4, 1, 1, 5, 5, 10, 20, 10, 10, 30, 30, 10, 5, 20, 30, 20, 5, 1, 5, 10, 10, 5, 1, 1, 6, 6, 15, 30, 15, 20, 60, 60, 20, 15, 60, 90, 60, 15, 6, 30, 60, 60, 30, 6, 1, 6, 15, 20, 15, 6, 1

Offset: 0

Lior Manor

Greatest numbers in each 2D triangle form A022916 (multinomial coefficient n!/([n/3]![(n+1)/3]![(n+2)/3]!).) 2D triangle sums are powers of 3. - Gerald McGarvey, Aug 15 2004
T(n,j,k) is the number of lattice paths from (0,0,0) to (n,j,k) with steps (1,0,0), (1,1,0) and (1,1,1). - Dimitri Boscainos, Aug 16 2015
T(n,j,k) is the number of k-dimensional hyperfaces in an n-dimensional hypercube at an edge distance of j from a given vertex. For example, the number of 2D faces in a 3D cube touching a given vertex is T(3,0,2) = 3, and the number of 3D cube 1D edges at a separation of 1 edge from a given vertex is T(3,1,1) = 6. - Eitan Y. Levine, Jul 22 2023
The sums along vertical lines within each slice (when oriented as in the example) give A027907. See "vertical sums" link. - Eitan Y. Levine, May 17 2023
Numbers of ways to classify n circles black, red, or green; classified first by how many circles there are altogether, then by how many are of each color. - J. Lowell, Nov 06 2024

			The first few slices of the tetrahedron (or pyramid) are:
  1
-----------------
   1
  1 1
-----------------
    1
   2 2
  1 2 1
-----------------
     1 .... Here is the third slice of the pyramid
    3 3
   3 6 3
  1 3 3 1
----------------
...

Marco Costantini: Metodo per elevare qualsiasi trinomio a qualsiasi potenza. Archimede, rivista per gli insegnanti e i cultori di matematiche pure e applicate, anno XXXVIII ottobre-dicembre 1986, pp. 205-209. [Vincenzo Librandi, Jul 19 2009]

Alois P. Heinz, Table of n, a(n) for n = 0..10659
Eitan Y. Levine, Vertical sums
Wikipedia, Pascal's pyramid
Index entries for sequences related to 3D arrays of numbers

Cf. A007318, A022916.
Entry [3, 2] in each slice gives A002378, entry [4, 3] in each slice gives A027480, entry [5, 2] in each slice gives A033488, entry [5, 3] in each slice gives A033487.
See A268240 for this read mod 2.
Cf. A013609 (row sums).

a046816 n = a046816_list !! n
a046816_list = concat $ concat $ iterate ([[1],[1,1]] *) [1]
instance Num a => Num [a] where
   fromInteger k = [fromInteger k]
   (p:ps) + (q:qs) = p + q : ps + qs
   ps + qs         = ps ++ qs
   (p:ps) * qs'@(q:qs) = p * q : ps * qs' + [p] * qs
    *                = []
-- Reinhard Zumkeller, Apr 02 2011

p:= proc(i, j, k) option remember;
      if k<0 or i<0 or i>k or j<0 or j>i then 0
    elif {i, j, k}={0} then 1
    else p(i, j, k-1) +p(i-1, j, k-1) +p(i-1, j-1, k-1)
      fi
    end:
seq(seq(seq(p(i, j, k), j=0..i), i=0..k), k=0..10);
#  Alois P. Heinz, Apr 03 2011

p[i_, j_, k_] := p[i, j, k] = Which[ k<0 || i<0 || i>k || j<0 || j>i, 0, {i, j, k} == {0, 0, 0}, 1, True, p[i, j, k-1] + p[i-1, j, k-1] + p[i-1, j-1, k-1]]; Table[p[i, j, k], {k, 0, 6}, {i, 0, k}, {j, 0, i}] // Flatten (* Jean-François Alcover, Dec 31 2012, translated from Alois P. Heinz's Maple program *)
(* or *)
Flatten[CoefficientList[CoefficientList[CoefficientList[Series[1/(1-x-x*y-x*y*z), {x, 0, 6}], x], y],z]] (* Georg Fischer, May 29 2019 *)

from math import comb as C
def T(n, j, k): return C(n, j) * C(n-j, k)
print([T(n, r-c, c) for n in range(7) for r in range(n+1) for c in range(r+1)]) # Michael S. Branicky, Dec 26 2024

Coefficients of x, y, z in (x+y+z)^n: Let T'(n; i,j,k) := T(n, j,k) where i = n-(j+k). Then T'(n+1; i,j,k) = T'(n; i-1,j,k)+T'(n; i,j-1,k)+T'(n; i,j,k-1), T'(n; i,j,-1) := 0, T'(n; i,j,k) is invariant under permutations of (i,j,k); T'(0, 0, 0)=1.
T'(n; i,j,k) = n!/(i!*j!*k!) and (x+y+z)^n = Sum_{i+j+k=n; 0 <= i,j,k <= n} T'(n; i,j,k)*x^i*y^j*z^k. Hence Sum_{i+j+k=n; 0 <= i,j,k <= n} T'(n; i,j,k) = 3^n. - Gregory Gerard Wojnar, Oct 08 2020
G.f.: 1/(1-x-x*y-x*y*z). - Georg Fischer, May 29 2019
T(n,j,k) = C(n,j) * C(n-j,k), where C(a,b) are the binomial coefficients, elements of A007318. In particular, T(n,j,0) = C(n,j). - Eitan Y. Levine, Jul 22 2023
(-1)^n * Sum_{i=ceiling(n/k)..n} (-1)^i * T(i*k,n-i,i) = k^n, for n,k > 0. - Eitan Y. Levine, Aug 31 2023

A013610 Triangle of coefficients in expansion of (1+3*x)^n.

Original entry on oeis.org

1, 1, 3, 1, 6, 9, 1, 9, 27, 27, 1, 12, 54, 108, 81, 1, 15, 90, 270, 405, 243, 1, 18, 135, 540, 1215, 1458, 729, 1, 21, 189, 945, 2835, 5103, 5103, 2187, 1, 24, 252, 1512, 5670, 13608, 20412, 17496, 6561, 1, 27, 324, 2268, 10206, 30618, 61236, 78732, 59049, 19683

Offset: 0

N. J. A. Sloane

T(n,k) is the number of lattice paths from (0,0) to (n,k) with steps (1,0) and three kinds of steps (1,1). The number of paths with steps (1,0) and s kinds of steps (1,1) corresponds to the expansion of (1+s*x)^n. - Joerg Arndt, Jul 01 2011
Rows of A027465 reversed. - Michael Somos, Feb 14 2002
T(n,k) equals the number of n-length words on {0,1,2,3} having n-k zeros. - Milan Janjic, Jul 24 2015
T(n-1,k-1) is the number of 3-compositions of n with zeros having k positive parts; see Hopkins & Ouvry reference.  - Brian Hopkins, Aug 16 2020

			Triangle begins
  1;
  1,    3;
  1,    6,    9;
  1,    9,   27,   27;
  1,   12,   54,  108,   81;
  1,   15,   90,  270,  405,  243;
  1,   18,  135,  540, 1215, 1458,  729;
  1,   21,  189,  945, 2835, 5103, 5103, 2187;

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
J. Goldman and J. Haglund, Generalized rook polynomials, J. Combin. Theory A91 (2000), 509-530, 1-rook coefficients on the 3xn board (all heights 3) with k rooks
Brian Hopkins and Stéphane Ouvry, Combinatorics of Multicompositions, arXiv:2008.04937 [math.CO], 2020.

Cf. A007318, A013609, A027465, etc.

Diagonals of the triangle: A000244 (k=n), A027471 (k=n-1), A027472 (k=n-2), A036216 (k=n-3), A036217 (k=n-4), A036219 (k=n-5), A036220 (k=n-6), A036221 (k=n-7), A036222 (k=n-8), A036223 (k=n-9), A172362 (k=n-10).

a013610 n k = a013610_tabl !! n !! k
a013610_row n = a013610_tabl !! n
a013610_tabl = iterate (\row ->
   zipWith (+) (map (* 1) (row ++ [0])) (map (* 3) ([0] ++ row))) [1]
-- Reinhard Zumkeller, May 26 2013

[3^k*Binomial(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, May 19 2021

T:= n-> (p-> seq(coeff(p, x, k), k=0..n))((1+3*x)^n):
seq(T(n), n=0..10);  # Alois P. Heinz, Jul 25 2015

t[n_, k_] := Binomial[n, k]*3^(n-k); Table[t[n, n-k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 05 2013 *)
BinomialROW[n_, k_, t_] := Sum[Binomial[n, k]*Binomial[k, j]*(-1)^(k - j)*t^j, {j, 0, k}]; Column[Table[BinomialROW[n, k, 4], {n, 0, 10}, {k, 0, n}], Center] (* Kolosov Petro, Jan 28 2019 *)
T[0, 0] := 1; T[n_, k_]/;0<=k<=n := T[n, k] = 3T[n-1, k-1]+T[n-1, k]; T[n_, k_] := 0; Flatten@Table[T[n, k], {n, 0, 7}, {k, 0, n}] (* Oliver Seipel, Jan 26 2025 *)

{T(n, k) = polcoeff((1 + 3*x)^n, k)}; /* Michael Somos, Feb 14 2002 */

/* same as in A092566 but use */
steps=[[1,0], [1,1], [1,1], [1,1]]; /* note triple [1,1] */
/* Joerg Arndt, Jul 01 2011 */

flatten([[3^k*binomial(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, May 19 2021

G.f.: 1 / (1 - x*(1+3*y)).
Row sums are 4^n. - Joerg Arndt, Jul 01 2011
T(n,k) = 3^k*C(n,k) = Sum_{i=n-k..n} C(i,n-k)*C(n,i)*2^(n-i). - Mircea Merca, Apr 28 2012
From Peter Bala, Dec 22 2014: (Start)
Riordan array ( 1/(1 - x), 3*x/(1 - x) ).
exp(3*x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(3*x)*(1 + 9*x + 27*x^2/2! + 27*x^3/3!) = 1 + 12*x + 90*x^2/2! + 540*x^3/3! + 2835*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), 3*x/(1 - x) ). (End)
T(n,k) = Sum_{j=0..k} (-1)^(k-j) * binomial(n,k) * binomial(k,j) * 4^j. - Kolosov Petro, Jan 28 2019
T(0,0)=1, T(n,k)=3*T(n-1,k-1)+T(n-1,k) for 0<=k<=n, T(n,k)=0 for k<0 or k>n. - Oliver Seipel, Feb 10 2025

A128099 Triangle read by rows: T(n,k) is the number of ways to tile a 3 X n rectangle with k pieces of 2 X 2 tiles and 3n-4k pieces of 1 X 1 tiles (0 <= k <= floor(n/2)).

Original entry on oeis.org

1, 1, 1, 2, 1, 4, 1, 6, 4, 1, 8, 12, 1, 10, 24, 8, 1, 12, 40, 32, 1, 14, 60, 80, 16, 1, 16, 84, 160, 80, 1, 18, 112, 280, 240, 32, 1, 20, 144, 448, 560, 192, 1, 22, 180, 672, 1120, 672, 64, 1, 24, 220, 960, 2016, 1792, 448, 1, 26, 264, 1320, 3360, 4032, 1792, 128, 1, 28

Offset: 0

Emeric Deutsch, Feb 18 2007

Row sums are the Jacobsthal numbers (A001045).
Apparently, T(n,k)/2^n equals the probability P that n will occur as a partial sum in a randomly-generated infinite sequence of 1s and 2s with n compositions (ordered partitions) into (n-2k) 1s and k 2s. Example: T(6,2)=24; P = 3/8 (24/2^6) that 6 will occur as a partial sum in the sequence with 2 (6-2*2) 1s and 2 2s. - Bob Selcoe, Jul 06 2013
From Johannes W. Meijer, Aug 28 2013: (Start)
The antidiagonal sums are A077949 and the backwards antidiagonal sums are A052947.
Moving the terms in each column of this triangle, see the example, upwards to row 0 gives the Pell-Jacobsthal triangle A013609 as a square array. (End)
The numbers in rows of the triangle are along "first layer" skew diagonals pointing top-right in center-justified triangle given in A013609 ((1+2*x)^n) and along (first layer) skew diagonals pointing top-left in center-justified triangle given in A038207 ((2+x)^n), see links. - Zagros Lalo, Jul 31 2018
If s(n) is the row sum at n, then the ratio s(n)/s(n-1) is approximately 2.000..., when n approaches infinity. - Zagros Lalo, Jul 31 2018
It appears that the rows of this array are the coefficients of the Jacobsthal polynomials (see MathWorld link). - Michel Marcus, Jun 15 2019

			Triangle starts:
  1;
  1;
  1,  2;
  1,  4;
  1,  6,  4;
  1,  8, 12;
  1, 10, 24,  8;
  1, 12, 40, 32;

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

G. C. Greubel, Table of n, a(n) for the first 100 rows, flattened
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.
Richard Fors, Independence Complexes of Certain Families of Graphs, Master's thesis in Mathematics at KTH, presented Aug 19 2011.
R. J. Mathar, Tiling n x m rectangles with 1 x 1 and s x s squares arXiv:1609.03964 [math.CO] (2016).
Zagros Lalo, First layer skew diagonals in center-justified triangle of coefficients in expansion of (1 + 2x)^n
Zagros Lalo, First layer skew diagonals in center-justified triangle of coefficients in expansion of (2 + x)^n
Eric Weisstein's World of Mathematics, Jacobsthal Polynomial

Cf. (Triangle sums) A001045, A095977, A077949, A052947, A113726, A052942, A077909.
Cf. (Similar triangles) A008315, A011973, A102541.
Cf. A013609, A038207
Cf. A207538

T := proc(n,k) if k<=n/2 then 2^k*binomial(n-k,k) else 0 fi end: for n from 0 to 16 do seq(T(n,k),k=0..floor(n/2)) od; # yields sequence in triangular form
T := proc(n, k) option remember: if k<0 or k > floor(n/2) then return(0) fi: if k = 0 then return(1) fi: 2*procname(n-2, k-1) + procname(n-1, k): end: seq(seq(T(n, k), k=0..floor(n/2)), n=0..13); # Johannes W. Meijer, Aug 28 2013

Table[2^k*Binomial[n - k, k] , {n,0,25}, {k,0,Floor[n/2]}] // Flatten  (* G. C. Greubel, Dec 28 2016 *)
t[0, 0] = 1; t[n_, k_] := t[n, k] = If[n < 0 || k < 0, 0, t[n - 1, k] + 2 t[n - 2, k - 1]]; Table[t[n, k], {n, 0, 15}, {k, 0, Floor[n/2]}] // Flatten (* Zagros Lalo, Jul 31 2018 *)

T(n, k) = 2^k*binomial(n-k,k) = 2^k*A011973(n,k).
G.f.: 1/(1-z-2*t*z^2).
Sum_{k=0..floor(n/2)} k*T(n,k) = A095977(n-1).
From Johannes W. Meijer, Aug 28 2013: (Start)
T(n, k) = 2*T(n-2, k-1) + T(n-1, k) with T(n, 0) = 1 and T(n, k) = 0 for k < 0 and k > floor(n/2).
T(n, k) = A013609(n-k, k), n >= 0 and 0 <= k <= floor(n/2). (End)

A039810 Matrix square of Stirling2 triangle A008277: 2-levels set partitions of [n] into k first-level subsets.

Original entry on oeis.org

1, 2, 1, 5, 6, 1, 15, 32, 12, 1, 52, 175, 110, 20, 1, 203, 1012, 945, 280, 30, 1, 877, 6230, 8092, 3465, 595, 42, 1, 4140, 40819, 70756, 40992, 10010, 1120, 56, 1, 21147, 283944, 638423, 479976, 156072, 24570, 1932, 72, 1, 115975, 2090424, 5971350, 5660615, 2350950, 487704, 53550, 3120, 90, 1

Offset: 1

Christian G. Bower, Feb 15 1999

This triangle groups certain generalized Stirling numbers of the second kind A000558, A000559, ... They can also be interpreted in terms of trees of height 3 with n leaves and constraints on the order of the root.
From Peter Bala, Jul 19 2014: (Start)
The (n,k)-th entry in this table gives the number of double partitions of the set [n] = {1,2,...,n} into k blocks. To form a double partition of [n] we first write [n] as a disjoint union X_1 U...U X_k of k nonempty subsets (blocks) X_i of [n]. Then each block X_i is further partitioned into sub-blocks to give a double partition. For instance, {1,2,4} U {3,5} is a partition of [5] into 2 blocks and {{1,4},{2}} U {{3},{5}} is a refinement of this partition to a double partition of [5] into 2 blocks (and 4 sub-blocks).
Compare the above interpretation for the (n,k)-th entry of this table with the interpretation of the (n,k)-th entry of A013609 (the square of Pascal's triangle but with the rows read in reverse order) as counting the pairs (X,Y) of subsets of [n] such that |Y| = k and X is contained in Y. (End)
Also the Bell transform of the shifted Bell numbers B(n+1) without column 0. For the definition of the Bell transform see A264428. - Peter Luschny, Jan 28 2016
T(n,k) is the number of partitions of an n-set into colored blocks, such that exactly k colors are used and the colors are introduced in increasing order. T(3,2) = 6: 1a|23b, 13a|2b, 12a|3b, 1a|2a|3b, 1a|2b|3a, 1a|2b|3b. - Alois P. Heinz, Aug 27 2019

			Triangle begins:
      k = 1    2    3    4    5          sum
  n
  1       1                                1
  2       2    1                           3
  3       5    6    1                     12
  4      15   32   12    1                60
  5      52  175  110   20    1          358
Matrix multiplication Stirling2 * Stirling2:
                  1  0  0  0
                  1  1  0  0
                  1  3  1  0
                  1  7  6  1
.
  1  0  0  0      1  0  0  0
  1  1  0  0      2  1  0  0
  1  3  1  0      5  6  1  0
  1  7  6  1     15 32 12  1
From _Peter Bala_, Jul 19 2014: (Start)
T(5,2) = 175: A 5-set can be partitioned into 2 blocks as either a union of a 3-set and a 2-set or as a union of a 4-set and a singleton set.
In the first case there are 10 ways of partitioning a 5-set into a 3-set and a 2-set. Each 3-set can be further partitioned into sub-blocks in Bell(3) = 5 ways and each 2-set can be further partitioned into sub-blocks in Bell(2) = 2 ways. So altogether we obtain 10*5*2 = 100 double partitions of this type.
In the second case, there are 5 ways of partitioning a 5-set into a 4-set and a 1-set. Each 4-set can be further partitioned in Bell(4) = 15 ways and each 1-set can be further partitioned in Bell(1) = 1 way. So altogether we obtain 5*15*1 = 75 double partitions of this type.
Hence, in total, T(5,2) = 100 + 75 = 175. (End)

Tilman Piesk, First 100 rows, flattened
A. Aboud, J.-P. Bultel, A. Chouria, J.-G. Luque, and O. Mallet, Bell polynomials in combinatorial Hopf algebras, arXiv preprint arXiv:1402.2960 [math.CO], 2014-2015.
T. Mansour, A. Munagi, and M. Shattuck, Recurrence Relations and Two-Dimensional Set Partitions , J. Int. Seq. 14 (2011) # 11.4.1.

Cf. A039811, A039814, A039813 (other products of Stirling matrices).
T(n, 1) = A000110(n) (first column) (Bell numbers).
T(n, 2) = A000558(n) 2-levels set partitions with 2 first-level classes.
T(n, n-1) = A002378(n-1) = n*(n-1) = 2*C(n,2) = set-partitions into (n-2) singletons and one of the two possible set partitions of [2].
Sum is A000258(n), 2-levels set partitions.
Another version with offset 0: A130191.
Horizontal mirror triangle is A046817.
T(2n,n) gives A321712.

# The function BellMatrix is defined in A264428.
# Adds (1,0,0,0, ..) as column 0.
BellMatrix(n -> combinat:-bell(n+1), 10); # Peter Luschny, Jan 28 2016

Flatten[Table[Sum[StirlingS2[n,i]*StirlingS2[i,k],{i,k,n}],{n,1,10},{k,1,n}]] (* Indranil Ghosh, Feb 22 2017 *)
rows = 10;
t = Table[BellB[n+1], {n, 0, rows}];
T[n_, k_] := BellY[n, k, t];
Table[T[n, k], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 22 2018, after Peter Luschny *)

T(n, k) = sum(j=0, n, stirling(n, j, 2)*stirling(j, k, 2)); \\ Seiichi Manyama, Feb 13 2022

S2 = A008277 (Stirling numbers of the second kind).
T = (S2)^2.
T(n,k) = Sum_{i=k..n} S2(n,i) * S2(i,k).
E.g.f. of k-th column: (exp(exp(x)-1)-1)^k/k!. [corrected by Seiichi Manyama, Feb 12 2022]
From Peter Bala, Jul 19 2014: (Start)
T(n,k) = Sum_{disjoint unions X_1 U...U X_k = [n]} Bell(|X_1|)*...*Bell(|X_k|), where Bell(n) = A000110(n).
Recurrence equation: T(n+1,k+1) = Sum_{j = k..n} Bell(n+1-j)*binomial(n,j)* T(j,k).
Row sums [1,3,12,60,358,...] = A000258. (End)

Definition and interpretation edited by Olivier Gérard, Jul 31 2011

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

Original entry on oeis.org

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

Offset: 1

Clark Kimberling, Feb 18 2012

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

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

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

Jean-Luc Baril and José Luis Ramírez, Descent distribution on Catalan words avoiding ordered pairs of Relations, arXiv:2302.12741 [math.CO], 2023.
S. Halici, On some Pell polynomials , Acta Universitatis Apulensis, No. 29/2012, pp. 105-112.
Zagros Lalo, First layer skew diagonals in center-justified triangle of coefficients in expansion of (1 + 2x)^n
Zagros Lalo, First layer skew diagonals in center-justified triangle of coefficients in expansion of (2 + x)^n

Cf. A028297, A207537, A133156, A038207, A099089.
Cf. A053117, A097610, A091894.
Cf. A013609, A038207.
Cf. A128099.

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

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

A001846 Centered 4-dimensional orthoplex numbers (crystal ball sequence for 4-dimensional cubic lattice).

Original entry on oeis.org

1, 9, 41, 129, 321, 681, 1289, 2241, 3649, 5641, 8361, 11969, 16641, 22569, 29961, 39041, 50049, 63241, 78889, 97281, 118721, 143529, 172041, 204609, 241601, 283401, 330409, 383041, 441729, 506921, 579081, 658689, 746241, 842249, 947241, 1061761, 1186369

Offset: 0

N. J. A. Sloane

a(n) is the number of points in the Z^4 lattice that are at distance at most n from the origin in the adjacency graph. - N. J. A. Sloane, Feb 19 2013
Number of nodes of degree 8 in virtual, optimal, chordal graphs of diameter d(G)=n. - S. Bujnowski & B. Dubalski (slawb(AT)atr.bydgoszcz.pl), Mar 07 2002
If Y_i (i=1,2,3,4) are 2-blocks of an (n+4)-set X then a(n-4) is the number of 8-subsets of X intersecting each Y_i (i=1,2,3,4). - Milan Janjic, Oct 28 2007
Equals binomial transform of [1, 8, 24, 32, 16, 0, 0, 0, ...] where (1, 8, 24, 32, 16) = row 4 of the Chebyshev triangle A013609. - Gary W. Adamson, Jul 19 2008
Comment from Ben Thurston, Feb 18 2013: In the plane, if you make a picture by taking one unit step in each of the basic 8 directions from a central dot, then from each of those going one unit step in each of the eight directions, ... (see illustration), it appears that the number of dots in the picture after n steps is equal to a(n). Response from N. J. A. Sloane, Feb 19 2013: This is correct, and follows from the fact that the Z-module Z[1,i,(+-1+i)/sqrt(2)] is essentially a copy of the Z^4 lattice.
a(n) = D(4,n) where D are the Delannoy numbers (A008288). As such, a(n) gives the number of grid paths from (0,0) to (4,n) using steps that move one unit north, east, or northeast. - Jack W Grahl, Feb 15 2021
The first comment above can be re-expressed and generalized as follows: a(n) is the number of points in Z^4 that are L1 (Manhattan) distance <= n from any given point. Equivalently, due to a symmetry that is easier to see in the Delannoy numbers array (A008288), as a special case of Dmitry Zaitsev's Dec 10 2015 comment on A008288, a(n) is the number of points in Z^n that are L1 (Manhattan) distance <= 4 from any given point. - Shel Kaphan, Jan 02 2023

			a(6)=1289: (2*6^4 + 4*6^3 + 10*6^2 + 8*6 + 3) / 3 = (2592 + 864 + 360 + 48 + 3) / 3 = 3867 / 3 = 1289.

L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81.
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).

T. D. Noe, Table of n, a(n) for n = 0..1000
D. Bump, K. Choi, P. Kurlberg, and J. Vaaler, A local Riemann hypothesis, I pages 16 and 17
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
Milan Janjic, Two Enumerative Functions. [Broken link; WayBackMachine archive]
G. Kreweras, Sur les hiérarchies de segments, Cahiers du Bureau Universitaire de Recherche Opérationnelle, Institut de Statistique, Université de Paris, #20 (1973).
G. Kreweras, Sur les hiérarchies de segments, Cahiers du Bureau Universitaire de Recherche Opérationnelle, Institut de Statistique, Université de Paris, #20 (1973). (Annotated scanned copy)
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
R. G. Stanton and D. D. Cowan, Note on a "square" functional equation, SIAM Rev., 12 (1970), 277-279.
Ben Thurston, Illustration of the first four clusters of points in two dimensions
Index entries for crystal ball sequences
Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5,1).

Cf. A005408, A001844, A001845, A001847, A059100, A013609.
First differences are A008412.
Cf. A240876.
Row/column 4 of A008288.

for n from 1 to k do eval((2*n^4+4*n^3+10*n^2+8*n+3)/3) od;
A001846:=-(z+1)**4/(z-1)**5; # conjectured (correctly) by Simon Plouffe in his 1992 dissertation

CoefficientList[Series[(-z^4-4 z^3-6 z^2-4 z-1)/(z-1)^5, {z, 0, 200}], z] (* Vladimir Joseph Stephan Orlovsky, Jun 19 2011 *)
Table[(((2 n + 4) n + 10) n + 8) n/3 + 1, {n, 0, 30}] (* Robert A. Russell, Jul 02 2025 *)

G.f.: (1+x)^4 /(1-x)^5.
a(n) = (2*n^4 + 4*n^3 + 10*n^2 + 8*n + 3)/3. - S. Bujnowski & B. Dubalski (slawb(AT)atr.bydgoszcz.pl), Mar 07 2002
From Jonathan Vos Post, Mar 15 2006: (Start)
a(n) = Sum_{i=0..n} A008412(i);
a(n) = Sum_{i=0..n} 8*i*(i^2 + 2)/3;
a(n) = Sum_{i=0..n} 8*i*(A059100(i))/3. (End)
a(n) = Sum_{k=0..min(4,n)} 2^k * binomial(4,k)* binomial(n,k). See Bump et al. - Tom Copeland, Sep 05 2014
E.g.f.: exp(x)*(3 + 24*x + 36*x^2 + 16*x^3 + 2*x^4)/3. - Stefano Spezia, Mar 14 2024
Sum_{n >= 1} (-1)^(n+1)/(n*a(n-1)*a(n)) = log(2) - 7/12 = log(2) - (1 - 1/2 + 1/3 - 1/4). - Peter Bala, Mar 23 2024

A001847 Crystal ball sequence for 5-dimensional cubic lattice.

Original entry on oeis.org

1, 11, 61, 231, 681, 1683, 3653, 7183, 13073, 22363, 36365, 56695, 85305, 124515, 177045, 246047, 335137, 448427, 590557, 766727, 982729, 1244979, 1560549, 1937199, 2383409, 2908411, 3522221, 4235671, 5060441, 6009091, 7095093, 8332863, 9737793, 11326283

Offset: 0

N. J. A. Sloane

Number of nodes degree 10 in virtual, optimal chordal graphs of diameter d(G)=n - S. Bujnowski & B. Dubalski (slawb(AT)atr.bydgoszcz.pl), Mar 07 2002
If Y_i (i=1,2,3,4,5) are 2-blocks of a (n+5)-set X then a(n-5) is the number of 10-subsets of X intersecting each Y_i (i=1,2,3,4,5). - Milan Janjic, Oct 28 2007
Equals binomial transform of [1, 10, 40, 80, 80, 32, 0, 0, 0, ...] where (1, 10, 40, 80, 80, 32) = row 5 of the Chebyshev triangle A013609. - Gary W. Adamson, Jul 19 2008
a(n) is the number of points in Z^5 that are L1 (Manhattan) distance <= n from any given point. Equivalently, due to a symmetry that is easier to see in the Delannoy numbers array (A008288), as a special case of Dmitry Zaitsev's Dec 10 2015 comment on A008288, a(n) is the number of points in Z^n that are L1 (Manhattan) distance <= 5 from any given point. - Shel Kaphan, Jan 02 2023

			a(5)=1683, (4*5^5 + 10*5^4 + 40*5^3 + 50*5^2 + 46*5 + 15)/15 = (12500 + 6250 + 5000 + 230 + 15)/15 = 25245/15 = 1683.

L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81.
E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 231.
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).

T. D. Noe, Table of n, a(n) for n = 0..1000
D. Bump, K. Choi, P. Kurlberg, and J. Vaaler, A local Riemann hypothesis, I pages 16 and 17
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
Milan Janjic, Two Enumerative Functions. [Broken link; WayBackMachine archive]
G. Kreweras, Sur les hiérarchies de segments, Cahiers du Bureau Universitaire de Recherche Opérationnelle, Institut de Statistique, Université de Paris, #20 (1973).
G. Kreweras, Sur les hiérarchies de segments, Cahiers du Bureau Universitaire de Recherche Opérationnelle, Institut de Statistique, Université de Paris, #20 (1973). (Annotated scanned copy)
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
R. G. Stanton and D. D. Cowan, Note on a "square" functional equation, SIAM Rev., 12 (1970), 277-279.
Index entries for crystal ball sequences
Index entries for linear recurrences with constant coefficients, signature (6,-15,20,-15,6,-1).

Cf. A005408, A001844, A001845, A001846, A013609.
Cf. A240876.
Row/column 5 of A008288.

for n from 1 to k do eval((4*n^5+10*n^4+40*n^3+50*n^2+46*n+15)/15) od;
A001847:=(z+1)**5/(z-1)**6; # conjectured by Simon Plouffe in his 1992 dissertation

CoefficientList[Series[(z+1)^5/(z-1)^6, {z, 0, 200}], z] (* Vladimir Joseph Stephan Orlovsky, Jun 19 2011 *)
Table[((((4 n + 10) n + 40) n + 50) n + 46) n/15 + 1, {n, 0, 30}] (* Robert A. Russell, Jul 02 2025 *)

G.f.: (1+x)^5 /(1-x)^6.
a(n) = (4*n^5+10*n^4+40*n^3+50*n^2+46*n+15)/15. - S. Bujnowski & B. Dubalski (slawb(AT)atr.bydgoszcz.pl), Mar 07 2002
a(n) = Sum_{k=0..min(5,n)} 2^k * binomial(5,k)* binomial(n,k). See Bump et al. - Tom Copeland, Sep 05 2014
E.g.f.: exp(x)*(15 + 150*x + 300*x^2 + 200*x^3 + 50*x^4 + 4*x^5)/15. - Stefano Spezia, Mar 17 2024
Sum_{n >= 1} (-1)^(n+1)/(n*a(n-1)*a(n)) = 47/60 - log(2) = (1 - 1/2 + 1/3 - 1/4 + 1/5) - log(2). - Peter Bala, Mar 23 2024

cp's OEIS Frontend

A001845 Centered octahedral numbers (crystal ball sequence for cubic lattice).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

A024023 a(n) = 3^n - 1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

PARI

PARI

Formula

A059304 a(n) = 2^n * (2*n)! / (n!)^2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

PARI

SageMath

Formula

A046816 Pascal's tetrahedron: entries in 3-dimensional version of Pascal triangle, or trinomial coefficients.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Python

Formula

A013610 Triangle of coefficients in expansion of (1+3*x)^n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

A128099 Triangle read by rows: T(n,k) is the number of ways to tile a 3 X n rectangle with k pieces of 2 X 2 tiles and 3n-4k pieces of 1 X 1 tiles (0 <= k <= floor(n/2)).