Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A005715
Coefficient of x^7 in expansion of (1+x+x^2)^n.
(Formerly M3632)
6
4, 30, 126, 393, 1016, 2304, 4740, 9042, 16236, 27742, 45474, 71955, 110448, 165104, 241128, 344964, 484500, 669294, 910822, 1222749, 1621224, 2125200, 2756780, 3541590, 4509180, 5693454, 7133130, 8872231, 10960608, 13454496
OFFSET
4,1
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).
LINKS
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Eric Weisstein's World of Mathematics, Trinomial Coefficient
FORMULA
a(n) = binomial(n, 4)*(n^3+27*n^2+116*n-120)/210, n >= 4.
G.f.: (x^4)*(x-2)*(x^2-2)/(1-x)^8. (Numerator polynomial is N3(7, x) from A063420).
a(n) = A027907(n, 7), n >= 4 (eighth column of trinomial coefficients).
a(n) = A111808(n,7) for n>6. - Reinhard Zumkeller, Aug 17 2005
a(n) = 8*a(n-1) -28*a(n-2) +56*a(n-3) -70*a(n-4) +56*a(n-5) -28*a(n-6) +8*a(n-7) -a(n-8). Vincenzo Librandi, Jun 16 2012
a(n) = 4*binomial(n,4) + 10*binomial(n,5) + 6*binomial(n,6) + binomial(n,7) (see our comment in A026729). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
a(n) = GegenbauerC(N, -n, -1/2) where N = 7 if 7<n else 2*n-7. - Peter Luschny, May 10 2016
MAPLE
A005715:=(z-2)*(z**2-2)/(z-1)**8; # Conjectured by Simon Plouffe in his 1992 dissertation.
A005715 := n -> GegenbauerC(`if`(7<n, 7, 2*n-7), -n, -1/2):
seq(simplify(A005715(n)), n=4..20); # Peter Luschny, May 10 2016
MATHEMATICA
CoefficientList[Series[(x-2)*(x^2-2)/(1-x)^8, {x, 0, 40}], x] (* Vincenzo Librandi, Jun 16 2012 *)
PROG
(Magma) I:=[4, 30, 126, 393, 1016, 2304, 4740, 9042]; [n le 8 select I[n] else 8*Self(n-1)-28*Self(n-2)+56*Self(n-3)-70*Self(n-4)+56*Self(n-5)-28*Self(n-6)+8*Self(n-7)-Self(n-8): n in [1..40]]; // Vincenzo Librandi, Jun 16 2012
(Magma) /* By definition: */ P<x>:=PolynomialRing(Integers()); [ Coefficients((1+x+x^2)^n)[8]: n in [4..33] ]; // Bruno Berselli, Jun 17 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
More terms from Vladeta Jovovic, Oct 02 2000
STATUS
approved