function [LUT, M] = step1_initialization(S, k)
  m = size(S, 1);
  LUT = zeros(1, m);
  M = zeros(k, size(S, 2));
  
  step = floor (m/k);
  
  M = S(1:step:end-1,:);
end
