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-9 of 9 results.

A000578 The cubes: a(n) = n^3.

Original entry on oeis.org

0, 1, 8, 27, 64, 125, 216, 343, 512, 729, 1000, 1331, 1728, 2197, 2744, 3375, 4096, 4913, 5832, 6859, 8000, 9261, 10648, 12167, 13824, 15625, 17576, 19683, 21952, 24389, 27000, 29791, 32768, 35937, 39304, 42875, 46656, 50653, 54872, 59319, 64000, 68921, 74088, 79507
Offset: 0

Views

Author

Keywords

Comments

a(n) is the sum of the next n odd numbers; i.e., group the odd numbers so that the n-th group contains n elements like this: (1), (3, 5), (7, 9, 11), (13, 15, 17, 19), (21, 23, 25, 27, 29), ...; then each group sum = n^3 = a(n). Also the median of each group = n^2 = mean. As the sum of first n odd numbers is n^2 this gives another proof of the fact that the n-th partial sum = (n(n + 1)/2)^2. - Amarnath Murthy, Sep 14 2002
Total number of triangles resulting from criss-crossing cevians within a triangle so that two of its sides are each n-partitioned. - Lekraj Beedassy, Jun 02 2004. See Propp and Propp-Gubin for a proof.
Also structured triakis tetrahedral numbers (vertex structure 7) (cf. A100175 = alternate vertex); structured tetragonal prism numbers (vertex structure 7) (cf. A100177 = structured prisms); structured hexagonal diamond numbers (vertex structure 7) (cf. A100178 = alternate vertex; A000447 = structured diamonds); and structured trigonal anti-diamond numbers (vertex structure 7) (cf. A100188 = structured anti-diamonds). Cf. A100145 for more on structured polyhedral numbers. - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Schlaefli symbol for this polyhedron: {4, 3}.
Least multiple of n such that every partial sum is a square. - Amarnath Murthy, Sep 09 2005
Draw a regular hexagon. Construct points on each side of the hexagon such that these points divide each side into equally sized segments (i.e., a midpoint on each side or two points on each side placed to divide each side into three equally sized segments or so on), do the same construction for every side of the hexagon so that each side is equally divided in the same way. Connect all such points to each other with lines that are parallel to at least one side of the polygon. The result is a triangular tiling of the hexagon and the creation of a number of smaller regular hexagons. The equation gives the total number of regular hexagons found where n = the number of points drawn + 1. For example, if 1 point is drawn on each side then n = 1 + 1 = 2 and a(n) = 2^3 = 8 so there are 8 regular hexagons in total. If 2 points are drawn on each side then n = 2 + 1 = 3 and a(n) = 3^3 = 27 so there are 27 regular hexagons in total. - Noah Priluck (npriluck(AT)gmail.com), May 02 2007
The solutions of the Diophantine equation: (X/Y)^2 - X*Y = 0 are of the form: (n^3, n) with n >= 1. The solutions of the Diophantine equation: (m^2)*(X/Y)^2k - XY = 0 are of the form: (m*n^(2k + 1), m*n^(2k - 1)) with m >= 1, k >= 1 and n >= 1. The solutions of the Diophantine equation: (m^2)*(X/Y)^(2k + 1) - XY = 0 are of the form: (m*n^(k + 1), m*n^k) with m >= 1, k >= 1 and n >= 1. - Mohamed Bouhamida, Oct 04 2007
Except for the first two terms, the sequence corresponds to the Wiener indices of C_{2n} i.e., the cycle on 2n vertices (n > 1). - K.V.Iyer, Mar 16 2009
Totally multiplicative sequence with a(p) = p^3 for prime p. - Jaroslav Krizek, Nov 01 2009
Sums of rows of the triangle in A176271, n > 0. - Reinhard Zumkeller, Apr 13 2010
One of the 5 Platonic polyhedral (tetrahedral, cube, octahedral, dodecahedral and icosahedral) numbers (cf. A053012). - Daniel Forgues, May 14 2010
Numbers n for which order of torsion subgroup t of the elliptic curve y^2 = x^3 - n is t = 2. - Artur Jasinski, Jun 30 2010
The sequence with the lengths of the Pisano periods mod k is 1, 2, 3, 4, 5, 6, 7, 8, 3, 10, 11, 12, 13, 14, 15, 16, 17, 6, 19, 20, ... for k >= 1, apparently multiplicative and derived from A000027 by dividing every ninth term through 3. Cubic variant of A186646. - R. J. Mathar, Mar 10 2011
The number of atoms in a bcc (body-centered cubic) rhombic hexahedron with n atoms along one edge is n^3 (T. P. Martin, Shells of atoms, eq. (8)). - Brigitte Stepanov, Jul 02 2011
The inverse binomial transform yields the (finite) 0, 1, 6, 6 (third row in A019538 and A131689). - R. J. Mathar, Jan 16 2013
Twice the area of a triangle with vertices at (0, 0), (t(n - 1), t(n)), and (t(n), t(n - 1)), where t = A000217 are triangular numbers. - J. M. Bergot, Jun 25 2013
If n > 0 is not congruent to 5 (mod 6) then A010888(a(n)) divides a(n). - Ivan N. Ianakiev, Oct 16 2013
For n > 2, a(n) = twice the area of a triangle with vertices at points (binomial(n,3),binomial(n+2,3)), (binomial(n+1,3),binomial(n+1,3)), and (binomial(n+2,3),binomial(n,3)). - J. M. Bergot, Jun 14 2014
Determinants of the spiral knots S(4,k,(1,1,-1)). a(k) = det(S(4,k,(1,1,-1))). - Ryan Stees, Dec 14 2014
One of the oldest-known examples of this sequence is shown in the Senkereh tablet, BM 92698, which displays the first 32 terms in cuneiform. - Charles R Greathouse IV, Jan 21 2015
From Bui Quang Tuan, Mar 31 2015: (Start)
We construct a number triangle from the integers 1, 2, 3, ... 2*n-1 as follows. The first column contains all the integers 1, 2, 3, ... 2*n-1. Each succeeding column is the same as the previous column but without the first and last items. The last column contains only n. The sum of all the numbers in the triangle is n^3.
Here is the example for n = 4, where 1 + 2*2 + 3*3 + 4*4 + 3*5 + 2*6 + 7 = 64 = a(4):
1
2 2
3 3 3
4 4 4 4
5 5 5
6 6
7
(End)
For n > 0, a(n) is the number of compositions of n+11 into n parts avoiding parts 2 and 3. - Milan Janjic, Jan 07 2016
Does not satisfy Benford's law [Ross, 2012]. - N. J. A. Sloane, Feb 08 2017
Number of inequivalent face colorings of the cube using at most n colors such that each color appears at least twice. - David Nacin, Feb 22 2017
Consider A = {a,b,c} a set with three distinct members. The number of subsets of A is 8, including {a,b,c} and the empty set. The number of subsets from each of those 8 subsets is 27. If the number of such iterations is n, then the total number of subsets is a(n-1). - Gregory L. Simay, Jul 27 2018
By Fermat's Last Theorem, these are the integers of the form x^k with the least possible value of k such that x^k = y^k + z^k never has a solution in positive integers x, y, z for that k. - Felix Fröhlich, Jul 27 2018

Examples

			For k=3, b(3) = 2 b(2) - b(1) = 4-1 = 3, so det(S(4,3,(1,1,-1))) = 3*3^2 = 27.
For n=3, a(3) = 3 + (3*0^2 + 3*0 + 3*1^2 + 3*1 + 3*2^2 + 3*2) = 27. - _Patrick J. McNab_, Mar 28 2016
		

References

  • Albert H. Beiler, Recreations in the theory of numbers, New York, Dover, (2nd ed.) 1966. See p. 191.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 43, 64, 81.
  • R. L. Graham, D. E. Knuth, and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 255; 2nd. ed., p. 269. Worpitzky's identity (6.37).
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.6 Figurate Numbers, p. 292.
  • T. Aaron Gulliver, "Sequences from cubes of integers", International Mathematical Journal, 4 (2003), no. 5, 439 - 445. See http://www.m-hikari.com/z2003.html for information about this journal. [I expanded the reference to make this easier to find. - N. J. A. Sloane, Feb 18 2019]
  • J. Propp and A. Propp-Gubin, "Counting Triangles in Triangles", Pi Mu Epsilon Journal (to appear).
  • 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).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 6-7.
  • D. Wells, You Are A Mathematician, pp. 238-241, Penguin Books 1995.

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.
For sums of cubes, cf. A000537 (partial sums), A003072, A003325, A024166, A024670, A101102 (fifth partial sums).
Cf. A001158 (inverse Möbius transform), A007412 (complement), A030078(n) (cubes of primes), A048766, A058645 (binomial transform), A065876, A101094, A101097.
Subsequence of A145784.
Cf. A260260 (comment). - Bruno Berselli, Jul 22 2015
Cf. A000292 (tetrahedral numbers), A005900 (octahedral numbers), A006566 (dodecahedral numbers), A006564 (icosahedral numbers).
Cf. A098737 (main diagonal).

Programs

  • Haskell
    a000578 = (^ 3)
    a000578_list = 0 : 1 : 8 : zipWith (+)
       (map (+ 6) a000578_list)
       (map (* 3) $ tail $ zipWith (-) (tail a000578_list) a000578_list)
    -- Reinhard Zumkeller, Sep 05 2015, May 24 2012, Oct 22 2011
    
  • Magma
    [ n^3 : n in [0..50] ]; // Wesley Ivan Hurt, Jun 14 2014
    
  • Magma
    I:=[0,1,8,27]; [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..45]]; // Vincenzo Librandi, Jul 05 2014
    
  • Maple
    A000578 := n->n^3;
    seq(A000578(n), n=0..50);
    isA000578 := proc(r)
        local p;
        if r = 0 or r =1 then
            true;
        else
            for p in ifactors(r)[2] do
                if op(2, p) mod 3 <> 0 then
                    return false;
                end if;
            end do:
            true ;
        end if;
    end proc: # R. J. Mathar, Oct 08 2013
  • Mathematica
    Table[n^3, {n, 0, 30}] (* Stefan Steinerberger, Apr 01 2006 *)
    CoefficientList[Series[x (1 + 4 x + x^2)/(1 - x)^4, {x, 0, 45}], x] (* Vincenzo Librandi, Jul 05 2014 *)
    Accumulate[Table[3n^2+3n+1,{n,0,20}]] (* or *) LinearRecurrence[{4,-6,4,-1},{1,8,27,64},20](* Harvey P. Dale, Aug 18 2018 *)
  • Maxima
    A000578(n):=n^3$
    makelist(A000578(n),n,0,30); /* Martin Ettl, Nov 03 2012 */
    
  • PARI
    A000578(n)=n^3 \\ M. F. Hasler, Apr 12 2008
    
  • PARI
    is(n)=ispower(n,3) \\ Charles R Greathouse IV, Feb 20 2012
    
  • Python
    A000578_list, m = [], [6, -6, 1, 0]
    for _ in range(10**2):
        A000578_list.append(m[-1])
        for i in range(3):
            m[i+1] += m[i] # Chai Wah Wu, Dec 15 2015
    
  • Scheme
    (define (A000578 n) (* n n n)) ;; Antti Karttunen, Oct 06 2017

Formula

a(n) = Sum_{i=0..n-1} A003215(i).
Multiplicative with a(p^e) = p^(3e). - David W. Wilson, Aug 01 2001
G.f.: x*(1+4*x+x^2)/(1-x)^4. - Simon Plouffe in his 1992 dissertation
Dirichlet generating function: zeta(s-3). - Franklin T. Adams-Watters, Sep 11 2005, Amarnath Murthy, Sep 09 2005
E.g.f.: (1+3*x+x^2)*x*exp(x). - Franklin T. Adams-Watters, Sep 11 2005 - Amarnath Murthy, Sep 09 2005
a(n) = Sum_{i=1..n} (Sum_{j=i..n+i-1} A002024(j,i)). - Reinhard Zumkeller, Jun 24 2007
a(n) = lcm(n, (n - 1)^2) - (n - 1)^2. E.g.: lcm(1, (1 - 1)^2) - (1 - 1)^2 = 0, lcm(2, (2 - 1)^2) - (2 - 1)^2 = 1, lcm(3, (3 - 1)^2) - (3 - 1)^2 = 8, ... - Mats Granvik, Sep 24 2007
Starting (1, 8, 27, 64, 125, ...), = binomial transform of [1, 7, 12, 6, 0, 0, 0, ...]. - Gary W. Adamson, Nov 21 2007
a(n) = A007531(n) + A000567(n). - Reinhard Zumkeller, Sep 18 2009
a(n) = binomial(n+2,3) + 4*binomial(n+1,3) + binomial(n,3). [Worpitzky's identity for cubes. See. e.g., Graham et al., eq. (6.37). - Wolfdieter Lang, Jul 17 2019]
a(n) = n + 6*binomial(n+1,3) = binomial(n,1)+6*binomial(n+1,3). - Ron Knott, Jun 10 2019
A010057(a(n)) = 1. - Reinhard Zumkeller, Oct 22 2011
a(n) = A000537(n) - A000537(n-1), difference between 2 squares of consecutive triangular numbers. - Pierre CAMI, Feb 20 2012
a(n) = A048395(n) - 2*A006002(n). - J. M. Bergot, Nov 25 2012
a(n) = 1 + 7*(n-1) + 6*(n-1)*(n-2) + (n-1)*(n-2)*(n-3). - Antonio Alberto Olivares, Apr 03 2013
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 6. - Ant King Apr 29 2013
a(n) = A000330(n) + Sum_{i=1..n-1} A014105(i), n >= 1. - Ivan N. Ianakiev, Sep 20 2013
a(k) = det(S(4,k,(1,1,-1))) = k*b(k)^2, where b(1)=1, b(2)=2, b(k) = 2*b(k-1) - b(k-2) = b(2)*b(k-1) - b(k-2). - Ryan Stees, Dec 14 2014
For n >= 1, a(n) = A152618(n-1) + A033996(n-1). - Bui Quang Tuan, Apr 01 2015
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Jon Tavasanis, Feb 21 2016
a(n) = n + Sum_{j=0..n-1} Sum_{k=1..2} binomial(3,k)*j^(3-k). - Patrick J. McNab, Mar 28 2016
a(n) = A000292(n-1) * 6 + n. - Zhandos Mambetaliyev, Nov 24 2016
a(n) = n*binomial(n+1, 2) + 2*binomial(n+1, 3) + binomial(n,3). - Tony Foster III, Nov 14 2017
From Amiram Eldar, Jul 02 2020: (Start)
Sum_{n>=1} 1/a(n) = zeta(3) (A002117).
Sum_{n>=1} (-1)^(n+1)/a(n) = 3*zeta(3)/4 (A197070). (End)
From Amiram Eldar, Jan 20 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = cosh(sqrt(3)*Pi/2)/Pi.
Product_{n>=2} (1 - 1/a(n)) = cosh(sqrt(3)*Pi/2)/(3*Pi). (End)
a(n) = Sum_{d|n} sigma_3(d)*mu(n/d) = Sum_{d|n} A001158(d)*A008683(n/d). Moebius transform of sigma_3(n). - Ridouane Oudra, Apr 15 2021

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

A140113 a(1)=1, a(n)=a(n-1)+n if n odd, a(n)=a(n-1)+ n^2 if n is even.

Original entry on oeis.org

1, 5, 8, 24, 29, 65, 72, 136, 145, 245, 256, 400, 413, 609, 624, 880, 897, 1221, 1240, 1640, 1661, 2145, 2168, 2744, 2769, 3445, 3472, 4256, 4285, 5185, 5216, 6240, 6273, 7429, 7464, 8760, 8797, 10241, 10280, 11880, 11921, 13685, 13728, 15664, 15709
Offset: 1

Views

Author

Artur Jasinski, May 08 2008

Keywords

Comments

One notices the powers 8, 256, 400, and 2744 (14^3) and wonders if the sum is ever again a power. [J. M. Bergot, Sep 07 2011]

Crossrefs

Cf. A136047.

Programs

  • Mathematica
    nxt[{n_,a_}]:={n+1,If[EvenQ[n],a+n+1,a+(n+1)^2]}; Transpose[ NestList[ nxt,{1,1},50]][[2]] (* or *) LinearRecurrence[{1,3,-3,-3,3,1,-1},{1,5,8,24,29,65,72},50] (* Harvey P. Dale, Jul 22 2014 *)
    CoefficientList[Series[(- x^4 + 4 x^3 + 4 x + 1)/(x^7 - x^6 - 3 x^5 + 3 x^4 + 3 x^3 - 3 x^2 - x + 1), {x, 0, 40}], x] (* Vincenzo Librandi, Jul 23 2014 *)
  • PARI
    print1(a=1);for(n=2,99,print1(", ",a+=n^(2-n%2))) \\ Charles R Greathouse IV, Jul 19 2011

Formula

O.g.f.: (-x^4 + 4*x^3 + 4*x + 1)/(x^7 - x^6 - 3*x^5 + 3*x^4 + 3*x^3 - 3*x^2 - x + 1). - Alexander R. Povolotsky, May 08 2008
a(2*n) = A185872(n,2); a(2*n-1) = A100178(n). - Franck Maminirina Ramaharo, Feb 26 2018

A059722 a(n) = n*(2*n^2 - 2*n + 1).

Original entry on oeis.org

0, 1, 10, 39, 100, 205, 366, 595, 904, 1305, 1810, 2431, 3180, 4069, 5110, 6315, 7696, 9265, 11034, 13015, 15220, 17661, 20350, 23299, 26520, 30025, 33826, 37935, 42364, 47125, 52230, 57691, 63520, 69729, 76330, 83335, 90756, 98605, 106894, 115635, 124840
Offset: 0

Views

Author

Henry Bottomley, Feb 07 2001

Keywords

Comments

Mean of the first four nonnegative powers of 2n+1, i.e., ((2n+1)^0 + (2n+1)^1 + (2n+1)^2 + (2n+1)^3)/4. E.g., a(2) = (1 + 3 + 9 + 27)/4 = 10.
Equatorial structured meta-diamond numbers, the n-th number from an equatorial structured n-gonal diamond number sequence. There are no 1- or 2-gonal diamonds, so 1 and (n+2) are used as the first and second terms since all the sequences begin as such. - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Starting with offset 1 = row sums of triangle A143803. - Gary W. Adamson, Sep 01 2008
Form an array from the antidiagonals containing the terms in A002061 to give antidiagonals 1; 3,3; 7,4,7; 13,8,8,13; 21,14,9,14,21; and so on. The difference between the sum of the terms in n+1 X n+1 matrices and those in n X n matrices is a(n) for n>0. - J. M. Bergot, Jul 08 2013
Sum of the numbers from (n-1)^2 to n^2. - Wesley Ivan Hurt, Sep 08 2014

Crossrefs

Cf. A000330, A005900, A081436, A100178, A100179, A059722: "equatorial" structured diamonds; A000447: "polar" structured meta-diamond; A006484 for other structured meta numbers; and A100145 for more on structured numbers. - James A. Record (james.record(AT)gmail.com), Nov 07 2004

Programs

Formula

a(n) = A053698(2*n-1)/4.
a(n) = Sum_{j=1..n} ((n+j-1)^2-j^2+1). - Zerinvary Lajos, Sep 13 2006
From R. J. Mathar, Sep 02 2008: (Start)
G.f.: x*(1 + x)*(1 + 5*x)/(1 - x)^4.
a(n) = A002414(n-1) + A002414(n).
a(n+1) - a(n) = A136392(n+1). (End)
a(n) = (A000290(n) + A000290(n+1)) * (A000217(n+1) - A000217(n)). - J. M. Bergot, Nov 02 2012
a(n) = n * A001844(n-1). - Doug Bell, Aug 18 2015
a(n) = A000217(n^2) - A000217(n^2-2*n). - Bruno Berselli, Jun 26 2018
E.g.f.: exp(x)*x*(1 + 4*x + 2*x^2). - Stefano Spezia, Jun 20 2021

Extensions

Edited with new definition by N. J. A. Sloane, Aug 29 2008

A100185 Structured meta-anti-prism numbers, the n-th number from a structured n-gonal anti-prism number sequence.

Original entry on oeis.org

1, 4, 19, 68, 185, 416, 819, 1464, 2433, 3820, 5731, 8284, 11609, 15848, 21155, 27696, 35649, 45204, 56563, 69940, 85561, 103664, 124499, 148328, 175425, 206076, 240579, 279244, 322393, 370360, 423491, 482144
Offset: 1

Views

Author

James A. Record (james.record(AT)gmail.com)

Keywords

Examples

			There are no 1- or 2-gonal anti-prisms, so 1 and (2n) are used as the first and second terms since all the sequences begin as such.
		

Crossrefs

Cf. A005900, A000447, A096000, A100178, A100157, A100185 - structured anti-prisms; A006484 for other structured meta numbers; and A100145 for more on structured numbers.

Programs

  • Magma
    [(1/6)*(3*n^4-8*n^3+9*n^2+2*n): n in [1..40]]; // Vincenzo Librandi, Aug 03 2011

Formula

a(n) = (1/6)*(3*n^4 - 8*n^3 + 9*n^2 + 2*n).
G.f.: x*(1 - x + 9*x^2 + 3*x^3)/(1-x)^5. [Colin Barker, Jun 08 2012]

A213841 Rectangular array: (row n) = b**c, where b(h) = 2*h-1, c(h) = 4*n-7+4*h, n>=1, h>=1, and ** = convolution.

Original entry on oeis.org

1, 8, 5, 29, 24, 9, 72, 65, 40, 13, 145, 136, 101, 56, 17, 256, 245, 200, 137, 72, 21, 413, 400, 345, 264, 173, 88, 25, 624, 609, 544, 445, 328, 209, 104, 29, 897, 880, 805, 688, 545, 392, 245, 120, 33, 1240, 1221, 1136, 1001
Offset: 1

Views

Author

Clark Kimberling, Jul 05 2012

Keywords

Comments

Principal diagonal: A213842.
Antidiagonal sums: A213843.
Row 1, (1,5,9,13,...)**(1,3,5,7,...): A100178.
Row 2, (1,5,9,13,...)**(3,5,7,9,...): (4*k^3 + 9*k^2 + 2*k)/3.
Row 3, (1,5,9,13,...)**(5,7,9,11,...): (4*k^3 + 21*k^2 + 2*k)/3.
For a guide to related arrays, see A212500.

Examples

			Northwest corner (the array is read by falling antidiagonals):
1....8....29....72....145
5....24...65....136...245
9....40...101...200...345
13...56...137...264...445
17...72...173...328...545
21...88...209...392...645
		

Crossrefs

Cf. A212500.

Programs

  • Mathematica
    b[n_]:=2n-1;c[n_]:=4n-3;
    t[n_,k_]:=Sum[b[k-i]c[n+i],{i,0,k-1}]
    TableForm[Table[t[n,k],{n,1,10},{k,1,10}]]
    Flatten[Table[t[n-k+1,k],{n,12},{k,n,1,-1}]]
    r[n_]:=Table[t[n,k],{k,1,60}] (* A213841 *)
    Table[t[n,n],{n,1,40}] (* A213842 *)
    s[n_]:=Sum[t[i,n+1-i],{i,1,n}]
    Table[s[n],{n,1,50}] (* A213843 *)

Formula

T(n,k) = 4*T(n,k-1)-6*T(n,k-2)+4*T(n,k-3)-T(n,k-4).
G.f. for row n: f(x)/g(x), where f(x) = x*(4*n-3 + 4*x - (4*n-7)*x^2) and g(x) = (1-x)^4.

A172482 a(n) = (1+n)*(9 + 11*n + 4*n^2)/3.

Original entry on oeis.org

3, 16, 47, 104, 195, 328, 511, 752, 1059, 1440, 1903, 2456, 3107, 3864, 4735, 5728, 6851, 8112, 9519, 11080, 12803, 14696, 16767, 19024, 21475, 24128, 26991, 30072, 33379, 36920, 40703, 44736, 49027, 53584, 58415, 63528, 68931, 74632, 80639, 86960, 93603
Offset: 0

Views

Author

Paul Curtz, Feb 04 2010

Keywords

Comments

One of the bisections of the left central column in the Janet table A172002.
Row 1 of the convolution array A213844. - Clark Kimberling, Jul 05 2012
With offset 2, this is 4*n^3/3 - 3*n^2 + 8*n/3 - 1, the number of divisions of a 2 X n board into 3 pieces where the rightmost squares separate. See Jacob Brown article. - Michel Marcus, Jun 29 2021

Crossrefs

Programs

Formula

a(n) = A131941(2n+2), where A100178(n) = A131941(2n-1).
a(n) = 4*a(n) - 6*a(n-2) + 4*a(n-3) - a(n-4).
a(n) mod 10 = 3, 6, 7, 4, 5, 8, 1, 2, 9, 0 (and repeat periodically).
G.f.: (x+3)*(1+x)/(x-1)^4.
E.g.f.: exp(x)*(9 + 39*x + 27*x^2 + 4*x^3)/3. - Stefano Spezia, Mar 02 2025

Extensions

Edited by R. J. Mathar, Feb 24 2010

A213838 Rectangular array: (row n) = b**c, where b(h) = 4*h-3, c(h) = 2*n-3+2*h, n>=1, h>=1, and ** = convolution.

Original entry on oeis.org

1, 8, 3, 29, 20, 5, 72, 59, 32, 7, 145, 128, 89, 44, 9, 256, 235, 184, 119, 56, 11, 413, 388, 325, 240, 149, 68, 13, 624, 595, 520, 415, 296, 179, 80, 15, 897, 864, 777, 652, 505, 352, 209, 92, 17, 1240, 1203, 1104, 959, 784
Offset: 1

Views

Author

Clark Kimberling, Jul 05 2012

Keywords

Comments

Principal diagonal: A213839.
Antidiagonal sums: A213840.
Row 1, (1,5,9,13,...)**(1,3,5,7,...): A100178.
Row 2, (1,5,9,13,...)**(3,5,7,9,...): (4*k^3 + 9*k^2 - 4*k)/3.
Row 3, (1,5,9,13,...)**(5,7,9,11,...): (4*k^3 + 21*k^2 - 10*k)/3.
For a guide to related arrays, see A212500.

Examples

			Northwest corner (the array is read by falling antidiagonals):
1....8....29....72....145
3....20...59....128...235
5....32...89....184...325
7....44...119...240...415
9....56...149...296...505
11...68...179...352...595
		

Crossrefs

Cf. A212500.

Programs

  • Mathematica
    b[n_]:=4n-3; c[n_]:=2n-1;
    t[n_,k_]:=Sum[b[k-i]c[n+i],{i,0,k-1}]
    TableForm[Table[t[n,k],{n,1,10},{k,1,10}]]
    Flatten[Table[t[n-k+1,k],{n,12},{k,n,1,-1}]]
    r[n_]:=Table[t[n,k],{k,1,60}] (* A213838 *)
    Table[t[n,n],{n,1,40}] (* A213839 *)
    s[n_]:=Sum[t[i,n+1-i],{i,1,n}]
    Table[s[n],{n,1,50}] (* A213840 *)

Formula

T(n,k) = 4*T(n,k-1)-6*T(n,k-2)+4*T(n,k-3)-T(n,k-4).
G.f. for row n: f(x)/g(x), where f(x) = x*(2*n-1 + 4*n*x - (6*n-9)*x^2) and g(x) = (1-x)^4.

A167471 Janet periodic table of the elements and structured hexagonal diamond numbers. a(n) = A166911(2n) + A166911(2n+1).

Original entry on oeis.org

16, 128, 464, 1152, 2320, 4096, 6608, 9984, 14352, 19840, 26576, 34688, 44304, 55552, 68560, 83456, 100368, 119424, 140752, 164480, 190736, 219648, 251344, 285952, 323600, 364416, 408528, 456064, 507152, 561920, 620496, 683008, 749584
Offset: 1

Views

Author

Paul Curtz, Nov 04 2009

Keywords

Comments

(A166911=3,13,39,89,171,293,). 5 formulas for a(n):fifth is a(n)=4a(n-1)-6a(n-2)+4a(n-3)-a(n-4), a(0)=16,a(1)=128,a(2)=464,a(3)=1152. See A166464 and its companion A166911 (offset 0). Note, for a(n) mod 10, period 5:repeat 6,8,4,2,0. Mathematically, a(n) (then A100178=1,8,29,72,145,) can be preceded by 0.

Programs

  • Magma
    [16*n*(2-3*n+4*n^2)/3: n in [1..40]]; // Vincenzo Librandi, Aug 03 2011
  • Mathematica
    Table[16*n*(2 - 3*n + 4*n^2)/3, {n,1,50}] (* G. C. Greubel, Jun 13 2016 *)

Formula

a(n) = 16*n*(2-3*n+4*n^2)/3 = 16*(2*n-3*n^2+4*n^3)/3 = 16*A100178(n).
Showing 1-9 of 9 results.