OFFSET
1,1
COMMENTS
n^2 = 3*(2*x+4)^2+16.
Essentially the same as A011943. - Chris Boyd, Apr 18 2015
LINKS
FORMULA
a(n+2) = 14*a(n+1)-a(n) for n>=2.
G.f.: x*(2-21*x+x^2)/(1-14*x+x^2). a(n)=7*A007655(n+2)-97*A007655(n+1), n>1. - R. J. Mathar, Sep 11 2008
EXAMPLE
The first examples are 2^3-0^3=2*2^2 ; 5^3-3^3=2*7^2 ; 57^3-55^3=2*97^2 ; 781^3-779^3=2*1351^2 ; 10865^3-10863^3=2*18817^2
MAPLE
2, seq(othopoly[T](n, 7), n=1..50); # Robert Israel, Apr 19 2015
MATHEMATICA
a[1]=2; a[2]=7; a[3]=97; a[n_] := a[n] = 14*a[n-1]-a[n-2]; Table[a[n], {n, 1, 17}] (* Jean-François Alcover, Dec 17 2013 *)
LinearRecurrence[{14, -1}, {2, 7, 97}, 20] (* Harvey P. Dale, Sep 26 2016 *)
PROG
(Magma) I:=[2, 7, 97]; [n le 3 select I[n] else 14*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Apr 19 2015
(PARI) Vec(x*(2-21*x+x^2)/(1-14*x+x^2) + O(x^30)) \\ Michel Marcus, Apr 19 2015
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Richard Choulet, Sep 08 2008
EXTENSIONS
More terms from Vincenzo Librandi, Apr 19 2015
STATUS
approved