移位运算

a ^ b

题目大意

在 $a$ 的 $b$ 次方对 $p$ 取模，其中 $1\le a,b,p\le 10^9$

思路

根据数学常识，每一个正整数都可以为一的表示为若干指数不重复的 $2$ 的次幂的和。$b=c_{k-1}{2}^{k-1}+c_{k-2}{2}^{k-2}+\dots \dots +c_0{2}^0$。$c_i$ 表示在二进制下 $b$ 的第 $i$ 位是 $0$ 还是 $1$。

时间复杂度

因为 $k=\lceil log(b+1)\rceil$ ,通过 $k$ 次递推出每一个乘积项，因此整个算法的时间复杂度为 $lo{g}_{2}b$。

代码

#include<bits/stdc++.h>
using namespace std;

typedef long long LL;

const int N = 1e5 + 10;
 
pair<string,int> a[N];
int m,n;

int main()
{
	int a,b,p;
	cin >> a >> b >> p;
	LL ans = 1 % p;
	
	for( ; b; b >>= 1){
		if(b & 1){
			ans = (LL)ans * a % p;// a ^ b 
		}
		a = (LL)a * a % p;
	}
	cout << ans << endl;
	return 0;
}

---

64位整数乘法

题目大意

求 $a$ 乘 $b$ 对 $p$ 的取模，其中 $1\le a,b,p\le 10^{18}$

解法 $1$

思路

类似于快速幂的思想，把整数 $b$ 用二进制表示，$b=c_{k-1}{2}^{k-1}+c_{k-2}{2}^{k-2}+\dots \dots +c_0{2}^0$，则 $a\ast b=a\ast (c_{k-1}{2}^{k-1}+c_{k-2}{2}^{k-2}+\dots \dots +c_0{2}^0)$。

代码

#include<bits/stdc++.h>
using namespace std;

typedef long long LL;

const int N = 1e5 + 10;
 
pair<string,int> a[N];
int m,n;

int main()
{
	LL a,b,p;
	cin >> a >> b >> p;
	LL ans = 0;
	for(; b; b >>= 1){
		if(b & 1){
			ans = (ans + a) % p;
		}
		a = a * 2 % p;
	}
	cout << ans << endl;
	return 0;
} 

解法 $2$

思路

利用 $a\ast b$ mod $p$ = $a\ast b$ - $\lfloor a\ast b/p\rfloor \ast p$。在 $a\ast b/p$ 中我们可以用浮点数来执行，因为浮点数的在十进制下的有效数字有 $18、19$ 位$(因为相除之后会变成小数点后的位数，足够存储)$,足够胜任。

代码

#include<bits/stdc++.h>
using namespace std;

typedef long long LL;

const int N = 1e5 + 10;
 
pair<string,int> a[N];
int m,n;

int main()
{
	LL a,b,p;
	cin >> a >> b >> p;
	a %= p,b %= p;
	LL c = (long double)a * b / p;
	LL ans = a * b - c * p;
	if(ans < 0) ans += p;
	else if(ans >= p)ans -= p;
	cout << ans << endl;
	return 0;
}