快速求幂取模

快速求幂取模

2016年02月06日 18:02:39

阅读数：1694

公式求幂→二分求幂→快速求幂→快速求幂取模 
等不急的可以直接下拉到最后看快速幂取模。

直接用C语言的库函数pow()（别忘了它的头文件#include<math.h>），似乎很简单，但是它的时间复杂度高达O(n)。 
显然，这很容易超时。 
于是有了下面的二分求幂（时间复杂度O(lgn)）

二分求幂的原理可以用下面这张图表示 

用递归来实现，虽然代码有点长，但是很好理解

int pow(int a,int n)//返回值是a的n次方
{
    if(n==0)//递归终止条件
        return 1;
    if(n==1)
        return a;
    int result=pow(a,n/2);//二分递归
    result=result*result;//这部分奇数偶数都一样
    if(n%2==1)//如果n是奇数，就要多乘一次
        result=result*a;
    return result;
}

• 1
• 2
• 3
• 4
• 5
• 6
• 7
• 8
• 9
• 10
• 11
• 12

用非递归，更加简洁

int pow(int a,int n)//返回值是a的n次方
{
    int result=1;
    while(n!=0)
    {
        if(n%2==1)//如果n是奇数
            result=result*a;//就要多乘一次
        a=a*a;
        n=n/2;//二分
    }
    return result;
}

• 1
• 2
• 3
• 4
• 5
• 6
• 7
• 8
• 9
• 10
• 11
• 12

快速幂顾名思义比二分幂又快一些， 
快速幂借助了强大的位运算，时间复杂度达到O(log₂N)。 
位运算http://baike.baidu.com/view/379209.htm 
原理其实不难理解，用的还是小学生的公式。 
例如 求a的11次方 
11的二进制是1011 
11 = 2³×1 + 2²×0 + 2¹×1 + 2º×1 
因此，我们将a¹¹转化为算  
用非递归的代码实现

int pow(int a,int n)//返回值是a的n次方
{
    int result=1,flag=a;
    while(n!=0)
    {
        if(n&1)//如果n是奇数，即n的二进制最末位为1时
            result=result*flag;
        flag=flag*flag;
        n=n>>1;//n的二进制右移一位，即n/2
    }
    return result;
}

• 1
• 2
• 3
• 4
• 5
• 6
• 7
• 8
• 9
• 10
• 11
• 12

当然还能用递归来实现，但是太复杂，我没学会…

刷题中让直接求幂的不多，求幂后取模的却不少，毕竟求幂结果太大了。 
水平所限，只会用二分幂取模，时间复杂度与二分幂一样O(lgn)。 
基本可以在各种比赛中顺利通过，也是目前比较常用的方法

原理同样很简单，都是小学学过的：积的取余等于取余的积取余 
接下来用代码实现

int pow(int a,int n,int b)//返回值是a的n次方对b取余后的值
{
    int result=1;
    a=a%b;//积的取余等于取余的积取余
    while(n>0)
    {
        if(n%2==1)
            result=result*a%b;//n是奇数的话就要多乘一次，原理和前面的二分求幂一样
        n=n/2;//二分
        a=a*a%b;//积的取余等于取余的积取余
    }
    return result;
}

• 1
• 2
• 3
• 4
• 5
• 6
• 7
• 8
• 9
• 10
• 11
• 12
• 13

影响计算机效率的是运算次数，而不是运算结果。 
所以前面几个算法都是通过增大运算结果，减少运算次数，提高计算机效率。