Result = a^b mod N
k = length(b);

Result = 1;
for i = k-1 to 0 do
  Result = (Result X Result) mod N;   // SQUARING

  if b[i] == 1
    Result = (Result X a) mod N;     // MULTIPLICATION
  end if

end for