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.

Previous Showing 11-20 of 83 results. Next

A100145 Structured great rhombicosidodecahedral numbers.

Original entry on oeis.org

1, 120, 579, 1600, 3405, 6216, 10255, 15744, 22905, 31960, 43131, 56640, 72709, 91560, 113415, 138496, 167025, 199224, 235315, 275520, 320061, 369160, 423039, 481920, 546025, 615576, 690795, 771904, 859125, 952680, 1052791, 1159680
Offset: 1

Views

Author

James A. Record (james.record(AT)gmail.com), Nov 07 2004

Keywords

Comments

Structured polyhedral numbers are a type of figurate polyhedral numbers. Structurate polyhedra differ from regular figurate polyhedra by having appropriate figurate polygonal faces at any iteration, i.e., a regular truncated octahedron, n=2, would have 7 points on its hexagonal faces, whereas a structured truncated octahedron, n=2, would have 6 points - just as a hexagon, n=2, would have. Like regular figurate polygons, structured polyhedra seem to originate at a vertex and since many polyhedra have different vertices (a pentagonal diamond has 2 "polar" vertices with 5 adjacent vertices and 5 "equatorial" vertices with 4 adjacent vertices), these polyhedra have multiple structured number sequences, dependent on the "vertex structures" which are each equal to the one vertex itself plus its adjacent vertices. For polystructurate polyhedra the notation, structured polyhedra (vertex structure x) is used to differentiate between alternate vertices, where VS stands for vertex structure.

Crossrefs

Cf. A051673, A100146 through A100156 - structured Archimedean solids; A100157 through A100175 - structured Catalan solids; A100147 - structured prisms; A000447 - structured diamonds; A100185 - structured anti-prisms; and A100188 - structured anti-diamonds.

Programs

Formula

a(n) = (1/6)*(222*n^3 - 312*n^2 + 96*n).
From Jaume Oliver Lafont, Sep 08 2009: (Start)
a(n) = (1+(n-1))*(1+22*(n-1)+37*(n-1)^2);
G.f.: x*(1+116*x+105*x^2)/(1-x)^4. (End)
E.g.f.: exp(x)*x*(1 + 59*x + 37*x^2). - Stefano Spezia, Jun 06 2025

Extensions

Corrected by T. D. Noe, Oct 25 2006

A006527 a(n) = (n^3 + 2*n)/3.

Original entry on oeis.org

0, 1, 4, 11, 24, 45, 76, 119, 176, 249, 340, 451, 584, 741, 924, 1135, 1376, 1649, 1956, 2299, 2680, 3101, 3564, 4071, 4624, 5225, 5876, 6579, 7336, 8149, 9020, 9951, 10944, 12001, 13124, 14315, 15576, 16909, 18316, 19799, 21360, 23001, 24724, 26531, 28424, 30405
Offset: 0

Views

Author

Keywords

Comments

Number of ways to color vertices (or edges) of a triangle using <= n colors, allowing only rotations.
Also: dot_product (1,2,...,n)*(2,3,...,n,1), n >= 0. - Clark Kimberling
Start from triacid and attach amino acids according to the reaction scheme that describes the reaction between the active sites. See the hyperlink below on chemistry. - Robert G. Wilson v, Aug 02 2002
Starting with offset 1 = row sums of triangle A158822 and binomial transform of (1, 3, 4, 2, 0, 0, 0, ...). - Gary W. Adamson, Mar 28 2009
One-ninth of sum of three consecutive cubes: a(n) = ((n-1)^3 + n^3 + (n+1)^3)/9. - Zak Seidov, Jul 22 2013
For n > 2, number of different cubes, formed after splitting a cube in color C_1, by parallel planes in the colors C_2, C_3, ..., C_n in three spatial dimensions (in the order of the colors from a fixed vertex). Generally, in a large hypercube n^d is f(n,d) = C(n+d-1, d) + C(n, d) different small hypercubes. See below for my formula a(n) = f(n,3). - Thomas Ordowski, Jun 15 2014
a(n) is a square for n = 1, 2 & 24; and for no other values up to 10^7 (see M. Gardner). - Michel Marcus, Sep 06 2015
Number of unit tetrahedra contained in an n-scale tetrahedron composed of a tetrahedral-octahedral honeycomb. - Jason Pruski, Aug 23 2017

References

  • M. Gardner, New Mathematical Diversions from Scientific American. Simon and Schuster, NY, 1966, p. 246.
  • S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see p. 483.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

(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.
Column 1 of triangle A094414. Row 6 of the array in A107735.
Cf. A000292 (unoriented), A000292(n-2) (chiral), A000290 (achiral) triangle colorings.
Row 2 of A324999 (simplex vertices and facets) and A327083 (simplex edges and ridges).

Programs

  • Haskell
    a006527 n = n * (n ^ 2 + 2) `div` 3  -- Reinhard Zumkeller, Jan 06 2014
  • Magma
    [(n^3 + 2*n)/3: n in [0..50]]; // Vincenzo Librandi, May 15 2011
    
  • Maple
    A006527:=z*(1+z**2)/(z-1)**4; # conjectured by Simon Plouffe in his 1992 dissertation
    with(combinat):seq(lcm(fibonacci(4,n),fibonacci(2,n))/3,n=0..42); # Zerinvary Lajos, Apr 20 2008
  • Mathematica
    Table[ (n^3 + 2*n)/3, {n, 0, 45} ]
    LinearRecurrence[{4,-6,4,-1},{0,1,4,11},46] (* or *) CoefficientList[ Series[(x+x^3)/(x-1)^4,{x,0,49}],x] (* Harvey P. Dale, Jun 13 2011 *)
  • PARI
    a(n)=n*(n^2+2)/3 \\ Charles R Greathouse IV, Jul 25 2011
    

Formula

a(0)=0, a(1)=1, a(2)=4, a(3)=11; for n > 3, a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Harvey P. Dale, Jun 13 2011
From Paul Barry, Mar 13 2003: (Start)
a(n) = 2*binomial(n+1, 3) + binomial(n, 1).
G.f.: x*(1+x^2)/(1-x)^4. (End)
a(n) = A000292(n) + A000292(n-2). - Alexander Adamchuk, May 20 2006
a(n) = n*A059100(n)/3. - Lekraj Beedassy, Feb 06 2007
a(n) = A054602(n)/3. - Zerinvary Lajos, Apr 20 2008
a(n) = ( n + Sum_{i=1..n} A177342(i) )/(n+1), with n > 0. - Bruno Berselli, May 19 2010
a(n) = A002264(A000578(n) + A005843(n)). - Reinhard Zumkeller, Jun 16 2011
a(n) = binomial(n+2, 3) + binomial(n, 3). - Thomas Ordowski, Jun 15 2014
a(n) = A000292(n) - A000292(-n). - Bruno Berselli, Sep 22 2016
E.g.f.: (x/3)*(3 + 3*x + x^2)*exp(x). - G. C. Greubel, Sep 01 2017
From Robert A. Russell, Oct 20 2020: (Start)
a(n) = 1*C(n,1) + 2*C(n,2) + 2*C(n,3), where the coefficient of C(n,k) is the number of oriented triangle colorings using exactly k colors.
a(n) = 2*A000292(n) - A000290(n) = 2*A000292(n-2) + A000290(n). (End)
Sum_{n>0} 1/a(n) = 3*(2*gamma + polygamma(0, 1-i*sqrt(2)) + polygamma(0, 1+i*sqrt(2)))/4 = 1.45245201414472469745354677573358867... where i denotes the imaginary unit. - Stefano Spezia, Aug 31 2023

Extensions

More terms from Alexander Adamchuk, May 20 2006
Corrected and replaced 5th formula from Harvey P. Dale, Jun 13 2011
Deleted an erroneous comment. - N. J. A. Sloane, Dec 10 2018

A002577 Number of partitions of 2^n into powers of 2.

Original entry on oeis.org

1, 2, 4, 10, 36, 202, 1828, 27338, 692004, 30251722, 2320518948, 316359580362, 77477180493604, 34394869942983370, 27893897106768940836, 41603705003444309596874, 114788185359199234852802340, 588880400923055731115178072778, 5642645813427132737155703265972004
Offset: 0

Views

Author

Keywords

Comments

For given m, the general formula for t_m(n, k) and the corresponding tables T, computed as in the example, determine a family of related sequences (placed in the rows or in the columns of T). For example, the numbers from the second row of T, computed for given m and n > 2, are the (m+2)-gonal numbers. So the second row contains the first members of: A000290 (the square numbers) when m=2, A000326 (the pentagonal numbers) when m=3, and so on. But rows IV, V etc. of the given table are not represented in the OEIS till now. - Valentin Bakoev, Feb 25 2009; edited by M. F. Hasler, Feb 09 2014

Examples

			To compute t_2(6,1) we can use a table T, defined as T[i,j]= t_2(i,j), for i=1,2,...,6(=n), and j= 0,1,2,...,32(= k*m^{n-1}). It is: 1,2,3,4,5,6,7,8,9...,33; 1,4,9,16,25,36,49...,81; (so the second row contains the first members of A000290 -- the square numbers) 1,10,35,84,165,...,969; (so the third row contains the first members of A000447. The r-th tetrahedral number is given by formula r(r+1)(r+2)/6. This row (also A000447) contains the tetrahedral numbers, obtained for r=1,3,5,7,...) 1,36,201,656,1625; 1,202,1827; 1,1828; Column 1 contains the first 6 members of A002577. - _Valentin Bakoev_, Feb 25 2009
G.f. = 1 + 2*x + 4*x^2 + 10*x^3 + 36*x^4 + 202*x^5 + 1828*x^6 + ...
		

References

  • R. F. Churchhouse, Binary partitions, pp. 397-400 of A. O. L. Atkin and B. J. Birch, editors, Computers in Number Theory. Academic Press, NY, 1971.
  • Lawrence, Jim. "Dual-Antiprisms and Partitions of Powers of 2 into Powers of 2." Discrete & Computational Geometry, Vol. 16 (2019): 465-478. See page 466.
  • 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

a(n) = A000123(2^(n-1)) = A018818(2^n).
Column k=2 of A145515, diagonal of A152977. - Alois P. Heinz, Mar 25 2012
See also A002575, A002576.
A column of A125790.

Programs

  • Haskell
    import Data.MemoCombinators (memo2, list, integral)
    a002577 n = a002577_list !! n
    a002577_list = f [1] where
       f xs = (p' xs $ last xs) : f (1 : map (* 2) xs)
       p' = memo2 (list integral) integral p
       p  0 = 1; p []  = 0
       p ks'@(k:ks) m = if m < k then 0 else p' ks' (m - k) + p' ks m
    -- Reinhard Zumkeller, Nov 27 2015
  • Maple
    A002577 := proc(n) if n<=1 then n+1 else A000123(2^(n-1)); fi; end;
  • Mathematica
    $RecursionLimit = 10^5; (* b = A000123 *) b[0] = 1; b[n_?EvenQ] := b[n] = b[n-1] + b[n/2]; b[n_?OddQ] := b[n] = b[n-1] + b[(n-1)/2]; a[n_] := b[2^(n-1)]; a[0] = 1; Table[a[n], {n, 0, 17}] (* Jean-François Alcover, Nov 23 2011 *)
    a[ n_] := SeriesCoefficient[ 1 / Product[ 1 - x^2^k, {k, 0, n}], {x, 0, 2^n}]; (* Michael Somos, Apr 21 2014 *)
  • PARI
    a(n)=polcoeff(prod(j=0,n,1/(1-x^(2^j)+x*O(x^(2^n)))),2^n) \\ Paul D. Hanna
    

Formula

a(n) is about 0.9233*Sum_j {i=0, 1, 2, 3, ...} 2^(j*(2n-j-1)/2)/j!. - Henry Bottomley, Jul 23 2003
a(n) = A078121(n+1, 1). - Paul D. Hanna, Sep 13 2004
A002577(n)-1 = A125792(n). - Let m > 1, n > 0 and k >= 0. The general formula for the number of all partitions of k*m^n into powers of m is t_m(n, k)= k+1 if n=1, t_m(n, k)= 1 if k=0, and t_m(n, k)= t_m(n, k-1) + t_m(n-1, k*m) if n > 1 and k > 0. A002577 is obtained for m=2 and n=1,2,3,... - Valentin Bakoev, Feb 25 2009
a(n) = [x^(2^n)] 1/Product_{j>=0} (1-x^(2^j)). - Alois P. Heinz, Sep 27 2011

Extensions

Edited by M. F. Hasler, Feb 09 2014

A007588 Stella octangula numbers: a(n) = n*(2*n^2 - 1).

Original entry on oeis.org

0, 1, 14, 51, 124, 245, 426, 679, 1016, 1449, 1990, 2651, 3444, 4381, 5474, 6735, 8176, 9809, 11646, 13699, 15980, 18501, 21274, 24311, 27624, 31225, 35126, 39339, 43876, 48749, 53970, 59551, 65504, 71841, 78574, 85715, 93276, 101269, 109706, 118599, 127960
Offset: 0

Views

Author

Keywords

Comments

Also as a(n)=(1/6)*(12*n^3-6*n), n>0: structured hexagonal anti-diamond numbers (vertex structure 13) (Cf. A005915 = alternate vertex; A100188 = structured anti-diamonds; A100145 for more on structured numbers). - James A. Record (james.record(AT)gmail.com), Nov 07 2004
The only known square stella octangula number for n>1 is a(169) = 169*(2*169^2 - 1) = 9653449 = 3107^2. - Alexander Adamchuk, Jun 02 2008
Ljunggren proved that 9653449 = (13*239)^2 is the only square stella octangula number for n>1. See A229384 and the Wikipedia link. - Jonathan Sondow, Sep 30 2013
4*A007588 = A144138(ChebyshevU[3,n]). - Vladimir Joseph Stephan Orlovsky, Jun 30 2011
If A016813 is regarded as a regular triangle (with leading terms listed in A001844), a(n) provides the row sums of this triangle: 1, 5+9=14, 13+17+21=51 and so on. - J. M. Bergot, Jul 05 2013
Shares its digital root, A267017, with n*(n^2 + 1)/2 ("sum of the next n natural numbers" see A006003). - Peter M. Chema, Aug 28 2016

References

  • J. H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, p. 51.
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 140.
  • W. Ljunggren, Zur Theorie der Gleichung x^2 + 1 = Dy^4, Avh. Norske Vid. Akad. Oslo. I. 1942 (5): 27.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Backwards differences give star numbers A003154: A003154(n)=a(n)-a(n-1).
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. A001653 = Numbers n such that 2*n^2 - 1 is a square.
a(169) = (A229384(3)*A229384(4))^2.

Programs

Formula

G.f.: x*(1+10*x+x^2)/(1-x)^4.
a(n) = n*A056220(n).
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4), n>3. - Harvey P. Dale, Sep 16 2011
From Ilya Gutkovskiy, Jul 02 2016: (Start)
E.g.f.: x*(1 + 6*x + 2*x^2)*exp(x).
Dirichlet g.f.: 2*zeta(s-3) - zeta(s-1). (End)
a(n) = A004188(n) + A135503(n). - Miquel Cerda, Dec 25 2016
a(n) = A061317(n) - A005843(n) = A062392(n) - A062392(n-1). - J.S. Seneschal, Jul 01 2025

Extensions

In the formula given in the 1995 Encyclopedia of Integer Sequences, the second 2 should be an exponent.

A004068 Number of atoms in a decahedron with n shells.

Original entry on oeis.org

0, 1, 7, 23, 54, 105, 181, 287, 428, 609, 835, 1111, 1442, 1833, 2289, 2815, 3416, 4097, 4863, 5719, 6670, 7721, 8877, 10143, 11524, 13025, 14651, 16407, 18298, 20329, 22505, 24831, 27312, 29953, 32759, 35735, 38886, 42217, 45733, 49439
Offset: 0

Views

Author

Albert D. Rich (Albert_Rich(AT)msn.com)

Keywords

Comments

Also as a(n)=(n/6)*(5*n^2+1), n>0: structured pentagonal diamond numbers (vertex structure 6) (cf. A081436 = alternate vertex; A000447 = structured diamonds; A100145 for more on structured numbers). - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Number of atoms in decahedron with n shells, number = 5/6*(n^3) + 1/6*(n) (T. P. Martin, Shells of atoms, eq.(3)). - Brigitte Stepanov, Jul 02 2011
a(n+1) is the number of triples (w,x,y) having all terms in {0,...,n} and x+y >= w. - Clark Kimberling, Jun 14 2012
a(n) = Sum_{k=1..n} A215630(n,k) for n > 0. - Reinhard Zumkeller, Nov 11 2012
a(n) - a(n-2) = A010001(n-1), for n>1. - K. G. Stier, Dec 21 2012
a(n) is also a figurate number representing a cube of side n with a vertex cut off by a tetrahedron of side n-1. As such, a(n) = A000578(n) - A000292(n-1), n > 0. - Jean M. Morales, Aug 11 2013
The sequence starting with 1 is the third partial sum of (1, 4, 5, 5, 5, ...) and the binomial transform of (1, 6, 10, 5, 0, 0, 0, ...). - Gary W. Adamson, Sep 27 2015

Crossrefs

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

Programs

Formula

a(n) = 5*binomial(n + 1, 3) + binomial(n, 1).
a(n) = 5*n^3/6 + n/6.
a(n) = Sum_{i=0..n-1} A005891(i). - Xavier Acloque, Oct 08 2003
G.f.: x*(1+3*x+x^2) / (1-x)^4. - R. J. Mathar, Jun 05 2011
E.g.f.: (x/6)*(5x^2 + 15x + 6)*exp(x). - G. C. Greubel, Sep 27 2015
Sum_{n>0} 1/a(n) = 3*(2*gamma + polygamma(0, 1-i/sqrt(5)) + polygamma(0, 1+i/sqrt(5))) = 1.233988011257952852492845364799197179252... where i denotes the imaginary unit. - Stefano Spezia, Aug 31 2023

Extensions

Typo in definition corrected by Jean M. Morales, Aug 11 2013

A081436 Fifth subdiagonal in array of n-gonal numbers A081422.

Original entry on oeis.org

1, 7, 24, 58, 115, 201, 322, 484, 693, 955, 1276, 1662, 2119, 2653, 3270, 3976, 4777, 5679, 6688, 7810, 9051, 10417, 11914, 13548, 15325, 17251, 19332, 21574, 23983, 26565, 29326, 32272, 35409, 38743, 42280, 46026, 49987, 54169, 58578, 63220
Offset: 0

Views

Author

Paul Barry, Mar 21 2003

Keywords

Comments

One of a family of sequences with palindromic generators.
Also as A(n) = (1/6)*(6*n^3 - 3*n^2 + 3*n), n>0: structured pentagonal diamond numbers (vertex structure 5). (Cf. A004068 = alternate vertex; A000447 = structured diamonds; A100145 for more on structured numbers.) - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Sequence of the absolute values of the z^1 coefficients of the polynomials in the GF4 denominators of A156933. See A157705 for background information. - Johannes W. Meijer, Mar 07 2009
Row 1 of the convolution arrays A213831 and A213833. - Clark Kimberling, Jul 04 2012
Partial sums of A056109. - J. M. Bergot, Jun 22 2013
Number of ordered pairs of intersecting multisets of size 2, each chosen with repetition from {1,...,n}. - Robin Whitty, Feb 12 2014
Row sums of A244418. - L. Edson Jeffery, Jan 10 2015

Crossrefs

Programs

Formula

a(n) = (n+1)*(2*n^2 + 3*n + 2)/2.
G.f.: (1+x)*(1+2*x)/(1-x)^4. (Convolution of A005408 and A016777.)
a(n) = A110449(n, n-1), for n>1.
a(n) = (n+1)*T(n+1) + n*T(n), where T( ) are triangular numbers. Binomial transform of [1, 6, 11, 6, 0, 0, 0, ...]. - Gary W. Adamson, Dec 28 2007
E.g.f.: exp(x)*(2 + 12*x + 11*x^2 + 2*x^3)/2. - Stefano Spezia, Apr 13 2021
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Wesley Ivan Hurt, Apr 14 2021

Extensions

G.f. simplified and crossrefs added by Johannes W. Meijer, Mar 07 2009

A008956 Triangle of central factorial numbers |4^k t(2n+1,2n+1-2k)| read by rows (n>=0, k=0..n).

Original entry on oeis.org

1, 1, 1, 1, 10, 9, 1, 35, 259, 225, 1, 84, 1974, 12916, 11025, 1, 165, 8778, 172810, 1057221, 893025, 1, 286, 28743, 1234948, 21967231, 128816766, 108056025, 1, 455, 77077, 6092515, 230673443, 3841278805, 21878089479, 18261468225, 1, 680
Offset: 0

Views

Author

Keywords

Comments

The n-th row gives the coefficients in the expansion of Product_{i=0..n-1} (x+(2i+1)^2), highest powers first (see the discussion of central factorial numbers in A008955). - N. J. A. Sloane, Feb 01 2011
Descending row polynomials in x^2 evaluated at k generate odd coefficients of e.g.f. sin(arcsin(kt)/k): 1, x^2 - 1, 9x^4 - 10x^2 + 1, 225x^6 - 259x^4 + 34x^2 - 1, ... - Ralf Stephan, Jan 16 2005
From Johannes W. Meijer, Jun 18 2009: (Start)
We define (Pi/2)*Beta(n-1/2-z/2,n-1/2+z/2)/Beta(n-1/2,n-1/2) = (Pi/2)*Gamma(n-1/2-z/2)* Gamma(n-1/2+z/2)/Gamma(n-1/2)^2 = sum(BG2[2m,n]*z^(2m), m = 0..infinity) with Beta(z,w) the Beta function. Our definition leads to BG2[2m,1] = 2*beta(2m+1) and the recurrence relation BG2[2m,n] = BG2[2m,n-1] - BG2[2m-2,n-1]/(2*n-3)^2 for m = -2, -1, 0, 1, 2, .. and n = 2, 3, .. , with beta(m) = sum((-1)^k/(1+2*k)^m, k=0..infinity). We observe that beta(2m+1) = 0 for m = -1, -2, -3, .. .We found for the BG2[2*m,n] = sum((-1)^(k+n)*t2(n-1,k-1)* 2*beta(2*m-2*n+2*k+1),k=1..n)/((2*n-3)!!)^2 with the central factorial numbers t2(n,m) as defined above; see also the Maple program.
From the BG2 matrix and the closely related EG2 and ZG2 matrices, see A008955, we arrive at the LG2 matrix which is defined by LG2[2m-1,1] = 2*lambda(2*m) and the recurrence relation LG2[2*m-1,n] = LG2[2*m-3,n-1]/((2*n-3)*(2*n-1)) - (2*n-3)*LG2[2*m-1,n-1]/(2*n-1) for m = -2, -1, 0, 1, 2, .. and n = 2, 3, .. , with lambda(m) = (1-2^(-m))*zeta(m) with zeta(m) the Riemann zeta function. We found for the matrix coefficients LG2[2m-1,n] = sum((-1)^(k+1)* t2(n-1,k-1)*2*lambda(2*m-2*n+2*k)/((2*n-1)!!*(2*n-3)!!), k=1..n) and we see that the central factorial numbers t2(n,m) once again play a crucial role.
(End)

Examples

			Triangle begins:
[1]
[1, 1]
[1, 10, 9]
[1, 35, 259, 225]
[1, 84, 1974, 12916, 11025]
[1, 165, 8778, 172810, 1057221, 893025]
[1, 286, 28743, 1234948, 21967231, 128816766, 108056025]
[1, 455, 77077, 6092515, 230673443, 3841278805, 21878089479, 18261468225]
...
		

References

  • P. L. Butzer, M. Schmidt, E. L. Stark and L. Vogt, Central Factorial Numbers: Their main properties and some applications, Numerical Functional Analysis and Optimization, 10 (5&6), 419-488 (1989). [From Johannes W. Meijer, Jun 18 2009]
  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 217.

Crossrefs

Cf. A008958.
Columns include A000447, A001823. Right-hand columns include A001818, A001824, A001825. Cf. A008955.
Appears in A160480 (Beta triangle), A160487 (Lambda triangle), A160479 (ZL(n) sequence), A161736, A002197 and A002198. - Johannes W. Meijer, Jun 18 2009
Cf. A162443 (BG1 matrix) and A162448 (LG1 matrix). - Johannes W. Meijer, Jul 06 2009
Cf. A001147.

Programs

  • Haskell
    a008956 n k = a008956_tabl !! n !! k
    a008956_row n = a008956_tabl !! n
    a008956_tabl = [1] : f [1] 1 1 where
       f xs u t = ys : f ys v (t * v) where
         ys = zipWith (+) (xs ++ [t^2]) ([0] ++ map (* u^2) (init xs) ++ [0])
         v = u + 2
    -- Reinhard Zumkeller, Dec 24 2013
  • Maple
    f:=n->mul(x+(2*i+1)^2,i=0..n-1);
    for n from 0 to 12 do
    t1:=eval(f(n)); t1d:=degree(t1);
    t12:=y^t1d*subs(x=1/y,t1); t2:=seriestolist(series(t12,y,20));
    lprint(t2);
    od: # N. J. A. Sloane, Feb 01 2011
    A008956 := proc(n,k) local i ; mul( x+2*i-2*n-1,i=1..2*n) ; expand(%) ; coeftayl(%,x=0,2*(n-k)) ; abs(%) ; end: for n from 0 to 10 do for k from 0 to n do printf("%a,",A008956(n,k)) ; od: od: # R. J. Mathar, May 29 2009
    nmax:=7: for n from 0 to nmax do t2(n, 0):=1: t2(n, n):=(doublefactorial(2*n-1))^2 od: for n from 1 to nmax do for k from 1 to n-1 do t2(n, k) := (2*n-1)^2*t2(n-1, k-1)+t2(n-1, k) od: od: seq(seq(t2(n, k), k=0..n), n=0..nmax); # Johannes W. Meijer, Jun 18 2009, Revised Sep 16 2012
  • Mathematica
    t[, 0] = 1; t[n, n_] := t[n, n] = ((2*n-1)!!)^2; t[n_, k_] := t[n, k] = (2*n-1)^2*t[n-1, k-1] + t[n-1, k]; Table[t[n, k], {n, 0, 8}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jan 07 2014, after Johannes W. Meijer *)
  • PARI
    {T(n, k) = if( n<=0, k==0, (-1)^k * polcoeff( numerator( 2^(2*n -1) / sum(j=0, 2*n - 1, binomial( 2*n - 1, j) / (x + 2*n - 1 - 2*j))), 2*n - 2*k))}; /* Michael Somos, Feb 24 2003 */
    

Formula

Conjecture row sums: Sum_{k=0..n} T(n,k) = |A101927(n+1)|. - R. J. Mathar, May 29 2009
May be generated by the recurrence t2(n,k) = (2*n-1)^2*t2(n-1,k-1)+t2(n-1,k) with t2(n,0) = 1 and t2(n,n)=((2*n-1)!!)^2. - Johannes W. Meijer, Jun 18 2009

Extensions

More terms from Vladeta Jovovic, Apr 16 2000
Edited by N. J. A. Sloane, Feb 01 2011

A004188 a(n) = n*(3*n^2 - 1)/2.

Original entry on oeis.org

0, 1, 11, 39, 94, 185, 321, 511, 764, 1089, 1495, 1991, 2586, 3289, 4109, 5055, 6136, 7361, 8739, 10279, 11990, 13881, 15961, 18239, 20724, 23425, 26351, 29511, 32914, 36569, 40485, 44671, 49136, 53889, 58939, 64295, 69966, 75961
Offset: 0

Views

Author

Albert D. Rich (Albert_Rich(AT)msn.com)

Keywords

Comments

3-dimensional analog of centered polygonal numbers.
(1), (4+7), (10+13+16), (19+22+25+28), ... - Jon Perry, Sep 10 2004

References

  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 140.
  • T. P. Martin, Shells of atoms, Phys. Reports, 273 (1996), 199-241, eq. (11).

Crossrefs

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. A236770 (partial sums).

Programs

Formula

Partial sums of n-1 3-spaced triangular numbers, e.g., a(4) = t(1) + t(4) + t(7) = 1 + 10 + 28 = 39. - Jon Perry, Jul 23 2003
a(n) = C(2*n+1,3) + C(n+1,3), n >= 0. - Zerinvary Lajos, Jan 21 2007
a(n) = A000447(n) + A000292(n). - Zerinvary Lajos, Jan 21 2007
G.f.: x*(1+7*x+x^2) / (x-1)^4. - R. J. Mathar, Oct 08 2011
From Miquel Cerda, Dec 25 2016: (Start)
a(n) = A000578(n) + A135503(n).
a(n) = A007588(n) - A135503(n). (End)
E.g.f.: (x/2)*(2 + 9*x + 3*x^2)*exp(x). - G. C. Greubel, Sep 01 2017

A004466 a(n) = n*(5*n^2 - 2)/3.

Original entry on oeis.org

0, 1, 12, 43, 104, 205, 356, 567, 848, 1209, 1660, 2211, 2872, 3653, 4564, 5615, 6816, 8177, 9708, 11419, 13320, 15421, 17732, 20263, 23024, 26025, 29276, 32787, 36568, 40629, 44980, 49631, 54592
Offset: 0

Views

Author

Albert D. Rich (Albert_Rich(AT)msn.com)

Keywords

Comments

3-dimensional analog of centered polygonal numbers.
Also as a(n)=(1/6)*(10*n^3-4*n), n>0: structured pentagonal anti-diamond numbers (vertex structure 11) (Cf. A051673 = alternate vertex A100188 = structured anti-diamonds; A100145 for more on structured numbers). - James A. Record (james.record(AT)gmail.com), Nov 07 2004
a(n+1)-10*a(n) = (n+1)*(5*(n+1)^2-2)/3 - (10n(n+1)(n+2)/6) = n. The unit digits are 0,1,2,3,4,5,6,7,8,9,0,1,2,3,4,5,6,7,8,9,... . - Eric Desbiaux, Aug 18 2008

References

  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 140.

Crossrefs

Cf. A062786 (first differences), A264853 (partial sums).
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.

Programs

Formula

G.f.: x*(1+8*x+x^2)/(1-x)^4. - Colin Barker, Jan 08 2012
E.g.f.: (x/3)*(3 + 15*x + 5*x^2)*exp(x). - G. C. Greubel, Sep 01 2017

A063521 a(n) = n*(7*n^2-4)/3.

Original entry on oeis.org

0, 1, 16, 59, 144, 285, 496, 791, 1184, 1689, 2320, 3091, 4016, 5109, 6384, 7855, 9536, 11441, 13584, 15979, 18640, 21581, 24816, 28359, 32224, 36425, 40976, 45891, 51184, 56869, 62960, 69471, 76416, 83809, 91664, 99995, 108816, 118141, 127984, 138359, 149280, 160761
Offset: 0

Views

Author

N. J. A. Sloane, Aug 02 2001

Keywords

Comments

Also as a(n)=(1/6)*(14*n^3-8*n), n>0: structured heptagonal anti-diamond numbers (vertex structure 15) (Cf. A100186 = alternate vertex; A100188 = structured anti-diamonds; A100145 for more on structured numbers). - James A. Record (james.record(AT)gmail.com), Nov 07 2004

Crossrefs

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.

Programs

Formula

G.f.: x*(1+12*x+x^2)/(1-x)^4. - Colin Barker, Jan 10 2012
E.g.f.: (x/3)*(3 + 21*x + 7*x^2)*exp(x). - G. C. Greubel, Sep 01 2017
Previous Showing 11-20 of 83 results. Next