A119800 -id:A119800 - cp's OEIS frontend

A005899 Number of points on surface of octahedron; also coordination sequence for cubic lattice: a(0) = 1; for n > 0, a(n) = 4n^2 + 2.

1, 6, 18, 38, 66, 102, 146, 198, 258, 326, 402, 486, 578, 678, 786, 902, 1026, 1158, 1298, 1446, 1602, 1766, 1938, 2118, 2306, 2502, 2706, 2918, 3138, 3366, 3602, 3846, 4098, 4358, 4626, 4902, 5186, 5478, 5778, 6086, 6402, 6726, 7058, 7398, 7746, 8102, 8466

Offset: 0

N. J. A. Sloane

Also, the number of regions the plane can be cut into by two overlapping concave (2n)-gons. - Joshua Zucker, Nov 05 2002
If X is an n-set and Y_i (i=1,2,3) are mutually disjoint 2-subsets of X then a(n-5) is equal to the number of 5-subsets of X intersecting each Y_i (i=1,2,3). - Milan Janjic, Aug 26 2007
Binomial transform of a(n) is A055580(n). - Wesley Ivan Hurt, Apr 15 2014
The identity (4*n^2+2)^2 - (n^2+1)*(4*n)^2 = 4 can be written as a(n)^2 - A002522(n)*A008586(n)^2 = 4. - Vincenzo Librandi, Jun 15 2014
Also the least number of unit cubes required, at the n-th iteration, to surround a 3D solid built from unit cubes, in order to hide all its visible faces, starting with a unit cube. - R. J. Cano, Sep 29 2015
Also, coordination sequence for "tfs" 3D uniform tiling. - N. J. A. Sloane, Feb 10 2018
Also, the number of n-th order specular reflections arriving at a receiver point from an emitter point inside a cuboid with reflective faces. - Michael Schutte, Sep 18 2018

H. S. M. Coxeter, "Polyhedral numbers," in R. S. Cohen et al., editors, For Dirk Struik. Reidel, Dordrecht, 1974, pp. 25-35.
Gmelin Handbook of Inorg. and Organomet. Chem., 8th Ed., 1994, TYPIX search code (225) cF8
B. Grünbaum, Uniform tilings of 3-space, Geombinatorics, 4 (1994), 49-56. See tilings #16 and #22.
R. W. Marks and R. B. Fuller, The Dymaxion World of Buckminster Fuller. Anchor, NY, 1973, p. 46.
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
Barry Balof, Restricted tilings and bijections, J. Integer Seq. 15 (2012), no. 2, Article 12.2.3, 17 pp.
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
Pierre de la Harpe, On the prehistory of growth of groups, arXiv:2106.02499 [math.GR], 2021.
R. W. Grosse-Kunstleve, Coordination Sequences and Encyclopedia of Integer Sequences
R. W. Grosse-Kunstleve, G. O. Brunner and N. J. A. Sloane, Algebraic Description of Coordination Sequences and Exact Topological Densities for Zeolites, Acta Cryst., A52 (1996), pp. 879-889.
Milan Janjic, Two Enumerative Functions
Milan Janjić, On Restricted Ternary Words and Insets, arXiv:1905.04465 [math.CO], 2019.
M. O'Keeffe, Coordination sequences for lattices, Zeit. f. Krist., 210 (1995), 905-908.
M. O'Keeffe, Coordination sequences for lattices, Zeit. f. Krist., 210 (1995), 905-908. [Annotated scanned copy]
Carlos I. Perez-Sanchez, The Spectral Action on quivers, arXiv:2401.03705 [math.RT], 2024.
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
Reticular Chemistry Structure Resource (RCSR), The pcu tiling (or net)
Reticular Chemistry Structure Resource (RCSR), The tfs tiling (or net)
B. K. Teo and N. J. A. Sloane, Magic numbers in polygonal and polyhedral clusters, Inorgan. Chem. 24 (1985),4545-4558.
N. J. A. Sloane, Illustration of a(0)=1, a(1)=6, a(2)=18 (from Teo-Sloane 1985)
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Partial sums give A001845.
Column 2 * 2 of array A188645.
Cf. A001844, A002522, A008586, A010014, A055580, A195322, A206399.
Cf. A000384, A014105, A000217, A008574, A008412, A035597.
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 3 of A035607, A266213, A343599.
Column 3 of A113413, A119800, A122542.

[4*n^2 + 2 : n in [0..50]]; // Wesley Ivan Hurt, Oct 26 2015

A005899:=n->4*n^2 + 2; seq(A005899(n), n=0..50); # Wesley Ivan Hurt, Apr 15 2014

Join[{1},4Range[40]^2+2] (* or *) Join[{1},LinearRecurrence[{3,-3,1},{6,18,38},40]] (* Harvey P. Dale, Nov 08 2011 *)

Vec(((1+x)/(1-x))^3 + O(x^100)) \\ Altug Alkan, Oct 26 2015

G.f.: ((1+x)/(1-x))^3. - Simon Plouffe in his 1992 dissertation
Binomial transform of [1, 5, 7, 1, -1, 1, -1, 1, ...]. - Gary W. Adamson, Nov 02 2007
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), with a(0)=1, a(1)=6, a(2)=18, a(3)=38. - Harvey P. Dale, Nov 08 2011
Recurrence: n*a(n) = (n-2)*a(n-2) + 6*a(n-1), a(0)=1, a(1)=6. - Fung Lam, Apr 15 2014
For n > 0, a(n) = A001844(n-1) + A001844(n) = (n-1)^2 + 2n^2 + (n+1)^2. - Doug Bell, Aug 18 2015
For n > 0, a(n) = A010014(n) - A195322(n). - R. J. Cano, Sep 29 2015
For n > 0, a(n) = A000384(n+1) + A014105(n-1). - Bruce J. Nicholson, Oct 08 2017
a(n) = A008574(n) + A008574(n-1) + a(n-1). - Bruce J. Nicholson, Dec 18 2017
a(n) = 2*d*Hypergeometric2F1(1-d, 1-n, 2, 2) where d=3, n>0. - Shel Kaphan, Feb 16 2023
a(n) = A035597(n)*3/n, for n>0. - Shel Kaphan, Feb 26 2023
E.g.f.: exp(x)*(2 + 4*x + 4*x^2) - 1. - Stefano Spezia, Mar 08 2023
Sum_{n>=0} 1/a(n) = 3/4 + Pi *sqrt(2)*coth( Pi/sqrt 2)/8 = 1.31858... - R. J. Mathar, Apr 27 2024

A003506 Triangle of denominators in Leibniz's Harmonic Triangle a(n,k), n >= 1, 1 <= k <= n.

Original entry on oeis.org

1, 2, 2, 3, 6, 3, 4, 12, 12, 4, 5, 20, 30, 20, 5, 6, 30, 60, 60, 30, 6, 7, 42, 105, 140, 105, 42, 7, 8, 56, 168, 280, 280, 168, 56, 8, 9, 72, 252, 504, 630, 504, 252, 72, 9, 10, 90, 360, 840, 1260, 1260, 840, 360, 90, 10, 11, 110, 495, 1320, 2310, 2772, 2310, 1320, 495, 110, 11

Offset: 1

N. J. A. Sloane

Array 1/Beta(n,m) read by antidiagonals. - Michael Somos, Feb 05 2004
a(n,3) = A027480(n-2); a(n,4) = A033488(n-3). - Ross La Haye, Feb 13 2004
a(n,k) = total size of all of the elements of the family of k-size subsets of an n-element set. For example, a 2-element set, say, {1,2}, has 3 families of k-size subsets: one with 1 0-size element, one with 2 1-size elements and one with 1 2-size element; respectively, {{}}, {{1},{2}}, {{1,2}}. - Ross La Haye, Dec 31 2006
Second slice along the 1-2-plane in the cube a(m,n,o) = a(m-1,n,o) + a(m,n-1,o) + a(m,n,o-1) with a(1,0,0)=1 and a(m<>1=0,n>=0,0>=o)=0, for which the first slice is Pascal's triangle (slice read by antidiagonals). - Thomas Wieder, Aug 06 2006
Triangle, read by rows, given by [2,-1/2,1/2,0,0,0,0,0,0,...] DELTA [2,-1/2,1/2,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Oct 07 2007
This sequence * [1/1, 1/2, 1/3, ...] = (1, 3, 7, 15, 31, ...). - Gary W. Adamson, Nov 14 2007
n-th row = coefficients of first derivative of corresponding Pascal's triangle row. Example: x^4 + 4x^3 + 6x^2 + 4x + 1 becomes (4, 12, 12, 4). - Gary W. Adamson, Dec 27 2007
From Paul Curtz, Jun 03 2011: (Start)
Consider
   1     1/2   1/3    1/4    1/5
  -1/2  -1/6  -1/12  -1/20  -1/30
   1/3   1/12  1/30   1/60   1/105
  -1/4  -1/20 -1/60  -1/140 -1/280
   1/5   1/30  1/105  1/280  1/630
This is an autosequence (the inverse binomial transform is the sequence signed) of the second kind: the main diagonal is 2 times the first upper diagonal.
Note that 2, 12, 60, ... = A005430(n+1), Apery numbers = 2*A002457(n). (End)
From Louis Conover (for the 9th grade G1c mathematics class at the Chengdu Confucius International School), Mar 02 2015: (Start)
The i-th order differences of n^-1 appear in the (i+1)th row.
  1,    1/2,   1/3,   1/4,    1/5,    1/6,    1/7,     1/8, ...
  1/2,  1/6,  1/12,  1/20,   1/30,   1/42,   1/56,    1/72, ...
  1/3, 1/12,  1/30,  1/60,  1/105,  1/168,  1/252,   1/360, ...
  1/4, 1/20,  1/60, 1/140,  1/280,  1/504,  1/840,  1/1320, ...
  1/5, 1/30, 1/105, 1/280,  1/630, 1/1260, 1/2310,  1/3960, ...
  1/6, 1/42, 1/168, 1/504, 1/1260, 1/2772, 1/5544, 1/12012, ...
(End)
T(n,k) is the number of edges of distance k from a fixed vertex in the n-dimensional hypercube. - Simon Burton, Nov 04 2022

			The triangle begins:
  1;
  1/2, 1/2;
  1/3, 1/6, 1/3;
  1/4, 1/12, 1/12, 1/4;
  1/5, 1/20, 1/30, 1/20, 1/5;
  ...
The triangle of denominators begins:
   1
   2   2
   3   6   3
   4  12  12    4
   5  20  30   20    5
   6  30  60   60   30    6
   7  42 105  140  105   42    7
   8  56 168  280  280  168   56    8
   9  72 252  504  630  504  252   72   9
  10  90 360  840 1260 1260  840  360  90  10
  11 110 495 1320 2310 2772 2310 1320 495 110 11

A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, see 130.
B. A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8. English translation published by Fibonacci Association, Santa Clara Univ., Santa Clara, CA, 1993; see p. 38.
G. Boole, A Treatise On The Calculus of Finite Differences, Dover, 1960, p. 26.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 83, Problem 25.
M. Elkadi and B. Mourrain, Symbolic-numeric methods for solving polynomial equations and applications, Chap 3. of A. Dickenstein and I. Z. Emiris, eds., Solving Polynomial Equations, Springer, 2005, pp. 126-168. See p. 152.
D. Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, 35.

Reinhard Zumkeller, Rows n = 1..120 of triangle, flattened
F. S. Al-Kharousi, A. Umar, and M. M. Zubairu, On injective partial Catalan monoids, arXiv:2501.00285 [math.GR], 2024. See pp. 8-9.
Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids, English translation published by Fibonacci Association, Santa Clara Univ., Santa Clara, CA, 1993.
H. J. Brothers, Pascal's Prism: Supplementary Material.
D. Dumont, Matrices d'Euler-Seidel, Sem. Loth. Comb. B05c (1981) 59-78.
Yilmaz Simsek, Construction of a generalization of the Leibnitz numbers and their properties, arXiv:2011.13701 [math.NT], 2020.
Eric Weisstein's World of Mathematics, Leibniz Harmonic Triangle

Cf. A007622, A128064, A164555/A027642, A356546.
Cf. A094305, A121547, A121306, A119800, A002378, A007318.
Row sums are in A001787. Central column is A002457. Half-diagonal is in A090816. A116071, A215652.
Denominators of i-th order differences of n^-1 are given in: (1st) A002378, (2nd) A027480, (3rd) A033488, (4th) A174002, (5th) A253946. - Louis Conover, Mar 02 2015
Columns k >= 1 (offset 1): A000027, A002378, A027480, A033488, A174002, A253946(n+4), ..., with sum of reciprocals: infinity, 1, 1/2, 1/3, 1/4, 1/5, ..., respectively. - Wolfdieter Lang, Jul 20 2022

a003506 n k = a003506_tabl !! (n-1) !! (n-1)
a003506_row n = a003506_tabl !! (n-1)
a003506_tabl = scanl1 (\xs ys ->
   zipWith (+) (zipWith (+) ([0] ++ xs) (xs ++ [0])) ys) a007318_tabl
a003506_list = concat a003506_tabl
-- Reinhard Zumkeller, Nov 14 2013, Nov 17 2011

with(combstruct):for n from 0 to 11 do seq(m*count(Combination(n), size=m), m = 1 .. n) od; # Zerinvary Lajos, Apr 09 2008
A003506 := (n,k) -> k*binomial(n,k):
seq(print(seq(A003506(n,k),k=1..n)),n=1..7); # Peter Luschny, May 27 2011

L[n_, 1] := 1/n; L[n_, m_] := L[n, m] = L[n - 1, m - 1] - L[n, m - 1]; Take[ Flatten[ Table[ 1 / L[n, m], {n, 1, 12}, {m, 1, n}]], 66]
t[n_, m_] = Gamma[n]/(Gamma[n - m]*Gamma[m]); Table[Table[t[n, m], {m, 1, n - 1}], {n, 2, 12}]; Flatten[%] (* Roger L. Bagula and Gary W. Adamson, Sep 14 2008 *)
Table[k*Binomial[n,k],{n,1,7},{k,1,n}] (* Peter Luschny, May 27 2011 *)
t[n_, k_] := Denominator[n!*k!/(n+k+1)!]; Table[t[n-k, k] , {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Nov 28 2013 *)

A(i,j)=if(i<1||j<1,0,1/subst(intformal(x^(i-1)*(1-x)^(j-1)),x,1))

A(i,j)=if(i<1||j<1,0,1/sum(k=0,i-1,(-1)^k*binomial(i-1,k)/(j+k)))

{T(n, k) = (n + 1 - k) * binomial( n, k - 1)} /* Michael Somos, Feb 06 2011 */

T_row = lambda n: (n*(x+1)^(n-1)).list()
for n in (1..10): print(T_row(n)) # Peter Luschny, Feb 04 2017
# Assuming offset 0:
def A003506(n, k):
    return falling_factorial(n+1,n)//(factorial(k)*factorial(n-k))
for n in range(9): print([A003506(n, k) for k in range(n+1)]) # Peter Luschny, Aug 13 2022

a(n, 1) = 1/n; a(n, k) = a(n-1, k-1) - a(n, k-1) for k > 1.
Considering the integer values (rather than unit fractions): a(n, k) = k*C(n, k) = n*C(n-1, k-1) = a(n, k-1)*a(n-1, k-1)/(a(n, k-1) - a(n-1, k-1)) = a(n-1, k) + a(n-1, k-1)*k/(k-1) = (a(n-1, k) + a(n-1, k-1))*n/(n-1) = k*A007318(n, k) = n*A007318(n-1, k-1). Row sums of integers are n*2^(n-1) = A001787(n); row sums of the unit fractions are A003149(n-1)/A000142(n). - Henry Bottomley, Jul 22 2002
From Vladeta Jovovic, Nov 01 2003: (Start)
G.f.: x*y/(1-x-y*x)^2.
E.g.f.: x*y*exp(x+x*y). (End)
T(n,k) = n*binomial(n-1,k-1) = n*A007318(n-1,k-1). - Philippe Deléham, Aug 04 2006
Binomial transform of A128064(unsigned). - Gary W. Adamson, Aug 29 2007
From Roger L. Bagula and Gary W. Adamson, Sep 14 2008: (Start)
t(n,m) = Gamma(n)/(Gamma(n - m)*Gamma(m)).
f(s,n) = Integral_{x=0..oo} exp(-s*x)*x^n dx = Gamma(n)/s^n; t(n,m) = f(s,n)/(f(s,n-m)*f(s,m)) = Gamma(n)/(Gamma(n - m)*Gamma(m)); the powers of s cancel out. (End)
From Reinhard Zumkeller, Mar 05 2010: (Start)
T(n,5) = T(n,n-4) = A174002(n-4) for n > 4.
T(2*n,n) = T(2*n,n+1) = A005430(n). (End)
T(n,k) = 2*T(n-1,k) + 2*T(n-1,k-1) - T(n-2,k) - 2*T(n-2,k-1) - T(n-2,k-2), T(1,1) = 1 and, for n > 1, T(n,k) = 0 if k <= 1 or if k > n. - Philippe Deléham, Mar 17 2012
T(n,k) = Sum_{i=1..k} i*binomial(k,i)*binomial(n+1-k,k+1-i). - Mircea Merca, Apr 11 2012
If we include a main diagonal of zeros so that the array is in the form
  0
  1   0
  2   2   0
  3   6   3   0
  4  12  12   4   0
...
then we obtain the exponential Riordan array [x*exp(x),x], which factors as [x,x]*[exp(x),x] = A132440*A007318. This array is the infinitesimal generator for A116071. A signed version of the array is the infinitesimal generator for A215652. - Peter Bala, Sep 14 2012
a(n,k) = (n-1)!/((n-k)!(k-1)!) if k > n/2 and a(n,k) = (n-1)!/((n-k-1)!k!) otherwise. [Forms 'core' for Pascal's recurrence; gives common term of RHS of T(n,k) = T(n-1,k-1) + T(n-1,k)]. - Jon Perry, Oct 08 2013
Assuming offset 0: T(n, k) = FallingFactorial(n + 1, n) / (k! * (n - k)!). The counterpart using the rising factorial is A356546. - Peter Luschny, Aug 13 2022

Edited by N. J. A. Sloane, Oct 07 2007

A035607 Table a(d,m) of number of points of L1 norm m in cubic lattice Z^d, read by antidiagonals (d >= 1, m >= 0).

Original entry on oeis.org

1, 1, 2, 1, 4, 2, 1, 6, 8, 2, 1, 8, 18, 12, 2, 1, 10, 32, 38, 16, 2, 1, 12, 50, 88, 66, 20, 2, 1, 14, 72, 170, 192, 102, 24, 2, 1, 16, 98, 292, 450, 360, 146, 28, 2, 1, 18, 128, 462, 912, 1002, 608, 198, 32, 2, 1, 20, 162, 688, 1666, 2364, 1970, 952, 258, 36, 2, 1, 22, 200, 978, 2816

Offset: 0

N. J. A. Sloane

Table also gives coordination sequences of same lattices.
Rows sums are given by A001333. Rising and falling diagonals are the tribonacci numbers A000213, A001590. - Paul Barry, Feb 13 2003
a(d,m) also gives the number of ways to choose m squares from a 2 X (d-1) grid so that no two squares in the selection are (horizontally or vertically) adjacent. - Jacob A. Siehler, May 13 2006
Mirror image of triangle A113413. - Philippe Deléham, Oct 15 2006
The Ca1 sums lead to A126116 and the Ca2 sums lead to A070550, see A180662 for the definitions of these triangle sums. - Johannes W. Meijer, Aug 05 2011
A035607 is jointly generated with the Delannoy triangle A008288 as an array of coefficients of polynomials v(n,x):  initially, u(1,x) = v(1,x) = 1; for n > 1, u(n,x) = x*u(n-1,x) + v(n-1) and v(n,x) = 2*x*u(n-1,x) + v(n-1,x).  See the Mathematica section. - Clark Kimberling, Mar 05 2012
Also, the polynomial v(n,x) above is x + (x + 1)*f(n-1,x), where f(0,x) = 1. - Clark Kimberling, Oct 24 2014
Rows also give the coefficients of the independence polynomial of the n-ladder graph. - Eric W. Weisstein, Dec 29 2017
Considering both sequences as square arrays (offset by one row), the rows of A035607 are the first differences of the rows of A008288, and the rows of A008288 are the partial sums of the rows of A035607. - Shel Kaphan, Feb 23 2023
Considering only points with nonnegative coordinates, the number of points at L1 distance = m in d dimensions is the same as the number of ways of putting m indistinguishable balls into d distinguishable urns, binomial(m+d-1, d-1). This is one facet of the cross-polytope. Allowing for + and - coordinates, there are binomial(d,i)*2^i facets containing points with up to i nonzero coordinates. Eliminating double counting of points with any coordinates = 0, there are Sum_{i=1..d} (-1)^(d-i)*binomial(m+i-1,i-1)*binomial(d,i)*2^i points at distance m in d dimensions. One may avoid the alternating sum by using binomial(m-1,i-1) to count only the points per facet with exactly i nonzero coordinates, avoiding any double counting, but the result is the same. - Shel Kaphan, Mar 04 2023

			From _Clark Kimberling_, Oct 24 2014: (Start)
As a triangle of coefficients in polynomials v(n,x) in Comments, the first 6 rows are
  1
  1   2
  1   4   2
  1   6   8   2
  1   8  18  12   2
  1  10  32  38  16   2
  ... (End)
From _Shel Kaphan_, Mar 04 2023: (Start)
For d=3, m=4:
There are binomial(3,1)*2^1 = 6 facets (vertices) of binomial(4+1-1,1-1) = 1 point with <= one nonzero coordinate.
There are binomial(3,2)*2^2 = 12 facets (edges) of binomial(4+2-1,2-1) = 5 points with <= two nonzero coordinates.
There are binomial(3,3)*2^3 = 8 facets (faces) of binomial(4+3-1,3-1) = 15 points with <= three nonzero coordinates.
a(3,4) = 8*15 - 12*5 + 6*1 = 120 - 60 + 6 = 66. (End)

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
Bela Bajnok, Additive Combinatorics: A Menu of Research Problems, arXiv:1705.07444 [math.NT], May 2017. See Sect. 2.3.
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
M. Janjic and B. Petkovic, A Counting Function, arXiv preprint arXiv:1301.4550 [math.CO], 2013. - From _N. J. A. Sloane_, Feb 13 2013
Vaclav Kotesovec, Non-attacking chess pieces, 6ed, 2013, p. 392.
R. J. Mathar, Bivariate Generating Functions for Non-attacking Wazirs on Rectangular Boards (2024), Table 2.
Emanuele Munarini, Combinatorial properties of the antichains of a garland, Integers, 9 (2009), 353-374.
Joan Serra-Sagrista, Enumeration of lattice points in l_1 norm, Inf. Proc. Lett. 76 (1-2) (2000) 39-44.
J. Siehler, Adjacency-free selections from a 2xN grid (Mathematica notebook) [broken link]
Jacob A. Siehler, Selections Without Adjacency on a Rectangular Grid, arXiv:1409.3869 [math.CO], 2014.
Eric Weisstein's World of Mathematics, Independence Polynomial
Eric Weisstein's World of Mathematics, Ladder Graph

Other versions: A113413, A119800, A122542, A266213.
Cf. A008288, which has g.f. 1/(1-x-x*y-x^2*y).
Cf. A078057 (row sums), A050146 (central terms).
Cf. A050146.

a035607 n k = a035607_tabl !! n !! k
a035607_row n = a035607_tabl !! n
a035607_tabl = map fst $ iterate
   (\(us, vs) -> (vs, zipWith (+) ([0] ++ us ++ [0]) $
                      zipWith (+) ([0] ++ vs) (vs ++ [0]))) ([1], [1, 2])
-- Reinhard Zumkeller, Jul 20 2013

A035607 := proc(d,m) local j: add(binomial(floor((d-1+j)/2),d-m-1)*binomial(d-m-1, floor((d-1-j)/2)),j=0..d-1) end: seq(seq(A035607(d,m),m=0..d-1),d=1..11); # d=dimension, m=norm # Johannes W. Meijer, Aug 05 2011

u[1, x_] := 1; v[1, x_] := 1; z = 16;
u[n_, x_] := x*u[n - 1, x] + v[n - 1, x];
v[n_, x_] := 2 x*u[n - 1, x] + v[n - 1, x];
Table[Expand[u[n, x]], {n, 1, z/2}]
Table[Expand[v[n, x]], {n, 1, z/2}]
cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
TableForm[cu]
Flatten[%]    (* A008288 *)
Table[Expand[v[n, x]], {n, 1, z}]
cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
TableForm[cv]
Flatten[%]    (* A035607 *)
(* Clark Kimberling, Mar 09 2012 *)
Reverse /@ CoefficientList[CoefficientList[Series[(1 + x)/(1 - x - x y - x^2 y), {x, 0, 10}], x], y] // Flatten (* Eric W. Weisstein, Dec 29 2017 *)

T(n, k) = if (k==0, 1, sum(i=0, k-1, binomial(n-k,i+1)*binomial(k-1,i)*2^(i+1)));
tabl(nn) = for (n=1, nn, for (k=0, n-1, print1(T(n, k), ", ")); print); \\ as a triangle; Michel Marcus, Feb 27 2018

def A035607_row(n):
    @cached_function
    def prec(n, k):
        if k==n: return 1
        if k==0: return 0
        return prec(n-1,k-1)+2*sum(prec(n-i,k-1) for i in (2..n-k+1))
    return [prec(n, n-k) for k in (0..n-1)]
for n in (1..10): print(A035607_row(n)) # Peter Luschny, Mar 16 2016

From Johannes W. Meijer, Aug 05 2011: (Start)
f(d,m) = Sum_{j=0..d-1} binomial(floor((d-1+j)/2), d-m-1)*binomial(d-m-1, floor((d-1-j)/2)), d >= 1 and 0 <= m <= d-1.
f(d,m) = f(d-1,m-1) + f(d-1,m) + f(d-2,m-1) (d >= 3 and 1 <= m <= d-1) with f(d,0) = 1 (d >= 1) and f(d,d-1) = 2 (d>=2). (End)
From Roger Cuculière, Apr 10 2006: (Start)
The generating function G(x,y) of this double sequence is the sum of a(n,p)*x^n*y^p, n=1..oo, p=0..oo, which is G(x,y) = x*(1+y)/(1-x-y-x*y).
The horizontal generating function H_n(y), which generates the rows of the table: (1, 2, 2, 2, 2, ...), (1, 4, 8, 12, 16, ...), (1, 6, 18, 38, 66, ...), is the sum of a(n,p)*y^p, p=0..oo, for each fixed n. This is H_n(y) = ((1+y)^n)/((1-y)^n).
The vertical generating function V_p(x), which generates the columns of the table: (1, 1, 1, 1, 1, ...), (2, 4, 6, 8, 10, ...), (2, 8, 18, 32, 50, ...), is the sum of a(n,p)*x^n, n=1..oo, for each fixed p. This is V_p(x) = 2*((1+x)^(p-1))/((1-x)^(p+1)) for p >= 1 and V_0(x) = x/(1-x). (End)
G.f.: (1+x)/(1-x-x*y-x^2*y). - Vladeta Jovovic, Apr 02 2002 (But see previous lines!)
T(2*n,n) = A050146(n+1). - Reinhard Zumkeller, Jul 20 2013
Seen as a triangle read by rows: T(n,0) = 1, for n > 1: T(n,n-1) = 2, T(n,k) = T(n-1,k-1) + T(n-1,k) + T(n-2,k-1), 0 < k < n. - Reinhard Zumkeller, Jul 20 2013
Seen as a triangle T(n,k) with 0 <= k < n read by rows: T(n,0)=1 for n > 0 and T(n,k) = Sum_{i=0..k-1} binomial(n-k,i+1)*binomial(k-1,i)*2^(i+1) for k > 0. - Werner Schulte, Feb 22 2018
With p >= 1 and q >= 0, as a square array a(p,q) = T(p+q-1,q) = 2*p*Hypergeometric2F1[1-p, 1-q, 2, 2] for q >= 1. Consequently, a(p,q) = a(q,p)*p/q. - Shel Kaphan, Feb 14 2023
For n >= 1, T(2*n,n) = A002003(n), T(3*n,2*n) = A103885(n) and T(4*n,3*n) = A333715(n). - Peter Bala, Jun 15 2023

More terms from David W. Wilson

Maple program corrected and information added by Johannes W. Meijer, Aug 05 2011

A122542 Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 2, -1, 0, 0, 0, 0, 0, ...] DELTA [1, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 2, 4, 1, 0, 2, 8, 6, 1, 0, 2, 12, 18, 8, 1, 0, 2, 16, 38, 32, 10, 1, 0, 2, 20, 66, 88, 50, 12, 1, 0, 2, 24, 102, 192, 170, 72, 14, 1, 0, 2, 28, 146, 360, 450, 292, 98, 16, 1, 0, 2, 32, 198, 608, 1002, 912, 462, 128, 18, 1

Offset: 0

Philippe Deléham, Sep 19 2006, May 28 2007

Riordan array (1, x*(1+x)/(1-x)). Rising and falling diagonals are the tribonacci numbers A000213, A001590.

			Triangle begins:
  1;
  0, 1;
  0, 2,  1;
  0, 2,  4,   1;
  0, 2,  8,   6,   1;
  0, 2, 12,  18,   8,    1;
  0, 2, 16,  38,  32,   10,   1;
  0, 2, 20,  66,  88,   50,  12,   1;
  0, 2, 24, 102, 192,  170,  72,  14,   1;
  0, 2, 28, 146, 360,  450, 292,  98,  16,  1;
  0, 2, 32, 198, 608, 1002, 912, 462, 128, 18, 1;

Reinhard Zumkeller, Rows n = 0..120 of triangle, flattened
Bela Bajnok, Additive Combinatorics: A Menu of Research Problems, arXiv:1705.07444 [math.NT], May 2017. See Sect. 2.3.
Huyile Liang, Yanni Pei, and Yi Wang, Analytic combinatorics of coordination numbers of cubic lattices, arXiv:2302.11856 [math.CO], 2023. See p. 1.

Other versions: A035607, A113413, A119800, A266213.
Diagonals: A000012, A005843, A001105, A035597 - A035606.
Columns: A000007, A040000, A008575, A005899, A008412-A008416, A008418, A008420, A035706 - A035745.
Sums include: A000007, A001333 (row), A001590 (diagonal), A007483, A057077 (signed row), A078016 (signed diagonal), A086901, A091928, A104934, A122558, A122690.
Cf. A155161, A059283.

a122542 n k = a122542_tabl !! n !! k
a122542_row n = a122542_tabl !! n
a122542_tabl = map fst $ iterate
   (\(us, vs) -> (vs, zipWith (+) ([0] ++ us ++ [0]) $
                      zipWith (+) ([0] ++ vs) (vs ++ [0]))) ([1], [0, 1])
-- Reinhard Zumkeller, Jul 20 2013, Apr 17 2013

function T(n, k) // T = A122542
  if k eq 0 then return 0^n;
  elif k eq n then return 1;
  else return T(n-1,k) + T(n-1,k-1) + T(n-2,k-1);
  end if;
end function;
[T(n, k): k in [0..n], n in [0..12]]; // G. C. Greubel, Oct 27 2024

CoefficientList[#, y]& /@ CoefficientList[(1-x)/(1 - (1+y)x - y x^2) + O[x]^11, x] // Flatten (* Jean-François Alcover, Sep 09 2018 *)
(* Second program *)
T[n_, k_]:= T[n, k]= If[k==n, 1, If[k==0, 0, T[n-1,k-1] +T[n-1,k] +T[n-2,k- 1] ]]; (* T = A122542 *)
Table[T[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Oct 27 2024 *)

def A122542_row(n):
    @cached_function
    def prec(n, k):
        if k==n: return 1
        if k==0: return 0
        return prec(n-1,k-1)+2*sum(prec(n-i,k-1) for i in (2..n-k+1))
    return [prec(n, k) for k in (0..n)]
for n in (0..10): print(A122542_row(n)) # Peter Luschny, Mar 16 2016

Sum_{k=0..n} x^k*T(n,k) = A000007(n), A001333(n), A104934(n), A122558(n), A122690(n), A091928(n)  for x = 0, 1, 2, 3, 4, 5. - Philippe Deléham, Jan 25 2012
Sum_{k=0..n} 3^(n-k)*T(n,k) = A086901(n).
Sum_{k=0..n} 2^(n-k)*T(n,k) = A007483(n-1), n >= 1. - Philippe Deléham, Oct 08 2006
T(2*n, n) = A123164(n).
T(n, k) = T(n-1,k) + T(n-1,k-1) + T(n-2,k-1), n > 1. - Philippe Deléham, Jan 25 2012
G.f.: (1-x)/(1-(1+y)*x-y*x^2). - Philippe Deléham, Mar 02 2012
From G. C. Greubel, Oct 27 2024: (Start)
Sum_{k=0..n} (-1)^k*T(n, k) = A057077(n).
Sum_{k=0..floor(n/2)} T(n-k, k) = A001590(n+1).
Sum_{k=0..floor(n/2)} (-1)^k*T(n-k, k) = A078016(n). (End)

A266213 Square array A(n,r), the number of neighbors at a sharp Manhattan distance r in a finite n-hypercube lattice, read by upwards antidiagonals; A(n,r) = Sum_{k=0..min(n,r)} binomial(r-1,k-1)binomial(n,k) 2^k.

Original entry on oeis.org

1, 1, 0, 1, 2, 0, 1, 4, 2, 0, 1, 6, 8, 2, 0, 1, 8, 18, 12, 2, 0, 1, 10, 32, 38, 16, 2, 0, 1, 12, 50, 88, 66, 20, 2, 0, 1, 14, 72, 170, 192, 102, 24, 2, 0, 1, 16, 98, 292, 450, 360, 146, 28, 2, 0, 1, 18, 128, 462, 912, 1002, 608, 198, 32, 2, 0

Offset: 0

Dmitry Zaitsev, Dec 24 2015

In an n-dimensional hypercube lattice, the array A(n,r) gives the number of nodes situated at a Manhattan distance equal to r, counting the current node. When counting coordinate offsets for neighboring nodes, binomial(n,k) chooses k nonzero coordinates from n coordinates, binomial(r-1,k-1) partitions the number r as the sum of exactly k nonzero numbers, and 2^k counts combinations of signs for coordinate offsets; starting indexing from 0 adds 1, which counts the current node.
In cellular automata theory, the cell surrounding with Manhattan distance less than or equal to r is called the von Neumann neighborhood of radius r or the diamond-shaped neighborhood to distinguish it from other generalizations of the von Neumann neighborhood for radius r>1, for instance, as a neighborhood having a difference in the range from -r to r in exactly one coordinate (the "narrow" von Neumann neighborhood of radius r).
The square array of partial sums of A(n,r) on rows gives us the Delannoy numbers A008288, which correspond to the number of nodes in the diamond-shaped neighborhood of radius r. - Dmitry Zaitsev, Dec 24 2015
For n >= 2, the term A(n,r) gives the number of polyominoes of bounding box 2 x (r+n-1) of area (r + 2(n-1)). Let A'(n,k) be the table A(n,k) without the first two rows. The sum of the terms in the i-th anti-diagonal of A'(n,k) gives the i-th term of A034182. - Louis Marin, Dec 11 2024

			The array A(n, k) begins:
n \ k  0  1   2   3    4     5     6      7      8      9
---------------------------------------------------------
0:     1  0   0   0    0     0     0      0      0      0
1:     1  2   2   2    2     2     2      2      2      2
2:     1  4   8  12   16    20    24     28     32     36
3:     1  6  18  38   66   102   146    198    258    326
4:     1  8  32  88  192   360   608    952   1408   1992
5:     1 10  50 170  450  1002  1970   3530   5890   9290
6:     1 12  72 292  912  2364  5336  10836  20256  35436
7:     1 14  98 462 1666  4942 12642  28814  59906 115598
8:     1 16 128 688 2816  9424 27008  68464 157184 332688
9:     1 18 162 978 4482 16722 53154 148626 374274 864146
...
For instance, in a 5-hypercube lattice there are 170 nodes situated at a Manhattan distance of 3 for a chosen node.
The triangle T(m, r) begins:
m\r 0  1   2   3   4    5   6   7  8 9 10 ...
0:  1
1:  1  0
2:  1  2   0
3:  1  4   2   0
4:  1  6   8   2   0
5:  1  8  18  12   2    0
6:  1 10  32  38  16    2   0
7:  1 12  50  88  66   20   2   0
8:  1 14  72 170 192  102  24   2  0
9:  1 16  98 292 450  360 146  28  2 0
10: 1 18 128 462 912 1002 608 198 32 2  0
... Formatted by _Wolfdieter Lang_, Jan 31 2016

Vincenzo Librandi, Table of n, a(n) for n = 0..5150
Bela Bajnok, Additive Combinatorics: A Menu of Research Problems, arXiv:1705.07444 [math.NT], May 2017. See Sect. 2.3.
Dmitry Zaitsev, k-neighborhood for Cellular Automata, arXiv preprint arXiv:1605.08870 [cs.DM], 2016.
D. A. Zaitsev, A generalized neighborhood for cellular automata, Theoretical Computer Science, 666 (2017), 21-35.

Other versions: A035607, A113413, A119800, A122542.
Partial sums on rows of A give A008288.
Cf. A001333 (row sums of T). A057077 (alternating row sums of T). - Wolfdieter Lang, Jan 31 2016

# Prints the array by rows.
gf := n -> ((1 + x)/(1 - x))^n: ser := n -> series(gf(n), x, 40):
seq(lprint(seq(coeff(ser(n), x, k), k=0..6)), n=0..9); # Peter Luschny, Mar 20 2020

Table[Sum[Binomial[r - 1, k - 1] Binomial[n - r, k] 2^k, {k, 0, Min[n - r, r]}], {n, 0, 10}, {r, 0, n}] // Flatten (* Michael De Vlieger, Dec 24 2015 *)

from sympy import binomial
def T(n, r):
    if r==0: return 1
    return sum(binomial(r - 1, k - 1) * binomial(n - r, k) * 2**k for k in range(min(n - r, r) + 1))
for n in range(11): print([T(n, r) for r in range(n + 1)]) # Indranil Ghosh, May 23 2017

A(n, 0)=1, n>=0, A(0, r)=0, r>0.
A(n, r) = A(n, r-1) + A(n-1, r-1) + A(n-1, r).
A(n, r) = Sum_{k=0..min(n,r)} binomial(r-1,k-1)*binomial(n,k)*2^k.
Triangle T(m, r) = A(m-r, r),  n >= 0, 0 <= r <= n, otherwise 0. - Wolfdieter Lang, Jan 31 2016
A(n, k) = [x^k] ((1 + x)/(1 - x))^n. - Ilya Gutkovskiy, May 23 2017

A008412 Coordination sequence for 4-dimensional cubic lattice (points on surface of 4-dimensional cross-polytope).

Original entry on oeis.org

1, 8, 32, 88, 192, 360, 608, 952, 1408, 1992, 2720, 3608, 4672, 5928, 7392, 9080, 11008, 13192, 15648, 18392, 21440, 24808, 28512, 32568, 36992, 41800, 47008, 52632, 58688, 65192, 72160, 79608, 87552, 96008, 104992, 114520, 124608, 135272

Offset: 0

N. J. A. Sloane

Coordination sequence for 4-dimensional cyclotomic lattice Z[zeta_8].

If Y_i (i=1,2,3,4) are 2-blocks of a (n+4)-set X then a(n-3) is the number of 7-subsets of X intersecting each Y_i (i=1,2,3,4). - Milan Janjic, Oct 28 2007

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from T. D. Noe)
M. Beck and S. Hosten, Cyclotomic polytopes and growth series of cyclotomic lattices, arXiv:math/0508136 [math.CO], 2005-2006.
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
Pierre de la Harpe, On the prehistory of growth of groups, arXiv:2106.02499 [math.GR], 2021.
Milan Janjic, Two Enumerative Functions
Milan Janjić, On Restricted Ternary Words and Insets, arXiv:1905.04465 [math.CO], 2019.
Ross McPhedran, Numerical Investigations of the Keiper-Li Criterion for the Riemann Hypothesis, arXiv:2311.06294 [math.NT], 2023. See p. 6.
M. O'Keeffe, Coordination sequences for lattices, Zeit. f. Krist., 210 (1995), 905-908.
M. O'Keeffe, Coordination sequences for lattices, Zeit. f. Krist., 210 (1995), 905-908. [Annotated scanned copy]
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Cf. A001845, A005899, A006527, A035598.
First differences of A001846.
Row 4 of A035607, A266213.
Column 4 of A113413, A119800, A122542.

I:=[1,8,32,88,192]; [n le 5 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..50]]; // Vincenzo Librandi, Jan 15 2018

Maple
```
8/3*n^3+16/3*n;
```

CoefficientList[Series[((1+x)/(1-x))^4,{x,0,40}],x] (* or *)
LinearRecurrence[{4, -6, 4, -1}, {1, 8, 32, 88, 192}, 41] (* Harvey P. Dale, Jun 10 2011 *)
f[n_] := 8 n (n^2 + 2)/3; f[0] = 1; Array[f, 38, 0] (* or *)
g[n_] := 4n^2 +2; f[n_] := f[n-1] + g[n] + g[n -1]; f[0] = 1; f[1] = 8; Array[f, 38, 0] (* Robert G. Wilson v, Dec 27 2017 *)

a(n)=if(n,8*(n^2+2)*n/3,1) \\ Charles R Greathouse IV, Jun 10 2011

G.f.: ((1+x)/(1-x))^4.
a(n) = 8*n*(n^2+2)/3 for n>1.
a(n) = 8*A006527(n) for n>0.
a(n) = A005899(n) + A005899(n-1) + a(n-1). - Bruce J. Nicholson, Dec 17 2017
n*a(n) = 8*a(n-1) + (n-2)*a(n-2) for n > 1. - Seiichi Manyama, Jun 06 2018
a(n) = 2*d*Hypergeometric2F1(1-d, 1-n, 2, 2) where d=4, for n>=1. - Shel Kaphan, Feb 14 2023
a(n) = A035598(n)*4/n, for n>0. - Shel Kaphan, Feb 28 2023
E.g.f.: 1 + 8*exp(x)*x*(3 + 3*x + x^2)/3. - Stefano Spezia, Mar 14 2024

A113413 A Riordan array of coordination sequences.

Original entry on oeis.org

1, 2, 1, 2, 4, 1, 2, 8, 6, 1, 2, 12, 18, 8, 1, 2, 16, 38, 32, 10, 1, 2, 20, 66, 88, 50, 12, 1, 2, 24, 102, 192, 170, 72, 14, 1, 2, 28, 146, 360, 450, 292, 98, 16, 1, 2, 32, 198, 608, 1002, 912, 462, 128, 18, 1, 2, 36, 258, 952, 1970, 2364, 1666, 688, 162, 20, 1, 2, 40, 326

Offset: 0

Paul Barry, Oct 29 2005

Columns include A040000, A008574, A005899, A008412, A008413, A008414. Row sums are A078057(n)=A001333(n+1). Diagonal sums are A001590(n+3). Reverse of A035607. Signed version is A080246. Inverse is A080245.
For another version see A122542. - Philippe Deléham, Oct 15 2006
T(n,k) is the number of length n words on alphabet {0,1,2} with no two consecutive 1's and no two consecutive 2's and having exactly k 0's. - Geoffrey Critzer, Jun 11 2015
From Eric W. Weisstein, Feb 17 2016: (Start)
Triangle of coefficients (from low to high degree) of x^-n * vertex cover polynomial of the n-ladder graph P_2 \square p_n:
  Psi_{L_1}: x*(2 + x) -> {2, 1}
  Psi_{L_2}: x^2*(2 + 4 x + x^2) -> {2, 4, 1}
  Psi_{L_3}: x^3*(2 + 8 x + 6 x^2 + x^3) -> {2, 8, 6, 1}
(End)
Let c(n, k), n > 0, be multiplicative sequences for some fixed integer k >= 0 with c(p^e, k) = T(e+k, k) for prime p and e >= 0. Then we have Dirichlet g.f.: Sum_{n>0} c(n, k) / n^s = zeta(s)^(2*k+2) / zeta(2*s)^(k+1). Examples: For k = 0 see A034444 and for k = 1 see A322328. Dirichlet convolution of c(n, k) and lambda(n) is Dirichlet inverse of c(n, k). - Werner Schulte, Oct 31 2022

			Triangle begins
  1;
  2,  1;
  2,  4,  1;
  2,  8,  6,  1;
  2, 12, 18,  8,  1;
  2, 16, 38, 32, 10,  1;
  2, 20, 66, 88, 50, 12,  1;

G. C. Greubel, Table of n, a(n) for the first 50 rows
Bela Bajnok, Additive Combinatorics: A Menu of Research Problems, arXiv:1705.07444 [math.NT], May 2017. See Sect. 2.3.
Paul Barry, Riordan arrays, generalized Narayana triangles, and series reversion, Linear Algebra and its Applications, 491 (2016) 343-385.
P. Holub, M. Miller, H. Perez-Roses, and J. Ryan, Degree diameter problem on honeycomb networks, Disc. Appl. Math. 179 (2014) 139-151.
Milan Janjić, Words and Linear Recurrences, J. Int. Seq. 21 (2018), #18.1.4.
Huyile Liang, Yanni Pei, and Yi Wang, Analytic combinatorics of coordination numbers of cubic lattices, arXiv:2302.11856 [math.CO], 2023. See p. 2.
Mirka Miller, Hebert Perez-Roses, and Joe Ryan, The Maximum Degree-and-Diameter-Bounded Subgraph in the Mesh, arXiv preprint arXiv:1203.4069 [math.CO], 2012.

Other versions: A035607, A119800, A122542, A266213.

nn = 10; Map[Select[#, # > 0 &] &, CoefficientList[Series[1/(1 - 2 x/(1 + x) - y x), {x, 0, nn}], {x, y}]] // Grid (* Geoffrey Critzer, Jun 11 2015 *)
CoefficientList[CoefficientList[Series[1/(1 - 2 x/(1 + x) - y x), {x, 0, 10}, {y, 0, 10}], x], y] (* Eric W. Weisstein, Feb 17 2016 *)

T = lambda n,k : binomial(n, k)*hypergeometric([-k-1, k-n], [-n], -1).simplify_hypergeometric()
A113413 = lambda n,k : 1 if n==0 and k==0 else T(n, k)
for n in (0..12): print([A113413(n,k) for k in (0..n)]) # Peter Luschny, Sep 17 2014 and Mar 16 2016

# Alternatively:
def A113413_row(n):
    @cached_function
    def prec(n, k):
        if k==n: return 1
        if k==0: return 0
        return prec(n-1,k-1)+2*sum(prec(n-i,k-1) for i in (2..n-k+1))
    return [prec(n, k) for k in (1..n)]
for n in (1..10): print(A113413_row(n)) # Peter Luschny, Mar 16 2016

From Paul Barry, Nov 13 2005: (Start)
Riordan array ((1+x)/(1-x), x(1+x)/(1-x)).
T(n, k) = Sum_{i=0..n-k} C(k+1, i)*C(n-i, k).
T(n, k) = Sum_{j=0..n-k} C(k+j, j)*C(k+1, n-k-j).
T(n, k) = D(n, k) + D(n-1, k) where D(n, k) = Sum_{j=0..n-k} C(n-k, j)*C(k, j)*2^j = A008288(n, k).
T(n, k) = T(n-1, k) + T(n-1, k-1) + T(n-2, k-1).
T(n, k) = Sum_{j=0..n} C(floor((n+j)/2), k)*C(k, floor((n-j)/2)). (End)
T(n, k) = C(n, k)*hypergeometric([-k-1, k-n], [-n], -1). - Peter Luschny, Sep 17 2014
T(n, k) = (Sum_{i=2..k+2} A137513(k+2, i) * (n-k)^(i-2)) / (k!) for 0 <= k < n (conjectured). - Werner Schulte, Oct 31 2022

A035597 Number of points of L1 norm 3 in cubic lattice Z^n.

Original entry on oeis.org

0, 2, 12, 38, 88, 170, 292, 462, 688, 978, 1340, 1782, 2312, 2938, 3668, 4510, 5472, 6562, 7788, 9158, 10680, 12362, 14212, 16238, 18448, 20850, 23452, 26262, 29288, 32538, 36020, 39742, 43712, 47938, 52428, 57190, 62232, 67562

Offset: 0

N. J. A. Sloane

Sums of the first n terms > 0 of A001105 in palindromic arrangement. a(n) = Sum_{i=1 .. n} A001105(i) + Sum_{i=1 .. n-1} A001105(i), e.g. a(3) = 38 = 2 + 8 + 18 + 8 + 2; a(4) = 88 = 2 + 8 + 18 + 32 + 18 + 8 + 2. - Klaus Purath, Jun 19 2020

Apart from multiples of 3, all divisors of n are also divisors of a(n), i.e. if n is not divisible by 3, a(n) is divisible by n. All divisors d of a(n) for d !== 0 (mod) 3 are also divisors of a(abs(n-d)) and a(n+d). For all n congruent to 0,2,7 (mod 9) a(n) is divisible by 3. If n is divisible by 3^k, a(n) is divisible by 3^(k-1). - Klaus Purath, Jul 24 2020

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
J. H. Conway and N. J. A. Sloane, Low-dimensional lattices. VII. Coordination sequences, Proc. Roy. Soc. Lond. A 458 (1996) 2369-2389.
Milzn Janjic, On a class of polynomials with integer coefficients, JIS 11 (2008) 08.5.2.
Milan Janjic and B. Petkovic, A Counting Function, arXiv preprint arXiv:1301.4550 [math.CO], 2013. - _N. J. A. Sloane_, Feb 13 2013
Milan Janjic and B. Petkovic, A Counting Function Generalizing Binomial Coefficients and Some Other Classes of Integers, J. Int. Seq. 17 (2014) # 14.3.5.
R. J. Mathar, Bivariate generating functions for non-attacking Wazirs on rectangular boards, viXra:2404.0122 (2024) page 14
Joan Serra-Sagrista, Enumeration of lattice points in l_1 norm, Inf. Proc. Lett. 76 (1-2) (2000) 39-44.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Partial sums of A069894.
Cf. A001105, A005900, A084570, A097869.
Cf. A001844, A001845, A005899.
Column 3 of A035607, A266213, A343599.
Row 3 of A113413, A119800, A122542.

[(4*n^3 + 2*n)/3: n in [0..40]]; // Vincenzo Librandi, Sep 19 2011

f := proc(n,m) local i; sum( 2^i*binomial(n,i)*binomial(m-1,i-1),i=1..min(n,m)); end; # n=dimension, m=norm

Table[(4n^3+2n)/3,{n,0,40}] (* or *) LinearRecurrence[{4,-6,4,-1},{0,2,12,38},41] (* Harvey P. Dale, Sep 18 2011 *)

a(n) = (4*n^3 + 2*n)/3.
a(n) = 2*A005900(n). - R. J. Mathar, Dec 05 2009
a(0)=0, a(1)=2, a(2)=12, a(3)=38, a(n)=4*a(n-1)-6*a(n-2)+4*a(n-3)-a(n-4). G.f.: (2*x*(x+1)^2)/(x-1)^4. - Harvey P. Dale, Sep 18 2011
a(n) = -a(-n), a(n+1) = A097869(4n+3) = A084570(2n+1). - Bruno Berselli, Sep 20 2011
a(n) = 2*n*Hypergeometric2F1(1-n,1-k,2,2), where k=3. Also, a(n) = A001845(n) - A001844(n). - Shel Kaphan, Feb 26 2023
a(n) = A005899(n)*n/3. - Shel Kaphan, Feb 26 2023
a(n) = A006331(n)+A006331(n-1). - R. J. Mathar, Aug 12 2025

A094305 Triangle read by rows: T(n,k) = ((n+1)(n+2)/2) * binomial(n,k) (0 <= k <= n).

Original entry on oeis.org

1, 3, 3, 6, 12, 6, 10, 30, 30, 10, 15, 60, 90, 60, 15, 21, 105, 210, 210, 105, 21, 28, 168, 420, 560, 420, 168, 28, 36, 252, 756, 1260, 1260, 756, 252, 36, 45, 360, 1260, 2520, 3150, 2520, 1260, 360, 45, 55, 495, 1980, 4620, 6930, 6930, 4620, 1980, 495, 55, 66

Offset: 0

Amarnath Murthy, Apr 29 2004

Sum of all possible sums of k+1 numbers chosen from among the first n+1 numbers. Additive analog of triangle of Stirling numbers of first kind (A008275). - David Wasserman, Oct 04 2007
Third slice along the 1-2-plane in the cube a(m,n,o) = a(m-1,n,o)+a(m,n-1,o)+a(m,n,o-1) with a(1,0,0)=1 and a(m<>1=0,n>=0,0>=o)=0, for which the first slice is Pascal's triangle (slice read by antidiagonals). - Thomas Wieder, Aug 06 2006
Triangle T(n,k), 0<=k<=n, read by rows given by [3,-1,2/3,-1/6,1/2,0,0,0,0,0,0,...] DELTA [3,-1,2/3,-1/6,1/2,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Oct 07 2007
T(n,k) is the number of ordered triples of bit strings with n bits and exactly k 1's over all bits in the triple.  For example for n=1 we have (0,e,e),(e,0,e),(e,e,0),(1,e,e),(e,1,e),(e,e,1) where e is the empty string. - Geoffrey Critzer, Apr 06 2013
T(n,k) = A000217(n+1) * A007318(n,k), 0 <= k <= n. - Reinhard Zumkeller, Jul 30 2013

			Triangle begins:
  1
  3 3
  6 12 6
  10 30 30 10
  15 60 90 60 15
  21 105 210 210 105 21
  ...
The n-th row is the product of the n-th triangular number and the n-th row of Pascal's triangle. The fifth row is (15,60,90,60,15) or 15*{1,4,6,4,1}.

A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, identity 152.

Reinhard Zumkeller, Rows n = 0..100 of table, flattened
Mircea Merca, A Special Case of the Generalized Girard-Waring Formula, J. Integer Sequences, Vol. 15 (2012), Article 12.5.7.

For a closely related array that also includes a row and column of zeros see A129533.
Columns include A000217. Row sums are A001788. Cf. A094306.
Cf. A003506, A121547, A121306, A119800, A000217, A007318.

a094305 n k = a094305_tabl !! n !! k
a094305_row n = a094305_tabl !! n
a094305_tabl = zipWith (map . (*)) (tail a000217_list) a007318_tabl
-- Reinhard Zumkeller, Jul 30 2013

A094305:= proc(n,k) (n+1)*(n+2)/2 * binomial(n,k); end;

nn=10; f[list_]:=Select[list,#>0&];a=1/(1-x-y x); Map[f,CoefficientList[Series[a^3,{x,0,nn}],{x,y}]]//Grid
(* Geoffrey Critzer, Apr 06 2013 *)
Flatten[Table[((n+1)(n+2))/2 Binomial[n,k],{n,0,10},{k,0,n}]] (* Harvey P. Dale, Aug 31 2014 *)

T(n,k) = Sum_{i=1..k+1} (-1)^(i+1)*i^2*binomial(n+2,k+i+1)*binomial(n+2,k-i+1). - Mircea Merca, Apr 05 2012

O.g.f.: 1/(1 - x - y*x)^3. - Geoffrey Critzer, Apr 06 2013

Edited by Ralf Stephan, Feb 04 2005
Further comments from David Wasserman, Oct 04 2007
Further editing by N. J. A. Sloane, Oct 07 2007

A008413 Coordination sequence for 5-dimensional cubic lattice.

Original entry on oeis.org

1, 10, 50, 170, 450, 1002, 1970, 3530, 5890, 9290, 14002, 20330, 28610, 39210, 52530, 69002, 89090, 113290, 142130, 176170, 216002, 262250, 315570, 376650, 446210, 525002, 613810, 713450, 824770, 948650, 1086002, 1237770, 1404930

Offset: 0

N. J. A. Sloane

If Y_i (i=1,2,3,4,5) are 2-blocks of a (n+5)-set X then a(n-4) is the number of 9-subsets of X intersecting each Y_i (i=1,2,3,4,5). - Milan Janjic, Oct 28 2007

Seiichi Manyama, Table of n, a(n) for n = 0..10000
Milan Janjic, Two Enumerative Functions
Milan Janjić, On Restricted Ternary Words and Insets, arXiv:1905.04465 [math.CO], 2019.
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5,1).

Cf. A035599.
Row 5 of A035607, A266213.
Column 5 of A113413, A119800, A122542.

Maple
```
4/3*n^4+20/3*n^2+2;
```

LinearRecurrence[{5,-10,10,-5,1},{1,10,50,170,450,1002},40] (* Harvey P. Dale, May 02 2016 *)
{1}~Join~Table[4/3 n^4 + 20/3 n^2 + 2, {n, 32}] (* or *)
CoefficientList[Series[((1 + x)/(1 - x))^5, {x, 0, 32}], x] (* Michael De Vlieger, Oct 04 2016 *)

G.f.: ((1+x)/(1-x))^5.
a(n) = (4/3)*n^4 + (20/3)*n^2 + 2 for n > 0. - Michael De Vlieger, Oct 04 2016
n*a(n) = 10*a(n-1) + (n-2)*a(n-2) for n > 1. - Seiichi Manyama, Jun 06 2018
From Shel Kaphan, Mar 03 2023: (Start)
a(n) = 2*d*Hypergeometric2F1(1-d, 1-n, 2, 2) where d=5, for n>=1.
a(n) = A035599(n)*5/n, for n>0. (End)

cp's OEIS Frontend

A005899 Number of points on surface of octahedron; also coordination sequence for cubic lattice: a(0) = 1; for n > 0, a(n) = 4n^2 + 2.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A003506 Triangle of denominators in Leibniz's Harmonic Triangle a(n,k), n >= 1, 1 <= k <= n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

SageMath

Formula

Extensions

A035607 Table a(d,m) of number of points of L1 norm m in cubic lattice Z^d, read by antidiagonals (d >= 1, m >= 0).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Sage

Formula

Extensions

A122542 Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 2, -1, 0, 0, 0, 0, 0, ...] DELTA [1, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

Sage

Formula

A266213 Square array A(n,r), the number of neighbors at a sharp Manhattan distance r in a finite n-hypercube lattice, read by upwards antidiagonals; A(n,r) = Sum_{k=0..min(n,r)} binomial(r-1,k-1)*binomial(n,k)* 2^k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Python

Formula

A008412 Coordination sequence for 4-dimensional cubic lattice (points on surface of 4-dimensional cross-polytope).

Views

Author

A266213 Square array A(n,r), the number of neighbors at a sharp Manhattan distance r in a finite n-hypercube lattice, read by upwards antidiagonals; A(n,r) = Sum_{k=0..min(n,r)} binomial(r-1,k-1)binomial(n,k) 2^k.