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

Wiederhole
  Miller-Rabin-TestMIt (t), t+2,t+6, t+8
  Falls Vierling gefunden dann eintragen/ausgeben
 
  wiederhole
    erhöhe t um 30 
    erhöhe P_v0 um 30 
    falls P_v0 > 2*3*5*..prim_n then
      P_v0 := 1
  bis (teste mit (P_v0) AND(P_v0+2)AND (P_v0+6)AND (P_v0+8) mit p=2 bis prim_n auf teilerfremd = wahr)
bis JetztIstGenug					

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:

program primvier;
{$IFDEF FPC}
   {$MODE Delphi}
{$ENDIF}   

{
10^(n-1)+a, 10^(n-1)+a+2, 10^(n-1)+a+6 und 10^(n-1)+a+8
haben Abstand 30 und Beginn bei 1 (1,3,7,9 -> war wohl nichts..)
Deshalb kein Test auf Teilbarkeit 2,3,5 
}
{Primzahlen bis 100
2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97
}


const
  TestTeiler : array[0..11] of integer=
                (7,11,13,17,19,23,29,31,37,41,43,47);
//MAXANZAHL4 = TestTeiler[0]*TestTeiler[1]*TestTeiler[0]*.. )
  MaxAnzahl4_1 :DWORD = 7*11*13*17*19*23*29;
  MAXANZAHL4_2 :DWORD = 31*37*41*43*47; 

var 
  cnt: integer;
  
function istTeilbar(t: int64 ): boolean;
var
  i: integer;
begin
  result := true;

  For i := Low(Testteiler) to High(Testteiler) do
    begin
    result := result AND ((t MOD Testteiler[i]) <> 0);
    iF Not(result) then
      EXIT;
    end;  
end;
var
  i,j : DWord;
  p1 :int64;
Begin
 writeln(MaxAnzahl4_1);
 writeln(MaxAnzahl4_2);
  cnt := 0;
  p1 := 1;  
  For j := 1 to MaxAnzahl4_2 do
    For i := 1 to MAXANZAHL4_1 do
      begin
      IF istTeilbar(p1) AND istTeilbar(p1+2)AND istTeilbar(p1+6)AND istTeilbar(p1+8) then
        begin
        inc(cnt);
        IF cnt <> 0 then
          begin
          // Primzahl, Anzahl der möglichen Vierers, noch verbliebene 
          writeln(p1:20,i:10,cnt:10,cnt/i*100.0:10:3,'%');
          end;
       end;
      inc(p1,30); 
      IF p1 > 15000 then 
        halt;  
      end; 
end.
215656441
95041567
                 541        19         1     5.263%
                 571        20         2    10.000%
                1291        44         3     6.818%
                1621        55         4     7.273%
                2851        96         5     5.208%
                3181       107         6     5.607%
                3301       111         7     6.306%
                3691       124         8     6.452%
                4441       149         9     6.040%
                6241       209        10     4.785%
                6421       215        11     5.116%
                7171       240        12     5.000%
                8521       285        13     4.561%
                9181       307        14     4.560%
                9811       328        15     4.573%
                9901       331        16     4.834%
               10111       338        17     5.030%
               11461       383        18     4.700%
               11491       384        19     4.948%
               13261       443        20     4.515%
               14011       468        21     4.487%
               14101       471        22     4.671%
               14401       481        23     4.782%
               14821       495        24     4.848%
Drücken Sie eine beliebige Taste . . .					

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

100000000000000000000453777211
100000000000000000000456577981
100000000000000000000463041031
100000000000000000000463596061
100000000000000000000467236321
100000000000000000000468464821
100000000000000000000470324671
100000000000000000000472370071
Abbruch bei 100000000000000000000473190241