q = x^3; s = x^2 + 1; z = 40;
p[n_, x_] := Fibonacci[n, x];
Table[Expand[p[n, x]], {n, 1, 7}]
reduce[{p1_, q_, s_, x_}] :=
FixedPoint[(s PolynomialQuotient @@ #1 +
PolynomialRemainder @@ #1 &)[{#1, q, x}] &, p1]
t = Table[reduce[{p[n, x], q, s, x}], {n, 1, z}];
u1 = Table[Coefficient[Part[t, n], x, 0], {n, 1, z}]
u2 = Table[Coefficient[Part[t, n], x, 1], {n, 1, z}]
u3 = Table[Coefficient[Part[t, n], x, 2], {n, 1, z}]
LinearRecurrence[{1, 3, -1, -3, 1, 1}, {1, 0, 1, 1, 2, 5}, 40] (*
Harvey P. Dale, Nov 07 2021 *)