基于Python模拟实现n乘法的逆元计算器

 更新时间:2026年06月29日 08:55:10   作者:金銀銅鐵  
本文介绍了扩展欧几里得算法及其在计算模n乘法逆元中的应用,通过Pygame实现图形化界面,开发了一个模n乘法逆元计算器,支持大数运算,高效便捷

背景

在 [Python] 扩展欧几里得算法 一文中,我们已经接触了扩展欧几里得算法。借助扩展欧几里得算法,我们可以高效地计算模 n 乘法的逆元。本文会介绍相关的内容。

要点

  • 对自然数 a 和正整数 n 而言,ax1(modn) 有解的充分必要条件是gcd(a,n)=1
  • 在扩展欧几里得算法的基础上,我们可以高效地求出 a 的逆元(即上式中的 x
  • 再借助Pygame 实现图形化界面,我们就可以制作一个模 n 乘法的逆元计算器了

对自然数 a 和 正整数 n 而言,是否可以找到一个整数 x 使得 ax1(modn) 成立呢?

由于以下两个等式成立(其中 k是任意整数),我们只需要关心满足0a,x<n 条件的整数a,x 就够了。

  • ax(a+kn)x(modn)
  • axa(x+kn)(modn)

我们可以先用比较小的 n 试一试。

用较小的n做尝试

我们可以用暴力的方法,对较小的 n 做尝试 ⬇️

def show_markdown_result(n):
    print(f"#### $n={n}$ 时")
    print(f"| $a$ | 满足 $ax\\equiv 1 \\pmod{{{n}}}$ 和 $0\\le x\\lt {n}$ 的 $x$ |")
    print("| --- | --- |")
    for a in range(n):
        matched_result = []
        for x in range(n):
            if (a * x - 1) % n == 0:
                matched_result.append(x)
        if matched_result:
            joined_result = ','.join(map(str, matched_result))
            print(f"| ${a}$ | ${joined_result}$ |")
        else:
            print(f"| ${a}$ | (不存在) |")

for n in range(1, 10):
    show_markdown_result(n)
    print("\n")

请将以上代码保存为 show.py,用下方的命令可以运行 show.py

python3 show.py

运行结果中包含了9markdown 表格 ⬇️

n=1时

a满足ax1(mod1)0x<1x
00

n=2时

a满足ax1(mod2)0x<2x
0(不存在)
11

n=3时

a满足 ax1(mod3)0x<3x
0(不存在)
11
22

n=4时

a满足ax1(mod4)0x<4x
0(不存在)
11
2(不存在)
33

n=5时

a满足ax1(mod5)0x<5x
0(不存在)
11
23
32
44

n=6时

a满足ax1(mod6)0x<6x
0(不存在)
11
2(不存在)
3(不存在)
4(不存在)
55

n=7时

a满足ax1(mod7)0x<7x
0(不存在)
11
24
35
42
53
66

n=8时

a满足ax1(mod8)0x<8x
0(不存在)
11
2(不存在)
33
4(不存在)
55
6(不存在)
77

n=9时

a满足ax1(mod9)0x<9x
0(不存在)
11
25
3(不存在)
47
52
6(不存在)
74
88

一些规律

从上一小节的数据来看,似乎有以下规律成立(有待证明或者证伪)

  1. an 互质(即 gcd(a,n)=1)时,似乎总能找到满足ax1(modn)x
  2. an 互质(即 gcd(a,n)≠1)时,似乎总是无法找到满足ax1(modn)x
  3. 当我们限制 x 满足 0x<n 时,能让ax1(modn) 成立的 x 似乎最多只有 1 个。

规律一的证明

第一个规律,我们可以用扩展欧几里得算法来证明。当 a 和 nn互质(即 gcd(a,n)=1)时,借助扩展欧几里得算法,可以找到满足ax0+ny0=1 的整数 x0,y0。于是 ax01=ny0,所以ax01(modn)。令 x=x0 mod n,就找到了符合要求的 x

规律二的证明

第二个规律,可以用反证法证明。当 a 和 n互质(即 gcd(a,n)≠1)时,假设可以找到某个 x使得 ax1(modn)x 成立,那么存在某个整数 k 使得 ax+kn=1 成立。以下两个等式显然成立

  • gcd(a,n)∣ax
  • gcd(a,n)∣kn

那么 gcd(a,n)(ax+kn) 成立,即gcd(a,n)1 成立,这与 gcd(a,n)≠1 矛盾。第二个规律得证。有了规律一和规律二,我们就知道 ax1(modn) 有解的充分必要条件是 gcd(a,n)=1

规律三的证明

第三个规律,也可以证明。假设可以找到x1,x2 分别满足

  • ax11(modn)
  • ax21(modn)

我们在第一个等式两边乘以 x2,会得到

x2​ax1​≡x2​(modn)

我们在第二个等式两边乘以x1,会得到

x1​ax2​≡x1​(modn)

而模 n 的乘法满足 交换律和结合律(这里就不证明了),所以

  • x2​ax1​≡ax2​x1​(modn)
  • ax2​x1​≡ax1​x2​(modn)
  • ax1​x2​≡x1​ax2​(modn)

于是x1x2(modn)

因此,当我们限制 x 满足 0x<n 时,能让ax1(modn) 成立的 x 最多只有 1

用程序来处理

核心逻辑

[Python] 扩展欧几里得算法 一文中已经提供了扩展欧几里得算法的代码。我们在它的基础上,再加上计算逆元的方法即可

def extended_euclidean(a, b):
    if (a, b) == (0, 0):
        raise ValueError("a and b cannot be both 0")
    if b == 0:
        return (1, 0, a)
    x, y, g = extended_euclidean(b, a % b)
    return (y, x - a // b * y, g)

def calc_inverse_element(a, n):
    x, _, gcd = extended_euclidean(a, n)
    if gcd == 1:
        return x % n
    return None

添加 GUI 之后的完整程序

 我完成了和Pygame 相关的代码。完整的 Python 程序如下 

import pygame

def extended_euclidean(a, b):
    if (a, b) == (0, 0):
        raise ValueError("a and b cannot be both 0")
    if b == 0:
        return (1, 0, a)
    x, y, g = extended_euclidean(b, a % b)
    return (y, x - a // b * y, g)

def calc_inverse_element(a, n):
    x, _, gcd = extended_euclidean(a, n)
    if gcd == 1:
        return x % n
    return None

# ===================== Pygame 初始化 =====================
pygame.init()

WIDTH, HEIGHT = 700, 480
screen = pygame.display.set_mode((WIDTH, HEIGHT))
pygame.display.set_caption("Inverse Element Calculator")

# 颜色配置
WHITE = (255, 255, 255)
BLACK = (0, 0, 0)
GRAY = (220, 220, 220)
BLUE = (50, 100, 200)
RED = (200, 50, 50)
GREEN = (50, 180, 50)

# 字体
font_input = pygame.font.SysFont("Arial", 36)
font_text = pygame.font.SysFont("Arial", 32)
font_result = pygame.font.SysFont("Arial", 38)
font_equation = pygame.font.SysFont("Arial", 32)
font_btn = pygame.font.SysFont("Arial", 34)

# ===================== 布局:全部水平居中,对称美观 =====================
# 输入框组居中
INPUT_X = WIDTH // 2 - 110  # 输入框水平居中
INPUT_WIDTH = 220
INPUT_HEIGHT = 50

# 输入a
input_a_rect = pygame.Rect(INPUT_X, 110, INPUT_WIDTH, INPUT_HEIGHT)
input_a_text = ""
a_active = False

# 输入n
input_n_rect = pygame.Rect(INPUT_X, 180, INPUT_WIDTH, INPUT_HEIGHT)
input_n_text = ""
n_active = False

# 按钮居中
btn_rect = pygame.Rect(WIDTH // 2 - 150, 260, 300, 60)
btn_text = "Calculate"

# 结果
result_text = ""
equation_text = ""
result_color = BLACK

# ===================== 主循环 =====================
running = True
while running:
    screen.fill(WHITE)

    # -------------------- 文字与输入框(居中对齐)--------------------
    # Input a 标签(与输入框垂直居中)
    label_a = font_text.render("Input a:", True, BLACK)
    screen.blit(label_a, (INPUT_X - 120, 118))
    
    # Input n 标签(与输入框垂直居中)
    label_n = font_text.render("Input n:", True, BLACK)
    screen.blit(label_n, (INPUT_X - 120, 188))

    # 绘制输入框a
    pygame.draw.rect(screen, BLUE if a_active else GRAY, input_a_rect, 2)
    screen.blit(font_input.render(input_a_text, True, BLACK), (input_a_rect.x+10, input_a_rect.y+5))

    # 绘制输入框n
    pygame.draw.rect(screen, BLUE if n_active else GRAY, input_n_rect, 2)
    screen.blit(font_input.render(input_n_text, True, BLACK), (input_n_rect.x+10, input_n_rect.y+5))

    # 绘制按钮(居中)
    pygame.draw.rect(screen, BLUE, btn_rect, border_radius=8)
    btn_surface = font_btn.render(btn_text, True, WHITE)
    screen.blit(btn_surface, (btn_rect.centerx - btn_surface.get_width()//2, btn_rect.centery - btn_surface.get_height()//2))

    # 结果文字(居中显示,更美观)
    result_surface = font_result.render(result_text, True, result_color)
    screen.blit(result_surface, (WIDTH // 2 - result_surface.get_width()//2, 350))

    # 公式提示(底部居中)
    formula = font_equation.render(equation_text, True, BLACK)
    screen.blit(formula, (WIDTH // 2 - formula.get_width()//2, 420))

    # ===================== 事件处理 =====================
    for event in pygame.event.get():
        if event.type == pygame.QUIT:
            running = False

        # 鼠标点击
        if event.type == pygame.MOUSEBUTTONDOWN:
            if input_a_rect.collidepoint(event.pos):
                a_active = True
                n_active = False
            elif input_n_rect.collidepoint(event.pos):
                n_active = True
                a_active = False
            elif btn_rect.collidepoint(event.pos):
                result_text = ""
                try:
                    a = int(input_a_text)
                    n = int(input_n_text)
                    if a <= 0 or n <= 0:
                        result_text = "ERROR: Please input positive integers!"
                        result_color = RED
                    else:
                        inv = calc_inverse_element(a, n)
                        if inv is not None:
                            result_text = f"x: {inv}"
                            result_color = GREEN
                            equation_text = f"{a} * {inv} ≡ 1 (mod {n})"
                        else:
                            result_text = f"ERROR: {a} and {n} are not relatively prime, no inverse exists!"
                            result_color = RED
                except ValueError:
                    result_text = "ERROR: Please input valid integers!"
                    result_color = RED
            else:
                a_active = False
                n_active = False

        # 键盘输入
        if event.type == pygame.KEYDOWN:
            if a_active:
                if event.key == pygame.K_BACKSPACE:
                    input_a_text = input_a_text[:-1]
                elif event.unicode.isdigit():
                    input_a_text += event.unicode
            if n_active:
                if event.key == pygame.K_BACKSPACE:
                    input_n_text = input_n_text[:-1]
                elif event.unicode.isdigit():
                    input_n_text += event.unicode

    pygame.display.flip()

pygame.quit()

请将完整的程序保存为 inverse_element_calculator.py。使用如下的命令可以运行 inverse_element_calculator.py

python3 inverse_element_calculator.py

初始界面如下图所示 

此时需要用户输入 a 和 n。我输入了 3423,然后点击 Calculate 按钮,就可以计算 a 对应的逆元了

由于扩展欧几里得算法很高效,所以即使输入比较大 aaannn,也可以立刻算出结果,示例效果如下  (a=63245986,n=102334155)

到此这篇关于基于Python模拟实现n乘法的逆元计算器的文章就介绍到这了,更多相关Python逆元计算器内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家!

相关文章

  • python如何删除列为空的行

    python如何删除列为空的行

    在本篇文章里小编给大家整理的是关于python删除列为空的行方法,对此有需要的朋友们可以学习下。
    2020-07-07
  • Python实现高斯函数的三维显示方法

    Python实现高斯函数的三维显示方法

    今天小编就为大家分享一篇Python实现高斯函数的三维显示方法,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2018-12-12
  • Python使用tkinter模块实现推箱子游戏

    Python使用tkinter模块实现推箱子游戏

    这篇文章主要介绍了Python使用tkinter模块实现推箱子游戏,主要分享两点,第一就是这个程序的实现过程,第二点就是我在编写过程中的一些思考。本文通过实例代码给大家介绍的非常详细,需要的朋友参看下吧
    2019-10-10
  • python自动化发送邮件实例讲解

    python自动化发送邮件实例讲解

    在本篇文章里小编给大家分享了一篇关于python自动化发送邮件实例讲解内容,有兴趣的朋友们可以学习参考下。
    2021-01-01
  • Python调用C++程序的方法详解

    Python调用C++程序的方法详解

    这篇文章主要介绍了Python调用C++程序的方法,文中通过示例代码介绍的详细,相信对大家具有一定的参考借鉴价值,需要的朋友们下面来一起看看吧。
    2017-01-01
  • Python Flask框架开发之运用SocketIO实现WebSSH方法详解

    Python Flask框架开发之运用SocketIO实现WebSSH方法详解

    Socket.IO本是一个面向实时web应用的JavaScript库,现在已成为拥有众多语言支持的Web即时通讯应用的框架。这篇文章主要介绍了Python 运用SocketIO实现WebSSH方法
    2022-10-10
  • python语法之通过value找key问题

    python语法之通过value找key问题

    这篇文章主要介绍了python语法之通过value找key问题,具有很好的参考价值,希望对大家有所帮助。如有错误或未考虑完全的地方,望不吝赐教
    2022-11-11
  • Python编程之微信推送模板消息功能示例

    Python编程之微信推送模板消息功能示例

    这篇文章主要介绍了Python编程之微信推送模板消息功能,结合实例形式分析了Python微信推送消息接口的调用相关操作技巧,需要的朋友可以参考下
    2017-08-08
  • python3射线法判断点是否在多边形内

    python3射线法判断点是否在多边形内

    这篇文章主要为大家详细介绍了python3射线法判断点是否在多边形内,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2019-06-06
  • pandas去除重复值的实战

    pandas去除重复值的实战

    本文主要介绍了pandas去除重复值的实战,文中通过示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2021-09-09

最新评论