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:
36:
37:
38:
39:
40:
41:
42:
43:
44:
45:
46:
47:
48:
49:
50:
51:
52:
53:
54:
55:
56:
57:
58:
59:
60:
61:
62:
63:
64:
65:
66:
67:
68:
69:
70:
71:
72:
73:
74:
75:
76:
77:
78:
79:
80:
81:
82:
83:
84:
85:
86:
87:
88:
89:
90:

function IsPrime(N: Cardinal): Boolean; register; 
// test if N is prime, do some small Strong Pseudo Prime test in certain bounds 
// not PIC safe !! 
// very fast, about 1200 clockcycles 
// copyright by Hagen Reddmann, don't remove this line 
asm 
       TEST  EAX,1            // Odd(N) ?? 
       JNZ   @@1 
       CMP   EAX,2            // N == 2 ?? 
       SETE  AL 
       RET 
@@1:   CMP   EAX,73           // N  < 73 ?? 
       JB    @@C 
       JE    @@E              // N == 73 ?? 
       PUSH  ESI 
       PUSH  EDI 
       PUSH  EBX 
       PUSH  EBP 
       PUSH  EAX              // save N as Param for @@5 
       LEA   EBP,[EAX - 1]    // M == N -1, Exponent 
       MOV   ECX,32           // calc remaining Bits of M and shift M' 
       MOV   ESI,EBP 
@@2:   DEC   ECX 
       ADD   ESI,ESI 
       JNC   @@2 
       PUSH  ECX              // save Bits as Param for @@5 
       PUSH  ESI              // save M' as Param for @@5 
       CMP   EAX,08A8D7Fh     // N >= 9080191 ?? 
       JAE   @@3 
// now if (N < 9080191) and SPP(31, N) and SPP(73, N) then N is prime 
       MOV   EAX,31 
       CALL  @@5              // 31^((N-1)(2^s)) mod N 
       JC    @@4 
       MOV   EAX,73           // 73^((N-1)(2^s)) mod N 
       PUSH  OFFSET @@4 
       JMP   @@5 
// now if (N < 4759123141) and SPP(2, N) and SPP(7, N) and SPP(61, N) then N is prime 
@@3:   MOV   EAX,2 
       CALL  @@5 
       JC    @@4 
       MOV   EAX,7 
       CALL  @@5 
       JC    @@4 
       MOV   EAX,61 
       CALL  @@5 
@@4:   SETNC AL 
       ADD   ESP,4 * 3 
       POP   EBP 
       POP   EBX 
       POP   EDI 
       POP   ESI 
       RET 
// do a Strong Pseudo Prime Test 
@@5:   MOV   EBX,[ESP + 12]   // N on stack 
       MOV   ECX,[ESP +  8]   // remaining Bits 
       MOV   ESI,[ESP +  4]   // M' 
       MOV   EDI,EAX          // T = b, temp. Base 
@@6:   DEC   ECX 
       MUL   EAX 
       DIV   EBX 
       MOV   EAX,EDX 
       ADD   ESI,ESI 
       JNC   @@7 
       MUL   EDI 
       DIV   EBX 
       AND   ESI,ESI 
       MOV   EAX,EDX 
@@7:   JNZ   @@6 
       CMP   EAX,1            // b^((N -1)(2^s)) mod N ==  1 mod N ?? 
       JE    @@A 
@@8:   CMP   EAX,EBP          // b^((N -1)(2^s)) mod N == -1 mod N ?? , EBP = N -1 
       JE    @@A 
       DEC   ECX              // second part to 2^s 
       JNG   @@9 
       MUL   EAX 
       DIV   EBX 
       CMP   EDX,1 
       MOV   EAX,EDX 
       JNE   @@8 
@@9:   STC 
@@A:   RET 
@@B:   DB    3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71 
@@C:   MOV   EDX,OFFSET @@B 
       MOV   ECX,18 
@@D:   CMP   AL,[EDX + ECX] 
       JE    @@E 
       DEC   ECX 
       JNL   @@D 
@@E:   SETE  AL 
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:
36:
37:
38:
39:
40:
41:
42:
43:
44:
45:
46:
47:
48:
49:
50:
51:
52:
53:
54:
55:
56:
57:
58:
59:
60:
61:
62:
63:
64:
65:
66:
67:
68:
69:
70:
71:
72:

procedure TMainForm.StartButClick(Sender: TObject);
var
  Start, Ende: Real;
  Z: PCardinal;
  X, Y, LPrim: Cardinal;  // Hier
  PrimF: TStringList;
  CountP: String;
begin
  try
    Von := StrToInt(VonEdit.Text);
    Bis := StrToInt(BisEdit.Text);

    if Bis > 10 then
      SetLength(PrimS, Trunc(0.4*Bis-(Bis/4)))
    else
      SetLength(PrimS, 4);

    LPrim := 0;
    Z := @PrimS[0];
    Y := 0;  // Hier

    if (Von >= 0) AND (Bis >= 0) AND (Von < Bis) then
    begin
      Start := GetTickCount;

      for X := Von to Bis do
      begin
        if Prim(X) then
        begin
          Z^ := X;
          Inc(Z);
          Inc(Y);  // Hier
          LPrim := X;
        end;
        if X mod 20000 = 0 then
        begin
          Application.ProcessMessages;
          PrimLab.Caption := 'Aktuelle Primzahl: ' + IntToStr(LPrim);    // Akt. Primzahl anzeigen
        end;
      end;

      Ende := GetTickCount;
      DauerLab.Caption := 'Diese Überprüfung hat ' + FloatToStr((Ende - Start) / 1000) + ' Sekunden gedauert.';

      CountP := IntToStr(Y);
      PrimLab.Caption := 'Speichern... (0 / ' + CountP + ')';

      PrimF := TStringList.Create;

      for X := 0 to Length(PrimS)-1 do
      begin
        if PrimS[X] = 0 then
          Break;
        if X mod 1000 = 0 then   // Hier
        begin
          Application.ProcessMessages;
          PrimLab.Caption := 'Speichern... (' + IntToStr(X+1) + ' / ' + CountP + ')';
        end;
        PrimF.Add(IntToStr(PrimS[X]));
      end;

      PrimF.SaveToFile('Paket.prim');
      PrimF.Free;

      PrimLab.Caption := 'Aktuelle Primzahl: ' + IntToStr(LPrim);
    end
    else
      ShowMessage('Ungültige Eingabe(n)!');
  except
    ShowMessage('Es ist ein Fehler aufgetreten!');
  end;
end;					

1:
2:
3:
4:
5:
6:

    //MaxTest := Trunc(Sqrt(CurrNum*1.0));
    Asm
        FILD    QWORD PTR [CurrNum]
        FSQRT
        FISTP   QWORD PTR [MaxTest]
    End;					

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

    Asm
        FLDCW   WORD PTR [@@TruncCW]
        JMP     @@Skip
        @@TruncCW:
        DW      $1772   //Immer abwärts runden
        @@Skip:
    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:

procedure SavePrimes(MaxPrime: Cardinal; const FileName: String = '');    
var
  isPrime: Array of Boolean;
  i, j: Cardinal;
  SL: TStringList;
begin
  SetLength(isPrime, MaxPrime+1);
  FillChar(isPrime[2], Length(isPrime)-2, 1);
  for i:=2 to Trunc(Sqrt(MaxPrime)) do
    if isPrime[i] then begin
      j:=i*2;
      while j<=MaxPrime do begin
        isPrime[j]:=False;
        Inc(j, i);
      end;
    end;

  if FileName<>'' then begin
    SL:=TStringList.Create;
    try
      for i:=2 To MaxPrime do
        if isPrime[i] then
          SL.Add(IntToStr(i));
      SL.SaveToFile(FileName);
    finally
      FreeAndNil(SL);
    end;
  end;
end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:

    FillChar(PrimeGrid, SizeOf(PrimeGrid), #1);
    PrimeGrid[1] := False;
    for CurrPrimeCandidate := 2 to MAX_PRIME do
    begin
      if PrimeGrid[CurrPrimeCandidate] then
      begin
        NotAPrime := CurrPrimeCandidate * 2;
        while NotAPrime <= MAX_PRIME do
        begin
          PrimeGrid[NotAPrime] := False;
          Inc(NotAPrime, CurrPrimeCandidate);
        end;
      end;
    end;					

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

Function IsPrim(zahl : Integer): boolean;
var
i: integer;
begin
  result := true;
  If zahl = 1 then
  begin
    result := false;
    exit;
  end;
  For i := 2 to (zahl div 2) do
  begin
    If ((zahl mod i) = 0) then
    begin
      result := false;
      exit;
    end;
  end;
end;					

1:

For i := 2 to (zahl div 2) do					

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:

procedure SavePrimes(MinPrime, MaxPrime: Cardinal; const FileName: String = '');
var
  isPrime: Array of Cardinal; 
  i, j: Cardinal;
  SL: TStringList;
begin
  SetLength(isPrime, MaxPrime div 32);
  FillChar(isPrime[0], Length(isPrime)*SizeOf(Cardinal), MAXBYTE); // alle bits auf 1 setzten
  isPrime[0]:=$FFFFFFFC; // bit 0 und 1 ist keine primzahl

  for i:=2 to Trunc(Sqrt(MaxPrime)) do
    if isPrime[i div 32] shr (i mod 32) and 1=1 then begin
      j:=i*2;
      while j<=MaxPrime do begin
        isPrime[j div 32]:=isPrime[j div 32] and not (1 shl (j mod 32));
        Inc(j, i);
      end;
    end;

  if FileName<>'' then begin
    SL:=TStringList.Create;
    try
      for i:=MinPrime To MaxPrime do
        if isPrime[i div 32] shr (i mod 32) and 1=1 then
          SL.Add(IntToStr(i));
      SL.SaveToFile(FileName);
    finally
      FreeAndNil(SL);
    end;
  end;
end;					

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

function tform1.Prim(Zahl: LongInt): Boolean;
var
  Teiler, Wurzel: LongInt;
  p1prim:PTprim;
begin
  Result := True;
  p1prim:=startprim;
  if (not odd(Zahl)) OR (Zahl <= 5) then
  begin
    if (Zahl <> 2) AND (Zahl <> 3) AND (Zahl <> 5) then
      Result := False;
    Exit;
  end;
  Wurzel := Trunc(sqrt(Zahl));
  while (p1prim^.prim <= Wurzel) AND (Result) do
  begin
    if Zahl mod p1prim^.prim = 0 then
      Result := False;
    p1prim:=p1prim^.nextprim;
  end;
end;					

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

 if prim(zutesten) then
 begin
  new(P1prim);
  p1prim^.nextprim:=nil;
  lastprim^.nextprim:=p1prim;
  p1prim^.prim:=zutesten;
  lastprim:=p1prim;
  inc(zutesten)
 end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:

procedure TForm1.FormCreate(Sender: TObject);
var
a,b,c:integer;
p1prim:PTprim;
begin
 new(startprim);
 new(P1prim);
 startprim^.prim:=2;
 startprim^.nextprim:=P1prim;
 p1prim^.prim:=3;
 p1prim^.nextprim:=nil;
 lastprim:=p1prim;
end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:

procedure TForm1.FormClose(Sender: TObject; var Action: TCloseAction);
var
p1prim,p2prim:PTprim;
begin
p1prim:=startprim;
while p1prim^.nextprim<>nil do
begin
p2prim:=p1prim;
p1prim:=p1prim^.nextprim;
dispose(p2prim);
end;
dispose(p1prim);
end;					

1:
2:
3:
4:
5:

 
  PTPrim=^TPrim;
  TPrim =record
   prim:integer;
   nextprim:PTPrim;