快速幂

1.a^b可以写成

a*a*a*a....a(b个a相乘)

而b又可以二进制分解为 b=c*2^k-1+c*2^k-2.....+c*2^0(c为该数在二进制下对应数字)

所以 用b>>1可以舍去b的最低位,而b&1可以求出b的最低位

所以

#include<cstdio>
using namespace std;

int pow(int a,int b,int p){
    int ans=1%p;
    for (;b;b>>=1){
        if(b&1) ans=ans*a%p;
        a=a*a%p;
    }
    return ans;
}

int main(){
    int a,b,p;
    scanf("%d%d%d",&a,&b,&p);
    printf("%d\n",pow(a,b,p));
return 0;
}

2.64位整数乘法

求a*b%p

#include<cstdio>
using namespace std;

long long llmod(long long a,long long b,long long p){
    long long ans=0;
    for (;b;b>>=1){
        if(b&1) ans=(ans+a)%p;
        a=a*2%p;
    }
    return ans;
}

int main(){
    int a,b,p;
    scanf("%d%d%d",&a,&b,&p);
    printf("%d\n",llmod(a,b,p));
return 0;
}

 

posted @ 2018-04-15 10:54  lmjer  阅读(111)  评论(0编辑  收藏  举报