费马小定理和欧拉定理作业

CINTA 作业三

---

1. 设 $p=23$ 和 $a=3$，使用费马小定理计算 $a^{2019}\mathrm{mod}p$ ？
   
   由题意可知：
   $a^{2019}\mathrm{mod}p⟺3^{22\ast 91+17}\equiv 3^{17}(\mathrm{mod}23)$
   
   又有：$3^{17}\equiv 3^{-5}(\mathrm{mod}23)$
   
   设 $c=3$ 则有：$c{c}^{-1}\equiv 1(mod23)$
   
   由扩展欧几里得定理得：$c^{-1}=8$
   
   又因为 $gcd(8,23)=1$ ，可以得到：
   原式为：$8^5\mathrm{mod}23=16$
   
   

2. 使用费马小定理求解同余方程 $x^{50}\equiv 2(\mathrm{mod}17)$。
   
   由题意得：
   $x^{50}\mathrm{mod}17⟺x^{3\ast 16+2}\equiv x^2(\mathrm{mod}17)$
   
   所以原式为：$x^2\equiv 2(\mathrm{mod}17)$
   
   那么：$x^2-2=k\ast 17(k\in Z)$
   
   所以当k为2时：$x=\pm 6$
   
   

3. 请证明 13 整除 $2^{70}+3^{70}$
   
   要证：$13|(2^{70}+3^{70})$
   
   只需证：$2^{70}+3^{70}\equiv 0(\mathrm{mod}13)$
   
   因为 $gcd(2,13)=gcd(3,13)=1$：
   
   所以原式可以等价为：$(2^{70}\mathrm{mod}13)+(3^{70}\mathrm{mod}13)$
   
   再由费马小定理得原式为：$(2^{-2}\mathrm{mod}13)+(3^{-2}\mathrm{mod}13)$
   
   由拓展欧几里得定理得 2 和 3 的乘法逆元：分别为 -6 和 -4
   
   那么原式又变作：$(36\mathrm{mod}13)+(16\mathrm{mod}13)$
   
   由同模数加法得原式为：$52\mathrm{mod}13=0$
   
   综上所述： 13 整除 $2^{70}+3^{70}$
   
   

4. 使用欧拉定理计算 $2^{100000}\mathrm{mod}55$
   
   由原式可知：
   $p=55,p$ 不为素数
   
   但 $p=5\ast 11$ 二者均为素数
   
   那么，$\varphi (5)=5-1=4$ ，而 $\varphi (11)=11-1=10$
   
   再有 $gcd(5,11)=1$
   
   由欧拉 $Phi$ 公式：
   $\varphi (55)=\varphi (5)\ast \varphi (11)=40$
   
   因为 $gcd(2,55)=1$ 得：$2^{100000}\equiv 2^{100000\mathrm{mod}\varphi (55)}(\mathrm{mod}55)$
   
   所以原式为：$2^{100000}\equiv 2^0(\mathrm{mod}55)=1$
   
   

5. 手动计算 $7^{1000}$ 的最后两个数位等于什么？
   
   由题意知：所求数为 $7^{1000}\mathrm{mod}100$
   
   据欧拉定理：$\varphi (100)=\varphi (4)\ast \varphi (25)=40$
   
   因为 $gcd(7,100)=1$ 得：$7^{1000}\equiv 7^{1000\mathrm{mod}\varphi (100)}(\mathrm{mod}100)$
   
   所以结果为：$7^{1000}\mathrm{mod}=1$ ，最后两个数位为： 01
   
   

6. 编写 Python 语言程序完成欧拉 $Phi$ 函数的计算，即输入正整数 $n$，计算并返回 $\varphi (n)$。

import math


# 朴素的欧拉函数

def classic_Euler(n):
    num = 1
    for i in range(2, n):
        if math.gcd(i, n) == 1:
            num += 1
    return num


# 优化欧拉函数

def obb_judge(n):
    li = []
    for num in range(2, n + 1):
        flag = 0
        for i in range(2, int(math.sqrt(num)) + 1):
            if num % i == 0:
                flag = 1
            # print('{} % {} == {}'.format(num, i, num % i))
        if flag == 0:
            li.append(num)
    return li


def new_Euler(n):
    p = obb_judge(n)
    if n in p:
        return n - 1
    else:
        ret = 0 + n
        for i in range(0, len(p)):
            if ret % p[i] == 0:
                n -= n // p[i]
                # print('n: ',n)
            while ret % p[i] == 0:
                ret //= p[i]
                # print('ret: ',ret)
        return n

2. 设 $p$ 是素数，计算 $(p-1)!\mathrm{mod}p$，并找出规律，写成定义，并给出证明。
   

# 测试计算使用代码
import math


def obb_judge(n):
    li = []
    for num in range(2, n + 1):
        flag = 0
        for i in range(2, int(math.sqrt(num)) + 1):
            if num % i == 0:
                flag = 1
            # print('{} % {} == {}'.format(num, i, num % i))
        if flag == 0:
            li.append(num)
    return li


def factorial(m):
    re = 1
    for i in range(1, m + 1):
        re *= i
    return re


def fac_to_mode(p):
    if p in obb_judge(p):
        return factorial(p - 1) % p
    else:
        return -1


if __name__ == '__main__':
    m = eval(input('p: '))
    for j in range(2, m):
        if fac_to_mode(j) != -1:
            print('#{}: {}'.format(j, fac_to_mode(j)))

【补充】这里发现一个令人疑惑的地方，直接使用 in range (2,int(math.sqrt(n))+1) 的条件会误判 5、7、11等等数为因子的倍数为素数，因此这里使用上一段代码的判定方式。

# 运行结果：
p: 128
#2: 1
#3: 2
#5: 4
#7: 6
#11: 10
#13: 12
#17: 16
#19: 18
#23: 22
...
#103: 102
#107: 106
#109: 108
#113: 112
#127: 126

这里我们发现对 $(p-1)!\mathrm{mod}p=p-1$ 于是原问题变为：求证：当 $p$ 为素数时， $(p-1)!\equiv (p-1)(\mathrm{mod}p)$ 且只当 $p$ 为素数时（威尔逊定理）

证明：
(1) $p$ 为素数时，总有 $(p-1)!\equiv (p-1)(\mathrm{mod}p)$

$(p-1)!\equiv (p-1)(\mathrm{mod}p)$ 与 $(p-1)!\equiv -1(\mathrm{mod}p)$ 等同
对于 $p=2$ 此结论总成立
对于 $p=3$ 此结论也显然成立

而其余质数 $p$ 满足 $p\mathrm{mod}2=1,p\ge 5$ ：
假设一数 $1<a\le (p-1)$ （因为 $a$ 总是在此区间内存在 [$p$ 为素数]）
那么： a ∈ M ( M = { x ∣ 2 ≤ x ≤ ( p − 2 ) } ) a\in M (M=\{x|2\leq x\leq (p-2)\}) a∈M(M={x∣2≤x≤(p−2)})
再设： N = { y ∣ y = x ∗ a , x ∈ M } N=\{y|y=x*a,x\in M\} N={y∣y=x∗a,x∈M}
那么显然可知：由 $M\to N$ 是双射的（当 $a$ 是常数时）
设方程 $x\ast a\equiv 1(\mathrm{mod}p)$ ，由以上结论可知：
我们只能找到一个 $x$ 与一个 $a$ 一一对应此方程：$x\ast a\equiv 1(\mathrm{mod}p)$
所以有以下所有方程：
$1\ast a_1\mathrm{mod}p=1$
$2\ast a_2\mathrm{mod}p=1$
…
$(p-2)\ast a_{p-2}\mathrm{mod}p=1$
又因 $p$ 为素数：
所以总有：$gcd(a_i,p)=1$
那么对于上述方程都有：$n\ast a_n\equiv n(\mathrm{mod}p)$
对于所有 $n$ 都有 $\gcd (n,p)=1$
所以由乘法原理：$(p-1)!\equiv (p-1)(\mathrm{mod}p)$ 成立
原式得证

(2) $(p-1)!\equiv (p-1)(\mathrm{mod}p)$ 成立时，$p$ 为素数
假设：$p$ 不为素数，那么必然有 $\exists a,b\le p$
那么：原式为 $(p-1)!=1\ast 2\ast 3\ast ...\ast a\ast ...\ast b\ast ...\ast (p-1)\mathrm{mod}p$
因此其中总存在：$\gcd (a\ast b,p)\ne 1$
所以总有：$(p-1)!\equiv 0(\mathrm{mod}p)$
与原式等同的 $(p-1)!\equiv -1(\mathrm{mod}p)$ 不符
所以 $p$ 必然为素数