A051740 -id:A051740 - cp's OEIS frontend

A051877 Partial sums of A051740.

1, 12, 57, 182, 462, 1008, 1974, 3564, 6039, 9724, 15015, 22386, 32396, 45696, 63036, 85272, 113373, 148428, 191653, 244398, 308154, 384560, 475410, 582660, 708435, 855036, 1024947, 1220842, 1445592, 1702272, 1994168, 2324784, 2697849, 3117324, 3587409

Offset: 0

Barry E. Williams, Dec 14 1999

Convolution of triangular numbers (A000217) and enneagonal numbers (A001106). - Bruno Berselli, Jul 21 2015

A. H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, pp. 194-196.
Herbert John Ryser, Combinatorial Mathematics, "The Carus Mathematical Monographs", No. 14, John Wiley and Sons, 1963, pp. 1-16.

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

Cf. A051740; A000217, A001106.

Cf. A093564 ((7, 1) Pascal, column m=5).

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

[(n+1)*(n+2)*(n+3)*(n+4)*(7*n+5)/120 : n in [0..40]]; // Wesley Ivan Hurt, May 02 2015

A051877:=n->binomial(n+4,4)*(7*n+5)/5: seq(A051877(n), n=0..40); # Wesley Ivan Hurt, May 02 2015

Table[(n+1)(n+2)(n+3)(n+4)(7n+5)/120, {n, 0, 40}] (* Vincenzo Librandi, May 03 2015 *)
LinearRecurrence[{6,-15,20,-15,6,-1},{1,12,57,182,462,1008},40] (* Harvey P. Dale, May 05 2022 *)

vector(40, n, (7*n-2)*binomial(n+3,4)/5) \\ G. C. Greubel, Aug 29 2019

[(7*n+5)*binomial(n+4,4)/5 for n in (0..40)] # G. C. Greubel, Aug 29 2019

a(n) = C(n+4, 4)*(7*n+5)/5.
G.f.: (1+6*x)/(1-x)^6.
From Wesley Ivan Hurt, May 02 2015: (Start)
a(n) = 6*a(n-1) -15*a(n-2) +20*a(n-3) -15*a(n-4) +6*a(n-5) -a(n-6).
a(n) = (n+1)*(n+2)*(n+3)*(n+4)*(7*n+5)/120. (End)
E.g.f.: (5! +1320*x +2040*x^2 +920*x^3 +145*x^4 +7*x^5)*exp(x)/5!

A080852 Square array of 4D pyramidal numbers, read by antidiagonals.

Original entry on oeis.org

1, 1, 4, 1, 5, 10, 1, 6, 15, 20, 1, 7, 20, 35, 35, 1, 8, 25, 50, 70, 56, 1, 9, 30, 65, 105, 126, 84, 1, 10, 35, 80, 140, 196, 210, 120, 1, 11, 40, 95, 175, 266, 336, 330, 165, 1, 12, 45, 110, 210, 336, 462, 540, 495, 220, 1, 13, 50, 125, 245, 406, 588, 750, 825, 715, 286

Offset: 0

Paul Barry, Feb 21 2003

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

			Array, n >= 0, k >= 0, begins
1 4 10 20  35  56 ...
1 5 15 35  70 126 ...
1 6 20 50 105 196 ...
1 7 25 65 140 266 ...
1 8 30 80 175 336 ...

Index to sequences related to pyramidal numbers

Rows include A000292, A000332, A002415, A001296, A002418, A002419, A051740, A051797.
Cf. A057145, A080851, A180266, A055796 (antidiagonal sums).
See A257200 for another version of the array.

vector(vector(poly_coeff(Taylor((1+kx)/(1-x)^5,x,11),x,n),n,0,11),k,-1,10)

VECTOR(VECTOR(comb(k+3,3)+comb(k+3,4)n, k, 0, 11), n, 0, 11)

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

T[n_, k_] := Binomial[k+3, 3] + Binomial[k+3, 4]n;
Table[T[n-k, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Jean-François Alcover, May 05 2023 *)

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

A093564 (7,1) Pascal triangle.

Original entry on oeis.org

1, 7, 1, 7, 8, 1, 7, 15, 9, 1, 7, 22, 24, 10, 1, 7, 29, 46, 34, 11, 1, 7, 36, 75, 80, 45, 12, 1, 7, 43, 111, 155, 125, 57, 13, 1, 7, 50, 154, 266, 280, 182, 70, 14, 1, 7, 57, 204, 420, 546, 462, 252, 84, 15, 1, 7, 64, 261, 624, 966, 1008, 714, 336, 99, 16, 1, 7, 71, 325, 885

Offset: 0

Wolfdieter Lang, Apr 22 2004

The array F(7;n,m) gives in the columns m>=1 the figurate numbers based on A016993, including the 9-gonal numbers A001106, (see the W. Lang link).
This is the seventh member, d=7, in the family of triangles of figurate numbers, called (d,1) Pascal triangles: A007318 (Pascal), A029653, A093560-3, for d=1..6.
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+6*z)/(1-(1+x)*z).
The SW-NE diagonals give A022097(n-1) = Sum_{k=0..ceiling((n-1)/2)} a(n-1-k,k), n >= 1, with n=0 value 6. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs.

			Triangle begins
  [1];
  [7,  1];
  [7,  8,  1];
  [7, 15,  9,  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: A000079(n+2), n>=1, 1 for n=0, alternating row sums are 1 for n=0, 6 for n=2 and 0 otherwise.
The column sequences give for m=1..9: A016993, A001106 (9-gonal), A007584, A051740, A051877, A050403, A027818, A034266, A055994.
Cf. A093565 (d=8).

a093564 n k = a093564_tabl !! n !! k
a093564_row n = a093564_tabl !! n
a093564_tabl = [1] : iterate
               (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [7, 1]
-- Reinhard Zumkeller, Sep 01 2014

N:= 20: # to get the first N rows
T:=Matrix(N,N):
T[1,1]:= 1:
for m from 2 to N do
T[m,1]:= 7:
T[m,2..m]:= T[m-1,1..m-1] + T[m-1,2..m];
od:
for m from 1 to N do
convert(T[m,1..m],list)
od; # Robert Israel, Dec 28 2014

a(n, m)=F(7;n-m, m) for 0<= m <= n, otherwise 0, with F(7;0, 0)=1, F(7;n, 0)=7 if n>=1 and F(7;n, m):=(7*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)=7 if n>=1; a(n, m)= a(n-1, m) + a(n-1, m-1).
G.f. column m (without leading zeros): (1+6*x)/(1-x)^(m+1), m>=0.
T(n, k) = C(n, k) + 6*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)*(7 + 15*x + 9*x^2/2! + x^3/3!) = 7 + 22*x + 46*x^2/2! + 80*x^3/3! + 125*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

A220212 Convolution of natural numbers (A000027) with tetradecagonal numbers (A051866).

Original entry on oeis.org

0, 1, 16, 70, 200, 455, 896, 1596, 2640, 4125, 6160, 8866, 12376, 16835, 22400, 29240, 37536, 47481, 59280, 73150, 89320, 108031, 129536, 154100, 182000, 213525, 248976, 288666, 332920, 382075, 436480, 496496, 562496, 634865, 714000, 800310, 894216, 996151

Offset: 0

Bruno Berselli, Dec 08 2012

Partial sums of A172073.

Apart from 0, all terms are in A135021: a(n) = A135021(A034856(n+1)) with n>0.

Bruno Berselli, 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. A135021, A172073.
Cf. convolution of the natural numbers (A000027) with the k-gonal numbers (* means "except 0"):
k= 2 (A000027 ): A000292;
k= 3 (A000217 ): A000332 (after the third term);
k= 4 (A000290 ): A002415 (after the first term);
k= 5 (A000326 ): A001296;
k= 6 (A000384*): A002417;
k= 7 (A000566 ): A002418;
k= 8 (A000567*): A002419;
k= 9 (A001106*): A051740;
k=10 (A001107*): A051797;
k=11 (A051682*): A051798;
k=12 (A051624*): A051799;
k=13 (A051865*): A055268.
Cf. similar sequences with formula n*(n+1)*(n+2)*(k*n-k+2)/12 listed in A264850.

A051866:=func; [&+[(n-k+1)*A051866(k): k in [0..n]]: n in [0..37]];

I:=[0,1,16,70,200]; [n le 5 select I[n] else 5*Self(n-1)-10*Self(n-2)+10*Self(n-3)-5*Self(n-4)+Self(n-5): n in [1..50]]; // Vincenzo Librandi, Aug 18 2013

A051866[k_] := k (6 k - 5); Table[Sum[(n - k + 1) A051866[k], {k, 0, n}], {n, 0, 37}]
CoefficientList[Series[x (1 + 11 x) / (1 - x)^5, {x, 0, 40}], x] (* Vincenzo Librandi, Aug 18 2013 *)

G.f.: x*(1+11*x)/(1-x)^5.
a(n) = n*(n+1)*(n+2)*(3*n-2)/6.
From Amiram Eldar, Feb 15 2022: (Start)
Sum_{n>=1} 1/a(n) = 3*(3*sqrt(3)*Pi + 27*log(3) - 17)/80.
Sum_{n>=1} (-1)^(n+1)/a(n) = 3*(6*sqrt(3)*Pi - 64*log(2) + 37)/80. (End)

A241765 a(n) = n(n + 1)(n + 2)(3n + 17)/24.

Original entry on oeis.org

0, 5, 23, 65, 145, 280, 490, 798, 1230, 1815, 2585, 3575, 4823, 6370, 8260, 10540, 13260, 16473, 20235, 24605, 29645, 35420, 41998, 49450, 57850, 67275, 77805, 89523, 102515, 116870, 132680, 150040, 169048, 189805, 212415, 236985, 263625, 292448

Offset: 0

Bruno Berselli, Apr 28 2014

Equivalently, Sum_{i=0..n} (i+4)*A000217(i).
Sequences of the type Sum_{i=0..n} (i+k)*A000217(i):
k = 0,  A001296: 0,  1,  7, 25,  65, 140, 266, 462, ...
k = 1,  A000914: 0,  2, 11, 35,  85, 175, 322, 546, ...
k = 2,  A050534: 0,  3, 15, 45, 105, 210, 378, 630, ... (deleting two 0)
k = 3,  A215862: 0,  4, 19, 55, 125, 245, 434, 714, ...
k = 4,     a(n): 0,  5, 23, 65, 145, 280, 490, 798, ...
k = 5,  A239568: 0,  6, 27, 75, 165, 315, 546, 882, ...
Antidiagonal sums (without 0) give A034263: 1, 9, 39, 119, 294, ...
Diagonal: 1, 11, 45, 125, 280, 546, ... is A051740.
Also: k = -1 gives A050534 deleting a 0; k = -2 gives 0 followed by A059302.
After 0, partial sums of A212343 and third column of A118788.
This sequence is even related to A005286 by a(n) = n*A005286(n) - Sum_{i=0..n-1} A005286(i).

			a(7) = 4*0 + 5*1 + 6*3 + 7*6 + 8*10 + 9*15 + 10*21 + 11*28 = 798.

Bruno Berselli, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5,1).

Cf. A000217, A005286, A118788, A212343, A227342.

Cf. similar sequences A000914, A001296, A050534, A059302, A215862, A239568 (see table in Comments lines).

/* By first comment: */ k:=4; A000217:=func; [&+[(i+k)*A000217(i): i in [0..n]]: n in [0..40]];

Maple

A241765:=n->n*(n + 1)*(n + 2)*(3*n + 17)/24; seq(A241765(n), n=0..40); # Wesley Ivan Hurt, May 09 2014

Mathematica

Table[n (n + 1) (n + 2) (3 n + 17)/24, {n, 0, 40}] (* or *) LinearRecurrence[{5, -10, 10, -5, 1}, {0, 5, 23, 65, 145}, 40] CoefficientList[Series[x (5 - 2 x)/(1 - x)^5, {x, 0, 40}], x] (* Vincenzo Librandi, May 09 2014 *)

Maxima

makelist(coeff(taylor(x*(5-2*x)/(1-x)^5, x, 0, n), x, n), n, 0, 40);

PARI

a(n)=n*(n+1)*(n+2)*(3*n+17)/24 \\ Charles R Greathouse IV, Oct 07 2015

PARI

x='x+O('x^99); concat(0, Vec(x*(5-2*x)/(1-x)^5)) \\ Altug Alkan, Apr 10 2016

Sage

[n*(n+1)*(n+2)*(3*n+17)/24 for n in (0..40)]

Formula

G.f.: x*(5 - 2*x)/(1 - x)^5.

a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5).

a(n) = A227342(A055998(n+1)).

a(n) = Sum_{j=0..n+2} (-1)^(n-j)*binomial(-j,-n-2)*S1(j,n), S1 Stirling cycle numbers A132393. - Peter Luschny, Apr 10 2016

cp's OEIS Frontend

A051877 Partial sums of A051740.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

A080852 Square array of 4D pyramidal numbers, read by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Derive

Derive

Maple

Mathematica

Formula

A093564 (7,1) Pascal triangle.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Formula

A220212 Convolution of natural numbers (A000027) with tetradecagonal numbers (A051866).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Magma

Mathematica

Formula

A241765 a(n) = n*(n + 1)*(n + 2)*(3*n + 17)/24.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Maxima

PARI

PARI

Sage

Formula

A134081 Triangle T(n, k) = binomial(n, k)*((2*k+1)*(n-k) +k+1)/(k+1), read by rows.

Views

Author

Keywords

Examples

Links

Crossrefs

A241765 a(n) = n(n + 1)(n + 2)(3n + 17)/24.

A134081 Triangle T(n, k) = binomial(n, k)((2k+1)*(n-k) +k+1)/(k+1), read by rows.