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

A001296 4-dimensional pyramidal numbers: a(n) = (3*n+1)*binomial(n+2, 3)/4. Also Stirling2(n+2, n).

Original entry on oeis.org

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
Offset: 0

Views

Author

Keywords

Comments

Permutations avoiding 12-3 that contain the pattern 31-2 exactly once.
Kekulé numbers for certain benzenoids. - Emeric Deutsch, Nov 18 2005
Partial sums of A002411. - Jonathan Vos Post, Mar 16 2006
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
Convolution of A000027 with A000326. - Bruno Berselli, Dec 06 2012
This sequence is related to A000292 by a(n) = n*A000292(n) - Sum_{i=0..n-1} A000292(i) for n>0. - Bruno Berselli, Nov 23 2017
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

Examples

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

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

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.
Cf. similar sequences listed in A241765 and A254142.
Cf. A000914.

Programs

  • Magma
    /* A000027 convolved with A000326: */ A000326:=func; [&+[(n-i+1)*A000326(i): i in [0..n]]: n in [0..40]]; // Bruno Berselli, Dec 06 2012
    
  • Magma
    [(3*n+1)*Binomial(n+2,3)/4: n in [0..40]]; // Vincenzo Librandi, Jul 30 2014
  • 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[StirlingS2[n+2, n], {n, 0, 40}] (* Jean-François Alcover, Jun 24 2015 *)
    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 *)
  • 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
    

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
a(n) = Sum_{j=0..n} j*A000217(j). - Jon Perry, Jul 28 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) = A153978(n)/2. - J. M. Bergot, Aug 09 2013
a(n) = A002817(n) + A000292(n-1). - J. M. Bergot, Aug 29 2013; [corrected by Cyril Damamme, Feb 26 2018]
a(n) = A000914(n+1) - 2 * A000330(n+1). - Antal Pinter, Dec 31 2015
a(n) = A080852(3,n-1). - R. J. Mathar, Jul 28 2016
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
a(n) = A000914(-2-n) for all n in Z. - Michael Somos, Sep 04 2017
a(n) = A000292(n) + A050534(n+1). - Cyril Damamme, Feb 26 2018
From Amiram Eldar, Jul 02 2020: (Start)
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)

A000914 Stirling numbers of the first kind: s(n+2, n).

Original entry on oeis.org

0, 2, 11, 35, 85, 175, 322, 546, 870, 1320, 1925, 2717, 3731, 5005, 6580, 8500, 10812, 13566, 16815, 20615, 25025, 30107, 35926, 42550, 50050, 58500, 67977, 78561, 90335, 103385, 117800, 133672, 151096, 170170, 190995, 213675, 238317, 265031
Offset: 0

Views

Author

Keywords

Comments

Sum of product of unordered pairs of numbers from {1..n+1}.
Number of edges of a complete k-partite graph of order k*(k+1)/2 (A000217), K_1,2,3,...,k. - Roberto E. Martinez II, Oct 18 2001
This sequence holds the x^(n-2) coefficient of the characteristic polynomial of the N X N matrix A formed by MAX(i,j), where i is the row index and j is the column index of element A[i][j], 1 <= i,j <= N. Here N >= 2. - Paul Max Payton, Sep 06 2005
The sequence contains the partial sums of A006002, which represent the areas beneath lines created by the triangular numbers plotted (t(1),t(2)) connected to (t(2),t(3)) then (t(3),t(4))...(t(n-1),t(n)) and the x-axis. - J. M. Bergot, May 05 2012
Number of functions f from [n+2] to [n+2] with f(x)=x for exactly n elements x of [n+2] and f(x)>x for exactly two elements x of [n+2]. To prove this, let the two elements of [n+2] with a larger image be labeled i and j. Note both i and j must be less than n+2. Then there are (n+2-i) choices for f(i) and (n+2-j) choices for f(j). Summing the product of the number of choices over all sets {i,j} gives us "Sum of product of unordered pairs of numbers from {1..n+1}" in the first line of the Comments Section. See the example in the Example Section below. - Dennis P. Walsh, Sep 06 2017
Zhu Shijie gives in his Magnus Opus "Jade Mirror of the Four Unknowns" the problem: "Apples are piled in the form of a triangular pyramid. The top apple is worth 2 and the price of the whole is 1320. Each apple in one layer costs 1 less than an apple in the next layer below." We find the solution 9 to this problem in this sequence 1320 = a(9). Zhu Shijie gave the solution polynomial: "Let the element tian be the number of apples in a side of the base. From the statement we have 31680 for the negative shi, 10 for the positive fang, 21 for the positive first lian, 14 for the positive last lian, and 3 for the positive yu." This translates into the polynomial equation: 3*x^4 + 14*x^3 + 21*x^2 + 10*x - 31680 = 0. - Thomas Scheuerle, Feb 10 2025

Examples

			Examples include E(K_1,2,3) = s(2+2,2) = 11 and E(K_1,2,3,4,5) = s(4+2,4) = 85, where E is the function that counts edges of graphs.
For n=2 the a(2)=11 functions f:[4]->[4] with exactly two f(x)=x and two f(x)>x are given by the 11 image vectors of form <f(1),f(2),f(3),f(4)> that follow: <1,3,4,4>, <1,4,4,4>, <2,2,4,4>, <3,2,4,4>, <4,2,4,4>, <2,3,3,4>, <2,4,3,4>, <3,3,3,4>, <3,4,3,4>, <4,3,3,4>, and <4,4,3,4>. - _Dennis P. Walsh_, Sep 06 2017
		

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. 833.
  • George E. Andrews, Number Theory, Dover Publications, New York, 1971, p. 4.
  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 227, #16.
  • F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 226.
  • H. S. Hall and S. R. Knight, Higher Algebra, Fourth Edition, Macmillan, 1891, p. 518.
  • Zhu Shijie, Jade Mirror of the Four Unknowns (Siyuan yujian), Book III Guo Duo Die Gang (Piles of Fruit), Problem number 1, 1303.
  • 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).

Crossrefs

Cf. similar sequences listed in A241765.
Cf. A001296.
Cf. A006325(n+1) (Zhu Shijie's problem number 2 uses a pyramid with square base).

Programs

  • Haskell
    a000914 n = a000914_list !! n
    a000914_list = scanl1 (+) a006002_list
    -- Reinhard Zumkeller, Mar 25 2014
    
  • Magma
    [StirlingFirst(n+2, n): n in [0..40]]; // Vincenzo Librandi, May 28 2019
  • Maple
    A000914 := n -> 1/24*(n+1)*n*(n+2)*(3*n+5);
    A000914 := proc(n)
        combinat[stirling1](n+2,n) ;
    end proc: # R. J. Mathar, May 19 2016
  • Mathematica
    Table[StirlingS1[n+2,n],{n,0,40}] (* Harvey P. Dale, Aug 24 2011 *)
    a[ n_] := n (n + 1) (n + 2) (3 n + 5) / 24; (* Michael Somos, Sep 04 2017 *)
  • PARI
    a(n)=sum(i=1,n+1,sum(j=1,n+1,i*j*(i
    				
  • PARI
    a(n)=sum(i=1,n+1,sum(j=1,i-1,i*j)) \\ Charles R Greathouse IV, Apr 07 2015
    
  • PARI
    a(n) = binomial(n+2, 3)*(3*n+5)/4 \\ Charles R Greathouse IV, Apr 07 2015
    
  • Sage
    [stirling_number1(n+2, n) for n in range(41)] # Zerinvary Lajos, Mar 14 2009
    

Formula

a(n) = binomial(n+2, 3)*(3*n+5)/4 = (n+1)*n*(n+2)*(3*n+5)/24.
E.g.f.: exp(x)*x*(48 + 84*x + 32*x^2 + 3*x^3)/24.
G.f.: (2*x+x^2)/(1-x)^5. - Simon Plouffe in his 1992 dissertation.
a(n) = Sum_{i=1..n} i*(i+1)^2/2. - Jon Perry, Jul 31 2003
a(n) = A052149(n+1)/2. - J. M. Bergot, Jun 02 2012
-(3*n+2)*(n-1)*a(n) + (n+2)*(3*n+5)*a(n-1) = 0. - R. J. Mathar, Apr 30 2015
a(n) = a(n-1) + (n+1)*binomial(n+1,2) for n >= 1. - Dennis P. Walsh, Sep 21 2015
a(n) = A001296(-2-n) for all n in Z. - Michael Somos, Sep 04 2017
From Amiram Eldar, Jan 10 2022: (Start)
Sum_{n>=1} 1/a(n) = 162*log(3)/5 - 18*sqrt(3)*Pi/5 - 384/25.
Sum_{n>=1} (-1)^(n+1)/a(n) = 36*sqrt(3)*Pi/5 - 96*log(2)/5 - 636/25. (End)
a(n) = 3*A000332(n+3) - A000292(n). - Yasser Arath Chavez Reyes, Apr 03 2024

Extensions

More terms from Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Jan 17 2000
Comments from Michael Somos, Jan 29 2000
Erroneous duplicate of the polynomial formula removed by R. J. Mathar, Sep 15 2009

A050534 Tritriangular numbers: a(n) = binomial(binomial(n,2),2) = n*(n+1)*(n-1)*(n-2)/8.

Original entry on oeis.org

0, 0, 0, 3, 15, 45, 105, 210, 378, 630, 990, 1485, 2145, 3003, 4095, 5460, 7140, 9180, 11628, 14535, 17955, 21945, 26565, 31878, 37950, 44850, 52650, 61425, 71253, 82215, 94395, 107880, 122760, 139128, 157080, 176715, 198135, 221445, 246753, 274170, 303810, 335790
Offset: 0

Views

Author

Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Dec 29 1999

Keywords

Comments

"There are n straight lines in a plane, no two of which are parallel and no three of which are concurrent. Their points of intersection being joined, show that the number of new lines drawn is (1/8)n(n-1)(n-2)(n-3)." (Schmall, 1915).
Several different versions of this sequence are possible, beginning with either one, two or three 0's.
If Y is a 3-subset of an n-set X then, for n>=6, a(n-4) is the number of (n-6)-subsets of X which have exactly one element in common with Y. - Milan Janjic, Dec 28 2007
Number of distinct ways to select 2 pairs of objects from a set of n+1 objects, when order doesn't matter. For example, with n = 3 (4 objects), the 3 possibilities are (12)(34), (13)(24), and (14)(23). - Brian Parsonnet, Jan 03 2012
Partial sums of A027480. - J. M. Bergot, Jul 09 2013
For the set {1,2,...,n}, the sum of the 2 smallest elements of all subsets with 3 elements is a(n) (see Bulut et al. link). - Serhat Bulut, Jan 20 2015
a(n) is also the number of subgroups of S_{n+1} (the symmetric group on n+1 elements) that are isomorphic to D_4 (the dihedral group of order 8). - Geoffrey Critzer, Sep 13 2015
a(n) is the coefficient of x1^(n-3)*x2^2 in exponential Bell polynomial B_{n+1}(x1,x2,...) (number of ways to select 2 pairs among n+1 objects, see above), hence its link with A000292 and A001296 (see formula). - Cyril Damamme, Feb 26 2018
Also the number of 4-cycles in the complete graph K_{n+1}. - Eric W. Weisstein, Mar 13 2018
Number of chiral pairs of colorings of the 4 edges or vertices of a square using n or fewer colors. Each member of a chiral pair is a reflection, but not a rotation, of the other. - Robert A. Russell, Oct 20 2020

Examples

			For a(3)=3, the chiral pairs of square colorings are AABC-AACB, ABBC-ACBB, and ABCC-ACCB. - _Robert A. Russell_, Oct 20 2020
		

References

  • Arthur T. Benjamin and Jennifer Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 154.
  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, Problem 1, page 72.
  • Richard P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 5.5, case k=2.

Crossrefs

Cf. A000217, A000332, A033487, A107394, A034827, A210569, Second column of triangle A001498.
Cf. similar sequences listed in A241765.
Cf. (square colorings) A006528 (oriented), A002817 (unoriented), A002411 (achiral),
Row 2 of A325006 (orthoplex facets, orthotope vertices) and A337409 (orthotope edges, orthoplex ridges).
Row 4 of A293496 (cycles of n colors using k or fewer colors).

Programs

  • GAP
    List([0..40],n->3*Binomial(n+1,4)); # Muniru A Asiru, Mar 20 2018
  • Magma
    [3*Binomial(n+1, 4): n in [0..40]]; // Vincenzo Librandi, Feb 14 2015
    
  • Maple
    [seq(binomial(n+1,4)*3,n=0..40)]; # Zerinvary Lajos, Jul 18 2006
  • Mathematica
    Table[Binomial[Binomial[n, 2], 2], {n, 0, 50}] (* Stefan Steinerberger, Apr 08 2006 *)
    LinearRecurrence[{5, -10, 10, -5, 1}, {0, 0, 0, 3, 15}, 40] (* Harvey P. Dale, Dec 14 2011 *)
    (* Start from Eric W. Weisstein, Mar 13 2018 *)
    Binomial[Binomial[Range[0, 20], 2], 2]
    Nest[Binomial[#, 2] &, Range[0, 20], 2]
    Nest[PolygonalNumber[# - 1] &, Range[0, 20], 2]
    CoefficientList[Series[3 x^3/(1 - x)^5, {x, 0, 20}], x]
    (* End *)
  • PARI
    a(n)=n*(n+1)*(n-1)*(n-2)/8 \\ Charles R Greathouse IV, Nov 20 2012
    
  • PARI
    x='x+O('x^100); concat([0, 0, 0], Vec(3*x^3/(1-x)^5)) \\ Altug Alkan, Nov 01 2015
    
  • Sage
    [(binomial(binomial(n,2),2)) for n in range(0, 39)] # Zerinvary Lajos, Nov 30 2009
    

Formula

a(n) = 3*binomial(n+1, 4) = 3*A000332(n+1).
From Vladeta Jovovic, May 03 2002: (Start)
Recurrence: a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5).
G.f.: 3*x^3 / (1-x)^5. (End)
a(n+1) = T(T(n)) - T(n); a(n+2) = T(T(n)+n) where T is A000217. - Jon Perry, Jun 11 2003
a(n+1) = T(n)^2 - T(T(n)) where T is A000217. - Jon Perry, Jul 23 2003
a(n) = T(T(n-1)-1) where T is A000217. - Jon E. Schoenfield, Dec 14 2014
a(n) = 3*C(n, 4) + 3*C(n, 3), for n>3.
From Alexander Adamchuk, Apr 11 2006: (Start)
a(n) = (1/2)*Sum_{k=1..n} k*(k-1)*(k-2).
a(n) = A033487(n-2)/2, n>1.
a(n) = C(n-1,2)*C(n+1,2)/2, n>2. (End)
a(n) = A052762(n+1)/8. - Zerinvary Lajos, Apr 26 2007
a(n) = (4x^4 - 4x^3 - x^2 + x)/2 where x = floor(n/2)*(-1)^n for n >= 0. - William A. Tedeschi, Aug 24 2010
E.g.f.: x^3*exp(x)*(4+x)/8. - Robert Israel, Nov 01 2015
a(n) = Sum_{k=1..n} Sum_{i=1..k} (n-i-1)*(n-k). - Wesley Ivan Hurt, Sep 12 2017
a(n) = A001296(n-1) - A000292(n-1). - Cyril Damamme, Feb 26 2018
Sum_{n>=3} 1/a(n) = 4/9. - Vaclav Kotesovec, May 01 2018
a(n) = A006528(n) - A002817(n) = (A006528(n) - A002411(n)) / 2 = A002817(n) - A002411(n). - Robert A. Russell, Oct 20 2020
Sum_{n>=3} (-1)^(n+1)/a(n) = 32*log(2)/3 - 64/9. - Amiram Eldar, Jan 09 2022
a(n) = Sum_{k=1..2} (-1)^(k+1)*binomial(n,2-k)*binomial(n,2+k). - Gerry Martens, Oct 09 2022

Extensions

Additional comments from Antreas P. Hatzipolakis, May 03 2002

A005286 a(n) = (n + 3)*(n^2 + 6*n + 2)/6.

Original entry on oeis.org

1, 6, 15, 29, 49, 76, 111, 155, 209, 274, 351, 441, 545, 664, 799, 951, 1121, 1310, 1519, 1749, 2001, 2276, 2575, 2899, 3249, 3626, 4031, 4465, 4929, 5424, 5951, 6511, 7105, 7734, 8399, 9101, 9841, 10620, 11439, 12299, 13201, 14146, 15135, 16169, 17249
Offset: 0

Views

Author

Keywords

Comments

Number of permutations of [n+3] with three inversions. - Michael Somos, Jun 25 2002
This sequence is related to A241765 by A241765(n) = n*a(n) - Sum_{i=0..n-1} a(i), with A241765(0)=0. For example: A241765(4) = 4*49 - (29+15+6+1) = 145. - Bruno Berselli, Apr 29 2014
For n >= 2, a(n) is also the number of multiplications between two nonzero matrix elements involved in calculating the product of an (n+1) X (n+1) Hessenberg matrix and an (n+1) X (n+1) upper triangular matrix. The formula for n X n matrices is (n+2)(n^2+4n-3)/6 multiplications, n >= 3. - John M. Coffey, Jul 18 2016

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 255, #2, b(n,3).
  • R. K. Guy, personal communication.
  • E. Netto, Lehrbuch der Combinatorik. 2nd ed., Teubner, Leipzig, 1927, p. 96.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 1, 1999; see Exercise 1.30, p. 49.

Crossrefs

Programs

  • Mathematica
    Table[(n + 3) (n^2 + 6*n + 2)/6, {n, 0, 100}] (* Vladimir Joseph Stephan Orlovsky, Jul 16 2011 *)
    LinearRecurrence[{4,-6,4,-1},{1,6,15,29},50] (* Harvey P. Dale, Mar 07 2012 *)
    Table[Binomial[n, 3] + Binomial[n, 2] - n, {n, 3, 47}] (* or *)
    CoefficientList[Series[(1 + 2 x - 3 x^2 + x^3)/(1 - x)^4, {x, 0, 44}], x] (* Michael De Vlieger, Jul 09 2016 *)
  • PARI
    a(n)=n+=3; (n^3-7*n)/6 /* Michael Somos, May 12 2005 */

Formula

G.f.: (1+2*x-3*x^2+x^3)/(1-x)^4. - Simon Plouffe in his 1992 dissertation
a(-6-n) = -a(n). - Michael Somos, May 12 2005
a(n) = a(n-1) + A000096(n+1) = A005581(n+2) - 1. - Henry Bottomley, Oct 25 2001
(m^3-7*m)/6 for m >= 3 gives the same sequence. - N. J. A. Sloane, Jul 15 2011
a(0)=1, a(1)=6, a(2)=15, a(3)=29, a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Harvey P. Dale, Mar 07 2012
E.g.f.: (6 + 30*x + 12*x^2 + x^3)*exp(x)/6. - Ilya Gutkovskiy, Jul 09 2016

A215862 Number of simple labeled graphs on n+2 nodes with exactly n connected components that are trees or cycles.

Original entry on oeis.org

0, 4, 19, 55, 125, 245, 434, 714, 1110, 1650, 2365, 3289, 4459, 5915, 7700, 9860, 12444, 15504, 19095, 23275, 28105, 33649, 39974, 47150, 55250, 64350, 74529, 85869, 98455, 112375, 127720, 144584, 163064, 183260, 205275, 229215, 255189, 283309, 313690, 346450
Offset: 0

Views

Author

Alois P. Heinz, Aug 25 2012

Keywords

Comments

Partial sums of A077414. - Bruno Berselli, Jul 30 2015

Examples

			a(1) = 4:
.1-2.  .1-2.  .1-2.  .1 2.
.|/ .  .|. .  . / .  .|/ .
.3...  .3...  .3...  .3...
		

Crossrefs

A diagonal of A215861.
Regarding the sixth formula, see similar sequences listed in A241765.

Programs

  • Maple
    a:= n-> binomial(n+2,3)*(3*n+13)/4:
    seq(a(n), n=0..40);
  • Mathematica
    Table[Binomial[n+2,3] (3n+13)/4,{n,0,40}] (* or *) LinearRecurrence[ {5,-10,10,-5,1},{0,4,19,55,125},40] (* Harvey P. Dale, Sep 10 2012 *)

Formula

G.f.: (x-4)*x/(x-1)^5.
a(n) = C(n+2,3)*(3*n+13)/4.
a(n) = 5*a(n-1)- 10*a(n-2)+ 10*a(n-3) -5*a(n-4)+a(n-5), n>4. - Harvey P. Dale, Sep 10 2012
a(n) = (1/n!) * Sum_{j=0..n} C(n,j)*(-1)^(n-j)*j^(n+1)*(j-1). - Vladimir Kruchinin, Jun 06 2013
a(n) = 4*A000332(n+2) - A000332(n+1). - R. J. Mathar, Aug 12 2013
a(n) = Sum_{i=0..n} (3+i)*A000217(i). - Bruno Berselli, Apr 29 2014

A239568 Number of ways to place 2 points on a triangular grid of side n so that they are not adjacent.

Original entry on oeis.org

0, 6, 27, 75, 165, 315, 546, 882, 1350, 1980, 2805, 3861, 5187, 6825, 8820, 11220, 14076, 17442, 21375, 25935, 31185, 37191, 44022, 51750, 60450, 70200, 81081, 93177, 106575, 121365, 137640, 155496, 175032, 196350, 219555, 244755, 272061, 301587, 333450, 367770
Offset: 2

Views

Author

Heinrich Ludwig, Mar 22 2014

Keywords

Crossrefs

Cf. A239567, A032091, A239569 (3 points), A239570 (4 points), A239571 (5 points), A282998 (6 points).
Regarding the third formula, see similar sequences listed in A241765.

Programs

  • PARI
    concat(0, Vec(3*x^3*(x-2)/(x-1)^5 + O(x^100))) \\ Colin Barker, Mar 22 2014

Formula

a(n) = n*(n-1)*(n-2)*(n+5)/8.
G.f.: 3*x^3*(x-2) / (x-1)^5. - Colin Barker, Mar 22 2014
a(n) = Sum_{i=0..n} (i+5)*A000217(i). - Bruno Berselli, Apr 29 2014
a(n) = t(t(n,k),n) + n, where t(n,k) = n*(n+1)/2 + k*n and t(n,1) = A000096(n). - Bruno Berselli, Feb 28 2017

A212343 a(n) = (n+1)*(n-2)*(n-3)/2.

Original entry on oeis.org

0, 0, 5, 18, 42, 80, 135, 210, 308, 432, 585, 770, 990, 1248, 1547, 1890, 2280, 2720, 3213, 3762, 4370, 5040, 5775, 6578, 7452, 8400, 9425, 10530, 11718, 12992, 14355, 15810, 17360, 19008, 20757, 22610, 24570, 26640, 28823, 31122, 33540, 36080, 38745, 41538, 44462, 47520, 50715, 54050, 57528
Offset: 2

Views

Author

N. J. A. Sloane, May 09 2012

Keywords

Comments

Sequence of coefficients of x^1 in marked mesh pattern generating function Q_{n,132}^(0,3,0,0)(x).
Is this row 2 of the convolution array A213819? - Clark Kimberling, Jul 04 2012

Crossrefs

Partial sums are in A241765.
Cf. similar sequences of the type m*(m+1)*(m+k)/2 listed in A267370.
Cf. also A212342.

Programs

  • Mathematica
    QQ0[t, x] = (1 - (1-4*x*t)^(1/2)) / (2*x*t); QQ1[t, x] = 1/(1 - t*QQ0[t, x]); QQ2[t, x] = (1 + t*(QQ1[t, x] - QQ0[t, x]))/(1 - t*QQ0[t, x]); QQ3[t, x] = (1 + t*(QQ2[t, x] - QQ0[t, x] + t*(QQ1[t, x] - QQ0[t,  x])))/(1 - t*QQ0[t, x]); CoefficientList[Coefficient[Simplify[Series[QQ3[t, x], {t, 0, 35}]],x],t]  (* Robert Price, Jun 04 2012 *)
    LinearRecurrence[{4,-6,4,-1},{0,0,5,18},60] (* Harvey P. Dale, Mar 15 2018 *)
  • PARI
    Vec(-x^4*(2*x-5)/(x-1)^4 + O(x^100)) \\ Colin Barker, Jul 10 2015

Formula

For n>=4, a(n) = (n-3)*A212342(n-1).
a(n) = 4*a(n-1)-6*a(n-2)+4*a(n-3)-a(n-4) for n>7. - Colin Barker, Jul 10 2015
G.f.: -x^4*(2*x-5) / (x-1)^4. - Colin Barker, Jul 10 2015
From Amiram Eldar, Apr 03 2022: (Start)
Sum_{n>=4} 1/a(n) = 23/72.
Sum_{n>=4} (-1)^n/a(n) = 4*log(2)/3 - 55/72. (End)

Extensions

a(10)-a(35) from Robert Price, Jun 02 2012
Entry revised by N. J. A. Sloane, Sep 10 2016

A269951 Triangle read by rows, T(n,k) = Sum_{j=0..n} (-1)^(n-j)*C(-j,-n)*S1(j,k), S1 the Stirling cycle numbers A132393, for n>=0 and 0<=k<=n.

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 5, 5, 1, 0, 16, 23, 9, 1, 0, 65, 116, 65, 14, 1, 0, 326, 669, 470, 145, 20, 1, 0, 1957, 4429, 3634, 1415, 280, 27, 1, 0, 13700, 33375, 30681, 14084, 3535, 490, 35, 1, 0, 109601, 283072, 284066, 147532, 43939, 7756, 798, 44, 1
Offset: 0

Views

Author

Peter Luschny, Apr 10 2016

Keywords

Examples

			Triangle begins:
  1;
  0,   1;
  0,   2,   1;
  0,   5,   5,   1;
  0,  16,  23,   9,   1;
  0,  65, 116,  65,  14,  1;
  0, 326, 669, 470, 145, 20, 1;
		

Crossrefs

A000522 (col. 1), A073596 (col. 2), A000096 (diag. n-1), A241765 (diag. n-2).
A001339 (row sums), A137597 (unsigned matrix inverse).

Programs

  • Maple
    A269951 := (n,k) -> add((-1)^(n-j)*binomial(-j,-n)*abs(Stirling1(j,k)), j=0..n):
    seq(seq(A269951(n,k), k=0..n), n=0..9);
  • Mathematica
    Flatten[ Table[ Sum[(-1)^(n-j) Binomial[-j,-n] Abs[StirlingS1[j,k]], {j,0,n}], {n,0,9}, {k,0,n}]]

A059302 A diagonal of A008296.

Original entry on oeis.org

-1, -1, 5, 25, 70, 154, 294, 510, 825, 1265, 1859, 2639, 3640, 4900, 6460, 8364, 10659, 13395, 16625, 20405, 24794, 29854, 35650, 42250, 49725, 58149, 67599, 78155, 89900, 102920, 117304, 133144, 150535, 169575, 190365, 213009, 237614, 264290, 293150, 324310, 357889, 394009, 432795, 474375
Offset: 2

Views

Author

N. J. A. Sloane, Jan 26 2001

Keywords

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 139, b(n,n-2).

Crossrefs

Cf. similar sequences listed in A241765.

Programs

  • Maple
    with(combinat): for n from 3 to 100 do for k from n-2 to n-2 do printf(`%d,`,sum(binomial(l,k)*k^(l-k)*stirling1(n,l), l=k..n)) od: od:
  • Mathematica
    f[n_] := 3*n - 1; s1 = s2 = s3 = 0; lst = {}; Do[a = f[n]; s1 += a;
    s2 += s1; s3 += s2; AppendTo[lst, s3], {n, 0, 5!}]; lst (* Vladimir Joseph Stephan Orlovsky, Jun 27 2009 *)
    Drop[CoefficientList[Series[-x^2 (1 - 4*x)/(1 - x)^5, {x, 0, 50}], x], 2] (* Vincenzo Librandi, Mar 22 2012 *)
    Rest[Table[(n - 1) n (n + 1)(3 n - 10)/24, {n, 50}]] (* or *) LinearRecurrence[ {5, -10, 10, -5, 1},{-1, -1, 5, 25, 70}, 50] (* Harvey P. Dale, Jun 05 2012 *)
  • PARI
    x='x+O('x^99); Vec(x^2*(-1+4*x)/(1-x)^5) \\ Altug Alkan, Sep 13 2017

Formula

a(n) = (n-1)n(n+1)(3n-10)/24.
If we define f(n,i,a) = Sum_{k=0..n-i} binomial(n,k)*Stirling1(n-k,i)*Product_{j=0..k-1} (-a-j), then a(n-1) = f(n,n-2,-1), for n >= 3. - Milan Janjic, Dec 20 2008
G.f.: -x^2*(1-4*x)/(1-x)^5. - Colin Barker, Mar 21 2012
a(2)=-1, a(3)=-1, a(4)=5, a(5)=25, a(6)=70, a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). - Harvey P. Dale, Jun 05 2012
a(n) = Sum_{k=1..n} Sum_{i=1..k} (n-i)*(n-k-1). - Wesley Ivan Hurt, Sep 12 2017

Extensions

More terms from James Sellers, Jan 26 2001

A118788 Triangle where T(n,k) = n!/(n-k)!*[x^k] ( x/(2*x + log(1-x)) )^(n+1), for n>=k>=0, read by rows.

Original entry on oeis.org

1, 1, 1, 1, 3, 5, 1, 6, 23, 41, 1, 10, 65, 255, 469, 1, 15, 145, 930, 3679, 6889, 1, 21, 280, 2590, 16429, 65247, 123605, 1, 28, 490, 6090, 54789, 344694, 1371887, 2620169, 1, 36, 798, 12726, 151599, 1338330, 8367785, 33347535, 64074901, 1, 45, 1230, 24360
Offset: 0

Views

Author

Paul D. Hanna, Apr 29 2006

Keywords

Comments

Row sums are A118789, where Sum_{n>=0} A118789(n)*x^n/n! = exp( Sum_{n>=1} A032188(n)*x^n/n! ).
Main diagonal is A032188(n) = number of labeled series-reduced mobiles (circular rooted trees) with n leaves.
Secondary diagonal is A118790.

Examples

			Row sums e.g.f. equals the exponential of the diagonal e.g.f.:
1 + x + 2*x^2/2! + 9*x^3/3! + 71*x^4/4! +...+ A118789(n)*x^n/n! +...
= exp(x + x^2/2! + 5*x^3/3! + 41*x^4/4! +...+ A032188(n)*x^n/n! +...).
Triangle begins:
  1;
  1, 1;
  1, 3, 5;
  1, 6, 23, 41;
  1, 10, 65, 255, 469;
  1, 15, 145, 930, 3679, 6889;
  1, 21, 280, 2590, 16429, 65247, 123605;
  1, 28, 490, 6090, 54789, 344694, 1371887, 2620169;
  1, 36, 798, 12726, 151599, 1338330, 8367785, 33347535, 64074901;
  ...
Triangle is formed from powers of F(x) = x/(2*x + log(1-x)):
  F(x)^1 = (1) + 1/2*x + 7/12*x^2 + 17/24*x^3 + 629/720*x^4 +...
  F(x)^2 = (1 + x) + 17/12*x^2 + 2*x^3 + 671/240*x^4 +...
  F(x)^3 = (1 + 3/2*x + 5/2*x^2) + 4*x^3 + 1489/240*x^4 +...
  F(x)^4 = (1 + 6/3*x + 23/6*x^2 + 41/6*x^3) + 8351/720*x^4 +...
  F(x)^5 = (1 + 10/4*x + 65/12*x^2 + 255/24*x^3 + 469/24*x^4) +...
		

Crossrefs

Third column is A241765.

Programs

  • PARI
    {T(n,k)=local(x=X+X^2*O(X^(k+2)));n!/(n-k)!*polcoeff((x/(2*x+log(1-x)))^(n+1),k,X)}

Formula

Main diagonal has e.g.f.: series_reversion[2*x+log(1-x)].
Conjecture: T(n,k) = Sum_{j=0..k} binomial(n+j, n-k)*A269940(k, j) for 0 <= k <= n. - Mikhail Kurkov, Feb 17 2025
Showing 1-10 of 10 results.