学习:费马小定理 & 欧拉定理
费马小定理
描述
若\(p\)为素数,\(a\in Z\),则有\(a^p\equiv a\pmod p\)。如果\(p\nmid a\),则有\(a^{p-1}\equiv 1\pmod p\)。
证明
费马小定理的证法有很多,此处介绍3种
证法一
摘自:《初等数论》 冯志刚 著,有改动
此处用归纳法证明。
当\(a=1\)时,原命题显然成立。
设\(a=n\)时命题成立,即\(n^p\equiv n\pmod p\),故\(n^p-n\equiv 0\pmod p\)。
考虑二项式系数\(C^k_p=\frac{p!}{k!(p-k)!}\),若\(k\)不为\(p\)或\(0\),由于分子有质数\(p\),但分母不含\(p\),故分子的\(p\)能保留,不被约分而除去,即\(C^k_p\)恒为\(p\)的倍数。
故有
\]
所以对于任意\(a\in N_+\),有\(p\mid (a^p-a)\),故有\(a^p\equiv a\pmod p\)。
证法二
摘自: Wikipedia,有改动
与证法一类似,先证明当\(p\)是质数时\(C^k_p\)恒为\(p\)的倍数。
然后就可以得到
\]
当\(b=a-1\)时,可得\(a^p\equiv a\pmod p\)
证法三
摘自:Matrix67的博客
首先我们证明这样一个结论:如果\(p\)是一个素数的话,那么对任意一个小于\(p\)的正整数\(a\),\(a, 2a, 3a,\cdots, (p-1)a\)除以\(p\)的余数正好是一个\(1\)到\(p-1\)的排列。
假如结论不成立的话,那么就是说有两个小于\(p\)的正整数\(m\)和\(n\)使得\(na\)和\(ma\)除以\(p\)的余数相同。不妨假设\(n>m\),则\(p\)可以整除\(a(n-m)\)。由于\(p\)是素数,那么\(a\)和\(n-m\)中至少有一个含有因子\(p\)。这显然是不可能的,因为\(a\)和\(n-m\)都比\(p\)小。
故:
\]
也即:
\]
两边同时除以\((p-1)!\),就得到了我们的最终结论:
\]
应用
例题:LJJ算数
题目描述
LJJ刚上完了一节课!这节课是数学课!他知道了加减属于一级运算,乘除属于二级运算,幂则属于三级运算,而幂的优先级>乘除的优先级>加减的优先级(这是几年级的数学课)。但是,从上一套试卷+上一题中,我们知道了LJJ是一个总是突发奇想并且智商不够的人(也就是说他又想出一个问题给你咯)。他发明了一种四级运算,我们姑且用符号#来表示(找不到别的符号了)。我们知道\(a\times b=a+a+a+\cdots+a\)(加b次),\(a^b=a\times a\times a\times a\times \cdots \times a\)(乘b次),则\(a\#b=((((a^a)^a)^a)^ \cdots )^a\)(进行幂运算b次),自然,#的优先级比幂的优先级高。那么,LJJ就请你来帮他求\(a\#b \mod {1000000007}\)咯。
输入格式:
输入仅1行,即\(a,b\)。
输出格式:
输出仅1行,即\(a\#b \mod 1000000007\)。
输入输出样例
输入样例#1:
3 5
输出样例#1:
968803245
说明
首先说明,样例答案不mod其实是4.4342648824303776994824963061915e+38(来自出题人的恶意)
然后,数据范围:
对于20%的数据,\(a\le1000,b\le1000\)
对于50%的数据,\(a\le 10^{16},b\le 10000\)
对于100%的数据,\(a\le 10^{16},b\le 10^{16}\)
题解
由费马小定理,得
\]
由此递归式,易知
\]
于是用两次快速幂+取模即可。
代码
#include <cstdio>
typedef long long LL;
inline LL pow_mod(LL a, LL b, LL mod)//快速幂模板
{
a %= mod;
LL ans = 1;
for(; b; b >>= 1,a *= a,a %= mod)
if(b & 1)
ans = ans*a%mod;
return ans%mod;
}
const LL mod = 1e9+7;
int main()
{
LL a,b;
scanf("%lld%lld", &a, &b);
printf("%lld", pow_mod(a, pow_mod(a, b-1, mod-1), mod));
return 0;
}
求数论倒数
数论倒数也称为模倒数,或者模逆元。若\(ab\equiv 1\pmod p\),则称\(b\)是\(a\)的数论倒数,亦可写作\(a^{-1}\equiv b\pmod p\)。
若\(p\)是素数,则\(a^{p-1}=a\times a^{p-2}\equiv 1\pmod p\),故\(a^{-1}\equiv a^{p-2}\pmod p\)。
但要注意前提条件是\(p\)是素数。
欧拉定理
同余类、剩余系、欧拉函数
对于任意整数,我们将它模\(n\),结果必定为\(0\)到\(n-1\)的整数之一。那么,我们可以把所有模\(n\)后同余的数看成一个集合,那么我们就把整数分成了\(n\)个集合,我们把这样的模\(n\)同余的所有整数组成的集合(即上述中的任意一个集合)称为同余类,标记为\({\overline {a}}_{n}\),假若从上下文知道模\(n\),则也可标记为\([a]\)。
那么,模\(n\)的同余类有\(n\)个,我们如何表示它们呢?和并查集的思想一样:任取集合中的一个数作为集合的代表元素。我们所取的那个元素就被称为该同余类的代表数。显然,同余类中的每个元素都可以作为该同余类的代表数。
剩余系指的是模\(n\)同余类的代表数集合。一个完全剩余系(完系)指的是模\(n\)的全部同余类的代表数的集合。例如,模\(3\)有三个同余类\([0],[1],[2]\),其完系可以是\(\{9,12+1,15+2\}\)。如果该集合是由每个同余类的最小非负整数所组成,亦即$ {0,1,2,...,n-1}\(,则称该集合为模\)n\(的**非负最小完全剩余系**。模\)n\(完整余数系统中,与模\)n\(互质的代表数所构成的集合,称为模\)n$的简化剩余系(简系)。
与\(m\)互素的所有模\(m\)的同余类的个数记为\(\varphi(n)\),通常称为欧拉函数。显然,\(\varphi(n)\)等于\(1,2,\cdots,m\)中与\(m\)互素的个数。
若\((a,m)=1\),而\(a_1,a_2,\cdots,a_{\varphi(m)}\)构成模\(m\)的一个简系,则我们可以证明\(aa_1,aa_2,\cdots,aa_{\varphi(m)}\)也是模\(m\)的简系。
欧拉定理
若\((a,n)=1\),则\(a^{\varphi(n)}\equiv 1\pmod n\)。(这里\(n\in N_+,a\in Z\))
证明
设\(a_1,a_2,\cdots,a_{\varphi(n)}\)是模\(n\)的简系,则由\((a,m)=1\),可知\(aa_1,aa_2,\cdots,aa_{\varphi(n)}\)也是模\(n\)的简系。因此,$$\prod^{\varphi(n)}{i=1}a_i\equiv \prod^{\varphi(n)}{i=1}aa_i\pmod n,$$即$$m\mid (\prod{\varphi(n)}_{i=1}a_i)(a{\varphi(n)}-1).$$又因为\((a_i,m)=1\),故$$a^{\varphi(n)}\equiv 1\pmod n.$$
就先写到这儿吧,欧拉定理的应用以后再说。
参考资料
Wikipedia
主要参考了:
《初等数论》 冯志刚 著
主要参考了:
2.2 同余类与剩余系
2.3 费马小定理与欧拉定理
Matrix67的博客——数论部分第一节:素数与素性测试
学习:费马小定理 & 欧拉定理的更多相关文章
- 费马小定理&欧拉定理
在p是素数的情况下,对任意整数x都有xp≡x(mod p).这个定理被称作费马小定理其中如果x无法被p整除,我们有xp-1≡1(mod p).利用这条性质,在p是素数的情况下,就很容易求出一个数的逆元 ...
- 【初等数论】费马小定理&欧拉定理&扩展欧拉定理(暂不含证明)
(不会证明--以后再说) 费马小定理 对于任意\(a,p \in N_+\),有 \(a^{p-1} \equiv 1\pmod {p}\) 推论: \(a^{-1} \equiv a^{p-2} \ ...
- 【BZOJ】3398: [Usaco2009 Feb]Bullcow 牡牛和牝牛(排列组合+乘法逆元+欧拉定理/费马小定理)
http://www.lydsy.com/JudgeOnline/problem.php?id=3398 以下牡牛为a,牝牛为b. 学完排列计数后试着来写这题,“至少”一词可以给我们提示,我们可以枚举 ...
- 【poj 1284】Primitive Roots(数论--欧拉函数 求原根个数){费马小定理、欧拉定理}
题意:求奇质数 P 的原根个数.若 x 是 P 的原根,那么 x^k (k=1~p-1) 模 P 为1~p-1,且互不相同. (3≤ P<65536) 解法:有费马小定理:若 p 是质数,x^( ...
- hdu_4869(费马小定理+快速幂)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4869 Turn the pokers Time Limit: 2000/1000 MS (Java/O ...
- hdu-5667 Sequence(矩阵快速幂+费马小定理+快速幂)
题目链接: Sequence Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) ...
- hdu 4704 Sum【组合数学/费马小定理/大数取模】By cellur925
首先,我们珂以抽象出S函数的模型:把n拆成k个正整数,有多少种方案? 答案是C(n-1,k-1). 然后发现我们要求的是一段连续的函数值,仔细思考,并根据组合数的性质,我们珂以发现实际上答案就是在让求 ...
- 费马小定理 x
费马小定理(Fermat Theory) 是数论中的一个重要定理,其内容为: 假如p是质数,且gcd(a,p)=1,那么 a(p-1)≡1(mod p).即:假如a是整数,p是质数,且a,p互质(即两 ...
- hdu 4704 Sum (整数和分解+快速幂+费马小定理降幂)
题意: 给n(1<n<),求(s1+s2+s3+...+sn)mod(1e9+7).其中si表示n由i个数相加而成的种数,如n=4,则s1=1,s2=3. ...
随机推荐
- Linux 下配置 iSCSI 客户端
安装客户端软件 Redhat/Centos: yum install -y iscsi-initiator-utils Debian/Ubuntu: apt-get install open-iscs ...
- GCN(Graph Convolutional Network)的简单公式推导
第一步:从前一个隐藏层到后一个隐藏层,对结点进行特征变换 第二步:对第一步进行具体实现 第三步:对邻接矩阵进行归一化(行之和为1) 邻接矩阵A的归一化,可以通过度矩阵D来实现(即通过D^-1*A来实现 ...
- lucene字典实现原理(转)
原文:https://www.cnblogs.com/LBSer/p/4119841.html 1 lucene字典 使用lucene进行查询不可避免都会使用到其提供的字典功能,即根据给定的term找 ...
- 2019年湖南省大学生计算机程序设计竞赛 (HNCPC2019) 简要题解
2019年湖南省大学生计算机程序设计竞赛 (HNCPC2019) 简要题解 update10.01 突然发现叉姐把这场的题传到牛客上了,现在大家可以有地方提交了呢. 不知道该干什么所以就来水一篇题解 ...
- AGC035
Contest Page A 唯一会做的题/kk 题目相当于要求相邻三个的异或和为\(0\). 当我们放入了三个数\(a,b,c\)时,接下来的放入顺序显然一定是\(a,b,c,a,b,c,...\) ...
- Java学习:构造方法
构造方法: 构造方法是专门用来创建对象的方法,当我们通过关键字new来创建对象时,其实就是再调用构造函数. 格式: public 类名称(参数类型 参数名称){ 方法体 } 注意事项: 构造方法的名称 ...
- SQL 查询表外键_T-Sql 2016——级联删除外键查询
SELECT fk.name AS foreign_key_name, oSub.name AS table_name, SubCol.name AS table_column, oMain.name ...
- Java自学-异常处理 Throwable
Java Throwable类 步骤 1 : Throwable Throwable是类,Exception和Error都继承了该类 所以在捕捉的时候,也可以使用Throwable进行捕捉 如图: 异 ...
- jQuery 前端复选框 全选 反选 下拉菜单联动
jQuery 页面中复选框全选.反选.下拉联动(级联) <!DOCTYPE html> <html lang="en"> <head> < ...
- Python 之 计算psnr和ssim值
基于python版的PSNR和ssim值计算 总所周知,图像质量评价的常用指标有PSNR和SSIM等,本博文是基于python版的图像numpy的float64格式和uint8格式计算两种指标值(附代 ...