夜深人静写算法(三十二)- 费马小定理

文章目录

  • 一、前言
  • 二、费马小定理
    • 1、费马小定理定义
    • 2、费马小定理证明
  • 三、素数判定和伪素数
    • 1、素数判定
    • 2、伪素数
  • 四、费马小定理的应用
    • 1、二分快速幂降幂
    • 2、模 p 逆元
    • 3、Rabin-Miller 素数判定
  • 五、费马小定理相关题集整理

一、前言

  今天要讲的内容,是数论中一个非常著名的定理 —— 费马小定理。阅读之前请确保对 二分快速幂 和 欧拉函数 已经有了一定的了解。本文的内容较短,但是作为 R a b i n   M i l l e r Rabin \ Miller Rabin Miller 大素数判定 有着重要的意义,所以打算单独拿出一个章节来讲,内容较为简单,适合打算放弃学习数论的小伙伴重新找回学习的动力和勇气。
  数论不像其它一些有趣的算法,可能对你没有很强的学习动力,但是所有数论的内容都是成体系的,任何一个公式的推导,可能都需要很强的前置知识,那么我们就把这些前置知识都学完,相信到那时候,任何一个数论题在你面前都易如反掌。在这里插入图片描述

二、费马小定理

1、费马小定理定义

【定义1】对于任意素数 p p p,和正整数 a a a,且 a a a 不是 p p p 的倍数,则: a p − 1 ≡ 1 ( m o d   p ) a^{p-1} \equiv 1 (mod \ p) ap11(mod p)

2、费马小定理证明

  • 我们在学习欧拉函数的时候,已经知道了欧拉定理,如下:
  • a ϕ ( n ) ≡ 1 ( m o d   n ) a^{\phi(n)} \equiv 1 (mod \ n) aϕ(n)1(mod n)
  • n n n 为 素数 p p p 时,它的欧拉函数为 ϕ ( p ) = p − 1 \phi(p) = p-1 ϕ(p)=p1,将它带入欧拉定理,得到:
  • a p − 1 ≡ 1 ( m o d   p ) a^{p-1} \equiv 1 (mod \ p) ap11(mod p)
  • 费马小定理,得证。
  • 注:关于欧拉定理的详细证明过程,可以参考以下这一篇文章,有兴趣的读者建议本文阅读完毕再看:
  • 夜深人静写算法(三十一)- 欧拉函数

三、素数判定和伪素数

1、素数判定

  • 我们可以用费马小定理来做什么?
  • 一个比较直观的想法就是:可以随机找几个和 n n n 互素的 a a a,然后对它计算:
  • a n − 1   m o d   n a^{n-1} \ mod \ n an1 mod n
  • 如果结果都为 1,我们就可以认为 n n n 是一个素数。
  • 如果这个结论成立,那么素数判定的时间复杂度就变成了 O ( C l o g 2 n ) O(Clog_2n) O(Clog2n),其中 C C C 为常数,代表找 C C C a a a 来做判定试验, O ( l o g 2 n ) O(log_2n) O(log2n) 则为利用二分快速幂进行判定的时间复杂度。
  • 以上假设成立吗?
  • 答案是 否!

2、伪素数

  • 事实上,费马小定理给出的是关于素数判定的 必要非充分 条件。
  • 如果 p p p 是素数,则 a p − 1 ≡ 1 ( m o d   p ) a^{p-1} \equiv 1 (mod \ p) ap11(mod p);相反,如果 a p − 1 ≡ 1 ( m o d   p ) a^{p-1} \equiv 1 (mod \ p) ap11(mod p),则不能推导出 p p p 是素数。
  • 原因是存在一些数 q q q,对于所有和 q q q 互素的 a a a,都能满足 a q − 1 ≡ 1 ( m o d   q ) a^{q-1} \equiv 1 (mod \ q) aq11(mod q),这样的数,我们称它为伪素数。
  • 第一个伪素数是 341,由 萨鲁斯 在 1819 年提出。

四、费马小定理的应用

1、二分快速幂降幂

【例题1】 给出一个大整数 n ( 1 ≤ n ≤ 1 0 100000 ) n(1 \le n \le 10^{100000}) n(1n10100000),求: 2 n   m o d   1000000007 2^n \ mod \ 1000000007 2n mod 1000000007

  • 任何一个正整数都就可以表示成 n = k x + m n = kx + m n=kx+m 的形式,其中 x x x 为除数, k k k 为商, m m m 为余数,并且 0 ≤ m < x 0 \le m \lt x 0m<x;其中除数 x x x 可以是任意非零整数。
  • 为了公式看起来整洁,我们用 p p p 来代替素数 1000000007 1000000007 1000000007,并且令 x = p − 1 x=p-1 x=p1 那么根据费马小定理,有:
  • 2 p − 1   m o d   p = 1 2^{p-1} \ mod \ p = 1 2p1 mod p=1
  • 原式就可以表示成:
  • 2 n   m o d   p = 2 k x + m   m o d   p = 2 k ( p − 1 ) + m   m o d   p = ( 2 k ( p − 1 )   m o d   p ) ∗ ( 2 m   m o d   p )   m o d   p = 2 m   m o d   p \begin{aligned}2^n \ mod \ p &= 2^{kx+m} \ mod \ p \\ &= 2^{k(p-1)+m} \ mod \ p \\ &= (2^{k(p-1)} \ mod \ p) * (2^{m} \ mod \ p) \ mod \ p \\ &= 2^m \ mod \ p \end{aligned} 2n mod p=2kx+m mod p=2k(p1)+m mod p=(2k(p1) mod p)(2m mod p) mod p=2m mod p
  • 这里的 m = n   m o d   ( p − 1 ) m = n \ mod \ (p-1) m=n mod (p1),可以利用大数取余求解。求得的 m ∈ [ 0 , p − 1 ) m \in [0, p-1) m[0,p1),再利用二分快速幂求解上面的式子即可。

2、模 p 逆元

【例题2】给定素数 p p p 和 正整数 a a a,求满足 a x ≡ 1 ( m o d   p ) ax \equiv 1 (mod \ p) ax1(mod p) 的最小正整数 x,如果不存在返回 -1。

  • 这个问题实际上是求正整数 a a a 在 模 p p p 域上的逆元。
  • 首先,当 a a a p p p 的倍数时, a x ≡ 0 ( m o d   p ) ax \equiv 0 (mod \ p) ax0(mod p),所以一定不存在,直接返回 -1;
  • 否则,根据费马小定理,我们可以知道:
  • a p − 1 ≡ 1 ( m o d   p ) a^{p-1} \equiv 1 (mod \ p) ap11(mod p)
  • 则可以得到:
  • a × a p − 2 ≡ 1 ( m o d   p ) a \times a^{p-2} \equiv 1 (mod \ p) a×ap21(mod p)
  • 对比原式,就可以得到:
  • x = a p − 2   m o d   p x = a^{p-2} \ mod \ p x=ap2 mod p

3、Rabin-Miller 素数判定

  • 对于一个很大的数 n n n(例如十进制表示有 100 100 100 位),如果还是采用试除法进行判定,时间复杂度必定难以承受,目前比较稳定的大素数判定法是 拉宾-米勒( R a b i n   M i l l e r Rabin \ Miller Rabin Miller)素数判定。
  • 判定过程用到了费马小定理,具体算法的推导和应用过程将在后续章节再进行展开讲解。

  • 关于 费马小定理 的内容到这里就暂时结束了,主要还是运用这个定理,对模数是素数的幂式进行降幂操作,如果还有不懂的问题可以留言告诉作者或者添加作者的微信公众号。

  • 本文所有示例代码均可在以下 github 上找到:github.com/WhereIsHeroFrom/Code_Templates

五、费马小定理相关题集整理

题目链接难度解析
HDU 1395 2^x mod n = 1★☆☆☆☆费马小定理 简化版
HDU 4704 Sum★☆☆☆☆【例题1】杨辉三角 + 费马小定理降幂
PKU 3641 Pseudoprime numbers★☆☆☆☆素数判定 + 费马小定理降幂
PKU 1845 Sumdiv★★☆☆☆因子和 + 逆元 + 费马小定理降幂
HDU 4869 Turn the pokers ★★★☆☆组合计数+ 费马小定理
HDU 3307 Description has only two Sentences★★★☆☆公式推导 + 费马小定理
HDU 6755 Fibonacci Sum★★★★☆逆元 + 组合公式 + 费马小定理 + 二分快速幂
HDU 4959 Poor Akagi★★★★★费马小定理 + 卢卡斯数通项公式 + 矩阵二分快速幂

热门文章

暂无图片
编程学习 ·

Java输出数组的内容

Java输出数组的内容_一万个小时-CSDN博客_java打印数组内容1. 输出内容最常见的方式// List<String>类型的列表List<String> list new ArrayList<String>();list.add("First");list.add("Second");list.add("Third");list.ad…
暂无图片
编程学习 ·

母螳螂的“魅惑之术”

在它们对大蝗虫发起进攻的时候&#xff0c;我认认真真地观察了一次&#xff0c;因为它们突然像触电一样浑身痉挛起来&#xff0c;警觉地面对限前这个大家伙&#xff0c;然后放下自己优雅的身段和祈祷的双手&#xff0c;摆出了一个可怕的姿势。我被眼前的一幕吓到了&#xff0c;…
暂无图片
编程学习 ·

疯狂填词 mad_libs 第9章9.9.2

#win7 python3.7.0 import os,reos.chdir(d:\documents\program_language) file1open(.\疯狂填词_d9z9d2_r.txt) file2open(.\疯狂填词_d9z9d2_w.txt,w) words[ADJECTIVE,NOUN,VERB,NOUN] str1file1.read()#方法1 for word in words :word_replaceinput(fEnter a {word} :)str1…
暂无图片
编程学习 ·

HBASE 高可用

为了保证HBASE是高可用的,所依赖的HDFS和zookeeper也要是高可用的. 通过参数hbase.rootdir指定了连接到Hadoop的地址,mycluster表示为Hadoop的集群. HBASE本身的高可用很简单,只要在一个健康的集群其他节点通过命令 hbase-daemon.sh start master启动一个Hmaster进程,这个Hmast…
暂无图片
编程学习 ·

js事件操作语法

一、事件的绑定语法 语法形式1 事件监听 标签对象.addEventListener(click,function(){}); 语法形式2 on语法绑定 标签对象.onclick function(){} on语法是通过 等于赋值绑定的事件处理函数 , 等于赋值本质上执行的是覆盖赋值,后赋值的数据会覆盖之前存储的数据,也就是on…
暂无图片
编程学习 ·

Photoshop插件--晕影动态--选区--脚本开发--PS插件

文章目录1.插件界面2.关键代码2.1 选区2.2 动态晕影3.作者寄语PS是一款栅格图像编辑软件&#xff0c;具有许多强大的功能&#xff0c;本文演示如何通过脚本实现晕影动态和选区相关功能&#xff0c;展示从互联网收集而来的一个小插件&#xff0c;供大家学习交流&#xff0c;请勿…
暂无图片
编程学习 ·

vs LNK1104 无法打开文件“xxx.obj”

写在前面&#xff1a; 向大家推荐两本新书&#xff0c;《深度学习计算机视觉实战》和《学习OpenCV4&#xff1a;基于Python的算法实战》。 《深度学习计算机视觉实战》讲了计算机视觉理论基础&#xff0c;讲了案例项目&#xff0c;讲了模型部署&#xff0c;这些项目学会之后可以…
暂无图片
编程学习 ·

工业元宇宙的定义与实施路线图

工业元宇宙的定义与实施路线图 李正海 1 工业元宇宙 给大家做一个关于工业元宇宙的定义。对于工业&#xff0c;从设计的角度来讲&#xff0c;现在的设计人员已经做到了普遍的三维设计&#xff0c;但是进入元宇宙时代&#xff0c;就不仅仅只是三维设计了&#xff0c;我们的目…
暂无图片
编程学习 ·

【leectode 2022.1.15】完成一半题目

有 N 位扣友参加了微软与力扣举办了「以扣会友」线下活动。主办方提供了 2*N 道题目&#xff0c;整型数组 questions 中每个数字对应了每道题目所涉及的知识点类型。 若每位扣友选择不同的一题&#xff0c;请返回被选的 N 道题目至少包含多少种知识点类型。 示例 1&#xff1a…
暂无图片
编程学习 ·

js 面试题总结

一、js原型与原型链 1. prototype 每个函数都有一个prototype属性&#xff0c;被称为显示原型 2._ _proto_ _ 每个实例对象都会有_ _proto_ _属性,其被称为隐式原型 每一个实例对象的隐式原型_ _proto_ _属性指向自身构造函数的显式原型prototype 3. constructor 每个prot…
暂无图片
编程学习 ·

java练习代码

打印自定义行数的空心菱形练习代码如下 import java.util.Scanner; public class daYinLengXing{public static void main(String[] args) {System.out.println("请输入行数");Scanner myScanner new Scanner(System.in);int g myScanner.nextInt();int num g%2;//…
暂无图片
编程学习 ·

RocketMQ-什么是死信队列?怎么解决

目录 什么是死信队列 死信队列的特征 死信消息的处理 什么是死信队列 当一条消息初次消费失败&#xff0c;消息队列会自动进行消费重试&#xff1b;达到最大重试次数后&#xff0c;若消费依然失败&#xff0c;则表明消费者在正常情况下无法正确地消费该消息&#xff0c;此时…
暂无图片
编程学习 ·

项目 cg day04

第4章 lua、Canal实现广告缓存 学习目标 Lua介绍 Lua语法 输出、变量定义、数据类型、流程控制(if..)、循环操作、函数、表(数组)、模块OpenResty介绍(理解配置) 封装了Nginx&#xff0c;并且提供了Lua扩展&#xff0c;大大提升了Nginx对并发处理的能&#xff0c;10K-1000K Lu…
暂无图片
编程学习 ·

输出三角形

#include <stdio.h> int main() { int i,j; for(i0;i<5;i) { for(j0;j<i;j) { printf("*"); } printf("\n"); } }
暂无图片
编程学习 ·

stm32的BOOTLOADER学习1

序言 最近计划学习stm32的BOOTLOADER学习,把学习过程记录下来 因为现在网上STM32C8T6还是比较贵的,根据我的需求flash空间小一些也可以,所以我决定使用stm32c6t6.这个芯片的空间是32kb的。 #熟悉芯片内部的空间地址 1、flash ROM&#xff1a; 大小32KB&#xff0c;范围&#xf…
暂无图片
编程学习 ·

通过awk和shell来限制IP多次访问之学不会你打死我

学不会你打死我 今天我们用shell脚本&#xff0c;awk工具来分析日志来判断是否存在扫描器来进行破解网站密码——限制访问次数过多的IP地址&#xff0c;通过Iptables来进行限制。代码在末尾 首先我们要先查看日志的格式&#xff0c;分析出我们需要筛选的内容&#xff0c;日志…
暂无图片
编程学习 ·

Python - 如何像程序员一样思考

在为计算机编写程序之前&#xff0c;您必须学会如何像程序员一样思考。学习像程序员一样思考对任何学生都很有价值。以下步骤可帮助任何人学习编码并了解计算机科学的价值——即使他们不打算成为计算机科学家。 顾名思义&#xff0c;Python经常被想要学习编程的人用作第一语言…
暂无图片
编程学习 ·

蓝桥杯python-数字三角形

问题描述 虽然我前后用了三种做法&#xff0c;但是我发现只有“优化思路_1”可以通过蓝桥杯官网中的测评&#xff0c;但是如果用c/c的话&#xff0c;每个都通得过&#xff0c;足以可见python的效率之低&#xff08;但耐不住人家好用啊&#xff08;哭笑&#xff09;&#xff09…