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

{---------------------------------------------------------------------------}
procedure mp_readrev(const a: ansistring; var b: mp_int);
  {-a umgekehrt nach b einlesen, nur Ziffern  0..9, Quickhack, ohne Checks}
var
  k: integer;
  d: mp_digit;
begin
  mp_zero(b);
  for k:=length(a) downto 1 do begin
     case a[k] of
       '0'..'9':  begin
                    mp_mul_d(b,10,b);
                    mp_add_d(b,ord(a[k])-ord('0'),b);
                  end;
           else   exit;
     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:
36:

procedure mp_readrev(const a: ansistring; var b: mp_int);
  {-a umgekehrt nach b einlesen, ganz ohne Checks}
var
  i,
  k: integer;
  d: mp_digit;
begin
  mp_zero(b);
  i := a;
  k := a div  9; //wenn m_digit = int32 passt 1e9
  for k:=k downto 1 do 
  begin
    d :=      (ord(a[i  ])-ord('0'));
    d := 10*d+(ord(a[i-1])-ord('0'));
    d := 10*d+(ord(a[i-2])-ord('0'));
    d := 10*d+(ord(a[i-3])-ord('0'));
    d := 10*d+(ord(a[i-4])-ord('0'));
    d := 10*d+(ord(a[i-5])-ord('0'));
    d := 10*d+(ord(a[i-6])-ord('0'));
    d := 10*d+(ord(a[i-7])-ord('0'));
    d := 10*d+(ord(a[i-8])-ord('0'));
    mp_mul_d(b,1000000000,b);
    mp_ad_d(b,d,b);
    dec(i,9);
  end;
  //Die restlichen 
  k := 1;
  d := 0;
  For i := i downto 0 do
  begin
    d := 10*d + (ord(a[i])-ord('0'));
    k := k*10;
  end; 
  mp_mul_d(b,k,b);
  mp_ad_d(b,d,b);
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:

procedure mp_readrev2(const a: ansistring; var b: mp_int);
  {-a umgekehrt nach b einlesen, nur Ziffern  0..9, Quickhack, ohne Checks}
var
  k: integer;
  d,s: mp_digit;
const
  dm = MP_DIGIT_MAX div 10;
begin
  mp_zero(b);
  d := 0;
  s := 1;
  for k:=length(a) downto 1 do begin
     case a[k] of
       '0'..'9':  begin
                    if s>dm then begin
                      mp_mul_d(b,s,b);
                      mp_add_d(b,d,b);
                      d := 0;
                      s := 1;
                    end;
                    d := d*10 + ord(a[k])-ord('0');
                    s := s*10;
                  end;
           else   exit;
     end;
  end;
  if d>0 then begin
    mp_mul_d(b,s,b);
    mp_add_d(b,d,b)
  end;
end;					

1:
2:
3:
4:
5:

      mp_miller_rabin(wx,1,prim);
      while not prim do begin
        mp_add_d(wx,2,wx);
        mp_miller_rabin(wx,1,prim);
      end;					

1:
2:

      //Umdrehen Zeile 91
      k:=mp_decimal(wx);					

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:

{$IFDEF FPC}
  {$MODE DELPHI}
{$ENDIF}  

procedure StrAdd(var K: AnsiString;kleineZahl:integer);
var
  i,r,ziffer : integer;
begin
  IF KleineZahl <= 0 then
  begin
    IF K ='' then
      k := '0';
    Exit;
  end;  
    
  i := length(k);
  IF i = 0 then
    begin
    k := '0'; 
    i :=1;
    end;
    
  while (KleineZahl > 0) do
  begin
    r := KleineZahl div 10;
    Ziffer := kleineZahl - 10*r; // mod 10;
    KleineZahl := r;        // div 10
    
    Ziffer := Ord(k[i])+Ziffer;
    If Ziffer > Ord('9') then
    begin
      inc(KleineZahl);      //Carry
      dec(Ziffer,10);
    end;
    k[i] := chr(Ziffer); //
    dec(i);
    IF (KleineZahl<>0) AND (i = 0) then
      begin
      k := '0'+k;
      inc(i);
      end;
  end;
end;

var
  z : Ansistring;
  i: integer;
  
Begin
  z := '9999999999999444444445';
  For i := 1 to 11 do  
  begin 
    StrAdd(z,555555555);
    writeln(i:4,'  ',z);
  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:

    k := edit1.text;
    //ungerade Zahl
    mp_mod_int(wx,2,rest);
    if rest=0 then 
    begin
      mp_add_d(wx,1,wx);
      StrAdd(k,1);
    end;
    //Folgensuche
    repeat
      mp_miller_rabin(wx,1,prim);
      while not prim do begin
        mp_add_d(wx,2,wx);
        StrAdd(k,2);
        mp_miller_rabin(wx,1,prim);
      end;
      mp_copy(wx,wx3);
      q:=0;

      //Umdrehen
//      k:=mp_decimal(wx);
       mp_readrev(k,??);
 ...					

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

Test of MP library V1.22.03 (31/32 bit)   (c) W.Ehrhardt 2006
nextsievex: ......Time [s]: 6.851
mp_Prime_Next_Prime: ......Time [s]: 12.428
mp_nexxxt: ......Time [s]: 8.281
mp_nextprime: ......Time [s]: 7.797
nextprime: ......Time [s]: 6.329
nextsieve: ......Time [s]: 6.939
nextsievew2: ......Time [s]: 6.932
nextsievew3: ......Time [s]: 6.932
nextsievew4: ......Time [s]: 6.038
nextsievew6: ......Time [s]: 6.040