OFFSET
1,2
COMMENTS
a(n) <= 10^(n/2)/sqrt(2). - Robert G. Wilson v, Jan 08 2005
LINKS
Ravi P. Agarwal, Pythagorean Triples Before and after Pythagoras, Computation 2020, 8(3), 62, p. 14.
Eric Weisstein's World of Mathematics, Twin Pythagorean Triple
EXAMPLE
For a(2)=7 because {3, 4, 5}, {5, 12, 13}, {7, 24, 25}, {20, 21, 29}, {9, 40, 41}, {11, 60, 61}, {13, 84, 85}.
MATHEMATICA
f[n_] := Round[ Sqrt[ (10^n - 1)/2] - 1]; a[n_] := (MatrixPower[{{1, 2, 2}, {2, 1, 2}, {2, 2, 3}}, n].{{1}, {0}, {1}})[[3, 1]]; t = Table[ a[n], {n, 2, 35}]; g[n_] := Length[ Select[t, # < 10^n &]]; Table[f[n] + g[n], {n, 25}] (* Robert G. Wilson v, Jan 08 2005 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Eric W. Weisstein, Dec 20 2004
EXTENSIONS
a(5) and up from Robert G. Wilson v, Jan 08 2005
STATUS
approved