cp's OEIS Frontend

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);