取余运算

问题描述

输入b，p，k的值，求b^p mod k的值。其中b，p，k*k为长整形数。

分析：参考以下公式将数据进行分解

如：
P=2 * P/2 + P%2，
如19=2 * 19/2 + 19%2=2*9+1

P 为偶数

b^p % k
=b^( 2 * P/2 )%k
=(b^( 2 * P/2 ))%k=(b^ P/2 * b^ P/2)%k
=(b^( 2 * P/2 ))%k=((b^ P/2)%k * (b^ P/2)%k)%k

temp=(b^ P/2) %k

b^p % k=((temp*temp)%k

P 为奇数

b^p % k
=b^( 2 * P/2 + 1 )%k
=(b^( 2 * P/2 )b)%k
=((b^ P/2 b^ P/2 )b)%k
=(b^( 2 * P/2 ) %kb %k)%k

temp=b^ P/2%k

b^p % k
=((temptemp)%k * b%k)%k
=(temptemp* b)%k

实现

#include<iostream>
#include<cstdio>
using namespace std;
int b,p,k,a;
int f(int p)   {
	if (p==0)
		return 1;
	int tmp;
	tmp=f(p/2)%k;
	tmp=(tmp*tmp) % k;
	if (p%2==1)
    tmp=(tmp*tmp * b) %k;
        return tmp;
}
int main(){
    cin>>b>>p>>k;
    int tmpb=b;
    b%=k;
    printf("%d^%d mod %d=%d\n",tmpb,p,k,f(p));
    return 0;
}