SICP Exercise 1.16.rb

# Iterative O(log n) exponentiation
def fast_expt(b, n)
  a = 1
  while n > 0
    if even?(n)
      n /= 2
      b *= b
    else
      n -= 1
      a *= b
    end
  end
  a
end