OFFSET
0,3
LINKS
Felix Huber, Table of n, a(n) for n=0..666
Felix Huber, Illustration of the term a(4) = 9
EXAMPLE
See the linked illustration for the term a(4) = 9.
MAPLE
A372915:=proc(n)
local p, q, g, h, u, v, x, y, L, M;
L:=[];
for g from 2 to n do
h:=2*n/g;
if type(h, integer) then
for x to n do
M:=[g, sqrt(x^2+h^2), sqrt((g-x)^2+h^2)];
M:=sort(M);
if not member(M, L) then
L:=[op(L), M];
fi;
od;
fi;
od;
for p to n do
for q from 1 to p do
g:=sqrt(p^2+q^2);
h:=2*n/g;
u:=h/g*q;
v:=q+h/g*p;
for x from max(1, ceil(p/q*(v-n)+u)) to min(n, floor(p/q*v+u)) do
y:=q/p*(u-x)+v;
if type(y, integer) and x <> p and y <> q then
M:=[g, sqrt(x^2+(y-q)^2), sqrt((x-p)^2+y^2)];
M:=sort(M);
if not member(M, L) then
L:=[op(L), M];
fi;
fi;
od;
od;
od;
return numelems(L);
end proc;
seq(A372915(n), n=0..53);
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Huber, Jun 02 2024
STATUS
approved