OFFSET
1,2
COMMENTS
Complement of A010814.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 100: # for terms <= N
S:= {$1..N}:
for x from 1 while 2*x^2 < N do
for y from 1 to x-1 do
p:= 2*x*(x+y);
if p > N then break fi;
if igcd(x, y) = 1 and (x::odd or y::odd) then
S:= S minus {seq(i*p, i=1..N/p)}
fi
od od:
sort(convert(S, list)); # Robert Israel, Mar 29 2026
MATHEMATICA
seq[max_] := Module[{s = {}, a, b, c}, Do[If[IntegerQ[c = Sqrt[a^2 + b^2]], AppendTo[s, a + b + c]], {a, max}, {b, Floor@Sqrt[a], a}]; Complement[Range[max], s]]; seq[100] (* Amiram Eldar, Oct 24 2021 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Oct 23 2021
STATUS
approved
