aa = {}; uu = 682 + 61 * Sqrt[125]; Do[vv = Expand[uu^(2 * n - 1)]; tt = ((-1)^n vv[[1]] + 57)/125; xx = (5^5 * tt^2 - 3000 * tt + 719); yy = Round[N[Sqrt[xx^3], 1000]]; dd = xx^3 - yy^2; AppendTo[aa, dd], {n, 1, 10}]; aa
dd = {-297, 548147655, -1019827620252441}; a0 = dd[[1]]; a1 = dd[[2]]; a2 = dd[[3]]; Do[a = a0 + 1860497 * a1 - 1860497 * a2; a0 = a1; a1 = a2; a2 = a; AppendTo[aa, a], {n, 1, 10}]; aa
(* Third one after Lucas numbers formula *)
Table[27/125 (-5 + (-1)^n ((-1)^(n + 1) 6 + LucasL[15 (-1 + 2 n)])), {n, 10}] (*
Artur Jasinski, Nov 18 2011*)