快速幂与快速乘算法C++代码

快速幂与快速乘算法，参考博客

#include <iostream>
using namespace std;
long long qpow(long long a, long long n, long long mod){	//快速幂，计算a^n % mod
    long long re = 1;
    while(n){
        if(n & 1)		//判断n的最后一位是否为1(n&1：n的末尾&1)
            re = (re * a) % mod;
        n >>= 1;		//舍去n的最后一位(n/=2) 
        a = (a * a) % mod;	//将a平方
    }
    return re;
}
long long qmul(long long a, long long n, long long mod){	//快速乘，计算a*n % mod
    long long re = 0;
    while(n){
        if(n & 1) 
			re = (re + a) % mod;
        n >>= 1;
		a = (a + a) % mod;
    }
    return re;
}
int main(){
	long long a,n,mod;
	cin>>a>>n>>mod;
	cout<<qpow(a,n,mod)<<" "<<qmul(a,n,mod);
	return 0;
}