Оптимизация кода

program e_1;
function NOK(n1,n2:LONGINT):LONGINT; //Функция нахождения НОК
var m,nk:LONGINT;
begin
if n1>n2 then m:=n1 else m:=n2;
nk:=m-1;
repeat
 inc(nk);
until ((nk mod n1)=0) and ((nk mod n2)=0);
NOK:=nk;
end;
function NOD(A, B:LONGINT):LONGINT; //Функция нахождения НОД
begin
 while a<>b do begin
 if a>b then a:=a-b else b:=b-a
 end;
 NOD := a
end;
var //переменные
input, HOD : array of LONGINT; // массивы
a, b, c, i, i2, i3, w, OTV, N, NHOD : LONGINT; // a, b, c - числа входящие в тройку для нахождения НОД. N - Кол-во символов первого массива, NHOD - Кол-во символов втрого массива.
begin
NHOD := 0; // обнуляем счетчик для второго массива
readln (N); // читаем кол-во символов для первого массива
SetLength (input, (N+1)); // установливаем длинну первого массива
for i := 1 to N do
read (input[i]); // записаваем входные данные
for i := 1 to N do   // перебор всех возможных троек чисел
begin                // перебор всех возможных троек чисел
a := input[i];       // перебор всех возможных троек чисел
for i2 := 1 to N do  // перебор всех возможных троек чисел
begin                // перебор всех возможных троек чисел
b := input[i2];      // перебор всех возможных троек чисел
for i3 := 1 to N do  // перебор всех возможных троек чисел
begin                // перебор всех возможных троек чисел
c := input[i3];      // перебор всех возможных троек чисел
if (i=i2) or (i2=i3) or (i=i2) then // убираем повторы
else
begin
NHOD := NHOD+1; // считаем длинну массива для НОД
SetLength (HOD, (NHOD+1)); // устанавливаем длинну массива
HOD[NHOD] := NOD(a, NOD(b,c)); // записываем в массив НОД троек чисел
end;
end;
end;
end;
OTV := HOD[1]; // обнуляем ответ
for w:=2 to NHOD do OTV := NOK (OTV, HOD[w]); // ищем НОК все чисел из второго массива
OTV := OTV mod (1000000000+7); // делим ответ на 10^9+7
writeln (OTV);
end.

function GCD(u, v: longint): longint;
var
  k: longint;
begin
  if (u > 1) and (v > 1)
    then
      begin
        k := 0;
        while not odd(u or v) do
          begin
            inc(k);
            u := u shr 1;
            v := v shr 1
          end;
        while v > 0 do
          begin
            while not odd(u) do u := u shr 1;
            while not odd(v) do v := v shr 1;
            if u > v
              then u := u - v
              else v := v - u;
          end;
        GCD := u shl k
      end
    else GCD := 1
end;
 
var
  n, i, j, k, LCM, t: longint;
  a: array of longint;
 
begin
  readln(n);
  setlength(a, n);
  for i := 0 to n - 1 do read(a[i]);
  readln;
  LCM := 1;
  for i := 0 to n - 3 do
    for j := i + 1 to n - 2 do
      begin
        t := GCD(a[i], a[j]);
        for k := j + 1 to n - 1 do
          begin
            t := GCD(t, a[k]);
            LCM := LCM * t div GCD(LCM, t)
          end
      end;
  write(LCM mod 1000000007)
end.