|
|
A219753
|
|
Expansion of x^4*(1-2*x+x^4)/((1+x)*(1-2*x)^2*(1-x-x^2)).
|
|
2
|
|
|
0, 0, 0, 0, 1, 2, 6, 13, 31, 68, 151, 327, 706, 1509, 3211, 6796, 14331, 30115, 63106, 131905, 275111, 572676, 1190031, 2469039, 5115426, 10584557, 21874931, 45158908, 93131731, 191884523, 395000226, 812445097, 1669750831, 3429176564, 7037627911, 14433683991
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,6
|
|
LINKS
|
|
|
FORMULA
|
G.f.: x^4*(1-2*x+x^4)/((1+x)*(1-2*x)^2*(1-x-x^2)) (half of A219752).
a(n) = (2^(n-5)*(3*n+16)+4*(-1)^n)/9+((1-sqrt(5))^(n-3)-(1+sqrt(5))^(n-3))/(2^(n-3)*sqrt(5)) with n>3, a(0)=a(1)=a(2)=a(3)=0. [Bruno Berselli, Nov 29 2012]
|
|
MAPLE
|
if n < 4 then
0 ;
else
128*(-1)^n+2^n*(16+3*n)-288*A000045(n-3) ;
%/288 ;
end if ;
end proc:
|
|
MATHEMATICA
|
CoefficientList[Series[x^4 (1 - 2 x + x^4)/((1 + x) (1 - 2 x)^2 (1 - x - x^2)), {x, 0, 35}], x] (* Bruno Berselli, Nov 30 2012 *)
LinearRecurrence[{4, -2, -7, 4, 4}, {0, 0, 0, 0, 1, 2, 6, 13, 31}, 40] (* Harvey P. Dale, Oct 05 2021 *)
|
|
PROG
|
(Maxima) makelist(coeff(taylor(x^4*(1-2*x+x^4)/((1+x)*(1-2*x)^2*(1-x-x^2)), x, 0, n), x, n), n, 0, 35); /* Bruno Berselli, Nov 29 2012 */
(Magma) I:=[0, 0, 0, 0, 1, 2, 6, 13, 31]; [n le 9 select I[n] else 4*Self(n-1) - 2*Self(n-2) - 7*Self(n-3) + 4*Self(n-4) + 4*Self(n-5): n in [1..40]]; // Vincenzo Librandi, Dec 14 2012
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|