#
# Usage:
#  factorial(n) = n*(n-1)*...*1
#  factorial(n, m) = n*(n-1)*...*m
def factorial (n, m=1)
  if m < n then
    k=(m+n)/2;
    factorial(k, m) * factorial(n, k+1);
  else
    m
  end
end