OFFSET
1,1
COMMENTS
Positive integers x in the solutions to 2*x^2-146*y^2-10512*y-254040 = 0.
LINKS
Colin Barker, Table of n, a(n) for n = 1..600
Index entries for linear recurrences with constant coefficients, signature (0,0,0,4562498,0,0,0,-1).
FORMULA
a(n) = 4562498*a(n-4)-a(n-8).
G.f.: -584*x*(x-1)*(7*x^6+47*x^5+2757*x^4+18250*x^3+2757*x^2+47*x+7) / ((x^4-2136*x^2-1)*(x^4+2136*x^2-1)).
EXAMPLE
4088 is in the sequence because 4088^2 = 16711744 = 442^2+443^2+...+514^2.
MATHEMATICA
LinearRecurrence[{0, 0, 0, 4562498, 0, 0, 0, -1}, {4088, 23360, 1582640, 9047912, 18642443912, 106578370640, 7220791811360, 41281080400088}, 40] (* Vincenzo Librandi, May 11 2015 *)
PROG
(PARI) Vec(-584*x*(x-1)*(7*x^6+47*x^5+2757*x^4+18250*x^3+2757*x^2+47*x+7) / ((x^4-2136*x^2-1)*(x^4+2136*x^2-1)) + O(x^100))
(Magma) I:=[4088, 23360, 1582640, 9047912, 18642443912, 106578370640, 7220791811360, 41281080400088]; [n le 8 select I[n] else 4562498*Self(n-4)-Self(n-8): n in [1..20]]; // Vincenzo Librandi, May 11 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Colin Barker, May 10 2015
STATUS
approved