(ACM数论)求乘法逆元的各种姿势~

首先,通过下面的式子来看看什么是乘法逆元~

x * n % P = 1,其中x和P为已知且互素,n未知(比如在 2 * n % 7 = 1 这个式子里,n就是乘法逆元)

弄懂什么是乘法逆元,来看看有什么姿势可以把它求出来吧~

姿势1.暴力(此姿势不要求P为素数)

没有什么问题是一个暴力解决不了的,如果有,那就两个(手动滑稽)

代码如下:

int n;
for(int i = 1 ;i <= P ;i++)
{
    if(x * i % P == 1)
    {
        n = i;
        break;
    }
}

虽然暴力简单粗暴,但是它的时间复杂度为O(P),只能在P较小的情况使用,局限性有点大。

姿势2.快速幂(此姿势要求P为素数)

这个姿势是由费马巨巨的费马小定理(x^(P-1) ≡ 1 mod P)得来的:

x^(P-1) ≡ 1 mod P 即 x^(P-1) % P = 1

已知 x * n % P = 1

所以有 x * n % P = x^(P-1) % P 即 n % P = x^(P-2) % P

到这里可以看出,逆元有无数个:n,n+P,n+2*P,…

一般情况下我们只要求那个最小的正整数,这样问题就转化为求 x^(P-2) % P

代码就不贴了,在网上找个有求模快速幂的模版,把x,P-2,P代进去就可以了

这个姿势的时间复杂度为O(log P),是个不错的姿势(手动滑稽)

姿势3.扩展欧几里德算法(此姿势不要求P为素数)

首先,我们将式子变形一下~

x * n % P = 1 等价于 n * x + k * P = 1,其中k为整数

因为 x 和 P 互素,所以 gcd (x ,P) = 1(最大公约数为1)

所以有 n * x + k * P = gcd (x ,P)

这样问题就转化为,求二元一次方程 n * x + k * P = gcd (x ,P) 的解

代码同样不贴了, 在网上找个扩展欧几里德的模版,把 x , P 代进去。

不过要注意,算出来的n有可能是负数,也有可能不是最小的。若需要得到最小的正整数,加上如下代码即可:

n = n % P;//处理不是最小的情况
while(n < 0)//处理负数的情况
{
    n =( n + P ) % P;
}

到这里可以得出结论:逆元的集合为{…,n - 2 * P,n - P,n , n + P,n + 2 * P,… }

这个姿势的时间复杂度为O(log P),是个不错的姿势(手动滑稽)

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值