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:

...
  for k := 0 to MaxPrime div PrimeBits do
  begin
    FillChar(Primes[0], PrimeLen, 0);
    for i := 0 to Trunc(Sqrt((k + 1) * PrimeBits) / 30) do
    begin
      for j := 0 to 7 do
      begin
        if PrimesLUT[i] AND (1 shl j) = 0 then
        begin
          s   := STEMPEL[j];
          Num := i * 30 + s;
          if k = 0 then
            Num2 := Num * Num
          else
            Num2 := Trunc(k * PrimeBits / Num) * Num + Num;

          mbit := Num2 mod 30;
          m    := (Num2 - mbit) div 30 - k * PrimeLen;
          while m < PrimeLen do
          begin
            primes[m] := Primes[m] or MODS[mbit];
            m    := m + i;
            mbit := mbit + s;
            if mbit > 29 then
            begin
              mbit := mbit - 30;
              inc(m);
            end;
          end;
        end;
      end;
    end;
  end;

  if not (FileName = '') then
  begin
    for k := 0 to MaxPrime div PrimeBits do
    begin
      for i := 0 to PrimeLen - 1 do
      begin
        for j := 0 to 7 do
        begin
          if k * PrimeBits + i * 30 + STEMPEL[j] > MaxPrime then
            Break;
          if not ((i = 0) AND (j = 0) AND (k = 0)) AND
                 (Primes[i] AND (1 shl j) = 0) then
           SList.Add(IntToStr(k * PrimeBits + i * 30 + STEMPEL[j]));
        end;
      end;
    end;

    SList.SaveToFile(FileName);
    SList.Free;
  end;
end;					

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

 
mbit = Num2 mod 30
mbit = Num² mod 30
mbit = (30i + s)² mod 30
mbit = (900i² + 60is + s²) mod 30 [(a + km) mod m = a mod m]
mbit = s² % 30					

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

m = (Num2 - mbit) div 30
m = (Num2 - s² mod 30) div 30
m = ( Num² - s² mod 30) div 30
m = ( (30i + s)² - s² mod 30) div 30
m = (900i² + 60is + s² - s² mod 30) div 30
m = 30i² + 2is + (s² - s² mod 30) div 30					

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:

procedure SavePrimes(MaxPrime: Cardinal; const FileName: string = '');
const
  CACHE = 64*1024;
  STEMPEL: array[0..7] of Byte = (1, 7, 11, 13, 17, 19, 23, 29);
  MODS: array[0..29] of Byte = (0, 1, 0, 0, 0, 0, 0, 2, 0, 0, 0, 4, 0, 8, 0, 0,
                                0, 16, 0, 32, 0, 0, 0, 64, 0, 0, 0, 0, 0, 128);
  mBitLUT: array[0..7] of Byte = (1, 19, 1, 19, 19, 1, 19, 1);
  mLUT: array[0..7] of Byte = (0, 1, 4, 5, 9, 12, 17, 28);
var
  Primes: array[0..CACHE-1] of Byte;
  PrimesLUT: array of Byte;
  i, j, k, PrimeLen, PrimeBits, Num, Num2, m, mbit, s: Cardinal;
  f: TextFile;
begin
  if FileName<>'' then begin
    AssignFile(f, FileName);
    ReWrite(f);
    WriteLn(f, '2'+#13#10+'3'+#13#10+'5');
  end;

  SetLength(PrimesLUT, Trunc(Sqrt(MaxPrime)/30)); // max 2184 Byte für 2^32 ;-)
  PrimesLUT[0]:=$01;
  PrimeLen:=Length(PrimesLUT);
  PrimeBits:=PrimeLen*30;
  for i:=0 to Trunc(Sqrt(PrimeBits)/30) do
    for j:=0 to 7 do
      if PrimesLUT[i] and (1 shl j)=0 then begin
        s:=STEMPEL[j];
        mbit:=mBitLUT[j];
        m:=30*(i*i) + 2*i*s + mLUT[j];
        while m<PrimeLen do begin
          PrimesLUT[m]:=PrimesLUT[m] or MODS[mbit];
          Inc(m, i);
          Inc(mbit, s);
          if mbit>29 then begin
            Dec(mbit, 30);
            Inc(m);
          end;
        end;
      end;

  PrimeLen:=Length(Primes);
  PrimeBits:=PrimeLen*30;
  for k:=0 to MaxPrime div PrimeBits do begin
    FillChar(Primes[0], PrimeLen, 0);
    for i:=0 to Trunc(Sqrt((k+1)*PrimeBits)/30) do
      for j:=0 to 7 do
        if PrimesLUT[i] and (1 shl j)=0 then begin
          s:=STEMPEL[j];
          if k=0 then begin
            mbit:=mbitLUT[j];
            m:=30*(i*i) + 2*i*s + mLUT[j];
          end else begin
            Num:=i*30+s;
            Num2:=Trunc(k*PrimeBits/Num)*Num+Num;
            mbit:=Num2 mod 30;
            m:=(Num2-mbit) div 30-k*PrimeLen;
          end;
          while m<PrimeLen do begin
            primes[m]:=Primes[m] or MODS[mbit];
            Inc(m, i);
            Inc(mbit, s);
            if mbit>29 then begin
              Dec(mbit, 30);
              Inc(m);
            end;
          end;
        end;
    if FileName<>'' then
      for i:=0 To PrimeLen-1 do
        for j:=0 to 7 do begin
          if k*PrimeBits+i*30+STEMPEL[j]>MaxPrime then
            Break;
          if not ((i=0) and (j=0) and (k=0)) and (Primes[i] and (1 shl j)=0) then
            WriteLn(f, IntToStr(k*PrimeBits+i*30+STEMPEL[j]));
        end;
  end;
  if FileName<>'' then
    CloseFile(f);
end;					

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

for j := 0 to 7 do
begin
  if k * PrimeBits + i * 30 + STEMPEL[j] > MaxPrime then
    Break;
  if not ((i = 0) AND (j = 0) AND (k = 0)) AND
         (Primes[i] AND (1 shl j) = 0) then
    SList.Add(IntToStr(k * PrimeBits + i * 30 + STEMPEL[j]));
end;					

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

for j := 1 to 7 do
begin
  if k * PrimeBits + i * 30 + STEMPEL[j] > MaxPrime then
    Break;
  if not (Primes[i] AND (1 shl j) = 0) then
    SList.Add(IntToStr(k * PrimeBits + i * 30 + STEMPEL[j]));
end;					

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

    if FileName<>'' then
      for i:=0 To PrimeLen-1 do
        for j:=0 to 7 do begin
          if k*PrimeBits+i*30+STEMPEL[j]>MaxPrime then
            Break;
          if not ((i=0) and (j=0) and (k=0)) and (Primes[i] and (1 shl j)=0) then
            WriteLn(f, IntToStr(k*PrimeBits+i*30+STEMPEL[j]));
        end;
  end;					

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

    if FileName<>'' then
      for i:=0 To PrimeLen-1 do
        for j:=0 to 7 do begin
          if k*PrimeBits+i*30+STEMPEL[j]>MaxPrime then
            Break;
          if (i or j or k or (Primes[i] and (1 shl j)) <> 0 then
            WriteLn(f, IntToStr(k*PrimeBits+i*30+STEMPEL[j]));
        end;
  end;					

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

var
    PrimPtr: PInteger;
    PrimSize: Integer;

{...}

PrimPtr := VirtualAlloc(65536*1024);

{...}

          if (i or j or k or (Primes[i] and (1 shl j)) <> 0 then
          Begin
            PrimPtr^ := k*PrimeBits+i*30+STEMPEL[j];
            Inc(PrimPtr);
            //Check for Buffer Overrun and Reallocate Memory if current size is not enought
          end;

{...}					

1:

(i or j or k or (Primes[i] and (1 shl j)) = 0)