monoid exponentiation

import Data.Monoid

power :: (Monoid a, Integral i) => a -> i -> a
power _ 0 = mempty
power a n | r==0 = b
          | otherwise = a `mappend` b
  where (q,r) = n `divMod` 2
        b = (a `mappend` a) `power` q