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

function  jacobi32(a,b: longint): integer;
  {-Compute the Jacobi/Legendre symbol (a|b), b: odd and > 2}

function  kronecker32(a,b: longint): integer;
  {-Compute the Kronecker symbol (a|b)}

function  mp_jacobi(const a,n: mp_int): integer;
  {-Compute the Jacobi/Legendre symbol (a|n), n: odd and > 2}

function  mp_jacobi_lm(a: longint; const n: mp_int): longint;
  {-Compute the Jacobi/Legendre symbol (a|n), n: odd and > 2}

function  mp_jacobi_ml(const a: mp_int; n: longint): longint;
  {-Compute the Jacobi/Legendre symbol (a|n), n: odd and > 2}

function  mp_kronecker(const a,n: mp_int): integer;
  {-Compute the Kronecker symbol (a|n)}                    

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:

function jac(x,y:integer):integer;
var m8, temp, res: integer;
begin
    if (y<3) or (not odd(y)) then
    begin
       jac:=0;
       exit
    end;
    res:=1;
    while true do
    begin
       x:=x mod y;
       if x=0 then begin jac:=0; exit end;
       while not odd(x) do
       begin
          x:=x div 2;
          m8:= y mod 8;
          if (m8 = 3) or (m8 = 5) then res := -res;
       end;
       if x = 1 then begin jac:=res; exit end;
       temp:= x;
       x := y;
       y := temp;
       if (x mod 4 = 3) and (y mod 4 = 3) then res:=-res;
    end;
end;