快速幂和快速乘

快速幂

转载：

算法学习笔记(4)：快速幂 - 知乎 (zhihu.com)

快速幂（Exponentiation by squaring，平方求幂）是一种简单而有效的小算法，它可以以O(log⁡n)的时间复杂度计算乘方。快速幂不仅本身非常常见，而且后续很多算法也都会用到快速幂。

让我们先来思考一个问题：7的10次方，怎样算比较快？

方法1：最朴素的想法，7*7=49，49*7=343，... 一步一步算，共进行了9次乘法。

这样算无疑太慢了，尤其对计算机的CPU而言，每次运算只乘上一个个位数，无疑太屈才了。这时我们想到，也许可以拆分问题。

方法2：先算7的5次方，即7*7*7*7*7，再算它的平方，共进行了5次乘法。

但这并不是最优解，因为对于“7的5次方”，我们仍然可以拆分问题。

方法3：先算7*7得49，则7的5次方为49*49*7，再算它的平方，共进行了4次乘法。

模仿这样的过程，我们得到一个在 O(log⁡n) 时间内计算出幂的算法，也就是快速幂。

1）递归快速幂

其中 odd是奇数，even是偶数

//递归快速幂
int qpow(int a, int n)
{
    if (n == 0)
        return 1;
    else if (n % 2 == 1)
        return qpow(a, n - 1) * a;
    else
    {
        int temp = qpow(a, n / 2);
        return temp * temp;
    }
}

 JAVA版本代码

    public double quickMul(double x, long N) {
        if (N == 0) {
            return 1.0;
        }
        double y = quickMul(x, N / 2);
        return N % 2 == 0 ? y * y : y * y * x;
    }

2）迭代快速幂

int qpow(int a, int n){
    int ans = 1;
    while(n){
        if(n&1)        //如果n的当前末位为1
            ans *= a;  //ans乘上当前的a
        a *= a;        //a自乘
        n >>= 1;       //n往右移一位
    }
    return ans;
}

JAVA版本代码 

    public double quickMul(double x, long N) {
        double ans = 1.0;
        // 贡献的初始值为 x
        double x_contribute = x;
        // 在对 N 进行二进制拆分的同时计算答案
        while (N > 0) {
            if (N % 2 == 1) {
                // 如果 N 二进制表示的最低位为 1，那么需要计入贡献
                ans *= x_contribute;
            }
            // 将贡献不断地平方
            x_contribute *= x_contribute;
            // 舍弃 N 二进制表示的最低位，这样我们每次只要判断最低位即可
            N /= 2;
        }
        return ans;
    }

 

最初ans为1，然后我们一位一位算：

1010的最后一位是0，所以a^1这一位不要。然后1010变为101，a变为a^2。

101的最后一位是1，所以a^2这一位是需要的，乘入ans。101变为10，a再自乘。

10的最后一位是0，跳过，右移，自乘。 

然后1的最后一位是1，ans再乘上a^8。循环结束，返回结果

快速乘

 

        在不用乘法的前提下，实现倍乘，这时候就可以用快速乘算法。它主要用加法和二进制的位运算来实现。同时，由于计算机底层设计的原因，做加法往往比乘法快的多，因此将乘法转换为加法计算将会大大提高（大数，比较小的数也没必要）乘法运算的速度，除此之外，当我们计算a*b%mod的时候，往往较大的数计算a＊b会超出long long int的范围，这个时候使用快速乘法方法也能解决上述问题．类似

两数相除 - 两数相除 - 力扣（LeetCode）

1）递归法

    public static int mul2(int a,int b){
        int result;
        if (b==0)
            return 0;
        else if (b%2==1){//奇数
            result=mul2(a,b-1)+a;
        }
        else {
            int temp=mul2(a,b/2);
            result=temp*2;
        }
        return result;
    }

2）迭代法

public int mul(int a,int b){
    int ans=0;          //定义一个返回值
    while(b>0){
        if(b&1==1){
             ans+=a;     //b最后一位为1时，结果可以加上a 
        }
        a+=a;         //a乘2
        b>>=1;        //b往右移一位
    }
    return ans;
}