pat每日刷题计划--day71

int的范围是-2^31~2^31-1

快速幂（所有可能溢出的地方都要mod）

递归写法

几点注意事项：1、temp在里面要先算出来一半再temp*temp,不然两次进循环复杂度增加

                         2、%m不要忘了写

#include<stdio.h>
#include<iostream>
#include<string.h>
using namespace std;
int binaryPow(int a,int b,int m)
{
    int temp;
    if(a==0) return 0;
    if(m==1) return 0;
    if(b==0) return 1;
    if(b==1) return a%m;
    a=a%m;
    if(b&1)
    {
        return a*binaryPow(a,b-1,m)%m;
    }
    else
        temp=binaryPow(a,b/2,m)%m;
        return temp*temp%m;
}
int main()
{
    int a,b,m;
    scanf("%d%d%d",&a,&b,&m);
    printf("%d^%d mod %d=%d",a,b,m,binaryPow(a,b,m));
    return 0;
}

View Code

迭代写法

注意a*a的位置也要modm 不然可能会溢出导致结果错误。

迭代中 eg:b=10  对应的二进制1001 这时候就是a^8*a^2

里面的a一直随着b走向高位每次都是a*a(注意是a*a不是a平方，因为二进制位数符号是16 8 4 2 这种的)

还有就是第一次是要先动ans再动a*a,否则会错位

#include<stdio.h>
#include<iostream>
#include<string.h>
using namespace std;
int binaryPow(int a,int b,int m)
{
    if(m==1) return 0;
    int ans=1;
    a=a%m;
    while(b>0)
    {

        if(b%2==1)
        {
            ans=ans*a%m;
        }
        a=a*a%m;
        b=b/2;
    }
    return ans;
}
int main()
{
    int a,b,m;
    scanf("%d%d%d",&a,&b,&m);
    printf("%d^%d mod %d=%d",a,b,m,binaryPow(a,b,m));
    return 0;
}

View Code

 实际操作差别不大，我觉得递归更方便一点