Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #26 Sep 08 2022 08:45:58
%S 1,2,2,7,16,44,113,298,778,2039,5336,13972,36577,95762,250706,656359,
%T 1718368,4498748,11777873,30834874,80726746,211345367,553309352,
%U 1448582692,3792438721,9928733474,25993761698,68052551623,178163893168,466439127884,1221153490481
%N Constant term in the reduction by (x^2->x+1) of the polynomial p(n,x) defined below at Comments.
%C The titular polynomial is defined by p(n,x) = x*p(n-1,x) +(x^2)*p(n-2,x), with p(0,x)=1, p(1,x)=1+x^2. For discussions of polynomial reduction, see A192232, A192744, and A192872.
%H Alois P. Heinz, <a href="/A192921/b192921.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (2,2,-1).
%F a(n) = 2*a(n-1) + 2*a(n-2) - a(n-3).
%F G.f.: (1-2*x)*(1+2*x) / ( (1+x)*(1-3*x+x^2) ). - _R. J. Mathar_, May 08 2014
%F a(n) = A059929(n-1) + 2*A059929(n-2). - _R. J. Mathar_, May 08 2014
%F a(n) = F(n-4)*F(n) + F(n-1)*F(n+2), where F(-4)=-3, F(-3)=2, F(-2)=-1, F(-1)=1. - _Bruno Berselli_, Nov 03 2015
%F a(n) = (2^(-n)*(-3*(-2)^n-(-4+sqrt(5))*(3+sqrt(5))^n+(3-sqrt(5))^n*(4+sqrt(5))))/5. - _Colin Barker_, Oct 01 2016
%e The coefficients in the polynomials p(n,x) are Fibonacci numbers. The first seven and their reductions:
%e ...
%e 1 -> 1
%e 1 + x^2 -> 2 + x
%e x + x^2 + x^3 -> 2 + 4*x
%e 2*x^2 + x^3 + 2*x^4 -> 7 + 10*x
%e 3*x^3 + 2*x^4 + 3*x^5 -> 16 + 27*x
%e 5*x^4 + 3*x^5 + 5*x^6 -> 44 + 70*x
%e 8*x^5 + 5*x^6 + 8*x^7 -> 113 + 184*x,
%e so that A192921=(1,2,2,7,16,44,113,...).
%p a:= n-> (<<0|1|0>, <0|0|1>, <-1|2|2>>^n. <<1,2,2>>)[1,1]:
%p seq(a(n), n=0..30); # _Alois P. Heinz_, Sep 28 2016
%t q = x^2; s = x + 1; z = 28;
%t p[0, x_] := 1; p[1, x_] := x^2 + 1;
%t p[n_, x_] := p[n - 1, x]*x + p[n - 2, x]*x^2;
%t Table[Expand[p[n, x]], {n, 0, 7}]
%t reduce[{p1_, q_, s_, x_}] :=
%t FixedPoint[(s PolynomialQuotient @@ #1 +
%t PolynomialRemainder @@ #1 &)[{#1, q, x}] &, p1]
%t t = Table[reduce[{p[n, x], q, s, x}], {n, 0, z}];
%t u1 = Table[Coefficient[Part[t, n], x, 0], {n, 1, z}]
%t (* A192921 *)
%t u2 = Table[Coefficient[Part[t, n], x, 1], {n, 1, z}]
%t (* A192879 *)
%t LinearRecurrence[{2,2,-1}, {1,2,2}, 30] (* _G. C. Greubel_, Feb 06 2019 *)
%o (PARI) a(n) = round((2^(-n)*(-3*(-2)^n-(-4+sqrt(5))*(3+sqrt(5))^n+(3-sqrt(5))^n*(4+sqrt(5))))/5) \\ _Colin Barker_, Oct 01 2016
%o (PARI) Vec(-(2*x-1)*(1+2*x)/((1+x)*(x^2-3*x+1)) + O(x^40)) \\ _Colin Barker_, Oct 01 2016
%o (PARI) {a(n) = fibonacci(n-2)^2 +fibonacci(n)*fibonacci(n+1)}; \\ _G. C. Greubel_, Feb 06 2019
%o (Magma) [Fibonacci(n-2)^2 + Fibonacci(n)*Fibonacci(n+1): n in [0..30]]; // _G. C. Greubel_, Feb 06 2019
%o (Sage) [fibonacci(n-2)^2 +fibonacci(n)*fibonacci(n+1) for n in range(30)] # _G. C. Greubel_, Feb 06 2019
%o (GAP) List([0..30], n -> Fibonacci(n-2)^2 +Fibonacci(n)*Fibonacci(n+1)); # _G. C. Greubel_, Feb 06 2019
%Y Cf. A192232, A192744, A192872.
%K nonn,easy
%O 0,2
%A _Clark Kimberling_, Jul 12 2011