login
A098890
a(n) = A001652(n)*A046090(n)+1 = A098602(n)+1.
2
1, 13, 421, 14281, 485113, 16479541, 559819261, 19017375313, 646030941361, 21946034630941, 745519146510613, 25325704946729881, 860328449042305321, 29225841562491651013, 992818284675673829101, 33726595837410418538401, 1145711440187278556476513
OFFSET
0,2
FORMULA
For n>0, a(n) = A001541(n)^2 + A001542(n)^2 - A098602(n-1).
G.f.: -(x^2-22*x+1) / ((x-1)*(x^2-34*x+1)). - Colin Barker, Jul 09 2015
a(n) = (5/8+1/16*(17+12*sqrt(2))^(-n)*(3-2*sqrt(2)+(3+2*sqrt(2))*(17+12*sqrt(2))^(2*n))). - Colin Barker, Mar 02 2016
EXAMPLE
a(3) = 14281 = 99^2+70^2-420.
PROG
(PARI) Vec(-(x^2-22*x+1)/((x-1)*(x^2-34*x+1)) + O(x^30)) \\ Colin Barker, Jul 09 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Charlie Marion, Nov 03 2004
EXTENSIONS
Name and a(0) corrected by Charlie Marion, Jul 09 2015
STATUS
approved