OFFSET
0,5
COMMENTS
Conjecture: this is the inverse Motzkin transform of A054393.
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..3435
Index entries for linear recurrences with constant coefficients, signature (2, -1, 2, -1, 1).
FORMULA
a(n)=1, n < 4. a(4)=2, a(5)=4. a(n) = 2*a(n-1) - a(n-2) + 2a(n-3) - a(n-4) + a(n-5), n > 5.
G.f.: 1 + x*(-1 + x + x^3)/(-1 + 2*x - x^2 + 2*x^3 - x^4 + x^5). - R. J. Mathar, May 26 2016
MAPLE
a:= n-> `if`(n=0, 1, (<<0|1|0|0|0>, <0|0|1|0|0>, <0|0|0|1|0>,
<0|0|0|0|1>, <1|-1|2|-1|2>>^n. <<0, 1, 1, 1, 2>>)[1$2]):
seq(a(n), n=0..40); # Alois P. Heinz, Nov 08 2016
MATHEMATICA
Join[{1}, LinearRecurrence[{2, -1, 2, -1, 1}, {1, 1, 1, 2, 4}, 50]] (* G. C. Greubel, Sep 01 2018 *)
PROG
(PARI) a(n) = {if (n==0, return(1)); if (n==1, return(1)); if (n==2, return(1)); if (n==3, return(1)); if (n==4, return(2)); if (n==5, return(4)); return (2*a(n-1)-a(n-2)+2*a(n-3)-a(n-4)+a(n-5)); } \\ Michel Marcus, Jul 23 2013
(PARI) x='x+O('x^50); Vec(1 + x*(-1+x+x^3)/(-1+2*x-x^2+2*x^3-x^4+x^5)) \\ G. C. Greubel, Sep 01 2018
(Magma) m:=50; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!(1 + x*(-1+x+x^3)/(-1+2*x-x^2+2*x^3-x^4+x^5))); // G. C. Greubel, Sep 01 2018
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
R. J. Mathar, Dec 11 2008
STATUS
approved