扩展欧几里得算法、求逆元


def exgcd(a, b):
    if b == 0:
        return a, 1, 0
    g, x2, y2 = exgcd(b, a % b)

    x1, y1 = y2, x2 - (a // b) * y2

    return g, x1, y1


# ax+by=m

def f(a, b, m):
    g, x1, y1 = exgcd(a, b)

    if m % g != 0:
        return None, None, None
    x0, y0 = x1 * m // g, y1 * m // g
    return g, x0, y0


def Inv(a, n):
    g, x, y = f(a, n, 1)
    if x is None:
        return None
    else:
        return (x % n + n) % n

n,m=map(int,input().split())
print(Inv(n,m))

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值