数论学习之乘法逆元
2016-11-09 21:10
232 查看
用法:用于除法取模
思路:扩欧
要求:b、p互质
设k为b的乘法逆元:
则在求解除法取模问题时:
有(a/b)%p =>(a*k)%p
当b很大时,用除法会出现精度问题。。so
乘法逆元:
如果b*k ≡ 1 (mod p)
则称k是b关于p的乘法逆元
我们可以通过求 b 关于 p 的乘法逆元 k,将 a 乘上 k 再模 p,即 (a * k) mod p。其结果与(a / b) mod p等价。
证:
因为 b * k ≡ 1
(mod p)
则有 b * k = p*
x+1
得到 k = (p * x
+ 1) / b
将 k 代入(a * k) mod p
得到:
(a * (p * x + 1) / b) mod p
=((a * p * x) / b + a / b) mod p
=[((a * p * x) / b) mod p +(a / b)] mod p
=[(p * (a * x) / b) mod p +(a / b)] mod p
=(0 + (a / b)) mod p
= (a/b) mod p
用欧几里得扩展求逆元要求 gcd(b, p) == 1
求乘法逆元可以用到欧几里得扩展:
void Euild(ll a, ll b, ll &x, ll &y) // x 是 a 关于 b 的乘法逆元
{
if(0 == b){
x = 1, y = 0;
return ;
}
Euild(b, a%b, x, y);
ll flag = x;
x = y;
y = flag - a/b * y;
}
思路:扩欧
要求:b、p互质
设k为b的乘法逆元:
则在求解除法取模问题时:
有(a/b)%p =>(a*k)%p
当b很大时,用除法会出现精度问题。。so
乘法逆元:
如果b*k ≡ 1 (mod p)
则称k是b关于p的乘法逆元
我们可以通过求 b 关于 p 的乘法逆元 k,将 a 乘上 k 再模 p,即 (a * k) mod p。其结果与(a / b) mod p等价。
证:
因为 b * k ≡ 1
(mod p)
则有 b * k = p*
x+1
得到 k = (p * x
+ 1) / b
将 k 代入(a * k) mod p
得到:
(a * (p * x + 1) / b) mod p
=((a * p * x) / b + a / b) mod p
=[((a * p * x) / b) mod p +(a / b)] mod p
=[(p * (a * x) / b) mod p +(a / b)] mod p
=(0 + (a / b)) mod p
= (a/b) mod p
用欧几里得扩展求逆元要求 gcd(b, p) == 1
求乘法逆元可以用到欧几里得扩展:
void Euild(ll a, ll b, ll &x, ll &y) // x 是 a 关于 b 的乘法逆元
{
if(0 == b){
x = 1, y = 0;
return ;
}
Euild(b, a%b, x, y);
ll flag = x;
x = y;
y = flag - a/b * y;
}
相关文章推荐
- 数论学习之乘法逆元
- 数论学习之乘法逆元
- 浅谈数论(四)扩展欧几里得算法与乘法逆元
- ACM学习历程—HDU5667 Sequence(数论 && 矩阵乘法 && 快速幂)
- 基础数论学习笔记-----------逆元【除法取模运算的关键】
- 数论继续学习6---数论倒数,又称逆元
- 数论学习_逆元意义及求法
- 浅谈数论(四)扩展欧几里得算法与乘法逆元
- 数论知识总结(乘法逆元,欧拉函数,线性筛,快速幂,快速乘等)
- poj 3696 The Luckiest number——数论 (快速幂取模,乘法取模模版)
- python 学习1阶段 编写简单乘法运算
- [转][数论]Farey Sequence深入学习
- 求逆元[数论]
- 数论学习笔记
- 数论学习之扩展欧几里德
- 求逆元的简单数论题
- csu1163扩展欧几里得之乘法逆元
- HDU 1576 A/B(数论简单题,求逆元)
- 数论剩余系学习
- Java学习笔记8:用键盘输入1~9,列出相对应的乘法口诀