计算乘法逆元

1.什么是乘法逆元

        众所周知,对于模运算     (a * b) % c = ((a % c) * (b % c)) % c

        除法不能应用其中,也就是说对于    (a / b) % c ,

        我们不能像乘法一样分配,这时我们可以利用逆元来解决这一难题,如果我们取得了b的逆元inv(b),上式可以转换为

        (a * inv(b)) % c = ((a % c) * (inv(b) % c)) % c

        那么什么是乘法逆元呢?

        首先我们依靠不怎么能看懂的"费马小定理"作为我们的理论支持

         可见   a * a^(p-2) = 1  ,1为幺元,所以可以得到数字a的乘法逆元为   a^(p-2)   ,需要注意这个算法有一大缺陷就是p必须是一个质数

2.乘法逆元怎么算

        简单来说乘法逆元的计算就是算出a^(p-2),a的p-2次方,但是一般来说算法题中p是很大的数,所以我们需要采取更高级的算法

        快速幂

        快速幂就是利用位运算快速求解的算法,比如要算2对于模100的乘法逆元,对于下面的代码,需要输入ksm(2,98,100);

//假如需要算a^(p-2),那么a=a,b=(p-2),p=p
inline int ksm(int a, int b, int p) {
	int ans = 1;
	while (b) {
		if (b & 1) { ans= ans* a % p; }
		a = a * a % p;
		b >>=1;
	}
	return ans;
}

        如果只是求次方问题的话可以去掉%和p,如求2的3次方,ksm(2,3);

inline int ksm(int a, int b) {
	int ans = 1;
	while (b) {
		if (b & 1)ans = ans * a;
		a = a * a;
		b >>= 1;
	}
	return ans;
}

        结束.

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论
26在模意义下的乘法元是指另一个整数x,使得26与模数取模后相乘等于1,即 $26x \equiv 1 \pmod{m}$,其中m为模数。 有多种方法可以计算乘法元,下面介绍两种常见的方法。 ### 扩展欧几里得算法 扩展欧几里得算法可以求解形如 $ax + by = \gcd(a,b)$ 的线性方程,其中a和b为任意整数。当a和模数m互质时,即 $\gcd(a,m) = 1$ 时,可以通过扩展欧几里得算法求出a在模m意义下的乘法元。 具体做法是,首先用欧几里得算法求出$\gcd(a,m)$,同时记录下每一步中的商和余数,然后从最后一步开始推,得到 $ax + my = \gcd(a,m)$ 的一组解,其中x即为a在模m意义下的乘法元。 以下是Python代码实现: ```python def ext_euclid(a, b): if b == 0: return a, 1, 0 else: d, x, y = ext_euclid(b, a % b) return d, y, x - (a // b) * y def mod_inv(a, m): d, x, y = ext_euclid(a, m) return x % m if d == 1 else None # 计算26在模37意义下的乘法元 inv = mod_inv(26, 37) print(inv) # 输出16 ``` ### 快速幂算法 当模数m为素数时,可以用费马小定理来计算乘法元。具体来说,如果a和m互质,那么 $a^{m-1} \equiv 1 \pmod{m}$,因此 $a \cdot a^{m-2} \equiv 1 \pmod{m}$,即$a^{m-2}$是a在模m意义下的乘法元。 对于一般的模数,可以先将模数分解为若干个不同的素数的乘积,然后对每个素数分别用费马小定理计算乘法元,最后利用中国剩余定理合并结果。但是这种方法比较繁琐,不如直接使用扩展欧几里得算法。 以下是Python代码实现: ```python def mod_pow(a, b, m): # 快速幂算法,计算a^b mod m res = 1 while b > 0: if b % 2 == 1: res = (res * a) % m a = (a * a) % m b //= 2 return res def mod_inv(a, m): # 计算a在模m意义下的乘法元 return mod_pow(a, m - 2, m) # 计算26在模37意义下的乘法元 inv = mod_inv(26, 37) print(inv) # 输出16 ``` 以上两种方法都可以求解26在模37意义下的乘法元,结果均为16。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

唏嘘南溪

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值