OFFSET
1,1
COMMENTS
If m is in sequence, so is any multiple of m. Primitive elements (terms which are not divisible by any previous term) are A354381.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
MAPLE
ishyp:= proc(n) local s; ormap(s -> s mod 4 = 1, numtheory:-factorset(n)) end proc:
filter:= proc(n) local s;
ormap(s -> ishyp(subs(s, x)) and ishyp(subs(s, y)), [isolve(x^2+y^2=n^2)])
end proc:
select(filter, [$1..1000]); # Robert Israel, Jan 10 2023
MATHEMATICA
ishyp[n_] := AnyTrue[FactorInteger[n][[All, 1]], Mod[#, 4] == 1&];
filter[n_] := AnyTrue[Solve[x^2 + y^2 == n^2, Integers], ishyp[x /. #] && ishyp[y /. #]&];
Select[Range[400], filter] (* Jean-François Alcover, May 11 2023, after Robert Israel *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Lamine Ngom, May 24 2022
STATUS
approved