login
A228208
y-values in the solution to x^2 - 20*y^2 = 176.
6
1, 2, 5, 7, 14, 19, 37, 50, 97, 131, 254, 343, 665, 898, 1741, 2351, 4558, 6155, 11933, 16114, 31241, 42187, 81790, 110447, 214129, 289154, 560597, 757015, 1467662, 1981891, 3842389, 5188658, 10059505, 13584083, 26336126, 35563591, 68948873, 93106690
OFFSET
1,2
COMMENTS
Also y-values in the solution of x^2-5*y^2=44.
FORMULA
G.f.: x*(x+1)*(x^2+x+1) / ((x^2-x-1)*(x^2+x-1)).
a(n) = 3*a(n-2)-a(n-4).
Let h(n) = hypergeom([(1 - n)/2, (n + 1) mod 2 - n/2], [1 - n], -4) then a(n) = h(n-1) + h(n) for n > 2. - Peter Luschny, Sep 03 2019
MATHEMATICA
CoefficientList[Series[(x + 1) (x^2 + x + 1) / ((x^2 - x - 1) (x^2 + x - 1)), {x, 0, 40}], x] (* Vincenzo Librandi, Aug 17 2013 *)
PROG
(PARI) Vec(x*(x+1)*(x^2+x+1)/((x^2-x-1)*(x^2+x-1)) + O(x^100))
(Magma) I:=[1, 2, 5, 7, 14]; [n le 4 select I[n] else 3*Self(n-2)-Self(n-4): n in [1..40]]; // Vincenzo Librandi, Aug 17 2013
CROSSREFS
Cf. A228207, A054486 (bisection).
Sequence in context: A022771 A132603 A022765 * A338231 A286573 A133511
KEYWORD
nonn,easy
AUTHOR
Colin Barker, Aug 16 2013
STATUS
approved