Поиск максимальной подматрицы

/Отсортировать n строк B в порядке убывания в списки list[i], т.о. в каждом списке list[i] будут находится индексы упорядоченных элементов строки i;
/Заведём множество V = {1, ..., n};
for i := 1 to n do begin
  for k := 1 to n do begin
    cand[k]:=first of list[k];
    d[k] := a[i, k] + b[k, cand[k]];
  end;
  / Упорядочим V как очередь с приоритетом по ключам d[1], ..., d[n];
  / Множество решений S - пустое;
  /* Фаза 1 : До критической точки */
  while |S| ≤ n − n log n do begin
    /Найти v с минимальным ключем в очереди;
  /Занести cand[v] в множество S;
    c[i, cand[v]] := d[v];
    /Пусть множество W = {w|cand[w] = cand[v]};
  for w in W do
      while cand[w] is in S do cand[w]:= next of list[w];
  /Отсортировать очередь W по новым ключам d[w] = a[i, w]+b[w, cand[w]];
  end;
  U := S;
  /* Фаза 2 : После критической точки */
  while |S| < n do begin
    /Найти v с минимальным ключём в очереди;
    if cand[v] is not in S then begin
      /Занести cand[v] в множество S;
      c[v, cand[v]] := d[v];
      /Пусть множество W = {w|cand[w] = cand[v]};
    end else W = {v};
    for w in W do
      cand[w]:=next of list[w];
    while cand[w] is in U do cand[w]:= next of list[w];
    /Отсортировать очередь W по ключам d[w] = a[i, w]+b[w, cand[w]];
  end;
end.

Uses Crt;
const
n=1000;
m=1000;
var
a: array [1..n, 1..m] of integer;
b: array [1..n, 1..m] of integer;
buf,c,d,e,f,k:integer;
i,j:byte;
begin
clrscr;
writeln('vvedite rasmetiri glavnih i pobosnih mtrx');
read(c,d,e,f);
writeln('glav mtrx');
for i:=1 to c do 
for j:=1 to d do
a[i,j]:=random(2);
for i:=1 to c do begin
for j:=1 to d do 
write(a[i,j]:2);
writeln;
end;
writeln('pobos mtrx');
for i:=1 to e do 
for j:=1 to f do
b[i,j]:=random(2);
for i:=1 to e do begin
for j:=1 to f do
write(b[i,j]:2);
writeln;
end;
buf:=0;
k:=1;
for i:=1 to c do begin
for j:=1 to d do 
if a[k,k]+b[i,j]<=1 then begin
buf:=1;
end;
k:=k+1;
end;
for i:=c downto 1 do begin
for j:=d downto 1 do 
if buf=1 then begin
a[k,k]:=a[k,k]+b[i,j];
end;
k:=k-1;
end;
writeln;
for i:=1 to c do begin
for j:=1 to d do 
write(a[i,j]:2);
writeln;
end;
end.