power

　　这个算法为SGI专属，并不在标准STL之列。用它来计算某数的n幂次方。这里所谓的n幂次是指自己对自己进行某种运算，达n次。运算类型可由外界指定；如果指定为乘法，就是乘幂。

　　版本1：乘幂

template <class T, class Integer>
inline T power(T x, Integer n)
{
    return power(x, n, multiplies<T>());//指定运算形式为乘法 
}

　　

　　版本2：幂次方。如果指定为乘法运算，则当n >= 0时返回 x^n。

　　注意，"MonoidOperation"必须满足结合律，但不需要满足交换律。

template <class T, class Integer, class MonoidOperation>
T power(T x, Integer n, MonoidOperation op)
{
    if( n == 0)
        return identity_element(op);//取出“证同元素”identity element 
    else //所谓“证同元素” 
    {
        while((n & 1) == 0)
        {
            n >>= 1;
            x = op(x, x);
        }
        T result = x;
        n >>=1;
        while ( n != 0)
        {
            x = op(x, x);
            if((n & 1) != 0)
                result = op(result, x);
            n >>= 1;
        }
        return result;
    }
}