【数学】乘法逆元进阶

2024-05-05 11:52
文章标签 进阶 数学 乘法 逆元

本文主要是介绍【数学】乘法逆元进阶,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

乘法逆元基础

可以使用扩展欧几里得算法或费马小定理求得。

乘法逆元进阶

有一种乘法逆元的线性递推算法:
显然 1 − 1 ≡ 1 ( m o d m ) 1^{-1} \equiv 1 \pmod m 111(modm)
对于 i − 1 i^{-1} i1,我们令
k = ⌊ m i ⌋ k = \lfloor \frac{m}{i} \rfloor k=im j = m m o d i j = m\bmod i j=mmodi,有 m = k i + j m= ki + j m=ki+j,即 k i + j ≡ 0 ( m o d m ) ki+j \equiv 0 \pmod m ki+j0(modm)
两边同时乘 i − 1 × j − 1 i^{-1} \times j^{-1} i1×j1 得: k j − 1 + i − 1 ≡ 0 ( m o d m ) kj^{-1}+i^{-1} \equiv 0 \pmod m kj1+i10(modm)
i − 1 ≡ − k j − 1 ( m o d m ) i^{-1} \equiv -kj^{-1} \pmod m i1kj1(modm)
代入 m = k i + j m= ki + j m=ki+j 得:

i − 1 ≡ − ⌊ m i ⌋ ( m m o d i ) − 1 ( m o d m ) i^{-1} \equiv -\lfloor\frac{m}{i}\rfloor (m \bmod i)^{-1} \pmod m i1im(mmodi)1(modm)
故利用迭代易知。

还有一个求 i ! − 1 i!^{-1} i!1 时的小技巧如下:
假设我们要求 i ∈ [ 1 , N ] i\in[1,N] i[1,N] i ! i! i! 在模 m m m 意义下的逆元。
不妨先使用费马小定理求出 N ! − 1 N!^{-1} N!1
接着我们发现 i ! − 1 ⋅ i ≡ ( i − 1 ) − 1 i!^{-1}\cdot i\equiv (i-1)^{-1} i!1i(i1)1
故可以简单地线性递推。

代码

  • i − 1 i^{-1} i1
inv[1]=1;
for (int i=2;i<=n;i++) inv[i]=(long long)(M-M/i)*inv[M%i]%M;
  • i ! − 1 i!^{-1} i!1
fac[0]=1;
for (int i=1;i<=n;i++) fac[i]=fac[i-1]*i%M;
inv[n]=pow(fac[n],M-2);
for (int i=K-1;i>=0;i--) inv[i]=inv[i+1]*(i+1)%M;

这篇关于【数学】乘法逆元进阶的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/961657

相关文章

Python进阶之Excel基本操作介绍

《Python进阶之Excel基本操作介绍》在现实中,很多工作都需要与数据打交道,Excel作为常用的数据处理工具,一直备受人们的青睐,本文主要为大家介绍了一些Python中Excel的基本操作,希望... 目录概述写入使用 xlwt使用 XlsxWriter读取修改概述在现实中,很多工作都需要与数据打交

使用C#代码计算数学表达式实例

《使用C#代码计算数学表达式实例》这段文字主要讲述了如何使用C#语言来计算数学表达式,该程序通过使用Dictionary保存变量,定义了运算符优先级,并实现了EvaluateExpression方法来... 目录C#代码计算数学表达式该方法很长,因此我将分段描述下面的代码片段显示了下一步以下代码显示该方法如

Spring Security 从入门到进阶系列教程

Spring Security 入门系列 《保护 Web 应用的安全》 《Spring-Security-入门(一):登录与退出》 《Spring-Security-入门(二):基于数据库验证》 《Spring-Security-入门(三):密码加密》 《Spring-Security-入门(四):自定义-Filter》 《Spring-Security-入门(五):在 Sprin

Java进阶13讲__第12讲_1/2

多线程、线程池 1.  线程概念 1.1  什么是线程 1.2  线程的好处 2.   创建线程的三种方式 注意事项 2.1  继承Thread类 2.1.1 认识  2.1.2  编码实现  package cn.hdc.oop10.Thread;import org.slf4j.Logger;import org.slf4j.LoggerFactory

hdu4828(卡特兰数+逆元)

这题的前几个数据分别为1,2,5,14,32......................然后确定这是个卡特兰数列 下面来介绍下卡特兰数,它的递推式为f[i+1] = f[i]*(4*n - 6)/n,其中f[2] = f[3] =1;f[4] = 2;f[5] = 14;f[6] = 32.................................. 但是这题的n太大了,所以要用到逆元,

hdu4869(逆元+求组合数)

//输入n,m,n表示翻牌的次数,m表示牌的数目,求经过n次操作后共有几种状态#include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<queue>#include<set>#include<map>#include<stdio.h>#include<stdlib.h>#includ

uva 10014 Simple calculations(数学推导)

直接按照题意来推导最后的结果就行了。 开始的时候只做到了第一个推导,第二次没有继续下去。 代码: #include<stdio.h>int main(){int T, n, i;double a, aa, sum, temp, ans;scanf("%d", &T);while(T--){scanf("%d", &n);scanf("%lf", &first);scanf

uva 10025 The ? 1 ? 2 ? ... ? n = k problem(数学)

题意是    ?  1  ?  2  ?  ...  ?  n = k 式子中给k,? 处可以填 + 也可以填 - ,问最小满足条件的n。 e.g k = 12  - 1 + 2 + 3 + 4 + 5 + 6 - 7 = 12 with n = 7。 先给证明,令 S(n) = 1 + 2 + 3 + 4 + 5 + .... + n 暴搜n,搜出当 S(n) >=

uva 11044 Searching for Nessy(小学数学)

题意是给出一个n*m的格子,求出里面有多少个不重合的九宫格。 (rows / 3) * (columns / 3) K.o 代码: #include <stdio.h>int main(){int ncase;scanf("%d", &ncase);while (ncase--){int rows, columns;scanf("%d%d", &rows, &col

【生成模型系列(初级)】嵌入(Embedding)方程——自然语言处理的数学灵魂【通俗理解】

【通俗理解】嵌入(Embedding)方程——自然语言处理的数学灵魂 关键词提炼 #嵌入方程 #自然语言处理 #词向量 #机器学习 #神经网络 #向量空间模型 #Siri #Google翻译 #AlexNet 第一节:嵌入方程的类比与核心概念【尽可能通俗】 嵌入方程可以被看作是自然语言处理中的“翻译机”,它将文本中的单词或短语转换成计算机能够理解的数学形式,即向量。 正如翻译机将一种语言