cp's OEIS Frontend

A122099 a(n) = -3*a(n-1) + a(n-3) for n>2, with a(0)=1, a(1)=1, a(2)=0.

%I A122099 #25 Sep 08 2022 08:45:28
%S A122099 1,1,0,1,-2,6,-17,49,-141,406,-1169,3366,-9692,27907,-80355,231373,
%T A122099 -666212,1918281,-5523470,15904198,-45794313,131859469,-379674209,
%U A122099 1093228314,-3147825473,9063802210,-26098178316,75146709475,-216376326215,623030800329,-1793945691512,5165460748321
%N A122099 a(n) = -3*a(n-1) + a(n-3) for n>2, with a(0)=1, a(1)=1, a(2)=0.
%H A122099 G. C. Greubel, <a href="/A122099/b122099.txt">Table of n, a(n) for n = 0..1000</a>
%H A122099 <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (-3,0,1).
%F A122099 G.f.: (1+4*x+3*x^2)/(1+3*x-x^3).
%F A122099 a(n) = (-1)^n*A122100(n). - _R. J. Mathar_, Sep 27 2014
%p A122099 seq(coeff(series((1+4*x+3*x^2)/(1+3*x-x^3), x, n+1), x, n), n = 0..40); # _G. C. Greubel_, Oct 02 2019
%t A122099 Transpose[NestList[{#[[2]],Last[#],First[#]-3Last[#]}&, {1,1,0},35]][[1]]  (* _Harvey P. Dale_, Mar 13 2011 *)
%t A122099 LinearRecurrence[{-3,0,1}, {1,1,0}, 40] (* _G. C. Greubel_, Oct 02 2019 *)
%o A122099 (PARI) Vec((1+4*x+3*x^2)/(1+3*x-x^3)+O(x^40)) \\ _Charles R Greathouse IV_, Jan 17 2012
%o A122099 (Magma) R<x>:=PowerSeriesRing(Integers(), 40); Coefficients(R!( (1+4*x+3*x^2)/(1+3*x-x^3) )); // _G. C. Greubel_, Oct 02 2019
%o A122099 (Sage)
%o A122099 def A122099_list(prec):
%o A122099     P.<x> = PowerSeriesRing(ZZ, prec)
%o A122099     return P((1+4*x+3*x^2)/(1+3*x-x^3)).list()
%o A122099 A122099_list(40) # _G. C. Greubel_, Oct 02 2019
%o A122099 (GAP) a:=[1,1,0];; for n in [4..40] do a[n]:=-3*a[n-1]+a[n-3]; od; a; # _G. C. Greubel_, Oct 02 2019
%Y A122099 Cf. A122100.
%K A122099 sign,easy,less
%O A122099 0,5
%A A122099 _Philippe Deléham_, Oct 18 2006