Напечатать все пифагоровы числа не привышающие 30 (алгоритм)

Перебираешь все взаимно простые m > n разной четности, так чтобы m*m было меньше 30. И формируешь тройки по формуле:

k * (m*m-n*n, 2*m*n, m*m + n*n)

Например так:

function gcd (a, b: integer) : integer;

begin

if b = 0 then gcd : = a

else if (a > b) and (b > 0) then gcd : = gcd (b, a mod b)

else gcd : = gcd (b, a) ;

end;

procedure Print (m, n:integer) ;

var i, a, b, c: integer;

begin

repeat

i : = i + 1;

a : = (m*m-n*n) * i;

b : = 2*m*n*i;

c : = (m*m+n*n) * i;

if c < 30 then

writeln (a, ' ', b, ' ', c)

else break;

until false;

end;

var m, n, s: integer;

begin

s : = Trunc (Sqrt (30)) ;

for m : = 1 to s do

for n : = m + 1 to s do

if (gcd (n, m) = 1) and ((m mod 2) (n mod 2)) then

Print (n, m)

end.