这里介绍一种分数取模的代码:(参考)
假设要求 (1/m) mod p
这里要引用小费马定理 a^p-1 mod p = 1 mod p (有兴趣的可以百度查下证明过程) ,这里对这个公式做点改动ap-1=ap-2 *a, 把 a 移到上述等式的右边有 a^p-2 mod p = a^-1 mod p , 那么这里,
a^-1 mod p 就是 我们要求的,这个值的结果恒等于 a^p-2 mod p ,
p是素数
long long fast_mod(long long a,long long b) {
long long r = 1;
a %= mod;
while (b) {
if (b & 1) r = (r*a) % mod;
a = (a*a) % mod;
b >>= 1;
}
//cout << r << endl;
return r;
}
欧拉定理:
这里先介绍欧拉函数
设m,n是互素的两个正整数,则
ϕ(mn)= ϕ(m) ϕ(n).
其中ϕ(m) , ϕ(n) 都是欧拉函数
`