OFFSET
0,3
COMMENTS
REFERENCES
T. Bier, Classifications of solutions of certain positive biquadratic division system, submitted May 12 2014.
T. Bier and O. Dira, Construction of integer sequences, submitted May 12 2014.
LINKS
Oboifeng Dira, Sequences solving division systems
Index entries for linear recurrences with constant coefficients, signature (4,0,-4,1).
FORMULA
a(n+1) = 4*a(n) - a(n-1) - p_n (n>0), where p_n = 0 if n is odd and p_n = 1 if n is even.
G.f.: (1 - 3*x - x^2 + 2*x^3)/((1 - x)*(1 + x)*(1 - 4*x + x^2)). [Bruno Berselli, May 22 2014]
a(n) = 4*a(n-1) - 4*a(n-3) + a(n-4) for n>3. [Bruno Berselli, May 22 2014]
a(n) = (3+(-1)^n+(4+sqrt(3))*(2-sqrt(3))^n+(4-sqrt(3))*(2+sqrt(3))^n)/12. [Bruno Berselli, May 25 2014]
EXAMPLE
a(2)=3, a(3)=10 satisfy: 10 divides 3^2+1 and 3 divides 10^2+10+1.
MAPLE
x0:=1: x1:=1: L:=[x0, x1]: for k from 1 to 30 do: if k mod 2 = 1 then z:=4*x1-x0: fi: if k mod 2 = 0 then z:=4*x1-x0-1: fi: L:=[op(L), z]: x0:=x1: x1:=z: od: print(L);
MATHEMATICA
Table[(3 + (-1)^n + (4 + Sqrt[3]) (2 - Sqrt[3])^n + (4 - Sqrt[3]) (2 + Sqrt[3])^n)/12, {n, 0, 30}] (* Bruno Berselli, May 25 2014 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Oboifeng Dira, May 21 2014
EXTENSIONS
Edited by Bruno Berselli, May 25 2014
STATUS
approved