Displaying 1-10 of 142 results found.
page
1
2
3
4
5
6
7
8
9
10
... 15
Pentagonal pyramidal numbers which are the sum of two other such numbers: A002411(k) = A002411(i)+ A002411(j) with i,j>0.
+20
3
37926, 279046, 2514726, 34797726, 71254566, 348145726, 566225926, 606245926, 828497488, 1206646551, 8750871976, 11481404326, 21078151200, 28120290876, 62963640936, 75128827176, 77009692800, 96868002726, 120282238710, 147061923750, 165998399400, 297022824126, 325443925926, 416397477888
COMMENTS
The corresponding indices are listed in A172437. This answers a question posed in A133459.
PROG
(PARI) for(n=1, 99999, for(m=1, n-1, isA002411( A002411(n)- A002411(m)) | next; print1( A002411(n)", "); break)) /* needs A002411() and isA240011() */
Numbers that are sums of consecutive pentagonal pyramidal numbers ( A002411).
+20
3
0, 1, 6, 7, 18, 24, 25, 40, 58, 64, 65, 75, 115, 126, 133, 139, 140, 196, 201, 241, 259, 265, 266, 288, 322, 397, 405, 437, 455, 461, 462, 484, 550, 610, 685, 693, 725, 726, 743, 749, 750, 889, 936, 955, 1015, 1090, 1130, 1148, 1154, 1155, 1183, 1243, 1276, 1439
Number of compositions (ordered partitions) of n into pentagonal pyramidal numbers ( A002411).
+20
3
1, 1, 1, 1, 1, 1, 2, 3, 4, 5, 6, 7, 9, 12, 16, 21, 27, 34, 44, 57, 74, 96, 124, 159, 205, 265, 343, 444, 574, 740, 954, 1231, 1590, 2055, 2656, 3430, 4428, 5716, 7380, 9531, 12312, 15902, 20536, 26518, 34242, 44218, 57106, 73751, 95245, 122999, 158837, 205117
FORMULA
G.f.: 1/(1 - Sum_{k>=1} x^(k^2*(k+1)/2)).
MATHEMATICA
nmax = 51; CoefficientList[Series[1/(1 - Sum[x^(k^2 (k + 1)/2), {k, 1, nmax}]), {x, 0, nmax}], x]
Indices of pentagonal pyramidal numbers which are the sum of two other such numbers: k such that A002411(k) = A002411(i)+ A002411(j) for some i,j>0.
+20
2
42, 82, 171, 411, 522, 886, 1042, 1066, 1183, 1341, 2596, 2842, 3480, 3831, 5012, 5316, 5360, 5786, 6219, 6650, 6924, 8406, 8666, 9408, 10707, 11735, 12590, 12891, 14422, 14646, 14826, 17351, 17702, 17757, 18882, 23210, 24108, 25127, 28175, 31980, 32400
COMMENTS
The numbers themselves are listed in A172425. This answers a question posed in A133459.
EXAMPLE
a(1)=42 because A002411(42) is the smallest term in that sequence which is the sum of two other (nonzero) terms of A002411.
PROG
(PARI) for(n=1, 99999, for(m=1, n-1, isA002411(p(n)-p(m)) & !print1(n", ") & break)) /* needs isA002411() and p() */
Triangle T(n,k) = A002411(k) read by rows.
+20
1
1, 1, 6, 1, 6, 18, 1, 6, 18, 40, 1, 6, 18, 40, 75, 1, 6, 18, 40, 75, 126, 1, 6, 18, 40, 75, 126, 196, 1, 6, 18, 40, 75, 126, 196, 288
FORMULA
Equals the matrix product A002260 * A127773 as infinite lower triangular matrices.
EXAMPLE
First few rows of the triangle are:
1;
1, 6;
1, 6, 18;
1, 6, 18, 40;
1, 6, 18, 40, 75;
...
Pentagonal numbers: a(n) = n*(3*n-1)/2.
(Formerly M3818 N1562)
+10
541
0, 1, 5, 12, 22, 35, 51, 70, 92, 117, 145, 176, 210, 247, 287, 330, 376, 425, 477, 532, 590, 651, 715, 782, 852, 925, 1001, 1080, 1162, 1247, 1335, 1426, 1520, 1617, 1717, 1820, 1926, 2035, 2147, 2262, 2380, 2501, 2625, 2752, 2882, 3015, 3151
COMMENTS
The average of the first n (n > 0) pentagonal numbers is the n-th triangular number. - Mario Catalani (mario.catalani(AT)unito.it), Apr 10 2003
a(n) is the sum of n integers starting from n, i.e., 1, 2 + 3, 3 + 4 + 5, 4 + 5 + 6 + 7, etc. - Jon Perry, Jan 15 2004 Partial sums of 1, 4, 7, 10, 13, 16, ... (1 mod 3), a(2k) = k(6k-1), a(2k-1) = (2k-1)(3k-2). - Jon Perry, Sep 10 2004 Starting with offset 1 = binomial transform of [1, 4, 3, 0, 0, 0, ...]. Also, A004736 * [1, 3, 3, 3, ...]. - Gary W. Adamson, Oct 25 2007 If Y is a 3-subset of an n-set X then, for n >= 4, a(n-3) is the number of 4-subsets of X having at least two elements in common with Y. - Milan Janjic, Nov 23 2007 Solutions to the duplication formula 2*a(n) = a(k) are given by the index pairs (n, k) = (5,7), (5577, 7887), (6435661, 9101399), etc. The indices are integer solutions to the pair of equations 2(6n-1)^2 = 1 + y^2, k = (1+y)/6, so these n can be generated from the subset of numbers [1+ A001653(i)]/6, any i, where these are integers, confined to the cases where the associated k=[1+ A002315(i)]/6 are also integers. - R. J. Mathar, Feb 01 2008 Even octagonal numbers divided by 8. - Omar E. Pol, Aug 18 2011 Sequence found by reading the line from 0, in the direction 0, 5, ... and the line from 1, in the direction 1, 12, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. - Omar E. Pol, Sep 08 2011 The hyper-Wiener index of the star-tree with n edges (see A196060, example). - Emeric Deutsch, Sep 30 2011 More generally the n-th k-gonal number is equal to n + (k-2)* A000217(n-1), n >= 1, k >= 3. In this case k = 5. - Omar E. Pol, Apr 06 2013 Note that both Euler's pentagonal theorem for the partition numbers and Euler's pentagonal theorem for the sum of divisors refer more exactly to the generalized pentagonal numbers, not this sequence. For more information see A001318, A175003, A238442. - Omar E. Pol, Mar 01 2014 The Fuss-Catalan numbers are Cat(d,k)= [1/(k*(d-1)+1)]*binomial(k*d,k) and enumerate the number of (d+1)-gon partitions of a (k*(d-1)+2)-gon (cf. Schuetz and Whieldon link). a(n)= Cat(n,3), so enumerates the number of (n+1)-gon partitions of a (3*(n-1)+2)-gon. Analogous sequences are A100157 (k=4) and A234043 (k=5). - Tom Copeland, Oct 05 2014 Binomial transform of (0, 1, 3, 0, 0, 0, ...) ( A169585 with offset 1) and second partial sum of (0, 1, 3, 3, 3, ...). - Gary W. Adamson, Oct 05 2015 For n > 0, a(n) is the number of compositions of n+8 into n parts avoiding parts 2 and 3. - Milan Janjic, Jan 07 2016 a(n) is also the number of edges in the Mycielskian of the complete graph K[n]. Indeed, K[n] has n vertices and n(n-1)/2 edges. Then its Mycielskian has n + 3n(n-1)/2 = n(3n-1)/2. See p. 205 of the West reference. - Emeric Deutsch, Nov 04 2016 Also the number of maximal cliques in the n-Andrásfai graph. - Eric W. Weisstein, Dec 01 2017 Coefficients in the hypergeometric series identity 1 - 5*(x - 1)/(2*x + 1) + 12*(x - 1)*(x - 2)/((2*x + 1)*(2*x + 2)) - 22*(x - 1)*(x - 2)*(x - 3)/((2*x + 1)*(2*x + 2)*(2*x + 3)) + ... = 0, valid for Re(x) > 1. Cf. A002412 and A002418. Column 2 of A103450. - Peter Bala, Mar 14 2019 A generalization of the Comment dated Apr 10 2003 follows. (k-3)* A000292(n-2) plus the average of the first n (2k-1)-gonal numbers is the n-th k-gonal number. - Charlie Marion, Nov 01 2020 a(n+1) is the number of Dyck paths of size (3,3n+1); i.e., the number of NE lattice paths from (0,0) to (3,3n+1) which stay above the line connecting these points. - Harry Richman, Jul 13 2021
REFERENCES
Tom M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, pages 2 and 311.
Raymond Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; p. 129.
Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6.
L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 1.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 284.
Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 64.
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).
André Weil, Number theory: an approach through history; from Hammurapi to Legendre, Birkhäuser, Boston, 1984; see p. 186.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, 1987, pp. 98-100.
Douglas B. West, Introduction to Graph Theory, 2nd ed., Prentice-Hall, NJ, 2001.
LINKS
Michel Waldschmidt, Continued fractions, Ecole de recherche CIMPA-Oujda, Théorie des Nombres et ses Applications, 18 - 29 mai 2015: Oujda (Maroc).
FORMULA
Product_{m > 0} (1 - q^m) = Sum_{k} (-1)^k*x^a(k). - Paul Barry, Jul 20 2003 G.f.: x*(1+2*x)/(1-x)^3.
E.g.f.: exp(x)*(x+3*x^2/2).
a(n) = n*(3*n-1)/2.
a(n) = binomial(3*n, 2)/3. - Paul Barry, Jul 20 2003 a(0) = 0, a(1) = 1; for n >= 2, a(n) = 2*a(n-1) - a(n-2) + 3. - Miklos Kristof, Mar 09 2005 a(n) = Sum_{k=1..n} (2*n - k). - Paul Barry, Aug 19 2005 a(n) = binomial(n+1, 2) + 2*binomial(n, 2).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0) = 0, a(1) = 1, a(2) = 5. - Jaume Oliver Lafont, Dec 02 2008 a(n) = a(n+9) - A016777(n)*9. (End) a(3*a(n) + 4*n + 1) = a(3*a(n) + 4*n) + a(3*n+1).
A generalization. Let {G_k(n)}_(n >= 0) be sequence of k-gonal numbers (k >= 3). Then the following identity holds: G_k((k-2)*G_k(n) + c(k-3)*n + 1) = G_k((k-2)*G_k(n) + c(k-3)*n) + G_k((k-2)*n + 1), where c = A000124. (End) Sum_{n>=1} (-1)^(n+1)/a(n) = 2*(sqrt(3)*Pi - 6*log(2))/3 = 0.85501000622865446... - Ilya Gutkovskiy, Jul 28 2016 In general, let P(k,n) be the n-th k-gonal number. Then P(k,m+n) = P(k,m) + (k-2)mn + P(k,n). - Charlie Marion, Apr 16 2017 Product_{n>=2} (1 - 1/a(n)) = 3/5. - Amiram Eldar, Jan 21 2021 (n+1)*(a(n^2) + a(n^2+1) + ... + a(n^2+n)) = n*(a(n^2+n+1) + ... + a(n^2+2n)). - Charlie Marion, Apr 28 2024 a(n) = Sum_{k = 0..3*n} (-1)^(n+k+1) * binomial(k, 2)*binomial(3*n+k-1, 2*k). - Peter Bala, Nov 04 2024
EXAMPLE
Illustration of initial terms:
.
. o
. o o
. o o o o
. o o o o o o
. o o o o o o o o o
. o o o o o o o o o o o
. o o o o o o o o o o o o o
. o o o o o o o o o o o o o o
. o o o o o o o o o o o o o o o
.
. 1 5 12 22 35
MATHEMATICA
LinearRecurrence[{3, -3, 1}, {0, 1, 5}, 61] (* Harvey P. Dale, Dec 27 2011 *) pentQ[n_] := IntegerQ[(1 + Sqrt[24 n + 1])/6]; pentQ[0] = True; Select[Range[0, 3200], pentQ@# &] (* Robert G. Wilson v, Mar 31 2014 *) Join[{0}, Accumulate[Range[1, 312, 3]]] (* Harvey P. Dale, Mar 26 2016 *) (* For Mathematica 10.4+ *) Table[PolygonalNumber[RegularPolygon[5], n], {n, 0, 46}] (* Arkadiusz Wesolowski, Aug 27 2016 *) CoefficientList[Series[x (-1 - 2 x)/(-1 + x)^3, {x, 0, 20}], x] (* Eric W. Weisstein, Dec 01 2017 *)
PROG
(PARI) a(n)=n*(3*n-1)/2
(PARI) vector(100, n, n--; binomial(3*n, 2)/3) \\ Altug Alkan, Oct 06 2015 (PARI) is_a000326 = my(s); n==0 || (issquare (24*n+1, &s) && s%6==5) \\ Hugo Pfoertner, Aug 03 2023 (Haskell)
(Python) # Intended to compute the initial segment of the sequence, not isolated terms.
def aList():
x, y = 1, 1
yield 0
while True:
yield x
x, y = x + y + 3, y + 3
CROSSREFS
The generalized pentagonal numbers b*n+3*n*(n-1)/2, for b = 1 through 12, form sequences A000326, A005449, A045943, A115067, A140090, A140091, A059845, A140672, A140673, A140674, A140675, A151542. Cf. A001318 (generalized pentagonal numbers), A049452, A033570, A010815, A034856, A051340, A004736, A033568, A049453, A002411 (partial sums), A033579. See A220083 for a list of numbers of the form n*P(s,n)-(n-1)*P(s,n-1), where P(s,n) is the n-th polygonal number with s sides. Cf. A240137: sum of n consecutive cubes starting from n^3. Cf. similar sequences listed in A022288.
KEYWORD
core,nonn,easy,nice,changed
Square pyramidal numbers: a(n) = 0^2 + 1^2 + 2^2 + ... + n^2 = n*(n+1)*(2*n+1)/6.
(Formerly M3844 N1574)
+10
507
0, 1, 5, 14, 30, 55, 91, 140, 204, 285, 385, 506, 650, 819, 1015, 1240, 1496, 1785, 2109, 2470, 2870, 3311, 3795, 4324, 4900, 5525, 6201, 6930, 7714, 8555, 9455, 10416, 11440, 12529, 13685, 14910, 16206, 17575, 19019, 20540, 22140, 23821, 25585, 27434, 29370
COMMENTS
The sequence contains exactly one square greater than 1, namely 4900 (according to Gardner). - Jud McCranie, Mar 19 2001, Mar 22 2007 [This is a result from Watson. - Charles R Greathouse IV, Jun 21 2013] [See A351830 for further related comments and references.] Number of rhombi in an n X n rhombus. - Matti De Craene (Matti.DeCraene(AT)rug.ac.be), May 14 2000
Number of acute triangles made from the vertices of a regular n-polygon when n is odd (cf. A007290). - Sen-Peng Eu, Apr 05 2001 Gives number of squares formed from an n X n square. In a 1 X 1 square, one is formed. In a 2 X 2 square, five squares are formed. In a 3 X 3 square, 14 squares are formed and so on. - Kristie Smith (kristie10spud(AT)hotmail.com), Apr 16 2002
a(n-1) = B_3(n)/3, where B_3(x) = x(x-1)(x-1/2) is the third Bernoulli polynomial. - Michael Somos, Mar 13 2004 Number of permutations avoiding 13-2 that contain the pattern 32-1 exactly once.
Since 3*r = (r+1) + r + (r-1) = T(r+1) - T(r-2), where T(r) = r-th triangular number r*(r+1)/2, we have 3*r^2 = r*(T(r+1) - T(r-2)) = f(r+1) - f(r-1) ... (i), where f(r) = (r-1)*T(r) = (r+1)*T(r-1). Summing over n, the right hand side of relation (i) telescopes to f(n+1) + f(n) = T(n)*((n+2) + (n-1)), whence the result Sum_{r=1..n} r^2 = n*(n+1)*(2*n+1)/6 immediately follows. - Lekraj Beedassy, Aug 06 2004 Also as a(n) = (1/6)*(2*n^3 + 3*n^2 + n), n > 0: structured trigonal diamond numbers (vertex structure 5) (cf. A006003 = alternate vertex; A000447 = structured diamonds; A100145 for more on structured numbers). - James A. Record (james.record(AT)gmail.com), Nov 07 2004 Number of triples of integers from {1, 2, ..., n} whose last component is greater than or equal to the others.
Sum of the first n positive squares. - Cino Hilliard, Jun 18 2007 Maximal number of cubes of side 1 in a right pyramid with a square base of side n and height n. - Pasquale CUTOLO (p.cutolo(AT)inwind.it), Jul 09 2007
If a 2-set Y and an (n-2)-set Z are disjoint subsets of an n-set X then a(n-3) is the number of 4-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007 We also have the identity 1 + (1+4) + (1+4+9) + ... + (1+4+9+16+ ... + n^2) = n(n+1)(n+2)(n+(n+1)+(n+2))/36; ... and in general the k-fold nested sum of squares can be expressed as n(n+1)...(n+k)(n+(n+1)+...+(n+k))/((k+2)!(k+1)/2). - Alexander R. Povolotsky, Nov 21 2007 The terms of this sequence are coefficients of the Engel expansion of the following converging sum: 1/(1^2) + (1/1^2)*(1/(1^2+2^2)) + (1/1^2)*(1/(1^2+2^2))*(1/(1^2+2^2+3^2)) + ... - Alexander R. Povolotsky, Dec 10 2007 Hankel transform of binomial(2*n-3, n-1) is -a(n). - Paul Barry, Feb 12 2008 Starting (1, 5, 14, 30, ...) = binomial transform of [1, 4, 5, 2, 0, 0, 0, ...]. - Gary W. Adamson, Jun 13 2008 Starting (1,5,14,30,...) = second partial sums of binomial transform of [1,2,0,0,0,...]. a(n) = Sum_{i=0..n} binomial(n+2,i+2)*b(i), where b(i)=1,2,0,0,0,... - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
Sequence of the absolute values of the z^1 coefficients of the polynomials in the GF1 denominators of A156921. See A157702 for background information. - Johannes W. Meijer, Mar 07 2009 The sequence is related to A000217 by a(n) = n* A000217(n) - Sum_{i=0..n-1} A000217(i) and this is the case d = 1 in the identity n^2*(d*n-d+2)/2 - Sum_{i=0..n-1} i*(d*i-d+2)/2 = n*(n+1)(2*d*n-2*d+3)/6, or also the case d = 0 in n^2*(n+2*d+1)/2 - Sum_{i=0..n-1} i*(i+2*d+1)/2 = n*(n+1)*(2*n+3*d+1)/6. - Bruno Berselli, Apr 21 2010, Apr 03 2012 a(n)/n = k^2 (k = integer) for n = 337; a(337) = 12814425, a(n)/n = 38025, k = 195, i.e., the number k = 195 is the quadratic mean (root mean square) of the first 337 positive integers. There are other such numbers -- see A084231 and A084232. - Jaroslav Krizek, May 23 2010 Also the number of moves to solve the "alternate coins game": given 2n+1 coins (n+1 Black, n White) set alternately in a row (BWBW...BWB) translate (not rotate) a pair of adjacent coins at a time (1 B and 1 W) so that at the end the arrangement shall be BBBBB..BW...WWWWW (Blacks separated by Whites). Isolated coins cannot be moved. - Carmine Suriano, Sep 10 2010 Using four consecutive numbers n, n+1, n+2, and n+3 take all possible pairs (n, n+1), (n, n+2), (n, n+3), (n+1, n+2), (n+1, n+3), (n+2, n+3) to create unreduced Pythagorean triangles. The sum of all six areas is 60*a(n+1).
Using three consecutive odd numbers j, k, m, (j+k+m)^3 - (j^3 + k^3 + m^3) equals 576*a(n) = 24^2*a(n) where n = (j+1)/2. (End)
For n > 0, the digital roots of this sequence A010888(a(n)) form the purely periodic 27-cycle {1, 5, 5, 3, 1, 1, 5, 6, 6, 7, 2, 2, 9, 7, 7, 2, 3, 3, 4, 8, 8, 6, 4, 4, 8, 9, 9}. For n > 0, the units' digits of this sequence A010879(a(n)) form the purely periodic 20-cycle {1, 5, 4, 0, 5, 1, 0, 4, 5, 5, 6, 0, 9, 5, 0, 6, 5, 9, 0, 0}. (End) Length of the Pisano period of this sequence mod n, n>=1: 1, 4, 9, 8, 5, 36, 7, 16, 27, 20, 11, 72, 13, 28, 45, 32, 17, 108, 19, 40, ... . - R. J. Mathar, Oct 17 2012 The number of intersections of diagonals in the interior of regular n-gon for odd n > 1 divided by n is a square pyramidal number; that is, A006561(2*n+1)/(2*n+1) = A000330(n-1) = (1/6)*n*(n-1)*(2*n-1). - Martin Renner, Mar 06 2013 A formula for the r-th successive summation of k^2, for k = 1 to n, is (2*n+r)*(n+r)!/((r+2)!*(n-1)!) (H. W. Gould). - Gary Detlefs, Jan 02 2014 The n-th square pyramidal number = the n-th triangular dipyramidal number (Johnson 12), which is the sum of the n-th + (n-1)-st tetrahedral numbers. E.g., the 3rd tetrahedral number is 10 = 1+3+6, the 2nd is 4 = 1+3. In triangular "dipyramidal form" these numbers can be written as 1+3+6+3+1 = 14. For "square pyramidal form", rebracket as 1+(1+3)+(3+6) = 14. - John F. Richardson, Mar 27 2014 Beukers and Top prove that no square pyramidal number > 1 equals a tetrahedral number A000292. - Jonathan Sondow, Jun 21 2014 Odd numbered entries are related to dissections of polygons through A100157. - Tom Copeland, Oct 05 2014 We construct a number triangle from the integers 1, 2, 3, ..., n as follows. The first column contains 2*n-1 integers 1. The second column contains 2*n-3 integers 2, ... The last column contains only one integer n. The sum of all the numbers in the triangle is a(n).
Here is an example with n = 5:
1
1 2
1 2 3
1 2 3 4
1 2 3 4 5
1 2 3 4
1 2 3
1 2
1
(End)
The Catalan number series A000108(n+3), offset 0, gives Hankel transform revealing the square pyramidal numbers starting at 5, A000330(n+2), offset 0 (empirical observation). - Tony Foster III, Sep 05 2016; see Dougherty et al. link p. 2. - Andrey Zabolotskiy, Oct 13 2016 Number of floating point additions in the factorization of an (n+1) X (n+1) real matrix by Gaussian elimination as e.g. implemented in LINPACK subroutines sgefa.f or dgefa.f. The number of multiplications is given by A007290. - Hugo Pfoertner, Mar 28 2018 The Jacobi polynomial P(n-1,-n+2,2,3) or equivalently the sum of dot products of vectors from the first n rows of Pascal's triangle ( A007318) with the up-diagonal Chebyshev T coefficient vector (1,3,2,0,...) ( A053120) or down-diagonal vector (1,-7,32,-120,400,...) ( A001794). a(5) = 1 + (1,1).(1,3) + (1,2,1).(1,3,2) + (1,3,3,1).(1,3,2,0) + (1,4,6,4,1).(1,3,2,0,0) = (1 + (1,1).(1,-7) + (1,2,1).(1,-7,32) + (1,3,3,1).(1,-7,32,-120) + (1,4,6,4,1).(1,-7,32,-120,400))*(-1)^(n-1) = 55. - Richard Turk, Jul 03 2018 Coefficients in the terminating series identity 1 - 5*n/(n + 4) + 14*n*(n - 1)/((n + 4)*(n + 5)) - 30*n*(n - 1)*(n - 2)/((n + 4)*(n + 5)*(n + 6)) + ... = 0 for n = 1,2,3,.... Cf. A002415 and A108674. - Peter Bala, Feb 12 2019 n divides a(n) iff n == +- 1 (mod 6) (see A007310). (See De Koninck reference.) Examples: a(11) = 506 = 11 * 46, and a(13) = 819 = 13 * 63. - Bernard Schott, Jan 10 2020 For n > 0, a(n) is the number of ternary words of length n+2 having 3 letters equal to 2 and 0 only occurring as the last letter. For example, for n=2, the length 4 words are 2221,2212,2122,1222,2220. - Milan Janjic, Jan 28 2020 Conjecture: Every integer can be represented as a sum of three generalized square pyramidal numbers. A related conjecture is given in A336205 corresponding to pentagonal case. A stronger version of these conjectures is that every integer can be expressed as a sum of three generalized r-gonal pyramidal numbers for all r >= 3. In here "generalized" means negative indices are included. - Altug Alkan, Jul 30 2020
REFERENCES
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 813.
A. H. Beiler, Recreations in the Theory of Numbers, Dover Publications, NY, 1964, p. 194.
A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 215,223.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 122, see #19 (3(1)), I(n); p. 155.
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.
S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (p.165).
J. M. De Koninck and A. Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problème 310, pp. 46-196, Ellipses, Paris, 2004.
E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93.
L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 2.
M. Gardner, Fractal Music, Hypercards and More, Freeman, NY, 1991, p. 293.
M. Holt, Math puzzles and games, Walker Publishing Company, 1977, p. 2 and p. 89.
Simon Singh, The Simpsons and Their Mathematical Secrets. London: Bloomsbury Publishing PLC (2013): 188.
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).
David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 126.
LINKS
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy]. Bruno Berselli, A description of the recursive method in Comments lines: website Matem@ticamente (in Italian). S. Butler and P. Karasik, A note on nested sums, J. Int. Seq. 13 (2010), 10.4.4, p=1 in first displayed equation page 4. Milan Janjic and B. Petkovic, A Counting Function, arXiv preprint arXiv:1301.4550 [math.CO], 2013. Eric Weisstein's World of Mathematics, PowerSum.
FORMULA
G.f.: x*(1+x)/(1-x)^4.
E.g.f.: (x + 3*x^2/2 + x^3/3)*exp(x).
a(n) = n*(n+1)*(2*n+1)/6 = binomial(n+2, 3) + binomial(n+1, 3).
Can be extended to Z with a(n) = -a(-1-n) for all n in Z.
a(n) = binomial(2*(n+1), 3)/4. - Paul Barry, Jul 19 2003 a(n) = (((n+1)^4 - n^4) - ((n+1)^2 - n^2))/12. - Xavier Acloque, Oct 16 2003
a(n) = sqrt(Sum_{j=1..n} Sum_{i=1..n} (i*j)^2).
a(n) = (Sum_{k=1..n} Sum_{j=1..n} Sum_{i=1..n} (i*j*k)^2)^(1/3). (End)
a(n) = Sum_{i=1..n} i*(2*n-2*i+1); sum of squares gives 1 + (1+3) + (1+3+5) + ... - Jon Perry, Dec 08 2004 Sum_{n>=1} 1/a(n) = 6*(3-4*log(2)); Sum_{n>=1} (-1)^(n+1)*1/a(n) = 6*(Pi-3). - Philippe Deléham, May 31 2005 Euler transform of length-2 sequence [ 5, -1 ]. - Michael Somos, Sep 04 2006 a(n) = binomial(n, 2) + 2*binomial(n, 3). - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009, corrected by M. F. Hasler, Jan 02 2024 a(n) = s(n+1, n)^2 - 2*s(n+1, n-1), where s(n, k) are Stirling numbers of the first kind, A048994. - Mircea Merca, Apr 03 2012 a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 2. - Ant King, Oct 17 2012 a(n) = Sum_{i=0..n-1} (n-i)*(2*i+1), with a(0) = 0. After 0, row sums of the triangle in A101447. - Bruno Berselli, Feb 10 2014 a(n) = (2*n^3 + 3*n^2 + n)/6, see Singh (2013). - Alonso del Arte, Feb 20 2015 Dirichlet g.f.: zeta(s-3)/3 + zeta(s-2)/2 + zeta(s-1)/6. - Ilya Gutkovskiy, Jun 26 2016 a(n) = binomial(n-1, 1) + binomial(n-1, 2) + binomial(n, 3) + binomial(n+1, 2) + binomial(n+1, 3). - Tony Foster III, Aug 24 2018 a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Nathan Fox, Dec 04 2019 Let T(n) = A000217(n), the n-th triangular number. Then a(n) = (T(n)+1)^2 + (T(n)+2)^2 + ... + (T(n)+n)^2 - (n+2)*T(n)^2. - Charlie Marion, Dec 31 2019 a(n) = 2*n - 1 - a(n-2) + 2*a(n-1). - Boštjan Gec, Nov 09 2023
EXAMPLE
G.f. = x + 5*x^2 + 14*x^3 + 30*x^4 + 55*x^5 + 91*x^6 + 140*x^7 + 204*x^8 + ...
MAPLE
A000330 := (1+z)/(z-1)^4; # Simon Plouffe (in his 1992 dissertation: generating function for sequence starting at a(1))
with(combstruct): ZL:=[st, {st=Prod(left, right), left=Set(U, card=r), right=Set(U, card=r), U=Sequence(Z, card>=1)}, unlabeled]: subs(r=1, stack): seq(count(subs(r=2, ZL), size=m*2), m=1..45) ; # Zerinvary Lajos, Jan 02 2008 a := n -> sum(k^2, k=1..n):seq(a(n), n=0...44); # Zerinvary Lajos, Jun 15 2008 nmax := 44; for n from 0 to nmax do fz(n) := product( (1-(2*m-1)*z)^(n+1-m) , m=1..n); c(n) := abs(coeff(fz(n), z, 1)); end do: a := n-> c(n): seq(a(n), n=0..nmax); # Johannes W. Meijer, Mar 07 2009
MATHEMATICA
Table[Binomial[w+2, 3] + Binomial[w+1, 3], {w, 0, 30}]
CoefficientList[Series[x(1+x)/(1-x)^4, {x, 0, 40}], x] (* Vincenzo Librandi, Jul 30 2014 *)
PROG
(PARI) {a(n) = n * (n+1) * (2*n+1) / 6};
(Haskell)
a000330 n = n * (n + 1) * (2 * n + 1) `div` 6
a000330_list = scanl1 (+) a000290_list
(Maxima) A000330(n):=binomial(n+2, 3)+binomial(n+1, 3)$ (Magma) [0] cat [((2*n+3)*Binomial(n+2, 2))/3: n in [0..40]]; // Vincenzo Librandi, Jul 30 2014 (Python) a=lambda n: (n*(n+1)*(2*n+1))//6 # Indranil Ghosh, Jan 04 2017 (Sage) [n*(n+1)*(2*n+1)/6 for n in (0..30)] # G. C. Greubel, Dec 31 2019 (GAP) List([0..30], n-> n*(n+1)*(2*n+1)/6); # G. C. Greubel, Dec 31 2019
CROSSREFS
Cf. A000217, A000292, A000537, A005408, A006003, A006331, A033994, A033999, A046092, A050409, A050446, A050447, A100157, A132124, A132112, A156921, A157702, A258708, A351830. Sums of 2 consecutive terms give A005900. Cf. A253903 (characteristic function). Cf. A034705 (differences of any two terms).
KEYWORD
nonn,easy,core,nice,changed
Hexagonal pyramidal numbers, or greengrocer's numbers.
(Formerly M4374 N1839)
+10
94
0, 1, 7, 22, 50, 95, 161, 252, 372, 525, 715, 946, 1222, 1547, 1925, 2360, 2856, 3417, 4047, 4750, 5530, 6391, 7337, 8372, 9500, 10725, 12051, 13482, 15022, 16675, 18445, 20336, 22352, 24497, 26775, 29190, 31746, 34447, 37297, 40300
COMMENTS
Binomial transform of (1, 6, 9, 4, 0, 0, 0, ...). - Gary W. Adamson, Oct 16 2007 a(n) is the sum of the maximum(m,n) over {(m,n): m,n in positive integers, m<=n}. - Geoffrey Critzer, Oct 11 2009 We obtain these numbers for d=2 in the identity n*(n*(d*n-d+2)/2)-sum(k*(d*k-d+2)/2, k=0..n-1) = n*(n+1)*(2*d*n-2*d+3)/6 (see Klaus Strassburger in Formula lines). - Bruno Berselli, Apr 21 2010, Nov 16 2010 q^a(n) is the Hankel transform of the q-Catalan numbers. - Paul Barry, Dec 15 2010 For n>0, the digital roots of this sequence A010888( A002412(n)) form the purely periodic 27-cycle {1,7,4,5,5,8,9,3,3,4,1,7,8,8,2,3,6,6,7,4,1,2,2,5,6,9,9}. For n>0, the units' digits of this sequence A010879( A002412(n)) form the purely periodic 20-cycle {1,7,2,0,5,1,2,2,5,5,6,2,7,5,0,6,7,7,0,0}. (End)
Number of orbits of Aut(Z^7) as function of the infinity norm (n+1) of the representative integer lattice point of the orbit, when the cardinality of the orbit is equal to 40320. - Philippe A.J.G. Chevalier, Dec 28 2015 Coefficients in the hypergeometric series identity 1 - 7*(x - 1)/(3*x + 1) + 22*(x - 1)*(x - 2)/((3*x + 1)*(3*x + 2)) - 50*(x - 1)*(x - 2)*(x - 3)/((3*x + 1)*(3*x + 2)*(3*x + 3)) + ... = 0, valid for Re(x) > 1. Cf. A000326 and A002418. Column 3 of A103450. - Peter Bala, Mar 14 2019
REFERENCES
A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 194.
E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93.
L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 2.
T. A. Gulliver, Sequences from Arrays of Integers, Int. Math. Journal, Vol. 1, No. 4, pp. 323-332, 2002.
I. Siap, Linear codes over F_2 + u*F_2 and their complete weight enumerators, in Codes and Designs (Ohio State, May 18, 2000), pp. 259-271. De Gruyter, 2002.
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).
LINKS
Bruno Berselli, A description of the transform in Comments lines: website Matem@ticamente (in Italian).
FORMULA
a(n) = n(n + 1)(4n - 1)/6.
G.f.: x*(1+3*x)/(1-x)^4. - Simon Plouffe in his 1992 dissertation. a(n) = n^3 - Sum_{i=1..n-1} i^2. - Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de)
Partial sums of n odd-indexed triangular numbers, e.g., a(3) = t(1)+t(3)+t(5) = 1+6+15 = 22. - Jon Perry, Jul 23 2003 a(n) = Sum_{i=0..n-1} (n - i)*(n + i). - Jon Perry, Sep 26 2004 a(n) = n* A000292(n) - (n-1)* A000292(n-1) = n*binomial((n+2),3) - (n-1)*binomial((n+1),3); e.g., a(5) = 95 = 5*35 - 4*20. - Gary W. Adamson, Dec 28 2007 a(n) = Sum_{i=0..n} (2i^2 + 3i + 1), for n >= 0 (Omits the leading 0). - William A. Tedeschi, Aug 25 2010 a(n) = 4*a(n-1)-6*a(n-2)+4*a(n-3)-a(n-4), with a(0)=0, a(1)=1, a(2)=7, a(3)=22. - Harvey P. Dale, Jul 16 2011 a(n) = sum a*b, where the summing is over all unordered partitions 2*n = a+b. - Vladimir Shevelev, May 11 2012 a(n) = a(n-1) + n*(2*n-1).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 4.
a(n) = binomial(n+2,3) + 3*binomial(n+1,3) = (4*n-1)*binomial(n+1,2)/3.
Sum_{n>=1} 1/a(n) = 6*(12*log(2)-2*Pi-1)/5 = 1.2414...
(End)
Sum_{n>=1} (-1)^(n+1)/a(n) = 6*(1 + 2*sqrt(2)*Pi - 2*(3+sqrt(2))*log(2) + 4*sqrt(2)*log(2-sqrt(2)))/5. - Amiram Eldar, Jan 04 2022
EXAMPLE
Let n=5, 2*n=10. Since 10 = 1 + 9 = 2 + 8 = 3 + 7 = 4 + 6 = 5 + 5, a(5) = 1*9 + 2*8 + 3*7 + 4*6 + 5*5 = 95. - Vladimir Shevelev, May 11 2012
MATHEMATICA
Figurate[ ngon_, rank_, dim_] := Binomial[rank + dim - 2, dim - 1] ((rank - 1)*(ngon - 2) + dim)/dim; Table[ Figurate[6, r, 3], {r, 0, 40}] (* Robert G. Wilson v, Aug 22 2010 *) Table[n(n+1)(4n-1)/6, {n, 0, 40}] (* or *) LinearRecurrence[{4, -6, 4, -1}, {0, 1, 7, 22}, 40] (* Harvey P. Dale, Jul 16 2011 *)
PROG
(PARI) v=vector(40, i, (i*(i+1))\2); s=0; print1(s", "); forstep(i=1, 40, 2, s+=v[i]; print1(s", "))
(GAP) List([0..40], n->n*(n+1)*(4*n-1)/6); # Muniru A Asiru, Mar 18 2019
CROSSREFS
a(n) = A093561(n+2, 3), (4, 1)-Pascal column. Cf. A220084 for a list of numbers of the form n*P(k,n)-(n-1)*P(k,n-1), where P(k,n) is the n-th k-gonal pyramidal number (see Adamson's formula). Cf. similar sequences listed in A237616.
0, 2, 12, 36, 80, 150, 252, 392, 576, 810, 1100, 1452, 1872, 2366, 2940, 3600, 4352, 5202, 6156, 7220, 8400, 9702, 11132, 12696, 14400, 16250, 18252, 20412, 22736, 25230, 27900, 30752, 33792, 37026, 40460, 44100, 47952, 52022, 56316, 60840
COMMENTS
(1) a(n) = sum of second string of n triangular numbers - sum of first n triangular numbers, or the 2n-th partial sum of triangular numbers ( A000217) - the n-th partial sum of triangular numbers ( A000217). The same for natural numbers gives squares. (2) a(n) = (n-th triangular number)*(the n-th even number) = n(n+1)/2 * (2n). - Amarnath Murthy, Nov 05 2002 Let M(n) be the n X n matrix m(i,j)=1/(i+j+x), let P(n,x) = (Product_{i=0..n-1} i!^2)/det(M(n)). Then P(n,x) is a polynomial with integer coefficients of degree n^2 and a(n) is the coefficient of x^(n^2-1). - Benoit Cloitre, Jan 15 2003 Y values of solutions of the equation: (X-Y)^3-X*Y=0. X values are a(n)=n*(n+1)^2 (see A045991) - Mohamed Bouhamida, May 09 2006 a(2d-1) is the number of self-avoiding walk of length 3 in the d-dimensional hypercubic lattice. - Michael Somos, Sep 06 2006 This sequence is related to A005449 by a(n) = n* A005449(n)-sum( A005449(i), i=0..n-1), and this is the case d=3 in the identity n^2*(d*n+d-2)/2 - Sum_{k=0..n-1} k*(d*k+d-2)/2 = n*(n+d)*(2*d*n+d-3)/6. - Bruno Berselli, Nov 18 2010 Using (n, n+1) to generate a primitive Pythagorean triangle, the sides will be 2*n+1, 2*(n^2+n), and 2*n^2+2*n+1. Inscribing the largest rectangle with integral sides will have sides of length n and n^2+n. Side n is collinear to side 2*n+1 of the triangle and side n^2+n is collinear to side 2*(n^2+n) of the triangle. The areas of theses rectangles are a(n). - J. M. Bergot, Sep 22 2011 Define a 'stable brick triangle' as:
-----
| c |
---------
| a | | b |
----------
with a, b, c > 0 and c <= a + b. This can be visualized as two bricks with a third brick on top. The third brick can only be as strong as a+b, otherwise the wall collapses - for example, (1,2,4) is unstable.
a(n) gives the number of stable brick triangles that can be formed if the two supporting bricks are 1 <= a <= n and 1 <= b <= n: a(n) = Sum_{a=1..n} Sum_{b=1..n} Sum_c 1 = n^3 + n^2 as given in the Adamchuk formula.
So for i=j=n=2 we have 4:
1 2 3 4
2 2 2 2 2 2 2 2
For example, n=2 gives 2 from [a=1,b=1], 3 from both [a=1,b=2] and [a=2,b=1] and 4 from [a=2,b=2] so a(2) = 2 + 3 + 3 + 4 = 12. (End)
Define the infinite square array m(n,k) by m(n,k) = (n-k)^2 if n >= k >= 0 and by m(n,k) = (k+n)*(k-n) if 0 <= n <= k. This contains A120070 below the diagonal. Then a(n) = Sum_{k=0..n} m(n,k) + Sum_{r=0..n} m(r,n), the "hook sum" of the terms to the left of m(n,n) and above m(n,n) with irrelevant (vanishing) terms on the diagonal. - J. M. Bergot, Aug 16 2013 a(n) is the sum of all pairs with repetition drawn from the set of odd numbers 2*n-3. This is similar to A027480 but using the odd integers instead. Example using n=3 gives the odd numbers 1,3,5: 1+1, 1+3, 1+5, 3+3, 3+5,5+5 having a total of 36=a(3). - J. M. Bergot, Apr 05 2016 a(n) is the first Zagreb index of the complete graph K[n+1]. The first Zagreb index of a simple connected graph is the sum of the squared degrees of its vertices. Alternately, it is the sum of the degree sums d(i)+d(j) over all edges ij of the graph. - Emeric Deutsch, Nov 07 2016 a(n-2) is the maximum sigma irregularity over all trees with n vertices. The extremal graphs are stars. (The sigma irregularity of a graph is the sum of squares of the differences between the degrees over all edges of the graph.) - Allan Bickle, Jun 14 2023
REFERENCES
L. B. W. Jolley, "Summation of Series", Dover Publications, 1961, pp. 50, 64.
LINKS
Bruno Berselli, A description of the recursive method in Comments lines: website Matem@ticamente (in Italian).
FORMULA
a(n) = Sum_{j=1..n} (Sum_{i=1..n} (i+j)), row sums of A126890 skipping numbers in the first column. - Alexander Adamchuk, Oct 12 2004 Sum_{n>0} 1/a(n) = (Pi^2 - 6)/6 = 0.6449340... [Jolley eq 272] - Gary W. Adamson, Dec 22 2006 Sum_{n>=1} (-1)^(n+1)/a(n) = 1 + Pi^2/12 - 2*log(2). - Amiram Eldar, Jul 04 2020
MATHEMATICA
LinearRecurrence[{4, -6, 4, -1}, {0, 2, 12, 36}, 40] (* Harvey P. Dale, Sep 13 2018 *)
PROG
(Haskell)
(GAP) List([0..40], n-> n^2*(n+1) ); # G. C. Greubel, Aug 10 2019
CROSSREFS
Cf. A000217, A000290, A000292, A000330, A000578, A002378, A002411, A002412, A002413, A005449, A013661, A022549, A027480, A045991, A049450, A120070, A126890, A195437, A245334.
4-dimensional pyramidal numbers: a(n) = (3*n+1)*binomial(n+2, 3)/4. Also Stirling2(n+2, n).
(Formerly M4385 N1845)
+10
59
0, 1, 7, 25, 65, 140, 266, 462, 750, 1155, 1705, 2431, 3367, 4550, 6020, 7820, 9996, 12597, 15675, 19285, 23485, 28336, 33902, 40250, 47450, 55575, 64701, 74907, 86275, 98890, 112840, 128216, 145112, 163625, 183855, 205905, 229881, 255892, 284050, 314470
COMMENTS
Permutations avoiding 12-3 that contain the pattern 31-2 exactly once.
If Y is a 3-subset of an n-set X then, for n>=6, a(n-5) is the number of 6-subsets of X having at least two elements in common with Y. - Milan Janjic, Nov 23 2007 Starting with 1 = binomial transform of [1, 6, 12, 10, 3, 0, 0, 0, ...]. Equals row sums of triangle A143037. - Gary W. Adamson, Jul 18 2008 Rephrasing the Perry formula of 2003: a(n) is the sum of all products of all two numbers less than or equal to n, including the squares. Example: for n=3 the sum of these products is 1*1 + 1*2 + 1*3 + 2*2 + 2*3 + 3*3 = 25. - J. M. Bergot, Jul 16 2011 Half of the partial sums of A011379. [Jolley, Summation of Series, Dover (1961), page 12 eq (66).] - R. J. Mathar, Oct 03 2011 Also the number of (w,x,y,z) with all terms in {1,...,n+1} and w < x >= y > z (see A211795). - Clark Kimberling, May 19 2012 a(n-2) is the maximum number of intersections made from the perpendicular bisectors of all pair combinations of n points. - Ian Tam, Dec 22 2020
REFERENCES
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 835.
A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 195.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 227, #16.
S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (see p. 166, Table 10.4/I/3).
F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 223.
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).
LINKS
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy]. P. Aluffi, Degrees of projections of rank loci, arXiv:1408.1702 [math.AG], 2014. ["After compiling the results of many explicit computations, we noticed that many of the numbers d_{n,r,S} appear in the existing literature in contexts far removed from the enumerative geometry of rank conditions; we owe this surprising (to us) observation to perusal of [Slo14]."] C. Krishnamachaki, The operator (xD)^n, J. Indian Math. Soc., 15 (1923),3-4. [Annotated scanned copy]
FORMULA
a(n) = n*(1+n)*(2+n)*(1+3*n)/24. - T. D. Noe, Jan 21 2008 G.f.: x*(1+2*x)/(1-x)^5. - Paul Barry, Jul 23 2003 E.g.f. with offset -1: exp(x)*(1*(x^2)/2! + 4*(x^3)/3! + 3*(x^4)/4!). For the coefficients [1, 4, 3] see triangle A112493. E.g.f. x*exp(x)*(24 + 60*x + 28*x^2 + 3*x^3)/24 (above e.g.f. differentiated).
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) + 3. - Kieren MacMillan, Sep 29 2008 a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). - Jaume Oliver Lafont, Nov 23 2008 O.g.f. is D^2(x/(1-x)) = D^3(x), where D is the operator x/(1-x)*d/dx. - Peter Bala, Jul 02 2012 a(n) = 1*(1+2+...+n) + 2*(2+3+...+n) + ... + n*n. For example, a(6) = 266 = 1(1+2+3+4+5+6) + 2*(2+3+4+5+6) + 3*(3+4+5+6) + 4*(4+5+6) + 5*(5+6) + 6*(6).- J. M. Bergot, Apr 20 2017 Sum_{n>=1} 1/a(n) = (6/5) * (47 - 3*sqrt(3)*Pi - 27*log(3)).
Sum_{n>=1} (-1)^(n+1)/a(n) = (6/5) * (16*log(2) + 6*sqrt(3)*Pi - 43). (End)
EXAMPLE
G.f. = x + 7*x^2 + 25*x^3 + 65*x^4 + 140*x^5 + 266*x^6 + 462*x^7 + 750*x^8 + 1155*x^9 + ...
MAPLE
A001296:=-(1+2*z)/(z-1)**5; # Simon Plouffe in his 1992 dissertation for sequence without the leading zero
MATHEMATICA
Table[n*(1+n)*(2+n)*(1+3*n)/24, {n, 0, 100}]
CoefficientList[Series[x (1 + 2 x)/(1 - x)^5, {x, 0, 40}], x] (* Vincenzo Librandi, Jul 30 2014 *) Table[ListCorrelate[Accumulate[Range[n]], Range[n]], {n, 0, 40}]//Flatten (* or *) LinearRecurrence[{5, -10, 10, -5, 1}, {0, 1, 7, 25, 65}, 40] (* Harvey P. Dale, Aug 14 2017 *)
PROG
(PARI) t(n)=n*(n+1)/2
for(i=1, 30, print1(", "sum(j=1, i, j*t(j))))
(PARI) {a(n) = n * (n+1) * (n+2) * (3*n+1) / 24}; /* Michael Somos, Sep 04 2017 */ (Sage) [stirling_number2(n+2, n) for n in range(0, 38)] # Zerinvary Lajos, Mar 14 2009 (Magma) [(3*n+1)*Binomial(n+2, 3)/4: n in [0..40]]; // Vincenzo Librandi, Jul 30 2014
CROSSREFS
a(n)=f(n, 2) where f is given in A034261. a(n)= A093560(n+3, 4), (3, 1)-Pascal column. Cf. A220212 for a list of sequences produced by the convolution of the natural numbers with the k-gonal numbers.
Search completed in 0.087 seconds
|