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

A050406 Partial sums of A051880.

Original entry on oeis.org

1, 16, 91, 336, 966, 2352, 5082, 10032, 18447, 32032, 53053, 84448, 129948, 194208, 282948, 403104, 562989, 772464, 1043119, 1388464, 1824130, 2368080, 3040830, 3865680, 4868955, 6080256, 7532721, 9263296
Offset: 0

Views

Author

Barry E. Williams, Dec 21 1999

Keywords

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, pp. 194-196.

Crossrefs

Cf. A051880.
Cf. A093645 ((10, 1) Pascal, column m=6).

Programs

  • GAP
    List([0..40], n-> Binomial(n+5,5)*(5*n+3)/3); # G. C. Greubel, Oct 30 2019
  • Magma
    [Binomial(n+5,5)*(5*n+3)/3: n in [0..40]]; // G. C. Greubel, Oct 30 2019
    
  • Maple
    seq(binomial(n+5,5)*(5*n+3)/3, n=0..40); # G. C. Greubel, Oct 30 2019
  • Mathematica
    Nest[Accumulate[#]&,Table[n(n+1)(10n-7)/6,{n,0,50}],3] (* Harvey P. Dale, Nov 13 2013 *)
  • PARI
    vector(41, n, binomial(n+4,5)*(5*n-2)/3) \\ G. C. Greubel, Oct 30 2019
    
  • Sage
    [binomial(n+5,5)*(5*n+3)/3 for n in (0..40)] # G. C. Greubel, Oct 30 2019
    

Formula

a(n) = C(n+5, 5)*(5*n + 3)/3.
G.f.: (1+9*x)/(1-x)^7.
E.g.f.: (360 +5400*x +10800*x^2 +6600*x^3 +1575*x^4 +153*x^5 +5*x^6) *exp(x)/360. - G. C. Greubel, Oct 30 2019

Extensions

Corrected by T. D. Noe, Nov 09 2006

A005585 5-dimensional pyramidal numbers: a(n) = n*(n+1)*(n+2)*(n+3)*(2n+3)/5!.

Original entry on oeis.org

1, 7, 27, 77, 182, 378, 714, 1254, 2079, 3289, 5005, 7371, 10556, 14756, 20196, 27132, 35853, 46683, 59983, 76153, 95634, 118910, 146510, 179010, 217035, 261261, 312417, 371287, 438712, 515592, 602888, 701624, 812889, 937839, 1077699, 1233765, 1407406
Offset: 1

Views

Author

Keywords

Comments

Convolution of triangular numbers (A000217) and squares (A000290) (n>=1). - Graeme McRae, Jun 07 2006
p^k divides a(p^k-3), a(p^k-2), a(p^k-1) and a(p^k) for prime p > 5 and integer k > 0. p^k divides a((p^k-3)/2) for prime p > 5 and integer k > 0. - Alexander Adamchuk, May 08 2007
If a 2-set Y and an (n-3)-set Z are disjoint subsets of an n-set X then a(n-5) is the number of 6-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 08 2007
5-dimensional square numbers, fourth partial sums of binomial transform of [1,2,0,0,0,...]. a(n) = Sum_{i=0..n} binomial(n+4, i+4)*b(i), where b(i)=[1,2,0,0,0,...]. - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
Antidiagonal sums of the convolution array A213550. - Clark Kimberling, Jun 17 2012
Binomial transform of (1, 6, 14, 16, 9, 2, 0, 0, 0, ...). - Gary W. Adamson, Jul 28 2015
2*a(n) is number of ways to place 4 queens on an (n+3) X (n+3) chessboard so that they diagonally attack each other exactly 6 times. The maximal possible attack number, p=binomial(k,2)=6 for k=4 queens, is achievable only when all queens are on the same diagonal. In graph-theory representation they thus form a corresponding complete graph. - Antal Pinter, Dec 27 2015
While adjusting for offsets, add A000389 to find the next in series A000389, A005585, A051836, A034263, A027800, A051843, A051877, A051878, A051879, A051880, A056118, A271567. (See Bruno Berselli's comments in A271567.) - Bruce J. Nicholson, Jun 21 2018
Coefficients in the terminating series identity 1 - 7*n/(n + 6) + 27*n*(n - 1)/((n + 6)*(n + 7)) - 77*n*(n - 1)*(n - 2)/((n + 6)*(n + 7)*(n + 8)) + ... = 0 for n = 1,2,3,.... Cf. A002415 and A040977. - Peter Bala, Feb 18 2019

Examples

			G.f. = x + 7*x^2 + 27*x^3 + 77*x^4 + 182*x^5 + 378*x^6 + 714*x^7 + 1254*x^8 + ... - _Michael Somos_, Jun 24 2018
		

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. 797.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n) = ((-1)^(n+1))*A053120(2*n+3, 5)/16, (1/16 of sixth unsigned column of Chebyshev T-triangle, zeros omitted).
Partial sums of A002415.
Cf. A006542, A040977, A047819, A111125 (third column).
Cf. a(n) = ((-1)^(n+1))*A084960(n+1, 2)/16 (compare with the first line). - Wolfdieter Lang, Aug 04 2014

Programs

  • Magma
    I:=[1, 7, 27, 77, 182, 378]; [n le 6 select I[n] else 6*Self(n-1)-15*Self(n-2)+20*Self(n-3)-15*Self(n-4)+6*Self(n-5)-Self(n-6): n in [1..40]]; // Vincenzo Librandi, Jun 09 2013
    
  • Maple
    [seq(binomial(n+2,6)-binomial(n,6), n=4..45)]; # Zerinvary Lajos, Jul 21 2006
    A005585:=(1+z)/(z-1)**6; # Simon Plouffe in his 1992 dissertation
  • Mathematica
    With[{c=5!},Table[n(n+1)(n+2)(n+3)(2n+3)/c,{n,40}]] (* or *) LinearRecurrence[ {6,-15,20,-15,6,-1},{1,7,27,77,182,378},40] (* Harvey P. Dale, Oct 04 2011 *)
    CoefficientList[Series[(1 + x) / (1 - x)^6, {x, 0, 50}], x] (* Vincenzo Librandi, Jun 09 2013 *)
  • PARI
    a(n)=binomial(n+3,4)*(2*n+3)/5 \\ Charles R Greathouse IV, Jul 28 2015

Formula

G.f.: x*(1+x)/(1-x)^6.
a(n) = 2*C(n+4, 5) - C(n+3, 4). - Paul Barry, Mar 04 2003
a(n) = C(n+3, 5) + C(n+4, 5). - Paul Barry, Mar 17 2003
a(n) = C(n+2, 6) - C(n, 6), n >= 4. - Zerinvary Lajos, Jul 21 2006
a(n) = Sum_{k=1..n} T(k)*T(k+1)/3, where T(n) = n(n+1)/2 is a triangular number. - Alexander Adamchuk, May 08 2007
a(n-1) = (1/4)*Sum_{1 <= x_1, x_2 <= n} |x_1*x_2*det V(x_1,x_2)| = (1/4)*Sum_{1 <= i,j <= n} i*j*|i-j|, where V(x_1,x_2) is the Vandermonde matrix of order 2. First differences of A040977. - Peter Bala, Sep 21 2007
a(n) = C(n+4,4) + 2*C(n+4,5). - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
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(1)=1, a(2)=7, a(3)=27, a(4)=77, a(5)=182, a(6)=378. - Harvey P. Dale, Oct 04 2011
a(n) = (1/6)*Sum_{i=1..n+1} (i*Sum_{k=1..i} (i-1)*k). - Wesley Ivan Hurt, Nov 19 2014
E.g.f.: x*(2*x^4 + 35*x^3 + 180*x^2 + 300*x + 120)*exp(x)/120. - Robert Israel, Nov 19 2014
a(n) = A000389(n+3) + A000389(n+4). - Bruce J. Nicholson, Jun 21 2018
a(n) = -a(-3-n) for all n in Z. - Michael Somos, Jun 24 2018
From Amiram Eldar, Jun 28 2020: (Start)
Sum_{n>=1} 1/a(n) = 40*(16*log(2) - 11)/3.
Sum_{n>=1} (-1)^(n+1)/a(n) = 20*(8*Pi - 25)/3. (End)
a(n) = A004302(n+1) - A207361(n+1). - J. M. Bergot, May 20 2022
a(n) = Sum_{i=0..n+1} Sum_{j=i..n+1} i*j*(j-i)/2. - Darío Clavijo, Oct 11 2023
a(n) = (A000538(n+1) - A000330(n+1))/12. - Yasser Arath Chavez Reyes, Feb 21 2024

A093645 (10,1) Pascal triangle.

Original entry on oeis.org

1, 10, 1, 10, 11, 1, 10, 21, 12, 1, 10, 31, 33, 13, 1, 10, 41, 64, 46, 14, 1, 10, 51, 105, 110, 60, 15, 1, 10, 61, 156, 215, 170, 75, 16, 1, 10, 71, 217, 371, 385, 245, 91, 17, 1, 10, 81, 288, 588, 756, 630, 336, 108, 18, 1, 10, 91, 369, 876, 1344, 1386, 966, 444, 126, 19, 1
Offset: 0

Views

Author

Wolfdieter Lang, Apr 22 2004

Keywords

Comments

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

Examples

			Triangle begins
   1;
  10,  1;
  10, 11,  1;
  10, 21, 12,  1;
  ...
		

References

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

Crossrefs

Row sums: 1 for n=0 and A005015(n-1), n >= 1, alternating row sums are 1 for n=0, 9 for n=2 and 0 otherwise.
The column sequences give for m=1..9: A017281, A051624 (12-gonal), A007587, A051799, A051880, A050406, A052254, A056125, A093646.

Programs

  • Haskell
    a093645 n k = a093645_tabl !! n !! k
    a093645_row n = a093645_tabl !! n
    a093645_tabl = [1] : iterate
                   (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [10, 1]
    -- Reinhard Zumkeller, Aug 31 2014
  • Mathematica
    t[0, 0] = 1; t[n_, k_] := Binomial[n, k] + 9*Binomial[n-1, k]; Table[t[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jul 05 2013, after Philippe Deléham *)

Formula

a(n, m) = F(10;n-m, m) for 0 <= m <= n, else 0, with F(10;0, 0)=1, F(10;n, 0)=10 if n >= 1 and F(10;n, m):=(10*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)=10 if n >= 1; a(n, m) = a(n-1, m) + a(n-1, m-1).
G.f. column m (without leading zeros): (1+9*x)/(1-x)^(m+1), m >= 0.
T(n, k) = C(n, k) + 9*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)*(10 + 21*x + 12*x^2/2! + x^3/3!) = 10 + 31*x + 64*x^2/2! + 110*x^3/3! + 170*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

A280880 Number T(n,k) of set partitions of [n] into exactly k blocks where sizes of distinct blocks are coprime; triangle T(n,k), n>=0, 0<=k<=n, read by rows.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 3, 1, 0, 1, 4, 6, 1, 0, 1, 15, 10, 10, 1, 0, 1, 6, 75, 20, 15, 1, 0, 1, 63, 21, 245, 35, 21, 1, 0, 1, 64, 476, 56, 630, 56, 28, 1, 0, 1, 171, 540, 2100, 126, 1386, 84, 36, 1, 0, 1, 130, 4185, 2640, 6930, 252, 2730, 120, 45, 1
Offset: 0

Views

Author

Alois P. Heinz, Jan 09 2017

Keywords

Examples

			T(5,1) = 1: 12345.
T(5,2) = 15: 1234|5, 1235|4, 123|45, 1245|3, 124|35, 125|34, 12|345, 1345|2, 134|25, 135|24, 13|245, 145|23, 14|235, 15|234, 1|2345.
T(5,3) = 10: 123|4|5, 124|3|5, 125|3|4, 134|2|5, 135|2|4, 1|234|5, 1|235|4, 145|2|3, 1|245|3, 1|2|345.
T(5,4) = 10: 12|3|4|5, 13|2|4|5, 1|23|4|5, 14|2|3|5, 1|24|3|5, 1|2|34|5, 15|2|3|4, 1|25|3|4, 1|2|35|4, 1|2|3|45.
T(5,5) = 1: 1|2|3|4|5.
Triangle T(n,k) begins:
  1;
  0, 1;
  0, 1,   1;
  0, 1,   3,    1;
  0, 1,   4,    6,    1;
  0, 1,  15,   10,   10,    1;
  0, 1,   6,   75,   20,   15,    1;
  0, 1,  63,   21,  245,   35,   21,    1;
  0, 1,  64,  476,   56,  630,   56,   28,   1;
  0, 1, 171,  540, 2100,  126, 1386,   84,  36,  1;
  0, 1, 130, 4185, 2640, 6930,  252, 2730, 120, 45, 1;
		

Crossrefs

T(n+k,n) for k=0-4 give: A000012, A000217, A000292, A051880(n-1) if n>0, A000389(n+4).
Row sums give A280275.
T(2n,n) gives A280889.

Programs

  • Maple
    with(numtheory):
    b:= proc(n, i, s) option remember; expand(
          `if`(n=0 or i=1, x^n, b(n, i-1, select(x->x<=i-1, s))+
          `if`(i>n or factorset(i) intersect s<>{}, 0, x*b(n-i, i-1,
          select(x->x<=i-1, s union factorset(i)))*binomial(n, i))))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(n$2, {})):
    seq(T(n), n=0..12);
  • Mathematica
    b[n_, i_, s_] := b[n, i, s] = Expand[If[n == 0 || i == 1, x^n, b[n, i - 1, Select[s, # <= i - 1 &]] + If[i > n || FactorInteger[i][[All, 1]]  ~Intersection~ s != {}, 0, x*b[n - i, i - 1, Select[ s ~Union~ FactorInteger[i][[All, 1]], # <= i - 1 &]]*Binomial[n, i]]]]; T[n_] := Function [p, Table[Coefficient[p, x, i], {i, 0, Exponent[p, x]}]][b[n, n, {}]]; Table[T[n], {n, 0, 12}] // Flatten (* Jean-François Alcover, Jan 20 2017, after Alois P. Heinz *)

A271567 Convolution of nonzero triangular numbers (A000217) and nonzero tetradecagonal numbers (A051866).

Original entry on oeis.org

1, 17, 87, 287, 742, 1638, 3234, 5874, 9999, 16159, 25025, 37401, 54236, 76636, 105876, 143412, 190893, 250173, 323323, 412643, 520674, 650210, 804310, 986310, 1199835, 1448811, 1737477, 2070397, 2452472, 2888952, 3385448, 3947944, 4582809, 5296809, 6097119
Offset: 0

Views

Author

Ilya Gutkovskiy, Apr 12 2016

Keywords

Comments

More generally, the ordinary generating function for the convolution of triangular numbers and k-gonal numbers is (1 + (k - 3)*x)/(1 - x)^6.

Crossrefs

Cf. similar sequences of the convolution of triangular numbers and k-gonal numbers: A005585 (k=4), A051836 (k=5), A034263 (k=6), A027800 (k=7), A051843 (k=8), A051877 (k=9), A051878 (k=10), A051879 (k=11), A051880 (k=12), A056118 (k=13), this sequence (k=14).

Programs

  • Magma
    /* From definition: */ P:=func; /*, where P(n, k) is the n-th k-gonal number, */ [&+[P(n+1-i, 3)*P(i, 14): i in [1..n]]: n in [1..40]]; // Bruno Berselli, Apr 18 2016
    
  • Magma
    [(n+1)*(n+2)*(n+3)*(n+4)*(12*n+5)/120: n in [0..40]]; // Bruno Berselli, Apr 18 2016
  • Mathematica
    LinearRecurrence[{6, -15, 20, -15, 6, -1}, {1, 17, 87, 287, 742, 1638}, 40]
    Table[(n + 1) (n + 2) (n + 3) (n + 4) (12 n + 5)/120, {n, 0, 40}]

Formula

O.g.f.: (1 + 11*x)/(1 - x)^6.
E.g.f.: (120 + 1920*x + 3240*x^2 + 1520*x^3 + 245*x^4 + 12*x^5)*exp(x)/120.
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)*(12*n + 5)/120.
Sum_{n>=0} 1/a(n) = 20*((15552*(6*log(2) + 3*log(3) + 2*sqrt(3)*log(2 - sqrt(3)) + (2 - sqrt(3))*Pi) - 29449)/531867) = 1.07654258697...

Extensions

Edited by Bruno Berselli, Apr 18 2016
Showing 1-5 of 5 results.