快速幂 快速乘法

#include <iostream>
#include <cstdio>
#include
#include <cmath>
#include <cstring>
#include<map>
using namespace std;

long long q_mul( long long a, long long b, long long mod ) //快速计算 (a*b) % mod
{
long long ans = 0;  // 初始化
while(b)                //根据b的每一位看加不加当前a
{
if(b & 1)           //如果当前位为1
{
b--;
ans =(ans+ a)%mod;   //ans+=a
}
b /= 2;                         //b向前移位
a = (a + a) % mod;          //更新a

}
return ans;
}

long long q_pow( long long a, long long b, long long mod ) //快速计算 (a^b) % mod
{
long long ans = 1; // 初始化
while(b)//根据b的每一位看乘不乘当前a
{
if(b & 1)   //如果当前位为1
{
ans = q_mul( ans, a, mod ); //ans*=a
}
b /= 2;                                     //b向前移位
a = q_mul( a, a, mod );         //更新a
}
return ans;
}

int main( )
{
long long a, b, n;
while(cin >> a >> b >> n)
{
cout << "a*b%n = " << q_mul( a, b, n ) << endl;
cout << "a^b%n = " << q_pow( a, b, n ) << endl;
}
return 0;
}