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.

Previous Showing 11-20 of 49 results. Next

A000574 Coefficient of x^5 in expansion of (1 + x + x^2)^n.

Original entry on oeis.org

3, 16, 51, 126, 266, 504, 882, 1452, 2277, 3432, 5005, 7098, 9828, 13328, 17748, 23256, 30039, 38304, 48279, 60214, 74382, 91080, 110630, 133380, 159705, 190008, 224721, 264306, 309256, 360096, 417384, 481712, 553707, 634032, 723387, 822510
Offset: 3

Views

Author

Keywords

Comments

If Y is a 3-subset of an n-set X then, for n>=7, a(n-4) is the number of 5-subsets of X having at most one element in common with Y. - Milan Janjic, Nov 23 2007

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 78.
  • 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

Column m=5 of (1, 3) Pascal triangle A095660.

Programs

  • Magma
    [3*Binomial(n+2,5)-2*Binomial(n+1,5): n in [3..50]]; // Vincenzo Librandi, Jun 10 2012
    
  • Maple
    A000574:=-(-3+2*z)/(z-1)**6; # conjectured by Simon Plouffe in his 1992 dissertation
    seq(3*binomial(n+2,5)-2*binomial(n+1,5),n=3..100); # Robert Israel, Aug 04 2015
    A000574 := n -> GegenbauerC(`if`(5A000574(n)), n=3..20); # Peter Luschny, May 10 2016
  • Mathematica
    CoefficientList[Series[(3-2*x)/(1-x)^6,{x,0,40}],x] (* Vincenzo Librandi, Jun 10 2012 *)
  • PARI
    x='x+O('x^50); Vec(x^3*(3-2*x)/(1-x)^6) \\ G. C. Greubel, Nov 22 2017

Formula

G.f.: x^3*(3-2*x)/(1-x)^6.
a(n) = 3*binomial(n+2,5) - 2*binomial(n+1,5).
a(n) = A111808(n,5) for n>4. - Reinhard Zumkeller, Aug 17 2005
a(n) = binomial(n+1, 4)*(n+12)/5 = 3*b(n-3)-2*b(n-4), with b(n)=binomial(n+5, 5); cf. A000389.
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). - Vincenzo Librandi, Jun 10 2012
a(n) = 3*binomial(n, 3) + 4*binomial(n, 4) + binomial(n, 5). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
a(n) = GegenbauerC(N, -n, -1/2) where N = 5 if 5Peter Luschny, May 10 2016
a(n) = Sum_{i=1..n-1} A000217(i)*A055998(n-1-i). - Bruno Berselli, Mar 05 2018
E.g.f.: exp(x)*x^3*(60 + 20*x + x^2)/120. - Stefano Spezia, Jul 09 2023

Extensions

More terms from Vladeta Jovovic, Oct 02 2000

A005582 a(n) = n*(n+1)*(n+2)*(n+7)/24.

Original entry on oeis.org

0, 2, 9, 25, 55, 105, 182, 294, 450, 660, 935, 1287, 1729, 2275, 2940, 3740, 4692, 5814, 7125, 8645, 10395, 12397, 14674, 17250, 20150, 23400, 27027, 31059, 35525, 40455, 45880, 51832, 58344, 65450, 73185, 81585, 90687, 100529, 111150, 122590, 134890
Offset: 0

Views

Author

Keywords

Comments

a(n) = number of Dyck (n+2)-paths with exactly 2 rows of peaks. A row of peaks is a maximal sequence of peaks all at the same height and 2 units apart. For example, UDUDUD ( = /\/\/\ ) contains exactly one row of peaks, as does UUUDDD, but UDUUDDUD has three and a(1)=2 counts UDUUDD, UUDDUD. - David Callan, Mar 02 2005
If X is an n-set and Y a fixed 2-subset of X then a(n-4) is equal to the number of (n-4)-subsets of X intersecting Y. - Milan Janjic, Jul 30 2007
Let I=I_n be the n X n identity matrix and P=P_n be the incidence matrix of the cycle (1,2,3,...,n). Then, for n>=7, a(n-7) is the number of (0,1) n X n matrices A<=P^(-1)+I+P having exactly two 1's in every row and column with perA=16. - Vladimir Shevelev, Apr 12 2010
Row 2 of the convolution array A213550. - Clark Kimberling, Jun 20 2012
a(n-1) = risefac(n, 4)/4! - risefac(n, 2)/2! is for n >= 1 also the number of independent components of a symmetric traceless tensor of rank 4 and dimension n. Here risefac is the rising factorial. - Wolfdieter Lang, Dec 10 2015
Consider the array formed by the second polygonal numbers of increasing rank:
A000217(-1-n): 0, 1, 3, 6, 10, 15, ...
A000270(-1-n): 1, 4, 9, 16, 25, 36, ...
A000326(-1-n): 2, 7, 15, 26, 40, 57, ...
A000384(-1-n): 3, 10, 21, 36, 55, 78, ...
Then the antidiagonal sums yield this sequence. - Michael Somos, Nov 23 2021

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), Table 22.7, p. 797.
  • Vladimir S. Shevelyov (Shevelev), Extension of the Moser class of four-line Latin rectangles, DAN Ukrainy, 3(1992),15-19. [From Vladimir Shevelev, Apr 12 2010]
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • A. M. Yaglom and I. M. Yaglom: Challenging Mathematical Problems with Elementary Solutions. Vol. I. Combinatorial Analysis and Probability Theory. New York: Dover Publications, Inc., 1987, p. 13, #51 (the case k=4) (First published: San Francisco: Holden-Day, Inc., 1964)

Crossrefs

Partial sums of A005581.

Programs

  • Maple
    [seq(binomial(n,4)+2*binomial(n,3), n=2..43)]; # Zerinvary Lajos, Jul 26 2006
    seq((n+4)*binomial(n,4)/n, n=3..43); # Zerinvary Lajos, Feb 28 2007
    A005582:=(-2+z)/(z-1)**5; # conjectured by Simon Plouffe in his 1992 dissertation
  • Mathematica
    Table[n(n+1)(n+2)(n+7)/24,{n,0,40}] (* Harvey P. Dale, Jun 01 2012 *)
  • PARI
    concat(0, Vec(x*(2-x)/(1-x)^5 + O(x^100))) \\ Altug Alkan, Dec 10 2015

Formula

a(n) = binomial(n+3, n-1) + binomial(n+2, n-1).
a(n) = binomial(n,4) + 2*binomial(n,3), n>=2. - Zerinvary Lajos, Jul 26 2006
From Colin Barker, Jan 28 2012: (Start)
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5).
G.f.: x*(2-x)/(1-x)^5. (End)
a(n) = Sum_{k=1..n} ( Sum_{i=1..k} i(n-k+2) ). - Wesley Ivan Hurt, Sep 26 2013
a(n+1) = A127672(8+n, n), n >= 0, with the Chebyshev C-polynomial coefficients A127672(n, k). See the Abramowitz-Stegun reference. - Wolfdieter Lang, Dec 10 2015
E.g.f.: (1/24)*x*(48 + 60*x + 16*x^2 + x^3)*exp(x). - G. C. Greubel, Jul 01 2017
Sum_{n>=1} 1/a(n) = 853/1225. - Amiram Eldar, Jan 02 2021
a(n) = A005587(-7-n) for all n in Z. - Michael Somos, Nov 23 2021

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Jun 01 2000

A005716 Coefficient of x^8 in expansion of (1+x+x^2)^n.

Original entry on oeis.org

1, 15, 90, 357, 1107, 2907, 6765, 14355, 28314, 52624, 93093, 157950, 258570, 410346, 633726, 955434, 1409895, 2040885, 2903428, 4065963, 5612805, 7646925, 10293075, 13701285, 18050760, 23554206, 30462615, 39070540, 49721892, 62816292, 78816012, 98253540
Offset: 4

Views

Author

Keywords

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 78.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    I:=[1, 15, 90, 357, 1107, 2907, 6765, 14355, 28314]; [n le 9 select I[n] else 9*Self(n-1)-36*Self(n-2)+84*Self(n-3)-126*Self(n-4)+126*Self(n-5)-84*Self(n-6)+36*Self(n-7)-9*Self(n-8)+Self(n-9): n in [1..40]]; // Vincenzo Librandi, Jun 16 2012
    
  • Magma
    /* By definition: */ P:=PolynomialRing(Integers()); [ Coefficients((1+x+x^2)^n)[9]: n in [4..32] ]; // Bruno Berselli, Jun 17 2012
  • Maple
    A005716:=-(6*z-9*z**2+3*z**3+1)/(z-1)**9; # Conjectured by Simon Plouffe in his 1992 dissertation.
    A005716 := n -> GegenbauerC(`if`(8A005716(n)), n=4..20); # Peter Luschny, May 10 2016
  • Mathematica
    CoefficientList[Series[(1+6*x-9*x^2+3*x^3)/(1-x)^9,{x,0,40}],x] (* Vincenzo Librandi, Jun 16 2012 *)

Formula

a(n) = binomial(n+1, 5)*(n^2+23*n-84)*(n+10)/336, n >= 4.
G.f.: (x^4)*(1+6*x-9*x^2+3*x^3)/(1-x)^9. (Numerator polynomial is N3(8, x) from A063420).
a(n) = A027907(n, 8), n >= 4 (ninth column of trinomial coefficients).
a(n) = A111808(n,8) for n>7. - Reinhard Zumkeller, Aug 17 2005
a(n) = 9*a(n-1) - 36*a(n-2) + 84*a(n-3) - 126*a(n-4) + 126*a(n-5) - 84*a(n-6) + 36*a(n-7) - 9*a(n-8) + a(n-9). Vincenzo Librandi, Jun 16 2012
a(n) = binomial(n,4) + 10*binomial(n,5) + 15*binomial(n,6) + 7*binomial(n,7) + binomial(n,8) (see our comment in A026729). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
a(n) = GegenbauerC(N, -n, -1/2) where N = 8 if 8Peter Luschny, May 10 2016

Extensions

More terms from Vladeta Jovovic, Oct 02 2000

A101986 Maximum sum of products of successive pairs in a permutation of order n+1.

Original entry on oeis.org

0, 2, 9, 23, 46, 80, 127, 189, 268, 366, 485, 627, 794, 988, 1211, 1465, 1752, 2074, 2433, 2831, 3270, 3752, 4279, 4853, 5476, 6150, 6877, 7659, 8498, 9396, 10355, 11377, 12464, 13618, 14841, 16135, 17502, 18944, 20463, 22061, 23740, 25502
Offset: 0

Views

Author

Eugene McDonnell (eemcd(AT)mac.com), Jan 29 2005

Keywords

Comments

1 3 5 4 2 is the 11th permutation, in lexical order. of order 5. Its reverse 2 4 5 3 1 is the 41st. The earliest permutation of order 6 is the 41st, 1 3 5 6 4 2. This pattern continues as far as I have looked, so its reversal 2 4 6 5 3 1 is the 191st and the earliest permutation of order 7 is the 191st, et cetera.
Comments from Dmitry Kamenetsky, Dec 15 2006: (Start)
This sequence is related to A026035, except here we take the maximum sum of products of successive pairs. Here is a method for generating such permutations. Start with two lists, the first has numbers 1 to n, while the second is empty.
Repeat the following operations until the first list is empty: 1. Move the smallest number of the first list to the leftmost available position in the second list. The move operation removes the original number from the first list. 2. Move the smallest number of the first list to the rightmost available position in the second list. For example when n=8, the permutation is 1, 3, 5, 7, 8, 6, 4, 2. (End)
Convolution of odd numbers and integers greater than 1. - Reinhard Zumkeller, Mar 30 2012
For n>0, a(n) is row 2 of the convolution array A213751. - Clark Kimberling, Jun 20 2012

Examples

			The permutations of order 5 with maximum sum of products is 1 3 5 4 2 and its reverse, since (1*3)+(3*5)+(5*4)+(4*2) is 46. All others are empirically less than 46. So a(4) = 46.
		

Crossrefs

Pairwise sums of A005581.

Programs

  • Haskell
    a101986 n = sum $ zipWith (*) [1,3..] (reverse [2..n+1])
    -- Reinhard Zumkeller, Mar 30 2012
  • J
    0 1 9 2 & p. % 6 & p. (A) NB. the polynomial P such that P(n) is a(n).
    NB. where 0 1 9 2 are the coefficients in ascending order of the numerator of a rational polynomial and 6 is the (constant) coefficient of its denominator. J's primitive function p. produces a polynomial with these coefficients. Division is indicated by % . Thus the J expression (A) is equivalent to the formula above.
    
  • Maple
    a:=n->add((n+j^2),j=1..n): seq(a(n),n=0..41); # Zerinvary Lajos, Jul 27 2006
  • Mathematica
    Table[(n + 9 n^2 + 2 n^3)/6, {n, 0, 41}] (* Robert G. Wilson v, Feb 04 2005 *)
  • PARI
    a(n)=n*(2*n^2+9*n+1)/6 \\ Charles R Greathouse IV, Jan 17 2012
    

Formula

a(n) = n*(2*n^2 + 9*n + 1)/6.
a(n+1) = a(n) + A008865(n+2); a(n) = A160805(n) - 4. [Reinhard Zumkeller, May 26 2009]
G.f.: x*(1+x)*(2-x)/(1-x)^4. - L. Edson Jeffery, Jan 17 2012
a(n) = 4*a(n-1)-6*a(n-2)+4*a(n-3)-a(n-4) for n>3, a(0)=0, a(1)=2, a(2)=9, a(3)=23. - L. Edson Jeffery, Jan 17 2012
a(n) = A000330(n) + A005449(n) - A000217(n). - Richard R. Forberg, Aug 07 2013
a(n) = 1 + sum( A008865(i), i=1..n+1 ). [Bruno Berselli, Jan 13 2015]
a(n) = A000290(n) + A000330(n). - J. M. Bergot, Apr 26 2018

Extensions

Edited by Bruno Berselli, Jan 13 2015
Name edited by Alois P. Heinz, Feb 02 2019

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

A005587 a(n) = n*(n+5)*(n+6)*(n+7)/24.

Original entry on oeis.org

0, 14, 42, 90, 165, 275, 429, 637, 910, 1260, 1700, 2244, 2907, 3705, 4655, 5775, 7084, 8602, 10350, 12350, 14625, 17199, 20097, 23345, 26970, 31000, 35464, 40392, 45815, 51765, 58275, 65379, 73112, 81510, 90610, 100450, 111069, 122507, 134805, 148005
Offset: 0

Views

Author

Keywords

Comments

a(n) = number of Standard Young Tableaux of shape (n+3,4). - David Callan, Aug 17 2004
a(n) = A214292(n+6,3). - Reinhard Zumkeller, Jul 12 2012
a(n) for n > 0 is the number of n-extended coalescent histories for a matching caterpillar gene tree and species tree with 5 leaves. - Noah A Rosenberg, Jun 16 2022

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

Crossrefs

Fifth diagonal of Catalan triangle A033184. Fifth column of Catalan triangle A009766.
Numerator polynomial 14 - 28x + 20x^2 - 5x^3 from fourth row of triangle A062991.

Programs

  • Magma
    [n*(n+5)*(n+6)*(n+7)/24: n in [0..40]]; // Vincenzo Librandi, Mar 20 2013
    
  • Maple
    A005587:=z*(-14+28*z-20*z**2+5*z**3)/(z-1)**5; # Simon Plouffe in his 1992 dissertation
    seq(numbperm(n,4)/24-numbperm(n,3)/6, n=7..46); # Zerinvary Lajos, May 20 2008
    a:=n->(sum(numbcomp(n,4), j=9..n)):seq(a(n)/4, n=8..47); # Zerinvary Lajos, Aug 26 2008
  • Mathematica
    Table[n (n + 5) (n + 6) (n + 7)/24, {n, 0, 60}] (* Vladimir Joseph Stephan Orlovsky, Jun 22 2011 *)
    LinearRecurrence[{5,-10,10,-5,1},{0,14,42,90,165},40] (* Harvey P. Dale, Aug 17 2017 *)
  • PARI
    x='x+O('x^50); concat([0], Vec((14 - 28*x + 20*x^2 - 5*x^3) / (1 - x)^5)) \\ G. C. Greubel, Jul 01 2017

Formula

G.f.: (14 - 28*x + 20*x^2 - 5*x^3) / (1 - x)^5.
a(n) = C(7+n, 4) - C(7+n, 3). - Zerinvary Lajos, Dec 09 2005
E.g.f.: (1/24)*x*(336 + 168*x + 24*x^2 + x^3)*exp(x). - G. C. Greubel, Jul 01 2017
From Amiram Eldar, Jun 28 2022: (Start)
Sum_{n>=1} 1/a(n) = 153/1225.
Sum_{n>=1} (-1)^(n+1)/a(n) = 288*log(2)/35 - 20759/3675. (End)
a(n) = A024191(n+1)-5. - R. J. Mathar, Nov 22 2024

Extensions

M4929 (this sequence) and M4930 were the same.
More terms from Matthew Conroy, Jan 16 2006
Plouffe Maple line edited by N. J. A. Sloane, May 13 2008

A051936 Truncated triangular numbers: a(n) = n*(n+1)/2 - 9.

Original entry on oeis.org

1, 6, 12, 19, 27, 36, 46, 57, 69, 82, 96, 111, 127, 144, 162, 181, 201, 222, 244, 267, 291, 316, 342, 369, 397, 426, 456, 487, 519, 552, 586, 621, 657, 694, 732, 771, 811, 852, 894, 937, 981, 1026, 1072, 1119, 1167, 1216, 1266, 1317, 1369, 1422, 1476
Offset: 4

Views

Author

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

Keywords

Comments

Equals binomial transform of [1, 5, 1, 0, 0, 0, ...]. - Gary W. Adamson, Apr 30 2008
Numbers m > 0 such that 8m+73 is a square. - Bruce J. Nicholson, Jul 29 2017

Examples

			Illustration of the initial terms:
                                                          .
                              .                         .   .
      .                     .   .                     o   o   o
    .   .                 o   o   o                 o   o   o   o
  .   o   .             .   o   o   .             .   o   o   o   .
.   .   .   .         .   .   o   .   .         .   .   o   o   .   .
----------------------------------------------------------------------
      1                       6                           12
----------------------------------------------------------------------
- _Bruno Berselli_, Oct 13 2016
		

Crossrefs

Cf. A000217.

Programs

Formula

G.f.: x^4*(-1-3*x+3*x^2) / (x-1)^3.
a(n) = n + a(n-1) for n>4, a(4)=1. - Vincenzo Librandi, Aug 06 2010
a(n) = 2*A000217(n-3) - A000217(n-6), with A000217(-2)=1, A000217(-1)=0. - Bruno Berselli, Oct 13 2016
Sum_{n>=4} 1/a(n) = 53/72 + 2*Pi*tan(sqrt(73)*Pi/2)/sqrt(73). - Amiram Eldar, Dec 13 2022

A342981 Triangle read by rows: T(n,k) is the number of rooted planar maps with n edges, k faces and no isthmuses, n >= 0, k = 1..n+1.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 1, 7, 5, 0, 1, 16, 37, 14, 0, 1, 30, 150, 176, 42, 0, 1, 50, 449, 1104, 794, 132, 0, 1, 77, 1113, 4795, 7077, 3473, 429, 0, 1, 112, 2422, 16456, 41850, 41504, 14893, 1430, 0, 1, 156, 4788, 47832, 189183, 319320, 228810, 63004, 4862
Offset: 0

Views

Author

Andrew Howroyd, Apr 02 2021

Keywords

Comments

The number of vertices is n + 2 - k.
For k >= 2, column k is a polynomial of degree 3*(k-2). This is because adding a face can increase the number of vertices whose degree is greater than two by at most two.
By duality, also the number of loopless rooted planar maps with n edges and k vertices.

Examples

			Triangle begins:
  1;
  0, 1;
  0, 1,   2;
  0, 1,   7,    5;
  0, 1,  16,   37,    14;
  0, 1,  30,  150,   176,    42;
  0, 1,  50,  449,  1104,   794,   132;
  0, 1,  77, 1113,  4795,  7077,  3473,   429;
  0, 1, 112, 2422, 16456, 41850, 41504, 14893, 1430;
  ...
		

Crossrefs

Columns k=3..4 are A005581, A006468.
Diagonals are A000108, A006419, A006420, A006421.
Row sums are A000260.

Programs

  • Mathematica
    G[m_, y_] := Sum[x^n*Sum[(n + k - 1)!*(2*n - k)!*y^k/(k!*(n + 1 - k)!*(2*k - 1)!*(2*n - 2*k + 1)!), {k, 1, n}], {n, 1, m}] + O[x]^m;
    H[n_] := With[{g = 1 + x*y + x*G[n - 1, y]}, Sqrt[InverseSeries[x/g^2 + O[x]^(n + 1), x]/x]];
    CoefficientList[#, y]& /@ CoefficientList[H[10], x] // Flatten (* Jean-François Alcover, Apr 15 2021, after Andrew Howroyd *)
  • PARI
    \\ here G(n, y) gives A082680 as g.f.
    G(n,y)={sum(n=1, n, x^n*sum(k=1, n, (n+k-1)!*(2*n-k)!*y^k/(k!*(n+1-k)!*(2*k-1)!*(2*n-2*k+1)!))) + O(x*x^n)}
    H(n)={my(g=1+x*y+x*G(n-1, y), v=Vec(sqrt(serreverse(x/g^2)/x))); vector(#v, n, Vecrev(v[n], n))}
    { my(T=H(8)); for(n=1, #T, print(T[n])) }

Formula

G.f. A(x,y) satisfies A(x) = G(x*A(x,y)^2, y) where G(x,y) = 1 + x*y + x*B(x,y) and B(x,y) is the g.f. of A082680.
A027836(n+1) = Sum_{k=1..n+1} k*T(n,k).
A002293(n) = Sum_{k=1..n+1} k*T(n,n+2-k).

A005583 Coefficients of Chebyshev polynomials.

Original entry on oeis.org

2, 11, 36, 91, 196, 378, 672, 1122, 1782, 2717, 4004, 5733, 8008, 10948, 14688, 19380, 25194, 32319, 40964, 51359, 63756, 78430, 95680, 115830, 139230, 166257, 197316, 232841, 273296, 319176, 371008, 429352, 494802, 567987, 649572, 740259, 840788
Offset: 1

Views

Author

Keywords

Comments

If X is an n-set and Y a fixed 2-subset of X then a(n-5) is equal to the number of (n-5)-subsets of X intersecting Y. - Milan Janjic, Jul 30 2007
a(n-1) = risefac(n,5)/5! - risefac(n,3)/3! is for n >= 1 also the number of independent components of a symmetric traceless tensor of rank 5 and dimension n. Here risefac is the rising factorial. Put a(0) = 0. - Wolfdieter Lang, Dec 10 2015

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), Table 22.7, p. 797.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Maple
    A005583:=-(-2+z)/(z-1)**6; # Simon Plouffe in his 1992 dissertation (this g.f. assumes offset 0)
  • PARI
    conv(u,v)=local(w); w=vector(length(u),i,sum(j=1,i,u[j]*v[i+1-j])); w;
    t(n)=n*(n+1)/2;
    u=vector(10,i,t(i));
    v=vector(10,i,t(i)-1);
    conv(u,v)
    
  • PARI
    a(n) = (1/120)*n*(n+9)*(n+3)*(n+2)*(n+1); \\ Joerg Arndt, Mar 05 2018

Formula

G.f.: x*(2-x)/(1-x)^6.
a(n) = binomial(n+4, n-1) + binomial(n+3, n-1) = (1/120)*n*(n+9)*(n+3)*(n+2)*(n+1).
a(n+1) = -A127672(10+n, n), n >= 0, with the coefficients of the Chebyshev C-polynomials A127672(n, k). - Wolfdieter Lang, Dec 10 2015
a(n) = Sum_{i=1..n} A000217(i)*A000096(n+1-i). - Bruno Berselli, Mar 05 2018
a(n) = binomial(n+3,5) + 2*binomial(n+3,4). - Yuchun Ji, May 23 2019
From Amiram Eldar, Feb 17 2023: (Start)
Sum_{n>=1} 1/a(n) = 40751/63504.
Sum_{n>=1} (-1)^(n+1)/a(n) = 1360*log(2)/63 - 922961/63504. (End)

Extensions

More terms from Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Dec 07 1999
More terms from Zerinvary Lajos, Jul 21 2006

A005714 Coefficient of x^6 in expansion of (1+x+x^2)^n.

Original entry on oeis.org

1, 10, 45, 141, 357, 784, 1554, 2850, 4917, 8074, 12727, 19383, 28665, 41328, 58276, 80580, 109497, 146490, 193249, 251713, 324093, 412896, 520950, 651430, 807885, 994266, 1214955, 1474795, 1779121, 2133792, 2545224, 3020424, 3567025, 4193322, 4908309, 5721717
Offset: 3

Views

Author

Keywords

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 78.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    I:=[1, 10, 45, 141, 357, 784, 1554]; [n le 7 select I[n] else 7*Self(n-1)-21*Self(n-2)+35*Self(n-3)-35*Self(n-4)+21*Self(n-5)-7*Self(n-6)+Self(n-7): n in [1..40]]; // Vincenzo Librandi, Jun 16 2012
    
  • Magma
    /* By definition: */ P:=PolynomialRing(Integers()); [ Coefficients((1+x+x^2)^n)[7]: n in [3..35] ]; // Bruno Berselli, Jun 17 2012
  • Maple
    A005714:=-(1+3*z-4*z**2+z**3)/(z-1)**7; # Conjectured by Simon Plouffe in his 1992 dissertation.
    A005714 := n -> GegenbauerC(`if`(6A005714(n)), n=3..20); # Peter Luschny, May 10 2016
  • Mathematica
    a[n_] := Coefficient[(1 + x + x^2)^n, x, 6]; Table[a[n], {n, 3, 35}]
    CoefficientList[Series[(1+3*x-4*x^2+x^3)/(1-x)^7,{x,0,40}],x] (* Vincenzo Librandi, Jun 16 2012 *)

Formula

a(n) = binomial(n, 3)*(n^3+18*n^2+17*n-120) /120.
G.f.: (x^3)*(1+3*x-4*x^2+x^3)/(1-x)^7. (Numerator polynomial is N3(6, x) from A063420).
a(n) = A027907(n, 6), n >= 3 (seventh column of trinomial coefficients).
a(n) = A111808(n,6) for n>5. - Reinhard Zumkeller, Aug 17 2005
a(n) = 7*a(n-1) -21*a(n-2) +35*a(n-3) -35*a(n-4) +21*a(n-5) -7*a(n-6) +a(n-7). Vincenzo Librandi, Jun 16 2012
a(n) = binomial(n,3) + 6*binomial(n,4) + 5*binomial(n,5) + binomial(n,6) (see our comment in A026729). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
a(n) = GegenbauerC(N, -n, -1/2) where N = 6 if 6Peter Luschny, May 10 2016
E.g.f.: exp(x)*x^3*(120 + 180*x + 30*x^2 + x^3)/720. - Stefano Spezia, Mar 28 2023

Extensions

More terms from Vladeta Jovovic, Oct 02 2000
Previous Showing 11-20 of 49 results. Next