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

A005900 Octahedral numbers: a(n) = n*(2*n^2 + 1)/3.

Original entry on oeis.org

0, 1, 6, 19, 44, 85, 146, 231, 344, 489, 670, 891, 1156, 1469, 1834, 2255, 2736, 3281, 3894, 4579, 5340, 6181, 7106, 8119, 9224, 10425, 11726, 13131, 14644, 16269, 18010, 19871, 21856, 23969, 26214, 28595, 31116, 33781, 36594, 39559, 42680
Offset: 0

Views

Author

Keywords

Comments

Series reversion of g.f.: A(x) is Sum_{n>0} - A066357(n)(-x)^n.
Partial sums of centered square numbers A001844. - Paul Barry, Jun 26 2003
Also as a(n) = (1/6)*(4n^3 + 2n), n>0: structured tetragonal diamond numbers (vertex structure 5) (cf. A000447 - structured diamonds); and structured trigonal anti-prism numbers (vertex structure 5) (cf. A100185 - structured anti-prisms). Cf. A100145 for more on structured polyhedral numbers. - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Schlaefli symbol for this polyhedron: {3,4}.
If X is an n-set and Y and Z are disjoint 2-subsets of X then a(n-4) is equal to the number of 5-subsets of X intersecting both Y and Z. - Milan Janjic, Aug 26 2007
Starting with 1 = binomial transform of [1, 5, 8, 4, 0, 0, 0, ...] where (1, 5, 8, 4) = row 3 of the Chebyshev triangle A081277. - Gary W. Adamson, Jul 19 2008
a(n) = largest coefficient of (1 + ... + x^(n-1))^4. - R. H. Hardin, Jul 23 2009
Convolution square root of (1 + 6x + 19x^3 + ...) = (1 + 3x + 5x^2 + 7x^3 + ...) = A005408(x). - Gary W. Adamson, Jul 27 2009
Starting with offset 1 = the triangular series convolved with [1, 3, 4, 4, 4, ...]. - Gary W. Adamson, Jul 28 2009
One of the 5 Platonic polyhedral (tetrahedral, cube, octahedral, dodecahedral, and icosahedral) numbers (cf. A053012). - Daniel Forgues, May 14 2010
Let b be any product of four different primes. Then the divisor lattice of b^n is of width a(n+1). - Jean Drabbe, Oct 13 2010
Arises in Bezdek's proof on contact numbers for congruent sphere packings (see preprint). - Jonathan Vos Post, Feb 08 2011
Euler transform of length 2 sequence [6, -2]. - Michael Somos, Mar 27 2011
a(n+1) is the number of 2 X 2 matrices with all terms in {0,1,...,n} and (sum of terms) = 2n. - Clark Kimberling, Mar 19 2012
a(n) is the number of semistandard Young tableaux over all partitions of 3 with maximal element <= n. - Alois P. Heinz, Mar 22 2012
Self convolution of the odd numbers. - Reinhard Zumkeller, Apr 04 2012
a(n) is the number of (w,x,y,z) with all terms in {1,...,n} and w+x=y+z; also the number of (w,x,y,z) with all terms in {0,...,n} and |w-x|<=y. - Clark Kimberling, Jun 02 2012
The sequence is the third partial sum of (0, 1, 3, 4, 4, 4, ...). - Gary W. Adamson, Sep 11 2015
a(n) is the number of join-irreducible elements in the Weyl group of type B_n with respect to the strong Bruhat order. - Rafael Mrden, Aug 26 2020
Number of unit octahedra contained in an n-scale octahedron composed of a tetrahedral-octahedral honeycomb. The number of unit tetrahedra in it is 8*A000292(n-1) = 4*(n^3 - n)/3. Also, the number of unit tetrahedra and unit octahedra contained in an n-scale tetrahedron composed of a tetrahedral-octahedral honeycomb is respectively A006527(n) = (n^3 + 2*n)/3 and A000292(n-1) = (n^3 - n)/6. - Jianing Song, Feb 24 2025

Examples

			G.f. = x + 6*x^2 + 19*x^3 + 44*x^4 + 85*x^5 + 146*x^6 + 231*x^7 + ...
		

References

  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 50.
  • H. S. M. Coxeter, Polyhedral numbers, pp. 25-35 of R. S. Cohen, J. J. Stachel and M. W. Wartofsky, eds., For Dirk Struik: Scientific, historical and political essays in honor of Dirk J. Struik, Reidel, Dordrecht, 1974.
  • 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 A001845. Cf. A001844.
1/12*t*(n^3-n)+n for t = 2, 4, 6, ... gives A004006, A006527, A006003, A005900, A004068, A000578, A004126, A000447, A004188, A004466, A004467, A007588, A062025, A063521, A063522, A063523.
Cf. A022521.
Cf. A081277.
Row n=3 of A210391. - Alois P. Heinz, Mar 22 2012
Cf. A005408.
Cf. A002061.
Cf. A000292 (tetrahedral numbers), A000578 (cubes), A006566 (dodecahedral numbers), A006564 (icosahedral numbers).
Similar sequence: A014820(n-1) (m=4), A069038 (m=5), A069039 (m=6), A099193(m=7), A099195 (m=8), A099196 (m=9), A099197 (m=10).

Programs

  • Haskell
    a005900 n = sum $ zipWith (*) odds $ reverse odds
                where odds = take n a005408_list
    a005900_list = scanl (+) 0 a001844_list
    -- Reinhard Zumkeller, Jun 16 2013, Apr 04 2012
    
  • Magma
    [n*(2*n^2+1)/3: n in [0..50]]; // Wesley Ivan Hurt, Sep 11 2015
    
  • Magma
    I:=[0,1,6,19]; [n le 4 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, Sep 12 2015
    
  • Maple
    al:=proc(s,n) binomial(n+s-1,s); end; be:=proc(d,n) local r; add( (-1)^r*binomial(d-1,r)*2^(d-1-r)*al(d-r,n), r=0..d-1); end; [seq(be(3,n), n=0..100)];
    A005900:=(z+1)**2/(z-1)**4; # Simon Plouffe in his 1992 dissertation
    with(combinat): seq(fibonacci(4,2*n)/12, n=0..40); # Zerinvary Lajos, Apr 21 2008
  • Mathematica
    Table[(2n^3+n)/3, {n,0,40}] (* or *) LinearRecurrence[{4,-6,4,-1}, {0,1,6,19},50] (* Harvey P. Dale, Oct 10 2013 *)
    CoefficientList[Series[x (1 + x)^2/(1 - x)^4, {x, 0, 45}], x] (* Vincenzo Librandi, Sep 12 2015 *)
  • Maxima
    makelist(n*(2*n^2+1)/3, n, 0, 20); /* Martin Ettl, Jan 07 2013 */
    
  • PARI
    {a(n) = n*(2*n^2+1)/3};
    
  • PARI
    concat([0],Vec(x*(1 + x)^2/(1 - x)^4 + O(x^50))) \\ Indranil Ghosh, Mar 16 2017
    
  • Python
    def a(n): return n*(2*n*n + 1)//3
    print([a(n) for n in range(41)]) # Michael S. Branicky, Sep 03 2021

Formula

a(n) = 1^2 + 2^2 + ... + (n-1)^2 + n^2 + (n-1)^2 + ... + 2^2 + 1^2. - Amarnath Murthy, May 28 2001
G.f.: x * (1 + x)^2 / (1 - x)^4. a(n) = -a(-n) = (2*n^3 + n) / 3.
a(n) = ( ((n+1)^5-n^5) - (n^5-(n-1)^5) )/30. - Xavier Acloque, Oct 17 2003
a(n) is the sum of the products pq, where p and q are both positive and odd and p + q = 2n, e.g., a(4) = 7*1 + 5*3 + 3*5 + 1*7 = 44. - Jon Perry, May 17 2005
a(n) = 4*binomial(n,3) + 4*binomial(n,2) + binomial(n,1). - Mitch Harris, Jul 06 2006
a(n) = binomial(n+2,3) + 2*binomial(n+1,3) + binomial(n,3), (this pair generalizes; see A014820, the 4-cross polytope numbers).
Sum_{n>=1} 1/a(n) = 3*gamma + 3*Psi((I*(1/2))*sqrt(2)) - (1/2)*(3*I)*Pi*coth((1/2)*Pi*sqrt(2)) - (1/2)*(3*I)*sqrt(2) = A175577, where I=sqrt(-1). - Stephen Crowley, Jul 14 2009
a(n) = A035597(n)/2. - J. M. Bergot, Jun 11 2012
a(n) = A000578(n) - 2*A000292(n-1) for n>0. - J. M. Bergot, Apr 05 2014
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4), n>3. - Wesley Ivan Hurt, Sep 11 2015
E.g.f.: (1/3)*x*(3 + 6*x + 2*x^2)*exp(x). - Ilya Gutkovskiy, Mar 16 2017
a(n) = (A002061(A002061(n+1)) - A002061(A002061(n)))/6. - Daniel Poveda Parrilla, Jun 10 2017
a(n) = 6*a(n-1)/(n-1) + a(n-2) for n > 1. - Seiichi Manyama, Jun 06 2018
Sum_{n >= 1} (-1)^(n+1)/(a(n)*a(n+1)) = 6*log(2) - 4 = 1/(6 + 2/(6 + 6/(6 + ... + n*(n-1)/(6 + ...)))). See A142983. - Peter Bala, Mar 06 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

A006331 a(n) = n*(n+1)*(2*n+1)/3.

Original entry on oeis.org

0, 2, 10, 28, 60, 110, 182, 280, 408, 570, 770, 1012, 1300, 1638, 2030, 2480, 2992, 3570, 4218, 4940, 5740, 6622, 7590, 8648, 9800, 11050, 12402, 13860, 15428, 17110, 18910, 20832, 22880, 25058, 27370, 29820, 32412, 35150, 38038, 41080, 44280
Offset: 0

Views

Author

Keywords

Comments

Triangles in rhombic matchstick arrangement of side n.
Maximum accumulated number of electrons at energy level n. - Scott A. Brown, Feb 28 2000
Let M_n denote the n X n matrix M_n(i,j)=i^2+j^2; then the characteristic polynomial of M_n is x^n - a(n)x^(n-1) - .... - Michael Somos, Nov 14 2002
Convolution of odds (A005408) and evens (A005843). - Graeme McRae, Jun 06 2006
a(n) is the number of non-monotonic functions with domain {0,1,2} and codomain {0,1,...,n}. - Dennis P. Walsh, Apr 25 2011
For any odd number 2n+1, find Sum_{aJ. M. Bergot, Jul 16 2011
a(n) gives the number of (n+1) X (n+1) symmetric (0,1)-matrices containing three ones (see [Cameron]). - L. Edson Jeffery, Feb 18 2012
a(n) is the number of 4-tuples (w,x,y,z) with all terms in {0,...,n} and |w - x| < y. - Clark Kimberling, Jun 02 2012
Partial sums of A001105. - Omar E. Pol, Jan 12 2013
Total number of square diagonals (of any size) in an n X n square grid. - Wesley Ivan Hurt, Mar 24 2015
Number of diagonal attacks of two queens on (n+1) X (n+1) chessboard. - Antal Pinter, Sep 20 2015
a(n) is the minimum value obtainable by partitioning either the set {x in the natural numbers | 1 <= x <= 2n} or the set {x in the natural numbers | 0 <= x <= 2n+1} into pairs, taking the product of all such pairs, and taking the sum of all such products. - Thomas Anton, Oct 21 2020
a(n) is the irregularity of the n-th power of a path of length at least 3*n. (The irregularity of a graph is the sum of the differences between the degrees over all edges of the graph.) - Allan Bickle, Jun 16 2023
a(n) is the maximum possible total number of inversions in all rows and all columns of a Latin square of order n+1. - Ivaylo Kortezov, Jun 28 2025

Examples

			For n=2, a(2)=10 since there are 10 non-monotonic functions f from {0,1,2} to {0,1,2}, namely, functions f = <f(1),f(2),f(3)> given by <0,1,0>, <0,2,0>, <0,2,1>, <1,0,1>, <1,0,2>, <1,2,0>, <1,2,1>, <2,0,1>, <2,0,2>, and <2,1,2>. - _Dennis P. Walsh_, Apr 25 2011
Let n=4, 2*n+1 = 9. Since 9 = 1+8 = 3+6 = 5+4 = 7+2, a(4) = 1*8 + 3*6 + 5*4 + 7*2 = 60. - _Vladimir Shevelev_, May 11 2012
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

A row of A132339.
Cf. A002378, A046092, A028896 (irregularities of maximal k-degenerate graphs).

Programs

  • Haskell
    a006331 n = sum $ zipWith (*) [2*n-1, 2*n-3 .. 1] [2, 4 ..]
    -- Reinhard Zumkeller, Feb 11 2012
  • Magma
    [n*(n+1)*(2*n+1)/3: n in [0..40]]; // Vincenzo Librandi, Aug 15 2011
    
  • Maple
    A006331 := proc(n)
        n*(n+1)*(2*n+1)/3 ;
    end proc:
    seq(A006331(n),n=0..80) ; # R. J. Mathar, Sep 27 2013
  • Mathematica
    Table[n(n+1)(2n+1)/3,{n,0,40}] (* or *) LinearRecurrence[{4,-6,4,-1},{0,2,10,28},50] (* Harvey P. Dale, Apr 12 2013 *)
  • PARI
    a(n)=if(n<0,0,n*(n+1)*(2*n+1)/3)
    

Formula

G.f.: 2*x*(1 + x)/(1 - x)^4. - Simon Plouffe (in his 1992 dissertation)
a(n) = 2*binomial(n+1,3) + 2*binomial(n+2,3).
a(n) = 2*A000330(n) = A002492(n)/2.
a(n) = Sum_{i=0..n} T(i,n-i), array T as in A048147. - N. J. A. Sloane, Dec 11 1999
From the formula for the sum of squares of positive integers 1^2 + 2^2 + 3^2 + ... + n^2 = n*(n+1)(2*n+1)/6, if we multiply both sides by 2 we get Sum_{k=0..n} 2*k^2 = n*(n+1)*(2*n+1)/3, which is an alternative formula for this sequence. - Mike Warburton, Sep 08 2007
10*a(n) = A016755(n) - A001845(n); since A016755 are odd cubes and A001845 centered octahedral numbers, 10*a(n) are the "odd cubes without their octahedral contents." - Damien Pras, Mar 19 2011
a(n) = sum(a*b), where the summing is over all unordered partitions 2*n+1=a+b. - Vladimir Shevelev, May 11 2012
a(n) = binomial(2*n+2, 3)/2. - Ronan Flatley, Dec 13 2012
a(n) = A000292(n) + A002411(n). - Omar E. Pol, Jan 11 2013
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) for n>3, with a(0)=0, a(1)=2, a(2)=10, a(3)=28. - Harvey P. Dale, Apr 12 2013
a(n) = A208532(n+1,2). - Philippe Deléham, Dec 05 2013
Sum_{n>0} 1/a(n) = 9 - 12*log(2). - Enrique Pérez Herrero, Dec 03 2014
a(n) = A000292(n-1) + (n+1)*A000217(n). - J. M. Bergot, Sep 02 2015
a(n) = 2*(A000332(n+3) - A000332(n+1)). - Antal Pinter, Sep 20 2015
From Bruno Berselli, May 17 2018: (Start)
a(n) = n*A002378(n) - Sum_{k=0..n-1} A002378(k) for n>0, a(0)=0. Also:
A163102(n) = n*a(n) - Sum_{k=0..n-1} a(k) for n>0, A163102(0)=0. (End)
a(n) = A005900(n) - A000290(n) = A096000(n) - A000578(n+1) = A000578(n+1) - A084980(n+1) = A000578(n+1) - A077415(n)-1 = A112524(n) + 1 = A188475(n) - 1 = A061317(n) - A100178(n) = A035597(n+1) - A006331(n+1). - Bruce J. Nicholson, Jun 24 2018
E.g.f.: (1/3)*exp(x)*x*(6 + 9*x + 2*x^2). - Stefano Spezia, Jan 05 2020
Sum_{n>=1} (-1)^(n+1)/a(n) = 3*Pi - 9. - Amiram Eldar, Jan 04 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)

A168380 Row sums of A168281.

Original entry on oeis.org

2, 4, 12, 20, 38, 56, 88, 120, 170, 220, 292, 364, 462, 560, 688, 816, 978, 1140, 1340, 1540, 1782, 2024, 2312, 2600, 2938, 3276, 3668, 4060, 4510, 4960, 5472, 5984, 6562, 7140, 7788, 8436, 9158, 9880, 10680, 11480, 12362, 13244, 14212, 15180, 16238, 17296, 18448, 19600, 20850, 22100
Offset: 1

Views

Author

Paul Curtz, Nov 24 2009

Keywords

Comments

The atomic numbers of the augmented alkaline earth group in Charles Janet's spiral periodic table are 0 and the first eight terms of this sequence (see Stewart reference). - Alonso del Arte, May 13 2011
Maximum number of 123 patterns in an alternating permutation of length n+3. - Lara Pudwell, Jun 09 2019

Examples

			From _Lara Pudwell_, Jun 09 2019: (Start)
a(1)=2. The alternating permutation of length 1+3=4 with the maximum number of copies of 123 is 1324.  The two copies are 124 and 134.
a(2)=4.  The alternating permutation of length 2+3=5 with the maximum number of copies of 123 is 13254.  The four copies are 124, 125, 134, and 135.
a(3)=12. The alternating permutation of length 3+3=6 with the maximum number of copies of 123 is 132546.  The twelve copies are 124, 125, 126, 134, 135, 136, 146, 156, 246, 256, 346, and 356. (End)
		

Programs

  • Magma
    [(n+1)*(3+2*n^2+4*n-3*(-1)^n)/12: n in [1..50] ]; // Vincenzo Librandi, Aug 06 2011
    
  • Mathematica
    LinearRecurrence[{2,1,-4,1,2,-1},{2, 4, 12, 20, 38, 56},50] (* G. C. Greubel, Jul 19 2016 *)
    Table[(n + 1) (3 + 2 n^2 + 4 n - 3 (-1)^n)/12, {n, 50}] (* Michael De Vlieger, Jul 20 2016 *)
  • PARI
    a(n)=([0,1,0,0,0,0; 0,0,1,0,0,0; 0,0,0,1,0,0; 0,0,0,0,1,0; 0,0,0,0,0,1; -1,2,1,-4,1,2]^(n-1)*[2;4;12;20;38;56])[1,1] \\ Charles R Greathouse IV, Jul 21 2016

Formula

a(n) = 2*A005993(n-1).
a(n) = (n+1)*(3 + 2*n^2 + 4*n - 3*(-1)^n)/12.
a(n+1) - a(n) = A093907(n) = A137583(n+1).
a(2n+1) = A035597(n+1), a(2n) = A002492(n).
a(n) = A099956(n-1), 2 <= n <= 7.
a(n) = 2*a(n-1) + a(n-2) - 4*a(n-3) + a(n-4) + 2*a(n-5) - a(n-6).
G.f.: 2*x*(1 + x^2) / ( (1+x)^2*(x-1)^4 ).
a(n) = A000292(n) + A027656(n-1). - Paul Curtz, Oct 26 2012
E.g.f.: (1/12)*(3*(x - 1) + (3 + 15*x + 12*x^2 + 2*x^3)*exp(2*x))*exp(-x). - G. C. Greubel, Jul 19 2016

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

A201780 Riordan array ((1-x)^2/(1-2x), x/(1-2x)).

Original entry on oeis.org

1, 0, 1, 1, 2, 1, 2, 5, 4, 1, 4, 12, 13, 6, 1, 8, 28, 38, 25, 8, 1, 16, 64, 104, 88, 41, 10, 1, 32, 144, 272, 280, 170, 61, 12, 1, 64, 320, 688, 832, 620, 292, 85, 14, 1, 128, 704, 1696, 2352, 2072, 1204, 462, 113, 16, 1
Offset: 0

Views

Author

Philippe Deléham, Dec 05 2011

Keywords

Comments

Diagonals ascending: 1, 0, 1, 1, 2, 2, 4, 5, 1, 8, 12, 4, ... (see A201509).

Examples

			Triangle begins:
  1;
  0,  1;
  1,  2,  1;
  2,  5,  4,  1;
  4, 12, 13,  6,  1;
  8, 28, 38, 25,  8,  1;
		

Crossrefs

Row sums: A052156

Programs

  • Mathematica
    CoefficientList[#, y]& /@ CoefficientList[(1-x)^2/(1-(y+2)*x) + O[x]^10, x] // Flatten (* Jean-François Alcover, Nov 03 2018 *)

Formula

T(n,k) = 2*T(n-1,k) + T(n-1,k-1) with T(0,0) = 0, T(1,0) = 0, T(2,0) = 0 and T(n,k)= 0 if k < 0 or if n < k.
Sum_{k=0..n} T(n,k)*x^k = A154955(n+1), A034008(n), A052156(n), A055841(n), A055842(n), A055846(n), A055270(n), A055847(n), A055995(n), A055996(n), A056002(n), A056116(n) for x = -1,0,1,2,3,4,5,6,7,8,9,10 respectively.
G.f.: (1-x)^2/(1-(y+2)*x).

A035598 Number of points of L1 norm 4 in cubic lattice Z^n.

Original entry on oeis.org

0, 2, 16, 66, 192, 450, 912, 1666, 2816, 4482, 6800, 9922, 14016, 19266, 25872, 34050, 44032, 56066, 70416, 87362, 107200, 130242, 156816, 187266, 221952, 261250, 305552, 355266, 410816, 472642, 541200, 616962, 700416, 792066
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [( 2*n^4 +4*n^2 )/3: n in [0..40]]; // Vincenzo Librandi, Apr 22 2012
  • Maple
    f := proc(d,m) local i; sum( 2^i*binomial(d,i)*binomial(m-1,i-1),i=1..min(d,m)); end; # n=dimension, m=norm
  • Mathematica
    CoefficientList[Series[2*x*(1+x)^3/(1-x)^5,{x,0,40}],x] (* Vincenzo Librandi, Apr 22 2012 *)
    LinearRecurrence[{5,-10,10,-5,1},{0,2,16,66,192},50] (* Harvey P. Dale, Dec 11 2019 *)
  • PARI
    a(n)=2*n^2*(n^2+2)/3 \\ Charles R Greathouse IV, Dec 07 2011
    

Formula

a(n) = 2*n^2*(n^2 + 2)/3. - Frank Ellermann, Mar 16 2002
G.f.: 2*x*(1+x)^3/(1-x)^5. - Colin Barker, Apr 15 2012
a(n) = 2*A014820(n-1). - R. J. Mathar, Dec 10 2013
a(n) = a(n-1) + A035597(n) + A035597(n-1). - Bruce J. Nicholson, Mar 11 2018
From Shel Kaphan, Feb 28 2023: (Start)
a(n) = 2*n*Hypergeometric2F1(1-n,1-k,2,2), where k=4.
a(n) = A001846(n) - A001845(n).
a(n) = A008412(n)*n/4. (End)
From Amiram Eldar, Mar 12 2023: (Start)
Sum_{n>=1} 1/a(n) = Pi^2/8 - 3*Pi*coth(sqrt(2)*Pi)/(8*sqrt(2)) + 3/16.
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/16 + 3*Pi*cosech(sqrt(2)*Pi)/(8*sqrt(2)) - 3/16. (End)
E.g.f.: 2*exp(x)*x*(3 + 9*x + 6*x^2 + x^3)/3. - Stefano Spezia, Mar 14 2024

A168234 A138100(n) + A168142(n).

Original entry on oeis.org

3, 3, 13, 13, 13, 13, 13, 13, 13, 13, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 171, 171, 171, 171, 171, 171, 171
Offset: 1

Views

Author

Paul Curtz, Nov 21 2009

Keywords

Comments

A138100 counts upwards in blocks of 2*k^2 numbers, restarting from 1,5,21,57,.. = A166464(k-1) = (2*k+1)*(2*k^2-4*k+3)/3, k>=1.
A168142 counts downwards in blocks of 2*k^2 numbers, restarting from 2*k^2, k>=1.
In consequence, the sequence here contains 2*k^2 copies of the number 1+2*k*(1+2*k^2)/3 = 1+A035597(k), k>=1,
where the sequence A035597 is a bisection of A168380.

Crossrefs

Extensions

Edited by R. J. Mathar, Feb 15 2010

A168547 a(n) = 1 - 2*n^2 + 4*n*(1 + 2*n^2)/3.

Original entry on oeis.org

1, 3, 17, 59, 145, 291, 513, 827, 1249, 1795, 2481, 3323, 4337, 5539, 6945, 8571, 10433, 12547, 14929, 17595, 20561, 23843, 27457, 31419, 35745, 40451, 45553, 51067, 57009, 63395, 70241, 77563, 85377, 93699, 102545, 111931, 121873, 132387, 143489, 155195
Offset: 0

Views

Author

Paul Curtz, Nov 29 2009

Keywords

Comments

Binomial transform of the quasi-finite sequence 1,2,12,16,0,... (0 continued).
A bisection of A168582.

Programs

Formula

a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4).
G.f.: (1 - x + 11*x^2 + 5*x^3)/(x-1)^4.
First differences: a(n+1) - a(n) = 2*A054569(n+1).
Second differences: a(n+2) - 2*a(n+1) + a(n) = 4*A004767(n).
Third differences: a(n+3) - 3*a(n+2) + 3*a(n+1) - a(n) = 16.
a(n+1) = A166464(n) + A035597(n+1).
a(n) = 1 - 2*n^2 + 4*A005900(n). - R. J. Mathar, Dec 05 2009
E.g.f.: (1/3)*(3 + 6*x + 18*x^2 + 8*x^3)*exp(x). - G. C. Greubel, Jul 26 2016

Extensions

Edited and extended by R. J. Mathar, Dec 05 2009
Showing 1-10 of 12 results. Next