A007585 -id:A007585 - cp's OEIS frontend

A051797 Partial sums of A007585.

1, 12, 50, 140, 315, 616, 1092, 1800, 2805, 4180, 6006, 8372, 11375, 15120, 19720, 25296, 31977, 39900, 49210, 60060, 72611, 87032, 103500, 122200, 143325, 167076, 193662, 223300, 256215, 292640, 332816, 376992, 425425, 478380, 536130

Offset: 0

Barry E. Williams, Dec 11 1999

a(n-1) is the n-th antidiagonal sum of the convolution array A213835. - Clark Kimberling, Jul 04 2012

Convolution of A000027 with A001107 (excluding 0). - Bruno Berselli, Dec 07 2012

Albert H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, pp. 194-196.
Murray R. Spiegel, Calculus of Finite Differences and Difference Equations, "Schaum's Outline Series", McGraw-Hill, 1971, pp. 10-20, 79-94.
Herbert John Ryser, Combinatorial Mathematics, "The Carus Mathematical Monographs", No. 14, John Wiley and Sons, 1963, pp. 1-8.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5,1).
Index to sequences related to pyramidal numbers.

Cf. A000027, A001107, A007585, A080851.
Cf. A093565 ((8, 1) Pascal, column m=4).
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers with the k-gonal numbers.

List([0..40], n-> (2*n+1)*Binomial(n+3,3)); # G. C. Greubel, Aug 30 2019

/* A000027 convolved with A001107 (excluding 0): */
A001107:=func; [&+[(n-i+1)*A001107(i): i in [1..n]]: n in [1..35]]; // Bruno Berselli, Dec 07 2012

[(2*n+1)*Binomial(n+3,3): n in [0..40]]; // G. C. Greubel, Aug 30 2019

seq((2*n+1)*binomial(n+3,3), n=0..40); # G. C. Greubel, Aug 30 2019

Table[(2*n+1)*Binomial[n+3,3], {n,0,40}] (*  Vladimir Joseph Stephan Orlovsky, Apr 19 2011, modified by G. C. Greubel, Aug 30 2019 *)

vector(40, n, (2*n-1)*binomial(n+2,3)) \\ G. C. Greubel, Aug 30 2019

[(2*n+1)*binomial(n+3,3) for n in (0..40)] # G. C. Greubel, Aug 30 2019

a(n) = binomial(n+3,3)*(2*n+1) = (n+1)*(n+2)*(n+3)*(2*n+1)/6.
G.f.: (1+7*x)/(1-x)^5.
a(n) = A080851(8,n). - R. J. Mathar, Jul 28 2016
E.g.f.: (6 + 66*x + 81*x^2 + 25*x^3 + 2*x^4)*exp(x)/6. - G. C. Greubel, Aug 30 2019
From Amiram Eldar, Feb 11 2022: (Start)
Sum_{n>=0} 1/a(n) = (32*log(2) - 11)/10.
Sum_{n>=0} (-1)^n/a(n) = (8*Pi - 56*log(2) + 23)/10. (End)

A344280 Numbers that are both 10-gonal numbers (A001107) and 10-gonal pyramidal numbers (A007585).

Original entry on oeis.org

0, 1, 175, 368050005576

Offset: 1

Seiichi Manyama, May 17 2021

Intersection of A001107 and A007585.

Numberphile, Cannon Ball Numbers.

Cf. A001107, A007585, A027568, A344376, A344377.

for(k=0, 1e4, if(ispolygonal(m=k*(k+1)*(8*k-5)/6, 10), print1(m", ")))

A001107 10-gonal (or decagonal) numbers: a(n) = n(4n-3).

Original entry on oeis.org

0, 1, 10, 27, 52, 85, 126, 175, 232, 297, 370, 451, 540, 637, 742, 855, 976, 1105, 1242, 1387, 1540, 1701, 1870, 2047, 2232, 2425, 2626, 2835, 3052, 3277, 3510, 3751, 4000, 4257, 4522, 4795, 5076, 5365, 5662, 5967, 6280, 6601, 6930, 7267, 7612, 7965, 8326

Offset: 0

N. J. A. Sloane

Write 0, 1, 2, ... in a square spiral, with 0 at the origin and 1 immediately below it; sequence gives numbers on the negative y-axis (see Example section).
Number of divisors of 48^(n-1) for n > 0. - J. Lowell, Aug 30 2008
a(n) is the Wiener index of the graph obtained by connecting two copies of the complete graph K_n by an edge (for n = 3, approximately: |>-<|). The Wiener index of a connected graph is the sum of the distances between all unordered pairs of vertices in the graph. - Emeric Deutsch, Sep 20 2010
This sequence does not contain any squares other than 0 and 1. See A188896. - T. D. Noe, Apr 13 2011
For n > 0: right edge of the triangle A033293. - Reinhard Zumkeller, Jan 18 2012
Sequence found by reading the line from 0, in the direction 0, 10, ... and the parallel line from 1, in the direction 1, 27, ..., in the square spiral whose vertices are the generalized decagonal numbers A074377. - Omar E. Pol, Jul 18 2012
Partial sums give A007585. - Omar E. Pol, Jan 15 2013
This is also a star pentagonal number: a(n) = A000326(n) + 5*A000217(n-1). - Luciano Ancora, Mar 28 2015
Also the number of undirected paths in the n-sunlet graph. - Eric W. Weisstein, Sep 07 2017
After 0, a(n) is the sum of 2*n consecutive integers starting from n-1. - Bruno Berselli, Jan 16 2018
Number of corona of an H0 hexagon with a T(n) triangle. - Craig Knecht, Dec 13 2024

			On a square lattice, place the nonnegative integers at lattice points forming a spiral as follows: place "0" at the origin; then move one step downward (i.e., in the negative y direction) and place "1" at the lattice point reached; then turn 90 degrees in either direction and place a "2" at the next lattice point; then make another 90-degree turn in the same direction and place a "3" at the lattice point; etc. The terms of the sequence will lie along the negative y-axis, as seen in the example below:
  99  64--65--66--67--68--69--70--71--72
   |   |                               |
  98  63  36--37--38--39--40--41--42  73
   |   |   |                       |   |
  97  62  35  16--17--18--19--20  43  74
   |   |   |   |               |   |   |
  96  61  34  15   4---5---6  21  44  75
   |   |   |   |   |       |   |   |   |
  95  60  33  14   3  *0*  7  22  45  76
   |   |   |   |   |   |   |   |   |   |
  94  59  32  13   2--*1*  8  23  46  77
   |   |   |   |           |   |   |   |
  93  58  31  12--11-*10*--9  24  47  78
   |   |   |                   |   |   |
  92  57  30--29--28-*27*-26--25  48  79
   |   |                           |   |
  91  56--55--54--53-*52*-51--50--49  80
   |                                   |
  90--89--88--87--86-*85*-84--83--82--81
[Edited by _Jon E. Schoenfield_, Jan 02 2017]

Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
Bruce C. Berndt, Ramanujan's Notebooks, Part II, Springer; see p. 23.
E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6.
S. M. Ellerstein, The square spiral, J. Recreational Mathematics 29 (#3, 1998) 188; 30 (#4, 1999-2000), 246-250.
R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 2nd ed., 1994, p. 99.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 0..1000
Soren Laing Aletheia-Zomlefer, Lenny Fukshansky, and Stephan Ramon Garcia, The Bateman-Horn Conjecture: Heuristics, History, and Applications, arXiv:1807.08899 [math.NT], 2018-2019. See 6.6.3 p. 33.
Emilio Apricena, A version of the Ulam spiral.
Yin Choi Cheng, Greedy Sidon sets for linear forms, J. Num. Theor. (2024).
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 344.
Craig Knecht, Corona of the H0 hexagon with a T(n) triangle.
Minh Nguyen, 2-adic Valuations of Square Spiral Sequences, Honors Thesis, Univ. of Southern Mississippi (2021).
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Leo Tavares, Illustration: Conjoined Hexagon/Square Pairs
Eric Weisstein's World of Mathematics, Barbell Graph.
Eric Weisstein's World of Mathematics, Decagonal Number.
Eric Weisstein's World of Mathematics, Graph Path.
Eric Weisstein's World of Mathematics, Sunlet Graph.
Index to sequences related to polygonal numbers
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A007585, A028994.
Cf. A093565 ((8, 1) Pascal, column m = 2). Partial sums of A017077.
Sequences from spirals: A001107 (this), A002939, A007742, A033951, A033952, A033953, A033954, A033989, A033990, A033991, A002943, A033996, A033988.
Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.
Cf. A003215.

[4*n^2-3*n : n in [0..50] ]; // Wesley Ivan Hurt, Jun 05 2014

A001107:=-(1+7*z)/(z-1)**3; # Simon Plouffe in his 1992 dissertation

LinearRecurrence[{3, -3, 1}, {0, 1, 10}, 60] (* Harvey P. Dale, May 08 2012 *)
Table[PolygonalNumber[RegularPolygon[10], n], {n, 0, 46}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
Table[4 n^2 - 3 n, {n, 0, 49}] (* Alonso del Arte, Jan 24 2017 *)
PolygonalNumber[10, Range[0, 20]] (* Eric W. Weisstein, Sep 07 2017 *)
LinearRecurrence[{3, -3, 1}, {1, 10, 27}, {0, 20}] (* Eric W. Weisstein, Sep 07 2017 *)

PARI
```
a(n)=4*n^2-3*n
```

a=lambda n: 4*n**2-3*n # Indranil Ghosh, Jan 01 2017
def aList(): # Intended to compute the initial segment of the sequence, not isolated terms.
     x, y = 1, 1
     yield 0
     while True:
         yield x
         x, y = x + y + 8, y + 8
A001107 = aList()
print([next(A001107) for i in range(49)]) # Peter Luschny, Aug 04 2019

a(n) = A033954(-n) = A074377(2*n-1).
a(n) = n + 8*A000217(n-1). - Floor van Lamoen, Oct 14 2005
G.f.: x*(1 + 7*x)/(1 - x)^3.
Partial sums of odd numbers 1 mod 8, i.e., 1, 1 + 9, 1 + 9 + 17, ... . - Jon Perry, Dec 18 2004
1^3 + 3^3*(n-1)/(n+1) + 5^3*((n-1)*(n-2))/((n+1)*(n+2)) + 7^3*((n-1)*(n-2)*(n-3))/((n+1)*(n+2)*(n+3)) + ... = n*(4*n-3) [Ramanujan]. - Neven Juric, Apr 15 2008
Starting (1, 10, 27, 52, ...), this is the binomial transform of [1, 9, 8, 0, 0, 0, ...]. - Gary W. Adamson, Apr 30 2008
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n>2, a(0)=0, a(1)=1, a(2)=10. - Jaume Oliver Lafont, Dec 02 2008
a(n) = 8*n + a(n-1) - 7 for n>0, a(0)=0. - Vincenzo Librandi, Jul 10 2010
a(n) = 8 + 2*a(n-1) - a(n-2). - Ant King, Sep 04 2011
a(n) = A118729(8*n). - Philippe Deléham, Mar 26 2013
a(8*a(n) + 29*n+1) = a(8*a(n) + 29*n) + a(8*n + 1). - Vladimir Shevelev, Jan 24 2014
Sum_{n >= 1} 1/a(n) = Pi/6 + log(2) = 1.216745956158244182494339352... = A244647. - Vaclav Kotesovec, Apr 27 2016
From Ilya Gutkovskiy, Aug 28 2016: (Start)
E.g.f.: x*(1 + 4*x)*exp(x).
Sum_{n >= 1} (-1)^(n+1)/a(n) = (sqrt(2)*Pi - 2*log(2) + 2*sqrt(2)*log(1 + sqrt(2)))/6 = 0.92491492293323294695... (End)
a(n) = A000217(3*n-2) - A000217(n-2). In general, if P(k,n) be the n-th k-gonal number and T(n) be the n-th triangular number, A000217(n), then P(T(k),n) = T((k-1)*n - (k-2)) - T(k-3)*T(n-2). - Charlie Marion, Sep 01 2020
Product_{n>=2} (1 - 1/a(n)) = 4/5. - Amiram Eldar, Jan 21 2021
a(n) = A003215(n-1) + A000290(n) - 1. - Leo Tavares, Jul 23 2022

A213500 Rectangular array T(n,k): (row n) = bc, where b(h) = h, c(h) = h + n - 1, n >= 1, h >= 1, and = convolution.

Original entry on oeis.org

1, 4, 2, 10, 7, 3, 20, 16, 10, 4, 35, 30, 22, 13, 5, 56, 50, 40, 28, 16, 6, 84, 77, 65, 50, 34, 19, 7, 120, 112, 98, 80, 60, 40, 22, 8, 165, 156, 140, 119, 95, 70, 46, 25, 9, 220, 210, 192, 168, 140, 110, 80, 52, 28, 10, 286, 275, 255, 228, 196, 161, 125, 90

Offset: 1

Clark Kimberling, Jun 14 2012

Principal diagonal: A002412.
Antidiagonal sums: A002415.
Row 1:  (1,2,3,...)**(1,2,3,...) = A000292.
Row 2:  (1,2,3,...)**(2,3,4,...) = A005581.
Row 3:  (1,2,3,...)**(3,4,5,...) = A006503.
Row 4:  (1,2,3,...)**(4,5,6,...) = A060488.
Row 5:  (1,2,3,...)**(5,6,7,...) = A096941.
Row 6:  (1,2,3,...)**(6,7,8,...) = A096957.
...
In general, the convolution of two infinite sequences is defined from the convolution of two n-tuples:  let X(n) = (x(1),...,x(n)) and Y(n)=(y(1),...,y(n)); then X(n)**Y(n) = x(1)*y(n)+x(2)*y(n-1)+...+x(n)*y(1); this sum is the n-th term in the convolution of infinite sequences:(x(1),...,x(n),...)**(y(1),...,y(n),...), for all n>=1.
...
In the following guide to related arrays and sequences, row n of each array T(n,k) is the convolution b**c of the sequences b(h) and c(h+n-1). The principal diagonal is given by T(n,n) and the n-th antidiagonal sum by S(n). In some cases, T(n,n) or S(n) differs in offset from the listed sequence.
b(h)........ c(h)........ T(n,k) .. T(n,n) .. S(n)
h .......... h .......... A213500 . A002412 . A002415
h .......... h^2 ........ A212891 . A213436 . A024166
h^2 ........ h .......... A213503 . A117066 . A033455
h^2 ........ h^2 ........ A213505 . A213546 . A213547
h .......... h*(h+1)/2 .. A213548 . A213549 . A051836
h*(h+1)/2 .. h .......... A213550 . A002418 . A005585
h*(h+1)/2 .. h*(h+1)/2 .. A213551 . A213552 . A051923
h .......... h^3 ........ A213553 . A213554 . A101089
h^3 ........ h .......... A213555 . A213556 . A213547
h^3 ........ h^3 ........ A213558 . A213559 . A213560
h^2 ........ h*(h+1)/2 .. A213561 . A213562 . A213563
h*(h+1)/2 .. h^2 ........ A213564 . A213565 . A101094
2^(h-1) .... h .......... A213568 . A213569 . A047520
2^(h-1) .... h^2 ........ A213573 . A213574 . A213575
h .......... Fibo(h) .... A213576 . A213577 . A213578
Fibo(h) .... h .......... A213579 . A213580 . A053808
Fibo(h) .... Fibo(h) .... A067418 . A027991 . A067988
Fibo(h+1) .. h .......... A213584 . A213585 . A213586
Fibo(n+1) .. Fibo(h+1) .. A213587 . A213588 . A213589
h^2 ........ Fibo(h) .... A213590 . A213504 . A213557
Fibo(h) .... h^2 ........ A213566 . A213567 . A213570
h .......... -1+2^h ..... A213571 . A213572 . A213581
-1+2^h ..... h .......... A213582 . A213583 . A156928
-1+2^h ..... -1+2^h ..... A213747 . A213748 . A213749
h .......... 2*h-1 ...... A213750 . A007585 . A002417
2*h-1 ...... h .......... A213751 . A051662 . A006325
2*h-1 ...... 2*h-1 ...... A213752 . A100157 . A071238
2*h-1 ...... -1+2^h ..... A213753 . A213754 . A213755
-1+2^h ..... 2*h-1 ...... A213756 . A213757 . A213758
2^(n-1) .... 2*h-1 ...... A213762 . A213763 . A213764
2*h-1 ...... Fibo(h) .... A213765 . A213766 . A213767
Fibo(h) .... 2*h-1 ...... A213768 . A213769 . A213770
Fibo(h+1) .. 2*h-1 ...... A213774 . A213775 . A213776
Fibo(h) .... Fibo(h+1) .. A213777 . A001870 . A152881
h .......... 1+[h/2] .... A213778 . A213779 . A213780
1+[h/2] .... h .......... A213781 . A213782 . A005712
1+[h/2] .... [(h+1)/2] .. A213783 . A213759 . A213760
h .......... 3*h-2 ...... A213761 . A172073 . A002419
3*h-2 ...... h .......... A213771 . A213772 . A132117
3*h-2 ...... 3*h-2 ...... A213773 . A214092 . A213818
h .......... 3*h-1 ...... A213819 . A213820 . A153978
3*h-1 ...... h .......... A213821 . A033431 . A176060
3*h-1 ...... 3*h-1 ...... A213822 . A213823 . A213824
3*h-1 ...... 3*h-2 ...... A213825 . A213826 . A213827
3*h-2 ...... 3*h-1 ...... A213828 . A213829 . A213830
2*h-1 ...... 3*h-2 ...... A213831 . A213832 . A212560
3*h-2 ...... 2*h-1 ...... A213833 . A130748 . A213834
h .......... 4*h-3 ...... A213835 . A172078 . A051797
4*h-3 ...... h .......... A213836 . A213837 . A071238
4*h-3 ...... 2*h-1 ...... A213838 . A213839 . A213840
2*h-1 ...... 4*h-3 ...... A213841 . A213842 . A213843
2*h-1 ...... 4*h-1 ...... A213844 . A213845 . A213846
4*h-1 ...... 2*h-1 ...... A213847 . A213848 . A180324
[(h+1)/2] .. [(h+1)/2] .. A213849 . A049778 . A213850
h .......... C(2*h-2,h-1) A213853
...
Suppose that u = (u(n)) and v = (v(n)) are sequences having generating functions U(x) and V(x), respectively. Then the convolution u**v has generating function U(x)*V(x). Accordingly, if u and v are homogeneous linear recurrence sequences, then every row of the convolution array T satisfies the same homogeneous linear recurrence equation, which can be easily obtained from the denominator of U(x)*V(x). Also, every column of T has the same homogeneous linear recurrence as v.

			Northwest corner (the array is read by southwest falling antidiagonals):
  1,  4, 10, 20,  35,  56,  84, ...
  2,  7, 16, 30,  50,  77, 112, ...
  3, 10, 22, 40,  65,  98, 140, ...
  4, 13, 28, 50,  80, 119, 168, ...
  5, 16, 34, 60,  95, 140, 196, ...
  6, 19, 40, 70, 110, 161, 224, ...
T(6,1) = (1)**(6) = 6;
T(6,2) = (1,2)**(6,7) = 1*7+2*6 = 19;
T(6,3) = (1,2,3)**(6,7,8) = 1*8+2*7+3*6 = 40.

Clark Kimberling, Antidiagonals n = 1..60, flattened

Cf. A000027.

b[n_] := n; c[n_] := n
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}]  (* A213500 *)

t(n,k) = sum(i=0, k - 1, (k - i) * (n + i));
tabl(nn) = {for(n=1, nn, for(k=1, n, print1(t(k,n - k + 1),", ");); print(););};
tabl(12) \\ Indranil Ghosh, Mar 26 2017

def t(n, k): return sum((k - i) * (n + i) for i in range(k))
for n in range(1, 13):
    print([t(k, n - k + 1) for k in range(1, n + 1)]) # Indranil Ghosh, Mar 26 2017

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

A080851 Square array of pyramidal numbers, read by antidiagonals.

Original entry on oeis.org

1, 1, 3, 1, 4, 6, 1, 5, 10, 10, 1, 6, 14, 20, 15, 1, 7, 18, 30, 35, 21, 1, 8, 22, 40, 55, 56, 28, 1, 9, 26, 50, 75, 91, 84, 36, 1, 10, 30, 60, 95, 126, 140, 120, 45, 1, 11, 34, 70, 115, 161, 196, 204, 165, 55, 1, 12, 38, 80, 135, 196, 252, 288, 285, 220, 66, 1, 13, 42, 90, 155, 231, 308, 372, 405, 385, 286, 78

Offset: 0

Paul Barry, Feb 21 2003

The first row contains the triangular numbers, which are really two-dimensional, but can be regarded as degenerate pyramidal numbers. - N. J. A. Sloane, Aug 28 2015

			Array begins (n>=0, k>=0):
1,  3,  6, 10,  15,  21,  28,  36,  45,   55, ... A000217
1,  4, 10, 20,  35,  56,  84, 120, 165,  220, ... A000292
1,  5, 14, 30,  55,  91, 140, 204, 285,  385, ... A000330
1,  6, 18, 40,  75, 126, 196, 288, 405,  550, ... A002411
1,  7, 22, 50,  95, 161, 252, 372, 525,  715, ... A002412
1,  8, 26, 60, 115, 196, 308, 456, 645,  880, ... A002413
1,  9, 30, 70, 135, 231, 364, 540, 765, 1045, ... A002414
1, 10, 34, 80, 155, 266, 420, 624, 885, 1210, ... A007584

Alois P. Heinz, Antidiagonals n = 0..140, flattened

Numerous sequences in the database are to be found in the array. Rows include the pyramidal numbers A000217, A000292, A000330, A002411, A002412, A002413, A002414, A007584, A007585, A007586.
Columns include or are closely related to A017029, A017113, A017017, A017101, A016777, A017305. Diagonals include A006325, A006484, A002417.
Cf. A057145, A027660 (antidiagonal sums).
See A257199 for another version of this array.

vector(vector(poly_coeff(Taylor((1+kx)/(1-x)^4,x,11),x,n),n,0,11),k,-1,10) VECTOR(VECTOR(comb(k+2,2)+comb(k+2,3)n, k, 0, 11), n, 0, 11)

A080851 := proc(n,k)
    binomial(k+3,3)+(n-1)*binomial(k+2,3) ;
end proc:
seq( seq(A080851(d-k,k),k=0..d),d=0..12) ; # R. J. Mathar, Oct 01 2021

pyramidalFigurative[ ngon_, rank_] := (3 rank^2 + rank^3 (ngon - 2) - rank (ngon - 5))/6; Table[ pyramidalFigurative[n-k-1, k], {n, 4, 15}, {k, n-3}] // Flatten (* Robert G. Wilson v, Sep 15 2015 *)

T(n, k) = binomial(k+3, 3) + (n-1)*binomial(k+2, 3), corrected Oct 01 2021.
T(n, k) = T(n-1, k) + C(k+2, 3) = T(n-1, k) + k*(k+1)*(k+2)/6.
G.f. for rows: (1 + n*x)/(1-x)^4, n>=-1.
T(n,k) = sum_{j=1..k+1} A057145(n+2,j). - R. J. Mathar, Jul 28 2016

A237616 a(n) = n(n + 1)(5*n - 4)/2.

Original entry on oeis.org

0, 1, 18, 66, 160, 315, 546, 868, 1296, 1845, 2530, 3366, 4368, 5551, 6930, 8520, 10336, 12393, 14706, 17290, 20160, 23331, 26818, 30636, 34800, 39325, 44226, 49518, 55216, 61335, 67890, 74896, 82368, 90321, 98770, 107730, 117216, 127243, 137826, 148980, 160720

Offset: 0

Bruno Berselli, Feb 10 2014

Also 17-gonal (or heptadecagonal) pyramidal numbers.

This sequence is related to A226489 by 2*a(n) = n*A226489(n) - Sum_{i=0..n-1} A226489(i).

			After 0, the sequence is provided by the row sums of the triangle:
   1;
   2,  16;
   3,  32,  31;
   4,  48,  62,  46;
   5,  64,  93,  92,  61;
   6,  80, 124, 138, 122,  76;
   7,  96, 155, 184, 183, 152,  91;
   8, 112, 186, 230, 244, 228, 182, 106;
   9, 128, 217, 276, 305, 304, 273, 212, 121;
  10, 144, 248, 322, 366, 380, 364, 318, 242, 136; etc.,
where (r = row index, c = column index):
T(r,r) = T(c,c) = 15*r-14 and T(r,c) = T(r-1,c)+T(r,r) = (r-c+1)*T(r,r), with r>=c>0.

E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93 (fifteenth row of the table).

Bruno Berselli, Table of n, a(n) for n = 0..1000
Eric Weisstein's World of Mathematics, Pyramidal Number.
Index to sequences related to pyramidal numbers.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Cf. A051869, A104728.

Cf. sequences with formula n*(n+1)*(k*n-k+3)/6: A000217 (k=0), A000292 (k=1), A000330 (k=2), A002411 (k=3), A002412 (k=4), A002413 (k=5), A002414 (k=6), A007584 (k=7), A007585 (k=8), A007586 (k=9), A007587 (k=10), A050441 (k=11), A172073 (k=12), A177890 (k=13), A172076 (k=14), this sequence (k=15), A172078(k=16), A237617 (k=17), A172082 (k=18), A237618 (k=19), A172117(k=20), A256718 (k=21), A256716 (k=22), A256645 (k=23), A256646(k=24), A256647 (k=25), A256648 (k=26), A256649 (k=27), A256650(k=28).

List([0..40], n-> n*(n+1)*(5*n-4)/2); # G. C. Greubel, Aug 30 2019

Magma
```
[n*(n+1)*(5*n-4)/2: n in [0..40]];
```

I:=[0,1,18,66]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..50]]; // Vincenzo Librandi, Feb 12 2014

seq(n*(n+1)*(5*n-4)/2, n=0..40); # G. C. Greubel, Aug 30 2019

Table[n(n+1)(5n-4)/2, {n, 0, 40}]
CoefficientList[Series[x (1+14x)/(1-x)^4, {x, 0, 40}], x] (* Vincenzo Librandi, Feb 12 2014 *)
LinearRecurrence[{4,-6,4,-1},{0,1,18,66},50] (* Harvey P. Dale, Jan 11 2015 *)

a(n)=n*(n+1)*(5*n-4)/2 \\ Charles R Greathouse IV, Sep 24 2015

[n*(n+1)*(5*n-4)/2 for n in (0..40)] # G. C. Greubel, Aug 30 2019

G.f.: x*(1 + 14*x)/(1 - x)^4.
For n>0, a(n) = Sum_{i=0..n-1} (n-i)*(15*i+1). More generally, the sequence with the closed form n*(n+1)*(k*n-k+3)/6 is also given by Sum_{i=0..n-1} (n-i)*(k*i+1) for n>0.
a(n) = A104728(A001844(n-1)) for n>0.
Sum_{n>=1} 1/a(n) = (2*sqrt(5*(5 + 2*sqrt(5)))*Pi + 10*sqrt(5)*arccoth(sqrt(5)) + 25*log(5) - 16)/72 = 1.086617842136293176... . - Vaclav Kotesovec, Dec 07 2016
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) for n >= 4. - Wesley Ivan Hurt, Dec 18 2020
E.g.f.: exp(x)*x*(2 + 16*x + 5*x^2)/2. - Elmo R. Oliveira, Aug 04 2025

A093565 (8,1) Pascal triangle.

Original entry on oeis.org

1, 8, 1, 8, 9, 1, 8, 17, 10, 1, 8, 25, 27, 11, 1, 8, 33, 52, 38, 12, 1, 8, 41, 85, 90, 50, 13, 1, 8, 49, 126, 175, 140, 63, 14, 1, 8, 57, 175, 301, 315, 203, 77, 15, 1, 8, 65, 232, 476, 616, 518, 280, 92, 16, 1, 8, 73, 297, 708, 1092, 1134, 798, 372, 108, 17, 1, 8, 81, 370, 1005

Offset: 0

Wolfdieter Lang, Apr 22 2004

The array F(8;n,m) gives in the columns m>=1 the figurate numbers based on A017077, including the decagonal numbers A001107,(see the W. Lang link).
This is the eighth member, d=8, in the family of triangles of figurate numbers, called (d,1) Pascal triangles: A007318 (Pascal), A029653, A093560-4, for d=1..7.
This is an example of a Riordan triangle (see A093560 for a comment and A053121 for a comment and the 1991 Shapiro et al. reference on the Riordan group). Therefore the o.g.f. for the row polynomials p(n,x):=Sum_{m=0..n} a(n,m)*x^m is G(z,x)=(1+7*z)/(1-(1+x)*z).
The SW-NE diagonals give A022098(n-1) = Sum_{k=0..ceiling((n-1)/2)} a(n-1-k,k), n >= 1, with n=0 value 7. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs.

			Triangle begins
  [1];
  [8,  1];
  [8,  9,  1];
  [8, 17, 10,  1];
  ...

Kurt Hawlitschek, Johann Faulhaber 1580-1635, Veroeffentlichung der Stadtbibliothek Ulm, Band 18, Ulm, Germany, 1995, Ch. 2.1.4. Figurierte Zahlen.
Ivo Schneider: Johannes Faulhaber 1580-1635, Birkhäuser, Basel, Boston, Berlin, 1993, ch.5, pp. 109-122.

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
W. Lang, First 10 rows and array of figurate numbers .

Row sums: A005010(n-1), n>=1, 1 for n=0, alternating row sums are 1 for n=0, 7 for n=2 and 0 else.
The column sequences give for m=1..9: A017077, A001107 (decagonal), A007585, A051797, A051878, A050404, A052226, A056001, A056122.
Cf. A093644 (d=9).

a093565 n k = a093565_tabl !! n !! k
a093565_row n = a093565_tabl !! n
a093565_tabl = [1] : iterate
               (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [8, 1]
-- Reinhard Zumkeller, Aug 31 2014

a(n, m)=F(8;n-m, m) for 0<= m <= n, otherwise 0, with F(8;0, 0)=1, F(8;n, 0)=8 if n>=1 and F(8;n, m):=(8*n+m)*binomial(n+m-1, m-1)/m if m>=1.
Recursion: a(n, m)=0 if m>n, a(0, 0)= 1; a(n, 0)=8 if n>=1; a(n, m)= a(n-1, m) + a(n-1, m-1).
G.f. column m (without leading zeros): (1+7*x)/(1-x)^(m+1), m>=0.
T(n, k) = C(n, k) + 7*C(n-1, k). - Philippe Deléham, Aug 28 2005
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(8 + 17*x + 10*x^2/2! + x^3/3!) = 8 + 25*x + 52*x^2/2! + 90*x^3/3! + 140*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ). - Peter Bala, Dec 22 2014

A076454 Sum of numbers that can be written as tn + u(n+1) for nonnegative integers t,u in exactly one way.

Original entry on oeis.org

1, 21, 102, 310, 735, 1491, 2716, 4572, 7245, 10945, 15906, 22386, 30667, 41055, 53880, 69496, 88281, 110637, 136990, 167790, 203511, 244651, 291732, 345300, 405925, 474201, 550746, 636202, 731235, 836535, 952816, 1080816, 1221297, 1375045, 1542870, 1725606, 1924111

Offset: 1

Floor van Lamoen, Oct 13 2002

This sequence is related to A007585 by a(n) = n*A007585(n) - Sum_{i=0..n-1} A007585(i). - Vincenzo Librandi, Aug 08 2010
In fact, this is the case d=4 in the identity n*(n*(n+1)*(2*d*n-2*d+3)/6) - Sum_{k=0..n-1} k*(k+1)*(2*d*k-2*d+3)/6 = n*(n+1)*(3*d*n^2-d*n+4*n-2*d+2)/12. - Bruno Berselli, Mar 01 2012
Bisection of A233329 (up to an offset). - L. Edson Jeffery, Jan 23 2014

Fred. Schuh, Vragen betreffende een onbepaalde vergelijking, Nieuw Tijdschrift voor Wiskunde, 52 (1964-1965) 193-198.

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
B. Berselli, A description of the recursive method in Comments lines: website Matem@ticamente (in Italian).
Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5,1).

Cf. A002417, A007585, A076455-A076459, A233329.

[n*(n+1)*(2*n^2-1)/2: n in [1..50]]; // Vincenzo Librandi, Dec 30 2013

Maple
```
seq(1/2*n*(n+1)*(2*n^2-1),n=1..40);
```

CoefficientList[Series[(1 + 16 x + 7 x^2)/(1 - x)^5, {x, 0, 50}], x] (* Vincenzo Librandi, Dec 30 2013 *)
LinearRecurrence[{5,-10,10,-5,1},{1,21,102,310,735},40] (* Harvey P. Dale, Jun 30 2023 *)

a(n) = n*(n+1)*(2*n^2-1)/2.
G.f.: x*(1+16*x+7*x^2)/(1-x)^5.
a(n) = 5*a(n-1)-10*a(n-2)+10*a(n-3)-5*a(n-4)+a(n-5), n>=6, with a(1)=1, a(2)=21, a(3)=102, a(4)=310, a(5)=735. - L. Edson Jeffery, Dec 30 2013

Comments rewritten from Bruno Berselli, Mar 01 2012

More terms from Vincenzo Librandi, Dec 30 2013

A132124 a(n) = n(n+1)(8*n + 1)/6.

Original entry on oeis.org

0, 3, 17, 50, 110, 205, 343, 532, 780, 1095, 1485, 1958, 2522, 3185, 3955, 4840, 5848, 6987, 8265, 9690, 11270, 13013, 14927, 17020, 19300, 21775, 24453, 27342, 30450, 33785, 37355, 41168, 45232, 49555, 54145, 59010, 64158, 69597, 75335, 81380, 87740, 94423

Offset: 0

Reinhard Zumkeller, Aug 12 2007

Convolution of the sequences (0,3,5,0,0,0,...) and (binomial(n+3, 3)), n >= 0. - Emeric Deutsch, Aug 30 2007

Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Cf. A000330, A033994, A132112, A050409.
Cf. A014105, A016061; A002412, A006331; A007585, A002378.
Row sums of A069011.

seq((1/6)*n*(n+1)*(8*n+1),n=0..40); # Emeric Deutsch, Aug 30 2007

a[n_] := n*(n + 1)*(8*n + 1)/6; Array[a, 42, 0] (* Amiram Eldar, May 20 2023 *)

a(n) = A132121(n,2) for n > 1.
G.f.: x*(3+5*x)/(1-x)^4. - Emeric Deutsch, Aug 30 2007
From Bruno Berselli, Nov 25 2010: (Start)
a(n) = n*A014105(n) - A016061(n-1), since A016061(n-1) = Sum_{k=0..n-1} A014105(k) (n > 0).
Also a(n) = A002412(n) + A006331(n) = A007585(n) + A002378(n). (End)
Sum_{n>=1} 1/a(n) = 54 - 24*(sqrt(2)+1)*Pi/7 - 24*(sqrt(2)+8)*log(2)/7 + 48*sqrt(2)*log(2-sqrt(2))/7. - Amiram Eldar, May 20 2023
E.g.f.: exp(x)*x*(18 + 33*x + 8*x^2)/6. - Stefano Spezia, Feb 21 2024

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

Original entry on oeis.org

1, 5, 3, 14, 11, 5, 30, 26, 17, 7, 55, 50, 38, 23, 9, 91, 85, 70, 50, 29, 11, 140, 133, 115, 90, 62, 35, 13, 204, 196, 175, 145, 110, 74, 41, 15, 285, 276, 252, 217, 175, 130, 86, 47, 17, 385, 375, 348, 308, 259, 205, 150, 98, 53, 19, 506, 495, 465, 420

Offset: 1

Clark Kimberling, Jun 20 2012

Principal diagonal:  A007585
Antidiagonal sums:  A002417
row 1,  (1,2,3,4,5,...)**(1,3,5,7,9,...): A000330
row 2,  (1,2,3,4,5,...)**(3,5,7,9,...): A051925
row 3,  (1,2,3,4,5,...)**(5,7,9,11,...): (2*k^3 + 15*k^2 + 13*k)/6
row 4,  (1,2,3,4,5,...)**(7,9,11,13,...): (2*k^3 + 21*k^2 + 19*k)/6
For a guide to related arrays, see A213500.

			Northwest corner (the array is read by falling antidiagonals):
1....5....14...30....55....91
3....11...26...50....85....133
5....17...38...70....115...175
7....23...50...90....145...217
9....29...62...110...175...259
11...35...74...130...205...301

Cf. A213500.

b[n_] := n; c[n_] := 2 n - 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}]  (* A213750 *)
d = Table[t[n, n], {n, 1, 40}] (* A007585 *)
s1 = Table[s[n], {n, 1, 50}] (* A002417 *)
FindLinearRecurrence[s1]
FindGeneratingFunction[s1, x]

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) = (2*n - 1) - (2*n - 3)*x and g(x) = (1 - x )^4.

cp's OEIS Frontend

A051797 Partial sums of A007585.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Magma

Magma

Maple

Mathematica

PARI

Sage

Formula

A344280 Numbers that are both 10-gonal numbers (A001107) and 10-gonal pyramidal numbers (A007585).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A001107 10-gonal (or decagonal) numbers: a(n) = n*(4*n-3).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Python

Formula

A213500 Rectangular array T(n,k): (row n) = b**c, where b(h) = h, c(h) = h + n - 1, n >= 1, h >= 1, and ** = convolution.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A080851 Square array of pyramidal numbers, read by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Derive

Maple

Mathematica

Formula

A237616 a(n) = n*(n + 1)*(5*n - 4)/2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

A001107 10-gonal (or decagonal) numbers: a(n) = n(4n-3).

A213500 Rectangular array T(n,k): (row n) = bc, where b(h) = h, c(h) = h + n - 1, n >= 1, h >= 1, and = convolution.

A237616 a(n) = n(n + 1)(5*n - 4)/2.

A076454 Sum of numbers that can be written as tn + u(n+1) for nonnegative integers t,u in exactly one way.

A132124 a(n) = n(n+1)(8*n + 1)/6.