OFFSET
1,2
COMMENTS
Note that 84847 = 7 * 17 * 23 * 31, the first four primes in A058529.
LINKS
Vincenzo Librandi and T. D. Noe, Table of n, a(n) for n = 1..1000 (first 159 terms from Vincenzo Librandi)
FORMULA
a(n) = a(n-1) + 6*a(n-81) - 6*a(n-82) - a(n-162) + a(n-163), where the 163 initial terms can be computed using the Mathematica program. The initial terms begin with 0 and end with 1696940.
MATHEMATICA
d = 84847; t = Select[Range[0, 68000], IntegerQ[Sqrt[#^2 + (#+d)^2]] &]
CROSSREFS
KEYWORD
nonn
AUTHOR
T. D. Noe, Feb 09 2012
STATUS
approved