OFFSET
1,1
LINKS
Robert Israel, Pythagorean triples and their sums of prime factors corresponding to the first 41 terms of the sequence.
Eric Weisstein's World of Mathematics, Pythagorean Triple.
EXAMPLE
a(1)=5 is in the sequence because (3,4,5) is a Pythagorean triple, whose sums of prime factors with repetition are again (3,4,5).
a(2)=289 is in the sequence because (161, 240, 289) is a Pythagorean triple whose sums of prime factors with repetition are (30, 16, 34), which again is a Pythagorean triple.
MAPLE
N:= 10^7: # for terms <= N
pyth:= [seq(seq([m^2-n^2, 2*m*n, m^2+n^2], n=1..min(m-1, floor(sqrt(N-m^2)))), m=1..floor(sqrt(N)))]:
filter:= t -> spf(t[3])^2-spf(t[1])^2-spf(t[2])^2=0:
R:= select(filter, pyth):
sort(map(t -> t[3], R));
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Nov 22 2020
STATUS
approved