费马小定理:任意整数a, a^p≡a (mod p),p是质数。

定理:1.对a/b≡a*x (mod n),则称x 为b的乘法逆元,记为b^-1。

2.当模数p为质数,且b与p互质(即b不为p的倍数)时,b的乘法逆元为b^(p-2)。

对2的简单证明:

b*b^-1≡1 (mod p)。①

由费马小定理,b^ (p-1)≡1 (mod p),bb^ -1b^ (p-1)≡1 (mod p),即 b*b^(p-2)≡ 1 (mod p)。对比①,所以得到2的结论。

所以,对等比数列和公式 a*(q^n-1)/(q-1),就可以用快速幂求得分母的逆元,然后乘上分子的模再取模。

对等比数列 1+a+a^ 2 +a^ 3+……+a^ n,结果对质数p取模,分两类情况讨论:

1.a-1与p不互质,那么上面的定理2就不能使用了。但是可以发现,a%p≡1 (mod p)。所以原式就变成了 n+1。
2.若 a-1与p不互质,则根据等比数列公式,原式就为1*[a^ (n+1)-1]%p *[(a-1)^(p-2)]%p。
这样,就在极短的时间内完成了等比数列的求和,主要的时间花费都集中在快速幂部分。

代码如下:

int qpow(int a,long long b,int mod)
{int ans=1;while(b>0){if(b&1)ans=(long long)ans*a%mod;a=(long long)a*a%mod;b=b>>1;}return ans;
}
int main()
{int a,n,p;cin>>a>>n>>p;int sum;if((a-1)%p==0){sum=n+1;}else{sum=qpow(a,(long long)n+1,p);sum=(sum-1+p)%p;sum=(long long)sum*qpow(a-1,p-2,p)%p;}cout<<sum;
}

费马小定理与等比数列的求和(MOD 质数)相关推荐

  1. 逆元+费马小定理+扩展欧几里得

    逆元:(即是逆元素)逆元素是指一个可以取消另一给定元素运算的元素. 在一个代数系统(S,*)中,存在单位元素e,如果对S内的元素a存在a^-1 * a = e,则将 a^-1称为a 的左逆元. 同理若 ...

  2. 数论杂谈(欧拉定理与费马小定理结论与应用)

    文章目录 欧拉定理: 欧拉定理性质: 扩展欧拉定理: 费马小定理: 指数循环节 费马大定理 逆元: 例题 原根 定义: 原根存在条件 例题 快速幂 代码 矩阵快速幂 原理: 代码: 欧拉定理: aφ( ...

  3. 循环小数与费马小定理

    循环小数与费马小定理 17/05/29 22:30:51 | Snakes 背景 题目出自之前亮灯问题.杨辉三角与Sierpinski三角形提及的生日题中的第三.四.五题. 题目 第三题 证明:对于任 ...

  4. [bzoj1951] [Sdoi2010]古代猪文 费马小定理+Lucas定理+CRT

    Description "在那山的那边海的那边有一群小肥猪.他们活泼又聪明,他们调皮又灵敏.他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心--" --选自猪王国民歌 很久 ...

  5. 费马小定理【模板例题】

    费马小定理: 如果p是一个质数,而整数a不是p的倍数, 则有a(p-1)≡1(mod p). 即:假如a是整数,p是质数,且a,p互质(即两者只有一个公约数1),那么a的(p-1)次方除以p的余数恒等 ...

  6. 【BZOJ】3398: [Usaco2009 Feb]Bullcow 牡牛和牝牛(排列组合+乘法逆元+欧拉定理/费马小定理)...

    http://www.lydsy.com/JudgeOnline/problem.php?id=3398 以下牡牛为a,牝牛为b. 学完排列计数后试着来写这题,"至少"一词可以给我 ...

  7. 快速幂-Jxc军训 (费马小定理+求概率)(2021-08-17)

    Jxc军训 题目 在文某路学车中学高一新生军训中,Jxc正站在太阳下站着军姿,对于这样的酷热的阳光,Jxc 表示非常不爽. Jxc将天空看做一个nn的矩阵,此时天上有m朵云,这些云会随机分布在m个不同 ...

  8. 同余定理 逆元 中国剩余定理 费马小定理

    同余定理 同余定理是数论中的重要概念.给定一个正整数m,如果两个整数a和b满足(a-b)能够被m整除,即(a-b)/m得到一个整数,那么就称整数a与b对模m同余,记作a≡b(mod m). 两个整数a ...

  9. 逆元的几种求法(扩展欧几里得,费马小定理或欧拉定理,特例,打表等)

    乘法逆元 对于缩系中的元素,每个数a均有唯一的与之对应的乘法逆元x,使得ax≡1(mod n) 一个数有逆元的充分必要条件是gcd(a,n)=1,此时逆元唯一存在 逆元的含义:模n意义下,1个数a如果 ...

最新文章

  1. 样本方差除以n-1而不是n的原因
  2. mysql 2013错误解决
  3. jaquery ready
  4. Invalid packaging for parent POM com.cloud:springcloud_demo:1.0-SNAPSHOT, must be “pom“ but is “jar“
  5. 任正非:华为鸿蒙将比安卓快 60%;小米已官方回应侵权偷图事件;博通正考虑收购软件公司,其中Tibco和赛门铁克为重点目标……...
  6. OpenSSL 宣布将开源许可证更改为 Apache 2.0
  7. 为什么使用HashMap需要重写hashcode和equals方法_不同时重写equals和hashCode又会怎样?听听过来人的经验...
  8. MQTT工作笔记0009---订阅主题和订阅确认
  9. 压力测试jmeter入门教程
  10. 已解决-Windows10没有windows照片查看器-Windows10打开照片是黑底的
  11. 朱晔的互联网架构实践心得S2E4:小议微服务的各种玩法(古典、SOA、传统、K8S、ServiceMesh)...
  12. mysql 多字段,多条件查找
  13. Excel2007将汉字转换为拼音首字母
  14. 深信服校园招聘c/c++软件开发A卷
  15. cesium两种方式鼠标移入移出实体显示提示框
  16. L3-020 至多删三个字符 (30 分)
  17. Win10右键菜单“发送到”添加/删除选项的方法
  18. 关闭休眠扩大c盘空间
  19. Appium移动自动化
  20. 第八章 对立统一——异步时钟同步化

热门文章

  1. 计算机专业pro配置,计算机专业,这种小米Pro配置可以吗?
  2. python TCP服务器v1.8 - PyQt5登录界面美化+淡入淡出
  3. 【网络】从日常开发说起,浅谈HTTP协议是做什么的
  4. (免费分享)基于springboot校园实习管理系统
  5. EC张星亮:2016年SaaS进入价值回归期
  6. 如何把视频转换成mp3格式?
  7. java使用国标方式取流,一种基于JAIN-SIP的国标28181平台分布式集群实现系统的制作方法...
  8. Collapsing margins
  9. Web前端开发 table标签相关用法和属性
  10. java中aes明文长度_AES密文与明文长度的关系