OFFSET
0,6
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (0,1,2).
FORMULA
a(n) = Sum_{m=1..n} Sum_{i=0..n-m} binomial(m+i-1,m-1)*Sum_{j=0..m} binomial(j,n-3*m+2*j-i)*2^(m-j)*binomial(m,j)*(-1)^(-n+3*m-j+i)). - Vladimir Kruchinin, May 12 2011
MAPLE
seq(coeff(series((1-x)/(-2*x^3-x^2+1), x, n+1), x, n), n = 0..50); # G. C. Greubel, Aug 04 2019
MATHEMATICA
CoefficientList[Series[(1-x)/(1-x^2-2*x^3), {x, 0, 50}], x] (* Wesley Ivan Hurt, Jan 24 2017 *)
PROG
(Maxima)
a(n):=sum(sum(binomial(m+i-1, m-1)*sum(binomial(j, n-3*m+2*j-i)*2^(m-j)*binomial(m, j)*(-1)^(-n+3*m-j+i), j, 0, m), i, 0, n-m), m, 1, n); /* Vladimir Kruchinin, May 12 2011 */
(PARI) Vec((1-x)/(1-x^2-2*x^3) + O(x^50)) \\ Felix Fröhlich, Jan 24 2017
(Magma) R<x>:=PowerSeriesRing(Integers(), 50); Coefficients(R!( (1-x)/(1-x^2-2*x^3) )); // G. C. Greubel, Aug 04 2019
(Sage) ((1-x)/(1-x^2-2*x^3)).series(x, 50).coefficients(x, sparse=False) # G. C. Greubel, Aug 04 2019
(GAP) a:=[1, -1, 1];; for n in [4..50] do a[n]:=a[n-2]+2*a[n-3]; od; a; # G. C. Greubel, Aug 04 2019
CROSSREFS
KEYWORD
sign,easy
AUTHOR
N. J. A. Sloane, Nov 17 2002
STATUS
approved