1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:

procedure TForm1.BerechnenClick(Sender: TObject);
var a,b,c:integer;
begin
  for a:=1 to 2500 do
  for b:=a to 2500 do
  for c:=b to 2500 do
  if a*a+b*b=c*c then
  if not((a mod 2=0) and (b mod 2=0) and (c mod 2=0)) then
  if not((a mod 3=0) and (b mod 3=0) and (c mod 3=0)) then
  if not((a mod 4=0) and (b mod 4=0) and (c mod 4=0)) then
  if not((a mod 5=0) and (b mod 5=0) and (c mod 5=0)) then
  if not((a mod 6=0) and (b mod 6=0) and (c mod 6=0)) then
  if not((a mod 7=0) and (b mod 7=0) and (c mod 7=0)) then
  if not((a mod 8=0) and (b mod 8=0) and (c mod 8=0)) then
  if not((a mod 9=0) and (b mod 9=0) and (c mod 9=0)) then
  if not((a mod 10=0) and (b mod 10=0) and (c mod 10=0)) then
  listbox1.items.add(inttostr(a)+' '+inttostr(b)+' '+inttostr(c));
end;					

1:

 if (ggT(a,b) = 1) and (ggT(a,c) = 1) and (ggt(b,c) = 1) then					

1:
2:
3:
4:
5:
6:
7:
8:

1. Erstelle Feld (Matrix), evtl. Shortint um Speicher zu sparen, der Dimension [1.. MaxZahl, 1..MaxZahl] und belege alle Werte mit "-1" (=Status des Zahlenpaars ist unbekannt)
2. für i=1 .. Maxzahl und j=i+1 .. Maxzahl:
 Wenn Feld[i,j]=-1  (also unbekannter Status des Paares (i,j) dann:
  Wenn sqrt (i*i+j*j) ganze Zahl dann:
   Feld [i,j] := 1 // Zahlenpaar (i,j) ist korrekte Lösung
   und weiterhin(wichtig!): Alle Feld  [i*x,j*x] mit i*x<=Maxzahl und j*x<Maxzahl auf 0 setzen. Bedeutet: Diese Zahlenpaare sind nur Vielfache der gefundenen Lösung, also für die weitere Berechnung überspringen

3. Ausgabe der Ergebnisse (gebe alle (i,j,i*i+j*j) aus, für die gilt: Feld[i,j]=1)					

1:

 if (ggT(a,b) = 1) and (ggT(a,c) = 1) and (ggt(b,c) = 1) then					

1:

 if (ggt(c,ggT(a,b)) = 1) then					

1:

if (ggt(a,b)=1) and (ggt(a,c)=1) and (ggt(b,c)=1) then					

1:

if ggt(ggt(a,b),c)=1 then					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:

procedure TForm1.BerechnenClick(Sender: TObject);
var a,b,c,cquadrat:integer;
function ggt(zahl1,zahl2:integer):integer;
var rest:integer;
begin
  repeat
  rest:=zahl1 mod zahl2;
  zahl1:=zahl2;
  zahl2:=rest;
  until rest=0;
  ggt:=zahl1;
end;
begin
  for a:=1 to 2500 do
  for b:=a to 2500 do
  begin
  cquadrat:=a*a + b*b;
  c:=trunc(sqrt(cquadrat));
  if (c<=2500) and (c * c = cquadrat) then
  begin
  if (ggt(a,b)=1) and (ggt(a,c)=1) and (ggt(b,c)=1) then
  listbox1.items.add(inttostr(a)+' '+inttostr(b)+' '+inttostr(c));
  end;
  end;
end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:

procedure TForm1.BerechnenClick(Sender: TObject);
var 
  a,b,c,cquadrat:integer;
  
  function ggt(zahl1,zahl2:integer):integer;
  var 
    rest:integer;
  begin
    repeat
      rest:=zahl1 mod zahl2;
      zahl1:=zahl2;
      zahl2:=rest;
    until rest=0;
    Result:=zahl1;
  end;
  
begin
  ListBox1.Items.BeginUpdate;
  try
    ListBox1.Items.Clear;
    for a:=1 to 2500 do
      for b:=a to 2500 do
      begin
        cquadrat:=sqr(a) + sqr(b);
        c:=trunc(sqrt(cquadrat));
        if (c<=2500) and (sqr(c) = cquadrat) then
        begin
          if (ggt(a,b)=1) and (ggt(a,c)=1) and (ggt(b,c)=1) then
            listbox1.items.add(Format('%d %d %d', [a, b, c]));
        end;
      end;
  finally
    ListBox1.Items.EndUpdate;
  end;
end;