Euclid算法

Euclid算法

Euclid除法

定理：（Elide除法）设a，b是两个整数，其中b>0,则存在唯一的整数q，r使得
$$a=qb+r0\le r<b$$

gcd

最大公因数，记a，b的最大公因数为$gcd(a,b)$或$(a,b)$。

定理：设a,b,r是三个不全为零的整数。如果
$$a=qb+r0\le r<b$$
其中q是整数，则$(a,b)=(b,r)$。

证明：
令$d=(a,b),d^{\prime }=(b,r)$,
则$d|a,d|b,d|a-qb=r$，
所以$d$是$d^{\prime }$的公因数即$d\le d^{\prime }$.
又$d^{\prime }|b,d^{\prime }|r,d^{\prime }|qb+r=a$,
所以$d^{\prime }$是$d$的公因数，所以$d^{\prime }\le d$,故$d=d^{\prime }$。证毕。

def gcd(a, b):
    while b != 0:
        r = a % b
        a = b
        b = r
    return a

最大公因数为$(a,b)=1$,表明$a,b$互素。

ext_gcd

扩展Euclid十分重要，可以用于求逆元。

给定$a,b$,则存在唯一的$s,t$使得$sa+tb=(a,b)$，若$(a,b)=1$，则$sa=1modb$,$s是a的逆元$.

def ext_gcd(a, b):
    s1, t1, r1 = 1, 0, a
    s2, t2, r2 = 0, 1, b
    while r2 != 0:
        q = r1 // r2
        temp1, temp2, temp3 = s1 - q * s2, t1 - q * t2, r1 - q * r2
        s1, t1, r1 = s2, t2, r2
        s2, t2, r2 = temp1, temp2, temp3
    return s1, t1, r1

invert

求$amodp的逆元$。

def invert(a, p):
    s1, t1, r1 = ext_gcd(a, p)
    if r1 != 1:
        print(str(s1) + "*" + str(a) + "+" + str(t1) + "*" + str(p) + "=" + str(r1))
        return None
    else:
        return s1 % p

最小公倍数

lcm

定理，$lcm(a,b)\ast gcd(a,b)=a\ast b$

def lcm(a, b):
    return a * b // gcd(a, b)