cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 19 results. Next

A287324 a(n) = A008412(n-1) + A008412(n-2) for n>1, a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 9, 40, 120, 280, 552, 968, 1560, 2360, 3400, 4712, 6328, 8280, 10600, 13320, 16472, 20088, 24200, 28840, 34040, 39832, 46248, 53320, 61080, 69560, 78792, 88808, 99640, 111320, 123880, 137352, 151768, 167160, 183560, 201000, 219512, 239128, 259880, 281800
Offset: 0

Views

Author

Leo James Borcherding, May 23 2017

Keywords

Comments

Let's iteratively apply the summation of two consecutive terms to A000292. It generates A000330, then A005900, then A001845, then A008412, then this sequence. Every sequence in this series starts with 1 followed by the sum of 1 and the next term in the previous sequence; because of that, for A008412 and this sequence, the initial term(s) are exceptions from the general formula.
From Leo James Borcherding, May 23 2017: (Start)
a(n) = f(9,n), where f(k,n) is the set of all series derived from the anchored series.
k = (All whole numbers (including negative values))
n = (All whole numbers >= 1)
The anchored series is f(0,n).
See the attached file for an in-depth explanation of the family of tetrahedron sequences that f(9,n) (this sequence) is a part of.
A Visual Representation of the summation process is as follows:
a.) f(7,n) + f(7,n-1) = f(8,n)
b.) f(8,n) + f(8,n-1) = f(9,n)
a.) b.)
1 + 0 = 1 1 + 0 = 1
7 + 1 = 8 8 + 1 = 9
25 + 7 = 32 32 + 8 = 40
63 + 25 = 88 88 + 32 = 120
129 + 63 = 192 192 + 88 = 280
231 + 129 = 360 360 + 192 = 552
377 + 231 = 608 608 + 360 = 968
575 + 377 = 952 952 + 608 = 1560
... iterate infinitely many times. (End)

References

  • William Dunham, Euler The Master of Us All, The Mathematical Association of America, 1999 p. 40.
  • Joseph and Frances Gies, Leonard of Pisa and the New Mathematics of the Middle Ages, Thomas Y. Crowell Company New York, 1969, p. 78.

Crossrefs

Programs

  • PARI
    concat(0, Vec(x*(x+1)^5/(x-1)^4 + O(x^30))) \\ Michel Marcus, May 24 2017

Formula

G.f.: x*(x + 1)^5 / (x - 1)^4.
a(n) = 8*(n - 1)*((n - 1)^2 + 2)/3 + 8*(n - 2)*((n - 2)^2 + 2)/3 = 8*(2*n - 3)*(n^2 - 3*n + 5)/3 for n>2, a(0)=0, a(1)=1, a(2)=9.
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) for n>6. - Colin Barker, Jun 05 2017

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

Original entry on oeis.org

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

Views

Author

Keywords

Comments

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

References

  • 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).

Crossrefs

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.
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.

Programs

Formula

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

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.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

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

References

  • 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).

Crossrefs

Partial sums give A001845.
Column 2 * 2 of array A188645.
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.

Programs

Formula

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

A058187 Expansion of (1+x)/(1-x^2)^4: duplicated tetrahedral numbers.

Original entry on oeis.org

1, 1, 4, 4, 10, 10, 20, 20, 35, 35, 56, 56, 84, 84, 120, 120, 165, 165, 220, 220, 286, 286, 364, 364, 455, 455, 560, 560, 680, 680, 816, 816, 969, 969, 1140, 1140, 1330, 1330, 1540, 1540, 1771, 1771, 2024, 2024, 2300, 2300, 2600, 2600, 2925, 2925, 3276, 3276
Offset: 0

Views

Author

Henry Bottomley, Nov 20 2000

Keywords

Comments

For n >= i, i = 6,7, a(n - i) is the number of incongruent two-color bracelets of n beads, i of which are black (cf. A005513, A032280), having a diameter of symmetry. The latter means the following: if we imagine (0,1)-beads as points (with the corresponding labels) dividing a circumference of a bracelet into n identical parts, then a diameter of symmetry is a diameter (connecting two beads or not) such that a 180-degree turn of one of two sets of points around it (obtained by splitting the circumference by this diameter) leads to the coincidence of the two sets (including their labels). - Vladimir Shevelev, May 03 2011
From Johannes W. Meijer, May 20 2011: (Start)
The Kn11, Kn12, Kn13, Fi1 and Ze1 triangle sums, see A180662 for their definitions, of the Connell-Pol triangle A159797 are linear sums of shifted versions of the duplicated tetrahedral numbers, e.g., Fi1(n) = a(n-1) + 5*a(n-2) + a(n-3) + 5*a(n-4).
The Kn11, Kn12, Kn13, Kn21, Kn22, Kn23, Fi1, Fi2, Ze1 and Ze2 triangle sums of the Connell sequence A001614 as a triangle are also linear sums of shifted versions of the sequence given above. (End)
The number of quadruples of integers [x, u, v, w] that satisfy x > u > v > w >= 0, n + 5 = x + u. - Michael Somos, Feb 09 2015
Also, this sequence is the fourth column in the triangle of the coefficients of the sum of two consecutive Fibonacci polynomials F(n+1, x) and F(n, x) (n>=0) in ascending powers of x. - Mohammad K. Azarian, Jul 18 2018

Crossrefs

Cf. A057884. Sum of 2 consecutive terms gives A006918, whose sum of 2 consecutive terms gives A002623, whose sum of 2 consecutive terms gives A000292, which is this sequence without the duplication. Continuing to sum 2 consecutive terms gives A000330, A005900, A001845, A008412 successively.

Programs

  • Haskell
    a058187 n = a058187_list !! n
    a058187_list = 1 : f 1 1 [1] where
       f x y zs = z : f (x + y) (1 - y) (z:zs) where
         z = sum $ zipWith (*) [1..x] [x,x-1..1]
    -- Reinhard Zumkeller, Dec 21 2011
    
  • Maple
    A058187:= proc(n) option remember; A058187(n):= binomial(floor(n/2)+3, 3) end: seq(A058187(n), n=0..51); # Johannes W. Meijer, May 20 2011
  • Mathematica
    a[n_]:= Length @ FindInstance[{x>u, u>v, v>w, w>=0, x+u==n+5}, {x, u, v, w}, Integers, 10^9]; (* Michael Somos, Feb 09 2015 *)
    With[{tetra=Binomial[Range[30]+2,3]},Riffle[tetra,tetra]] (* Harvey P. Dale, Mar 22 2015 *)
  • PARI
    {a(n) = binomial(n\2+3, 3)}; /* Michael Somos, Jun 07 2005 */
    
  • Sage
    [binomial((n//2)+3, 3) for n in (0..60)] # G. C. Greubel, Feb 18 2022

Formula

a(n) = A006918(n+1) - a(n-1).
a(2*n) = a(2*n+1) = A000292(n) = (n+1)*(n+2)*(n+3)/6.
a(n) = (2*n^3 + 21*n^2 + 67*n + 63)/96 + (n^2 + 7*n + 11)(-1)^n/32. - Paul Barry, Aug 19 2003
a(n) = A108299(n-3,n)*(-1)^floor(n/2) for n > 2. - Reinhard Zumkeller, Jun 01 2005
Euler transform of finite sequence [1, 3]. - Michael Somos, Jun 07 2005
G.f.: 1 / ((1 - x) * (1 - x^2)^3) = 1 / ((1 + x)^3 * (1 - x)^4). a(n) = -a(-7-n) for all n in Z.
a(n) = binomial(floor(n/2) + 3, 3). - Vladimir Shevelev, May 03 2011
a(-n) = -a(n-7); a(n) = A000292(A008619(n)). - Guenther Schrack, Sep 13 2018
Sum_{n>=0} 1/a(n) = 3. - Amiram Eldar, Aug 18 2022

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

Views

Author

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

Keywords

Comments

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

Examples

			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;
		

Crossrefs

Other versions: A035607, A113413, A119800, A266213.
Sums include: A000007, A001333 (row), A001590 (diagonal), A007483, A057077 (signed row), A078016 (signed diagonal), A086901, A091928, A104934, A122558, A122690.

Programs

  • Haskell
    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
    
  • Magma
    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
  • Mathematica
    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 *)
  • Sage
    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
    

Formula

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)

A144097 The 4-Schroeder numbers: a(n) = number of lattice paths (Schroeder paths) from (0,0) to (3n,n) with unit steps N=(0,1), E=(1,0) and D=(1,1) staying weakly above the line y = 3x.

Original entry on oeis.org

1, 2, 14, 134, 1482, 17818, 226214, 2984206, 40503890, 561957362, 7934063678, 113622696470, 1646501710362, 24098174350986, 355715715691350, 5289547733908510, 79163575684710818, 1191491384838325474
Offset: 0

Views

Author

Joachim Schroeder (schroderjd(AT)qwa.uovs.ac.za), Sep 10 2008

Keywords

Comments

a(n) is also the number of lattice path from (0,0) to (4n,0) with unit steps (1,3), (2,2) and (1,-1) staying weakly above the x-axis.
Also, the number of planar rooted trees with n non-leaf vertices such that each non-leaf vertex has either 3 or 4 children. - Cameron Marcott, Sep 18 2013
a(n) equals the number of ordered complete 4-ary trees with 3*n + 1 leaves, where the internal vertices come in two colors and such that each vertex and its rightmost child have different colors. See Drake, Example 1.6.9. - Peter Bala, Apr 30 2023

Examples

			a(2)=14, because
  01: NNNENNNE,
  02: NNDNNNE,
  03: NNNENND,
  04: NNDNND,
  05: NNNDNNE,
  06: NNNDND,
  07: NNNNENNE,
  08: NNNNEND,
  09: NNNNDNE,
  10: NNNNDD,
  11: NNNNNENE,
  12: NNNNNED,
  13: NNNNNDE,
  14: NNNNNNEE
are all the paths from (0,0) to (2,6) with steps N,E and D weakly above y=3x.
		

References

  • Sheng-Liang Yang and Mei-yang Jiang, The m-Schröder paths and m-Schröder numbers, Disc. Math. (2021) Vol. 344, Issue 2, 112209. doi:10.1016/j.disc.2020.112209. See Table 1.

Crossrefs

Cf. A027307 (the case y=2x), A008288 (Delannoy numbers), A008412 (4-dimensional coordination numbers).
This appears to equal 2*A243675. - N. J. A. Sloane, Mar 28 2021
The sequences listed in Yang-Jiang's Table 1 appear to be A006318, A001003, A027307, A034015, A144097, A243675, A260332, A243676. - N. J. A. Sloane, Mar 28 2021

Programs

  • Maple
    Schr:=proc(n,m,l)(n-l*m+1)/m*sum(2^v*binomial(m,v)*binomial(n,v-1),v=1..m) end proc; where n=3m and l=3, also
    Schr:=proc(n,m,l)(n-l*m+1)/(n+1)*sum(2^v*binomial(m-1,v-1)*binomial(n+1,v),v=0..m) end proc; where n=3m and l=3, also
    Schr:=proc(n,m,l)(n-l*m+1)/m*sum(binomial(m,v)*binomial(n+v,m-1),v=0..m) end proc; where n=3m and l=3, also
    Schr:=proc(n,m,l)(n-l*m+1)/(n+1)*sum(binomial(n+1,v)*binomial(m-1+v,n),v=0..n+1) end proc; where n=3m and l=3.
    # alternative Maple program:
    a:= proc(n) option remember; `if`(n<2, n+1,
          ((15610*n^5 -67123*n^4 +106824*n^3 -77633*n^2
           +25514*n-3000)*a(n-1) -(3*(n-2))*(3*n-4)*
           (3*n-5)*(35*n^2-28*n+5)*a(n-2)) / ((3*(3*n-1))
           *(3*n+1)*n*(35*n^2-98*n+68)))
        end:
    seq(a(n), n=0..20);  # Alois P. Heinz, May 26 2015
  • Mathematica
    d[n_, k_] := Binomial[n+k, k] Hypergeometric2F1[-k, -n, -n-k, -1]; a[0] = 1; a[n_] = d[3n, n] - 3d[3n+1, n-1] - 2d[3n, n-1]; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Feb 25 2017 *)
  • PARI
    {a(n) = sum(k=0, n, binomial(n, k) * binomial(3*n+k+1, n)/(3*n+k+1))} \\ Seiichi Manyama, Jul 25 2020
    
  • PARI
    {a(n) = if(n==0, 1, sum(k=1, n, 2^k*binomial(n, k) * binomial(3*n, k-1)/n))} \\ Seiichi Manyama, Jul 25 2020

Formula

G.f. A(z) satisfies A(z) = 1 + z(A(z)^3 + A(z)^4) a(n)= S_{3n+1}(n) - 3S_n(3n + 1), where S_a(b) are coordination numbers, i.e., the number of points in the a-dimensional cubic lattice Z^a having distance b in the L_1 norm.
Also a(n) = D(3n,n) - 3D(3n + 1,n-1) - 2D(3n,n-1), where D(a,b) are the Delannoy numbers, i.e., the number of paths with N, E and D steps from (0,0) to (a,b).
D-finite with recurrence 3*n*(3*n-1)*(3*n+1)*(35*n^2-98*n+68) *a(n) +(-15610*n^5+67123*n^4-106824*n^3+77633*n^2-25514*n+3000)*a(n-1) +3*(n-2) *(3*n-4) *(3*n-5) *(35*n^2-28*n+5) *a(n-2)=0. - R. J. Mathar, Sep 06 2016
From Seiichi Manyama, Jul 25 2020: (Start)
a(n) = Sum_{k=0..n} binomial(n,k) * binomial(3*n+k+1, n)/(3*n+k+1).
a(n) = (1/n) * Sum_{k=1..n} 2^k * binomial(n,k) * binomial(3*n,k-1) for n > 0. (End)
a(n) ~ sqrt(12160 + 3853*sqrt(10)) * 3^(3*n - 9/2) / (2*sqrt(5*Pi) * n^(3/2) * (223 - 70*sqrt(10))^(n - 1/2)). - Vaclav Kotesovec, Jul 31 2021
Series reversion of x*(1 - x^3)/(1 + x^3) = x + 2*x^4 + 14*x^7 + 134*x^10 + ... = Sum_{n >= 0} a(n)*x^(3*n+1). - Peter Bala, Apr 30 2023
From Peter Bala, Jun 16 2023: (Start)
The g.f. A(x) = 1 + 2*x + 14*x^2 + 134*x^3 + ... satisfies A(x)^3 = (1/x) * the series reversion of ((1 - x)/(1 + x))^3.
Define b(n) = [x^(3*n)] ( (1 + x)/(1 - x) )^n = (1/3) * [x^n] ((1 + x)/(1 - x))^(3*n) = A333715(n). Then A(x) = exp( Sum_{n >= 1} b(n)*x^n/n ).
a(n) = 2*hypergeom([1 - n, -3*n], [2], 2) for n >= 1. (End)
a(n) = (1/n) * Sum_{k=0..n-1} (-1)^k * 2^(n-k) * binomial(n,k) * binomial(4*n-k,n-1-k) for n > 0. - Seiichi Manyama, Aug 09 2023

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

Views

Author

Paul Barry, Oct 29 2005

Keywords

Comments

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

Examples

			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;
		

Crossrefs

Other versions: A035607, A119800, A122542, A266213.

Programs

  • Mathematica
    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 *)
  • Sage
    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
    
  • Sage
    # 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

Formula

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

A119800 Array of coordination sequences for cubic lattices (rows) and of numbers of L1 forms in cubic lattices (columns) (array read by antidiagonals).

Original entry on oeis.org

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

Views

Author

Thomas Wieder, Jul 30 2006, Aug 06 2006

Keywords

Examples

			The second row of the table is: 6, 18, 38, 66, 102, 146, 198, 258, 326, ... = A005899 = number of points on surface of octahedron.
The third column of the table is: 12, 38, 88, 170, 292, 462, 688, 978, 1340, ... = A035597 = number of points of L1 norm 3 in cubic lattice Z^n.
The first rows are: A008574, A005899, A008412, A008413, A008414, A008415, A008416, A008418, A008420.
The first columns are: A005843, A001105, A035597, A035598, A035599, A035600, A035601, A035602, A035603.
The main diagonal seems to be A050146.
Square array A(n,k) begins:
   4,   8,   12,   16,    20,    24,     28,     32,      36, ...
   6,  18,   38,   66,   102,   146,    198,    258,     326, ...
   8,  32,   88,  192,   360,   608,    952,   1408,    1992, ...
  10,  50,  170,  450,  1002,  1970,   3530,   5890,    9290, ...
  12,  72,  292,  912,  2364,  5336,  10836,  20256,   35436, ...
  14,  98,  462, 1666,  4942, 12642,  28814,  59906,  115598, ...
  16, 128,  688, 2816,  9424, 27008,  68464, 157184,  332688, ...
  18, 162,  978, 4482, 16722, 53154, 148626, 374274,  864146, ...
  20, 200, 1340, 6800, 28004, 97880, 299660, 822560, 2060980, ...
		

Crossrefs

Programs

  • Maple
    A:= proc(m, n)  option remember;
          `if`(n=0, 1, `if`(m=0, 2, A(m, n-1) +A(m-1, n) +A(m-1, n-1)))
        end:
    seq(seq(A(n, 1+d-n), n=1..d), d=1..10);  # Alois P. Heinz, Apr 21 2012
  • Mathematica
    A[m_, n_] := A[m, n] = If[n == 0, 1, If[m == 0, 2, A[m, n-1] + A[m-1, n] + A[m-1, n-1]]]; Table[Table[A[n, 1+d-n], {n, 1, d}], {d, 1, 10}] // Flatten (* Jean-François Alcover, Mar 09 2015, after Alois P. Heinz *)

Formula

A(m,n) = A(m,n-1) + A(m-1,n) + A(m-1,n-1), A(m,0)=1, A(0,0)=1, A(0,n)=2.

Extensions

Offset and typos corrected by Alois P. Heinz, Apr 21 2012

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

Views

Author

Keywords

Comments

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

Examples

			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.
		

References

  • 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).

Crossrefs

First differences are A008412.
Cf. A240876.
Row/column 4 of A008288.

Programs

  • Maple
    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
  • Mathematica
    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 *)

Formula

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

A019560 Coordination sequence for C_4 lattice.

Original entry on oeis.org

1, 32, 192, 608, 1408, 2720, 4672, 7392, 11008, 15648, 21440, 28512, 36992, 47008, 58688, 72160, 87552, 104992, 124608, 146528, 170880, 197792, 227392, 259808, 295168, 333600, 375232, 420192, 468608
Offset: 0

Views

Author

mbaake(AT)sunelc3.tphys.physik.uni-tuebingen.de (Michael Baake)

Keywords

Crossrefs

Cf. A103884 (row 4). For coordination sequences of other C_n lattices see A022144 (C_2), A010006 (C3), A019560 - A019564 (C_4 through C_8), A035746 - A035787 (C_9 through C_50).

Programs

  • Magma
    [1] cat [(32/3)*n*(1 + 2*n^2): n in [1..40]]; // Vincenzo Librandi, Apr 10 2017
  • Mathematica
    Join[{1}, Table[(32/3) n (1 + 2 n^2), {n, 30}]] (* Vincenzo Librandi, Apr 10 2017 *)

Formula

a(n) = (32/3)*n*(1 + 2*n^2) for n>0.
G.f.: (1 + 28*x + 70*x^2 + 28*x^3 + x^4)/(1 - x)^4.
G.f. for sequence with interpolated zeros: cosh(8*arctanh(x)) = 1/2*(((1 + x)/(1 - x))^4 + ((1 - x)/(1 + x))^4) = 1 + 32*x^2 + 192*x^4 + 608*x^6 + .... Cf. A057813. - Peter Bala, Apr 09 2017
a(n) = A008412(2*n). - Seiichi Manyama, Jun 08 2018
Showing 1-10 of 19 results. Next