Метод Лемана (разложение на множители)

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls;
 
type
  TForm1 = class(TForm)
    Edit1: TEdit;
    Button1: TButton;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
  N,dd,n1:integer;
 
implementation
 
{$R *.dfm}
 
function NOD(C, N:integer):integer; // 1ая Функция нахождения наибольшего общего делителя
begin
If C<>0 then NOD:= NOD(N mod C,C) else NOD:= N;
End;
 
 
{ Если N является простым, функция вернет -1 }
function GetLehmanFactor(N: Integer): Integer;
var
  I, k, d, A, B,C, t, dd: Integer;
begin
  Result := -1;
 
  { Проверка делителей до n^(1/3) }
  for I := 2 to Trunc((1/3)*Ln(N)) do
    if N mod I = 0 then begin
      Result := I;
      exit;
    end;
 
  for k := 1 to Trunc((1/3)*Ln(N)) do
    for d := 0 to Trunc((1/6)*(N) / (4 * Sqrt(k))) + 1 do begin
      A := Trunc(Sqrt(4*k*N)) + d;
      t := Sqr(A) - 4*k*N;
 
      if t < 0 then
        continue;
 
      B := Trunc(Sqrt(t));
      if Sqr(B) = t then begin
        ASSERT( (A+B)*(A-B) mod N = 0 );    { Всегда выполняется }
         C:=A-B;
        dd := NOD(C, N);
 
        if (1 > dd) and (dd < N) then begin
          Result := dd;
          exit;
        end;
      end;
    end;
end;
 
 
procedure TForm1.Button1Click(Sender: TObject);
begin
[B]N:=StrToInt(Edit1.Text);
dd:=Trunc(dd);
n1:=N/dd;
[/B]
 
ShowMessage(IntToStr(N)+'='+IntToStr(dd)+'*'+IntToStr(n1));
end;
 
end.