OFFSET
0,2
COMMENTS
a(0) = 1 is included by convention. This corresponds to the Pythagorean triple 1^2 + 0^2 = 1^2.
REFERENCES
Miguel Ángel Pérez García-Ortega, José Manuel Sánchez Muñoz and José Miguel Blanco Casado, El Libro de las Ternas Pitagóricas, Preprint 2024.
LINKS
Index entries for linear recurrences with constant coefficients, signature (41,-246,246,-41,1).
FORMULA
EXAMPLE
MATHEMATICA
s[n_]:=s[n]=Module[{a, b}, a=((1+Sqrt[2])^(2n+1)-(Sqrt[2]-1)^(2n+1))/2; b=(a^2-1)/2; {(a+2b+1)/2}]; semis={}; Do[semis=Join[semis, FullSimplify[s[n]]], {n, 0, 17}]; semis
PROG
(PARI) Vec((1 - 13*x - 41*x^2 + 21*x^3)/((1 - 34*x + x^2)*(1 - 6*x + x^2)*(1 - x)) + O(x^20)) \\ Andrew Howroyd, Oct 14 2024
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Miguel-Ángel Pérez García-Ortega, Oct 13 2024
STATUS
approved