Возведение матрицы в степень

[Jumbo]

Написал прогу, возводящую квадратную матрицу N*N в степень K и ответ по модулю S (1<=N<=100,1<=K<=10^9,1<=s<=1000). Но выдаёт неправильный ответ на 4 тесте. Забивал разные матрицы, возводил в первую степень, возводил матрицу из одного элемента - везде ответы сходились.

program power;

{$APPTYPE CONSOLE}

uses
  SysUtils;
type mas=array[1..100,1..100] of int64;
var t,res:mas;  f,f1:text; k,n,s,i,j:integer;

function mul(a,b:mas):mas;{перемножение матриц}
var i,j,d:integer; c:mas;
begin
 for i:=1 to n do
  for j:=1 to n do begin
   c[i,j]:=0;
   for d:=1 to n do
    c[i,j]:=(c[i,j]+ a[i,d]*b[d,j]) mod s ;
  end;

 mul:=c;
end;

procedure power_m(k: integer);{возведение матрицы t в степень k}
var i:integer;
begin
  for i:=1 to n do res[i,i]:=1;

  while (k > 0) do begin

    if (k and 1 = 1) then
      res := mul(res,t);
    if k<>1 then t := mul(t,t);
    k := k shr 1;

  end;

end;

begin
 reset(f,'power.in');
 rewrite(f1,'power.out');
 readln(f,n,k,s);

 for i:=1 to n do begin
  for j:=1 to n do read(f,t[i,j]);
  readln(f);
 end;
 close(f); {считан массив}

 power_m(k);

 for i:=1 to n do begin
  for j:=1 to n do write(f1,res[i,j] mod s,' ');
  writeln(f1);
 end;
 close(f1);
end.