码迷,mamicode.com
首页 > 其他好文 > 详细

power(乘幂)函数剖析

时间:2017-07-02 22:16:35      阅读:241      评论:0      收藏:0      [点我收藏+]

标签:ret   iat   交换   returns   开始   ext   意思   简洁   else   

近来学习STL,看到power函数的实现感觉挺有趣,记录一下。

1. 一般情况下,我自己要实现乘幂函数会这样实现:

int power(int x,size_t n)
{
    int result = 1;
    while (n--)
        result *= x;
    return result;
}

这样即使实现,这里的时间复杂度和n有关,时间复杂度为0(n)。

2. 看了stl源码实现是这样:

// Returns __x ** __n, where __n >= 0.  _Note that "multiplication"
// is required to be associative, but not necessarily commutative.  
//意思是multiplication要满足结合律,但不需要满足交换律

template <class _Tp, class _Integer, class _MonoidOperation>  
_Tp __power(_Tp __x, _Integer __n, _MonoidOperation __opr)  //这里第三个参数是一个二元操作函数
{
    if (__n == 0)
        return identity_element(__opr);           //返回1
    else 
    {
        while ((__n & 1) == 0)                   //如果n为偶数
        {
            __n >>= 1;
            __x = __opr(__x, __x);
        }

        _Tp __result = __x;
        __n >>= 1;
        while (__n != 0) 
        {
            __x = __opr(__x, __x);
            if ((__n & 1) != 0)                    //如果n为奇数
                __result = __opr(__result, __x);
            __n >>= 1;
        }
        return __result;
    }
}

template <class _Tp, class _Integer>
inline _Tp __power(_Tp __x, _Integer __n)
{
    return __power(__x, __n, multiplies<_Tp>());
}

// Alias for the internal name __power.  Note that power is an extension,
// not part of the C++ standard.

template <class _Tp, class _Integer, class _MonoidOperation>
inline _Tp power(_Tp __x, _Integer __n, _MonoidOperation __opr)   //也可以自定义,传进去一个函数
{
    return __power(__x, __n, __opr);
}

template <class _Tp, class _Integer>
inline _Tp power(_Tp __x, _Integer __n)   //默认情况下是乘幂
{
    return __power(__x, __n);
}

这里:当n为偶数时,X^n=(X^2)^(2/n),此时看2/n是否还是偶数,如果是则继续,否则,转向n是奇数的情况;

           当n为奇数时,X^n=X*X^(n-1);

例如计算5^6=(5^2)^3,也就是计算25^3。

这种情况和第三种情况类似,只是当n开始为偶数时,比第三种方法的效率更高,n为奇数时,和第三种方法时间复杂度一样。

3. 看了还有一种实现,和stl源码差不多,只是表示更简洁,和第二种情况相比效率会差点:

int power(int x, size_t n)
{
    if(n == 0)
        return 1;
    int result = 1;
    while(n)
    {
        if(n & 1)                    //n为奇数时
            result *= x;            
        n >>= 1;
        x *= x;                 
    }
}

如计算5^6是6化为二进制为0110,所以这里5^6=5^2*5^4,这里时间复杂度为0((log2n)+1)

power(乘幂)函数剖析

标签:ret   iat   交换   returns   开始   ext   意思   简洁   else   

原文地址:http://www.cnblogs.com/liuamin/p/7107202.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!