cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-9 of 9 results.

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

A103371 Number triangle T(n,k) = C(n,n-k)*C(n+1,n-k).

Original entry on oeis.org

1, 2, 1, 3, 6, 1, 4, 18, 12, 1, 5, 40, 60, 20, 1, 6, 75, 200, 150, 30, 1, 7, 126, 525, 700, 315, 42, 1, 8, 196, 1176, 2450, 1960, 588, 56, 1, 9, 288, 2352, 7056, 8820, 4704, 1008, 72, 1, 10, 405, 4320, 17640, 31752, 26460, 10080, 1620, 90, 1, 11, 550, 7425, 39600, 97020
Offset: 0

Views

Author

Paul Barry, Feb 03 2005

Keywords

Comments

Columns include A000027, A002411, A004302, A108647, A134287. Row sums are C(2n+1,n+1) or A001700.
T(n-1,k-1) is the number of ways to put n identical objects into k of altogether n distinguishable boxes. See the partition array A035206 from which this triangle arises after summing over all entries related to partitions with fixed part number k.
T(n, k) is also the number of order-preserving full transformations (of an n-chain) of height k (height(alpha) = |Im(alpha)|). - Abdullahi Umar, Oct 02 2008
The o.g.f. of the (n+1)-th diagonal is given by G(n, x) = (n+1)*Sum_{k=1..n} A001263(n, k)*x^(k-1) / (1 - x)^(2*n+1), for n >= 1 and for n = 0 it is G(0, x) = 1/(1-x). - Wolfdieter Lang, Jul 31 2017

Examples

			The triangle T(n, k) begins:
n\k  0   1    2     3     4     5     6    7  8 9 ...
0:   1
1:   2   1
2:   3   6    1
3:   4  18   12     1
4:   5  40   60    20     1
5:   6  75  200   150    30     1
6:   7 126  525   700   315    42     1
7:   8 196 1176  2450  1960   588    56    1
8:   9 288 2352  7056  8820  4704  1008   72  1
9:  10 405 4320 17640 31752 26460 10080 1620 90 1
...  reformatted. - _Wolfdieter Lang_, Jul 31 2017
From _R. J. Mathar_, Mar 29 2013: (Start)
The matrix inverse starts
       1;
      -2,       1;
       9,      -6,      1;
     -76,      54,    -12,      1;
    1055,    -760,    180,    -20,   1;
  -21906,   15825,  -3800,    450, -30,   1;
  636447, -460026, 110775, -13300, 945, -42, 1; (End)
O.g.f. of 4th diagonal [4, 40,200, ...] is G(3, x) = 4*(1 + 3*x + x^2)/(1 - x)^7, from the n = 3 row [1, 3, 1] of A001263. See a comment above. - _Wolfdieter Lang_, Jul 31 2017
		

Crossrefs

Cf. A007318, A000894 (central terms), A132813 (mirrored).

Programs

  • Haskell
    a103371 n k = a103371_tabl !! n !! k
    a103371_row n = a103371_tabl !! n
    a103371_tabl = map reverse a132813_tabl
    -- Reinhard Zumkeller, Apr 04 2014
    
  • Magma
    /* As triangle */ [[Binomial(n,n-k)*Binomial(n+1,n-k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Aug 01 2017
    
  • Maple
    A103371 := (n,k) -> binomial(n,k)^2*(n+1)/(k+1);
    seq(print(seq(A103371(n, k), k=0..n)), n=0..7); # Peter Luschny, Oct 19 2011
  • Mathematica
    Flatten[Table[Binomial[n,n-k]Binomial[n+1,n-k],{n,0,10},{k,0,n}]] (* Harvey P. Dale, May 26 2014 *)
    CoefficientList[Series[Series[E^(x(1+y))(BesselI[0,2*x*Sqrt[y]]+BesselI[1,2*x*Sqrt[y]]/Sqrt[y]),{x,0,8}],{y,0,8}],{x,y}]*Range[0,8]! (* Natalia L. Skirrow, Apr 14 2025 *)
  • Maxima
    create_list(binomial(n,k)*binomial(n+1,k+1),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */
    
  • PARI
    for(n=0,10, for(k=0,n, print1(binomial(n,k)*binomial(n+1,k+1), ", "))) \\ G. C. Greubel, Nov 09 2018

Formula

Number triangle T(n, k) = C(n, n-k)*C(n+1, n-k) = C(n, k)*C(n+1, k+1); Column k of this triangle has g.f. Sum_{j=0..k} (C(k, j)*C(k+1, j) * x^(k+j))/(1-x)^(2*k+2); coefficients of the numerators are the rows of the reverse triangle C(n, k)*C(n+1, k).
T(n,k) = C(n, k)*Sum_{j=0..(n-k)} C(n-j, k). - Paul Barry, Jan 12 2006
T(n,k) = (n+1-k)*N(n+1,k+1), with N(n,k):=A001263(n,k), the Narayana triangle (with offset [1,1]).
O.g.f.: ((1-(1-y)*x)/sqrt((1-(1+y)*x)^2-4*x^2*y) -1)/2, (from o.g.f. of A001263, Narayana triangle). - Wolfdieter Lang, Nov 13 2007
From Peter Bala, Jan 24 2008: (Start)
Matrix product of A007318 and A122899.
O.g.f. for row n: (1-x)^n*P(n,1,0,(1+x)/(1-x)) = 1/(2*x)*(1-x)^(n+1)*( Legendre_P(n+1,(1+x)/(1-x)) - Legendre_P(n,(1+x)/(1-x)) ), where P(n,a,b,x) denotes the Jacobi polynomial.
O.g.f. for column k: x^k/(1-x)^(k+2)*P(k,0,1,(1+x)/(1-x)). Compare with A008459. (End)
Let S(n,k) = binomial(2*n,n)^(k+1)*((n+1)^(k+1)-n^(k+1))/(n+1)^k. Then T(2*n,n) = S(n,1). (Cf. A194595, A197653, A197654). - Peter Luschny, Oct 20 2011
T(n,k) = A003056(n+1,k+1)*C(n,k)^2/(k+1). - Peter Luschny, Oct 29 2011
T(n,k) = A007318(n, k)*A135278(n, k), n >= k >= 0. - Wolfdieter Lang, Jul 31 2017
From Natalia L. Skirrow, Apr 14 2025: (Start)
T(n,k) = A008459(n,k) + n*N(n,k+1).
E.g.f.: e^(x*(1+y))*(I_0(2*x*sqrt(y)) + I_1(2*x*sqrt(y))/sqrt(y)), where I_n is the modified Bessel function of the first kind. (The I_0 contributes A008459(n,k), the I_1 contributes n*N(n,k+1))
O.g.f. for row n: (n+1)*2F1(-n,-n;2;y) = (n+1)*2F1(2+n,2+n;2;y)*(1-y)^(2*(n+1)) (by Euler's hypergeometric transformation); (n+1)*2F1(2+n,2+n;2;y) is the o.g.f. for row n of (k+n+1)!^2/(k!*(k+1)!*n!*(n+1)!), which is column n+1 of A132812.
O.g.f. for column k: 2F1(1+k,2+k;1;x)*x^k = 2F1(-k,-1-k;1;x)*x^k/(1-x)^(2+2*k). 2F1(-k,-1-k;1;x) is the kth row of A132813, the reflection of the kth row of this triangle.
O.g.f. for diagonal d (beginning at a(d,0)): (d+1)*x^d*2F1(d+1,d+2;2;x*y). 2F1(d+1,d+2;2;x) = 2F1(1-d,-d;2;x)/(1-x)^(2*d+1), numerator being the o.g.f. of row d of the Narayana triangle.
These respectively yield:
T(n,k) = Sum_{i=0..n+k} C(2*(n+1),i)*(-1)^i*A132812(n+1+k-i,n+1),
T(d+k,k) = Sum_{i=0..k} C(d-i+1+2*k,d-i)*T(k,k-i),
T(d+k,k) = Sum_{i=0..d} C(k-i + 2*d,k-i)*N(d,i+1)*(d+1).
E.g.f. for column k: 1F1(2+k;1;x)*x^k/k!.
E.g.f. for diagonal d: (d+1)*x^d*1F1(d+2;2;x*y)/d!. (End)

A088459 Triangle read by rows: T(n,k) represents the number of lozenge tilings of an (n,1,n)-hexagon which include the non-vertical tile above the main diagonal starting in position k+1.

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 1, 3, 6, 6, 3, 1, 1, 4, 12, 18, 18, 12, 4, 1, 1, 5, 20, 40, 60, 60, 40, 20, 5, 1, 1, 6, 30, 75, 150, 200, 200, 150, 75, 30, 6, 1, 1, 7, 42, 126, 315, 525, 700, 700, 525, 315, 126, 42, 7, 1, 1, 8, 56, 196, 588, 1176, 1960, 2450, 2450, 1960, 1176, 588, 196, 56, 8, 1
Offset: 1

Views

Author

Christopher Hanusa (chanusa(AT)washington.edu), Nov 14 2003

Keywords

Comments

Rows are of length 2, 4, 6, 8, 10, 12, ...
T(n,k)= number of symmetric Dyck paths of length 4n and having k peaks. Example: T(2,3)=2 because we have UU*DU*DU*DD and U*DUU*DDU*D, where U=(1,1), D=(1,-1) and * shows the peaks. - Emeric Deutsch, Feb 22 2004
T(n,k) is also the number of nodes at distance k from a specified node in the n-odd graph for k in 1..n-1. - Eric W. Weisstein, Mar 23 2018
T(n,k) seems to be the k-th Lie-Betti number of the star graph on n vertices. See A360571 for additional information and references. - Samuel J. Bevins, Feb 12 2023

Examples

			For example, the number of tilings of a 4,1,4 hexagon which includes the non-vertical tile above the main diagonal starting in position 3 is T(4,2)=12.
Triangle T(n, k) begins:
[1] 1,1,
[2] 1,2, 2,  1,
[3] 1,3, 6,  6,   3,   1,
[4] 1,4,12, 18,  18,  12,   4,   1,
[5] 1,5,20, 40,  60,  60,  40,  20,   5,   1,
[6] 1,6,30, 75, 150, 200, 200, 150,  75,  30,   6,  1,
[7] 1,7,42,126, 315, 525, 700, 700, 525, 315, 126, 42,    7,   1,
[8] 1,8,56,196, 588,1176,1960,2450,2450,1960,1176,588,  196,  56,  8, 1,
[9] 1,9,72,288,1008,2352,4704,7056,8820,8820,7056,4704,2352,1008,288,72,9,1
		

Crossrefs

Columns 0-5 are sequences A000012, A000027, A002378, A002411, A006011 and A004302.
Cf. A000984 (row sums).

Programs

  • Maple
    A088459 := proc(n,k)
        binomial(n,ceil(k/2))*binomial(n-1,floor(k/2)) ;
    end proc:
    seq(seq(A088459(n,k),k=0..2*n-1),n=1..10) ; # R. J. Mathar, Apr 02 2017
  • Mathematica
    Table[Binomial[n, Ceiling[k/2]] Binomial[n - 1, Floor[k/2]], {n, 10}, {k, 0, 2 n - 1}] // Flatten (* Eric W. Weisstein, Mar 23 2018 *)

Formula

T(n, k) = binomial(n, ceiling(k/2))* binomial(n-1, floor(k/2)), n>0 and k=0 to 2n-1.

Extensions

Edited and extended by Ray Chandler, Nov 17 2003

A132813 Triangle read by rows: A001263 * A127648 as infinite lower triangular matrices.

Original entry on oeis.org

1, 1, 2, 1, 6, 3, 1, 12, 18, 4, 1, 20, 60, 40, 5, 1, 30, 150, 200, 75, 6, 1, 42, 315, 700, 525, 126, 7, 1, 56, 588, 1960, 2450, 1176, 196, 8, 1, 72, 1008, 4704, 8820, 7056, 2352, 288, 9, 1, 90, 1620, 10080, 26460, 31752, 17640, 4320, 405, 10
Offset: 0

Views

Author

Gary W. Adamson, Sep 01 2007

Keywords

Comments

Also T(n,k) = binomial(n-1, k-1)*binomial(n, k-1), related to Narayana polynomials (see Sulanke reference). - Roger L. Bagula, Apr 09 2008
h-vector for cluster complex associated to the root system B_n. See p. 8, Athanasiadis and C. Savvidou. - Tom Copeland, Oct 19 2014

Examples

			First few rows of the triangle are:
  1;
  1,  2;
  1,  6,   3;
  1, 12,  18,   4;
  1, 20,  60,  40,   5;
  1, 30, 150, 200,  75,   6;
  1, 42, 315, 700, 525, 126, 7;
  ...
		

Crossrefs

Family of polynomials (see A062145): A008459 (c=1), this sequence (c=2), A062196 (c=3), A062145 (c=4), A062264 (c=5), A062190 (c=6).
Columns: A000012 (k=0), A002378 (k=1), A006011 (k=2), 4*A006542 (k=3), 5*A006857 (k=4), 6*A108679 (k=5), 7*A134288 (k=6), 8*A134289 (k=7), 9*A134290 (k=8), 10*A134291 (k=9).
Diagonals: A000027 (k=n), A002411 (k=n-1), A004302 (k=n-2), A108647 (k=n-3), A134287 (k=n-4).
Main diagonal: A000894.
Sums: (-1)^floor((n+1)/2)*A001405 (signed row), A001700 (row), A203611 (diagonal).
Cf. A103371 (mirrored).

Programs

  • GAP
    Flat(List([0..10],n->List([0..n], k->(k+1)*Binomial(n+1,k+1)*Binomial(n+1,k)/(n+1)))); # Muniru A Asiru, Feb 26 2019
    
  • Haskell
    a132813 n k = a132813_tabl !! n !! k
    a132813_row n = a132813_tabl !! n
    a132813_tabl = zipWith (zipWith (*)) a007318_tabl $ tail a007318_tabl
    -- Reinhard Zumkeller, Apr 04 2014
    
  • Magma
    /* triangle */ [[(k+1)*Binomial(n+1,k+1)*Binomial(n+1,k)/(n+1): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Oct 19 2014
    
  • Maple
    P := (n, x) -> hypergeom([1-n, -n], [1], x): for n from 1 to 9 do PolynomialTools:-CoefficientList(simplify(P(n,x)),x) od; # Peter Luschny, Nov 26 2014
  • Mathematica
    T[n_,k_]=Binomial[n-1,k-1]*Binomial[n,k-1]; Table[Table[T[n,k],{k,1,n}],{n,1,11}]; Flatten[%] (* Roger L. Bagula, Apr 09 2008 *)
    P[n_, x_] := HypergeometricPFQ[{1-n, -n}, {1}, x]; Table[CoefficientList[P[n, x], x], {n, 1, 10}] // Flatten (* Jean-François Alcover, Nov 27 2014, after Peter Luschny *)
  • PARI
    tabl(nn) = {for (n = 1, nn, for (k = 1, n, print1(binomial(n-1, k-1)*binomial(n, k-1) , ", ");););} \\ Michel Marcus, Feb 12 2014
    
  • SageMath
    def A132813(n,k): return binomial(n,k)*binomial(n+1,k)
    print(flatten([[A132813(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Mar 12 2025

Formula

T(n,k) = (k+1)*binomial(n+1,k+1)*binomial(n+1,k)/(n+1), n >= k >= 0.
From Roger L. Bagula, May 14 2010: (Start)
T(n, m) = coefficients(p(x,n)), where
p(x,n) = (1-x)^(2*n)*Sum_{k >= 0} binomial(k+n-1, k)*binomial(n+k, k)*x^k,
or p(x,n) = (1-x)^(2*n)*Hypergeometric2F1([n, n+1], [1], x). (End)
T(n,k) = binomial(n,k) * binomial(n+1,k). - Reinhard Zumkeller, Apr 04 2014
These are the coefficients of the polynomials Hypergeometric2F1([1-n,-n], [1], x). - Peter Luschny, Nov 26 2014
G.f.: A(x,y) = A281260(x,y)/(1-A281260(x,y))/x. - Vladimir Kruchinin, Oct 10 2020

A128629 A triangular array generated by moving Pascal sequences to prime positions and embedding new sequences at the nonprime locations. (cf. A007318 and A000040).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 6, 4, 1, 1, 4, 9, 10, 5, 1, 1, 6, 10, 16, 15, 6, 1, 1, 5, 18, 20, 25, 21, 7, 1, 1, 8, 15, 40, 35, 36, 28, 8, 1, 1, 9, 27, 35, 75, 56, 49, 36, 9, 1
Offset: 1

Views

Author

Alford Arnold, Mar 29 2007

Keywords

Comments

The array can be constructed by beginning with A007318 (Pascal's triangle) placing each diagonal on a prime row. The other rows are filled in by mapping the prime factorization of the row number to the known sequences on the prime rows and multiplying term by term.

Examples

			Row six begins 1 6 18 40 75 126 ... because rows two and three are
1 2 3 4 5 6 ...
1 3 6 10 15 21 ...
The array begins
1 1 1 1 1 1 1 1 1 A000012
1 2 3 4 5 6 7 8 9 A000027
1 3 6 10 15 21 28 36 45 A000217
1 4 9 16 25 36 49 64 81 A000290
1 4 10 20 35 56 84 120 165 A000292
1 6 18 40 75 126 196 288 405 A002411
1 5 15 35 70 126 210 330 495 A000332
1 8 27 64 125 216 343 512 729 A000578
1 9 36 100 225 441 784 1296 2025 A000537
1 8 30 80 175 336 588 960 1485 A002417
1 6 21 56 126 252 462 792 1287 A000389
1 12 54 160 375 756 1372 2304 3645 A019582
1 7 28 84 210 462 924 1716 3003 A000579
1 10 45 140 350 756 1470 2640 4455 A027800
1 12 60 200 525 1176 2352 4320 7425 A004302
1 16 81 256 625 1296 2401 4096 6561 A000583
1 8 36 120 330 792 1716 3432 6435 A000580
1 18 108 400 1125 2646 5488 10368 18225 A019584
1 9 45 165 495 1287 3003 6435 12870 A000581
1 16 90 320 875 2016 4116 7680 13365 A119771
1 15 90 350 1050 2646 5880 11880 22275 A001297
1 12 63 224 630 1512 3234 6336 11583 A027810
1 10 55 220 715 2002 5005 11440 24310 A000582
1 24 162 640 1875 4536 9604 18432 32805 A019583
1 16 100 400 1225 3136 7056 14400 27225 A001249
1 14 84 336 1050 2772 6468 13728 27027 A027818
1 27 216 1000 3375 9261 21952 46656 91125 A059827
1 20 135 560 1750 4536 10290 21120 40095 A085284
		

Crossrefs

Cf. A064553 (second diagonal), A080688 (second diagonal resorted).

Programs

  • Maple
    A128629 := proc(n,m) if n = 1 then 1; elif isprime(n) then p := numtheory[pi](n) ; binomial(p+m-1,p) ; else a := 1 ; for p in ifactors(n)[2] do a := a* procname(op(1,p),m)^ op(2,p) ; od: fi; end: # R. J. Mathar, Sep 09 2009

Extensions

A-number added to each row of the examples by R. J. Mathar, Sep 09 2009

A004282 a(n) = n*(n+1)^2*(n+2)^2/12.

Original entry on oeis.org

0, 3, 24, 100, 300, 735, 1568, 3024, 5400, 9075, 14520, 22308, 33124, 47775, 67200, 92480, 124848, 165699, 216600, 279300, 355740, 448063, 558624, 690000, 845000, 1026675, 1238328, 1483524, 1766100, 2090175
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

Formula

a(n) = C(2+n, 2)*C(2+n, 3) = A000217(n+1)*A000292(n). - Zerinvary Lajos, Jan 10 2006
a(n-1) = Sum_{1 <= x_1, x_2 <= n} x_1*(det V(x_1,x_2))^2 = Sum_{1 <= i,j <= n} i*(i-j)^2, where V(x_1,x_2) is the Vandermonde matrix of order 2. - Peter Bala, Sep 21 2007
G.f.: x*(3+6*x+x^2)/(1-x)^6. - Colin Barker, Feb 09 2012
a(n) = Sum_{k=0..n} Sum_{i=0..n} (n-i+1) * C(k+1,k-1). - Wesley Ivan Hurt, Sep 21 2017
a(n) = A004302(n+1) - A000537(n+1). - J. M. Bergot, Mar 28 2018
From Amiram Eldar, May 29 2022: (Start)
Sum_{n>=1} 1/a(n) = 30 - 3*Pi^2.
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/2 - 24*log(2) + 12. (End)

A028725 a(n) = floor(n/2) * floor((n-1)/2) * floor((n-2)/2) * floor((n-3)/2) * floor((n-4)/2) / 12.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 3, 12, 24, 60, 100, 200, 300, 525, 735, 1176, 1568, 2352, 3024, 4320, 5400, 7425, 9075, 12100, 14520, 18876, 22308, 28392, 33124, 41405, 47775, 58800, 67200, 81600, 92480, 110976, 124848, 148257, 165699, 194940, 216600, 252700, 279300
Offset: 0

Views

Author

Keywords

Crossrefs

Bisections: A004282, A004302.

Programs

  • Magma
    [(&*[Floor((n-j)/2):j in [0..4]])/12: n in [0..60]]; // G. C. Greubel, Apr 08 2022
    
  • Mathematica
    Table[(Times@@Floor/@(n/2-Range[0,4]/2))/12,{n,0,50}] (* or *) LinearRecurrence[ {1,5,-5,-10,10,10,-10,-5,5,1,-1}, {0,0,0,0,0,0,1,3,12,24,60}, 50] (* Harvey P. Dale, Jun 26 2012 *)
  • PARI
    concat([0,0,0,0,0,0], Vec(x^6*(x^4+2*x^3+4*x^2+2*x+1)/((x-1)^6*(x+1)^5) + O(x^100))) \\ Colin Barker, Mar 01 2015
    
  • SageMath
    [(1/768)*((-1)^n*(45 -65*n +38*n^2 -10*n^3 +n^4) -45 +193*n -230*n^2 +114*n^3 -25*n^4 +2*n^5) for n in (0..60)] # G. C. Greubel, Apr 08 2022

Formula

a(n) = a(n-1) + 5*a(n-2) - 5*a(n-3) - 10*a(n-4) + 10*a(n-5) + 10*a(n-6) - 10*a(n-7) - 5*a(n-8) + 5*a(n-9) + a(n-10) - a(n-11), with a(0)=0, a(1)=0, a(2)=0, a(3)=0, a(4)=0, a(5)=0, a(6)=1, a(7)=3, a(8)=12, a(9)=24, a(10)=60. - Harvey P. Dale, Jun 26 2012
G.f.: x^6*(1+2*x+4*x^2+2*x^3+x^4)/((1-x)^6*(1+x)^5). - Colin Barker, Mar 01 2015
From R. J. Mathar, Sep 23 2021: (Start)
a(2*n+1) = A004282(n-2).
a(2*n) = A004302(n-2).
a(n) = A028724(n)*A002620(n-4)/6. (End)
From G. C. Greubel, Apr 08 2022: (Start)
a(n) = (1/768)*((-1)^n*(45 -65*n +38*n^2 -10*n^3 +n^4) -45 +193*n -230*n^2 +114*n^3 -25*n^4 +2*n^5).
E.g.f.: (1/768)*((45 +36*x +15*x^2 +4*x^3 +x^4)*exp(-x) + (-45 +54*x -33*x^2 + 14*x^3 -5*x^4 +2*x^5)*exp(x)). (End)

A103252 Array A000292(n)*A000217(k), read by antidiagonals.

Original entry on oeis.org

1, 4, 3, 10, 12, 6, 20, 30, 24, 10, 35, 60, 60, 40, 15, 56, 105, 120, 100, 60, 21, 84, 168, 210, 200, 150, 84, 28, 120, 252, 336, 350, 300, 210, 112, 36, 165, 360, 504, 560, 525, 420, 280, 144, 45, 220, 495, 720, 840, 840, 735, 560, 360, 180, 55
Offset: 1

Views

Author

Gary W. Adamson, Mar 20 2005

Keywords

Examples

			Array begins
   1,   3,   6,  10,  15,  21,  28,   36,   45,   55, ...
   4,  12,  24,  40,  60,  84, 112,  144,  180,  220, ...
  10,  30,  60, 100, 150, 210, 280,  360,  450,  550, ...
  20,  60, 120, 200, 300, 420, 560,  720,  900, 1100, ...
  35, 105, 210, 350, 525, 735, 980, 1260, 1575, 1925, ...
  ...
		

Crossrefs

Cf. A000579 (antidiagonal sums).
Main diagonal gives A004302.

Programs

  • Mathematica
    A[n_,k_]:=Binomial[n+2,3]Binomial[k+1,2]; Table[A[n-k+1,k],{n,10},{k,n}]//Flatten (* Stefano Spezia, May 21 2023 *)

Formula

G.f.: x*y/((1 - x)^4*(1 - y)^3). - Stefano Spezia, May 21 2023

Extensions

More terms from Stefano Spezia, May 21 2023

A192832 Molecular topological indices of the lattice graphs.

Original entry on oeis.org

0, 48, 576, 2880, 9600, 25200, 56448, 112896, 207360, 356400, 580800, 906048, 1362816, 1987440, 2822400, 3916800, 5326848, 7116336, 9357120, 12129600, 15523200, 19636848, 24579456, 30470400, 37440000, 45630000, 55194048, 66298176, 79121280, 93855600
Offset: 1

Views

Author

Eric W. Weisstein, Jul 11 2011

Keywords

Comments

Lattice graphs are defined for n>=2; extended to n=1 using closed form.

Programs

  • GAP
    List([0..30], n -> 4*n^2*(n+1)*(n-1)^2); # G. C. Greubel, Jan 04 2019
  • Magma
    [4*n^2*(n+1)*(n-1)^2: n in [1..30]]; // G. C. Greubel, Jan 04 2019
    
  • Mathematica
    Table[4*n^2*(n+1)*(n-1)^2, {n,1,30}] (* G. C. Greubel, Jan 04 2019 *)
  • PARI
    vector(30, n, 4*n^2*(n+1)*(n-1)^2) \\ G. C. Greubel, Jan 04 2019
    
  • Sage
    [4*n^2*(n+1)*(n-1)^2 for n in (1..30)] # G. C. Greubel, Jan 04 2019
    

Formula

a(n) = 4*n^2*(n+1)*(n-1)^2.
a(n) = 48*A004302(n).
G.f.: 48*x^2*(1+6*x+3*x^2)/(1-x)^6. - Colin Barker, Aug 07 2012
E.g.f.: 4*x^2*(6 +18*x +9*x^2 +x^3)*exp(x). - G. C. Greubel, Jan 04 2019
Showing 1-9 of 9 results.