login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A077996 Expansion of (1-x)/(1-2*x-x^2-2*x^3). 2
1, 1, 3, 9, 23, 61, 163, 433, 1151, 3061, 8139, 21641, 57543, 153005, 406835, 1081761, 2876367, 7648165, 20336219, 54073337, 143779223, 382304221, 1016534339, 2702931345, 7187005471, 19110010965, 50812890091, 135109802089, 359252516199, 955240614669 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

LINKS

G. C. Greubel, Table of n, a(n) for n = 0..1000

Index entries for linear recurrences with constant coefficients, signature (2,1,2).

FORMULA

If p[1]=1, p[2]=2, p[i]=4, (i>2), and if A is Hessenberg matrix of order n defined by: A[i,j]=p[j-i+1], (i<=j), A[i,j]=-1, (i=j+1), and A[i,j]=0 otherwise. Then, for n>=1, a(n)=det A. - Milan Janjic, May 02 2010

a(n) = Sum_{k=1..n} Sum_{i=k..n} Sum_{j=0..k} binomial(j,-3*k+2*j+i) * 2^(k-j)*binomial(k,j)*binomial(n+k-i-1,k-1). - Vladimir Kruchinin, May 05 2011

MATHEMATICA

LinearRecurrence[{2, 1, 2}, {1, 1, 3}, 40] (* or *) CoefficientList[Series[(1 -x)/(1-2*x-x^2-2*x^3), {x, 0, 40}], x] (* G. C. Greubel, Jun 27 2019 *)

PROG

(Maxima)

a(n):=sum(sum((sum(binomial(j, -3*k+2*j+i)*2^(k-j)*binomial(k, j), j, 0, k) )*binomial(n+k-i-1, k-1), i, k, n), k, 1, n); /* Vladimir Kruchinin, May 05 2011 */

(PARI) Vec((1-x)/(1-2*x-x^2-2*x^3)+O(x^40)) \\ Charles R Greathouse IV, Sep 27 2012

(MAGMA) R<x>:=PowerSeriesRing(Integers(), 30); Coefficients(R!( (1-x)/(1-2*x-x^2-2*x^3) )); // G. C. Greubel, Jun 27 2019

(Sage) ((1-x)/(1-2*x-x^2-2*x^3)).series(x, 30).coefficients(x, sparse=False) # G. C. Greubel, Jun 27 2019

(GAP) a:=[1, 1, 3];; for n in [4..30] do a[n]:=2*a[n-1]+a[n-2]+2*a[n-3]; od; a; # G. C. Greubel, Jun 27 2019

CROSSREFS

Sequence in context: A026599 A061647 A207008 * A029852 A282732 A047085

Adjacent sequences:  A077993 A077994 A077995 * A077997 A077998 A077999

KEYWORD

nonn,easy

AUTHOR

N. J. A. Sloane, Nov 17 2002

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 17 04:03 EDT 2019. Contains 327119 sequences. (Running on oeis4.)