OFFSET
1,2
COMMENTS
The corresponding sequence for which 6n^2-6 is a square is A001079.
Positive values of x (or y) satisfying x^2 - 10xy + y^2 + 20 = 0. - Colin Barker, Feb 09 2014
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..200
Eric Weisstein's World of Mathematics, Pell Equation.
Index entries for linear recurrences with constant coefficients, signature (0,10,0,-1).
FORMULA
a(n) = 10*a(n-2)-a(n-4).
G.f.: x*(1-x)*(1+4*x+x^2)/(1-10*x^2+x^4). - Colin Barker, Jun 13 2012
a(2*n+1) = ((6+r)*(5+2*r)^n+(6-r)*(5+2*r)^n)/12, a(2*n+2) = ((18+7*r)*(5+2*r)^n+(18-7*r)*(5-2*r)^n)/12, where r=sqrt(6) and n>=0. - Paul Weisenhorn, Sep 01 2012
EXAMPLE
29 is a term of the sequence since 6*29^2 - 5 = 5041 = 71^2.
MATHEMATICA
Do[ If[ IntegerQ[ Sqrt[6n^2 - 5]], Print[n]], {n, 1, 3*10^7}]
a[1]=1; a[2]=3; a[3]=7; a[4]=29; a[n_] := a[n]=10a[n-2]-a[n-4]
CoefficientList[Series[(1 - x) (1 + 4 x + x^2)/(1 - 10 x^2 + x^4), {x, 0, 30}], x] (* Vincenzo Librandi, Feb 10 2014 *)
PROG
(Magma) I:=[1, 3, 7, 29]; [n le 4 select I[n] else 10*Self(n-2)-Self(n-4): n in [1..30]]; // Vincenzo Librandi, Feb 10 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
John W. Layman, Mar 24 2003
EXTENSIONS
Extended by Robert G. Wilson v, Apr 14 2003
STATUS
approved