[数学/数论] 普通求逆

网上的求逆方法普遍长这样

a[1]=1;
for(int i=2;i<=n;++i) a[i]=(mod-(mod/i))%mod*a[mod%i]%mod;

下面介绍一种简单且万能的

pr[1]=1;
for(int i=2;i<=n;++i) pr[i]=pr[i-1]i%mod;
b[n]=qpow(pr[n],mod-2);
for(int i=n-1;i;–i) b[i]=b[i+1](i+1)%mod;
b[1]=1;
for(int i=2;i<=n;++i) b[i]=b[i]*pr[i-1]%mod;

讲解开始

逆元的定义

在$(\mathrm{m}\mathrm{o}\mathrm{d}p)$意义下，$i\times i^{-1}\equiv 1$，则$i$与$i^{-1}$互为逆元

使用逆元的意义：为了保证精度，维护整数问题（即数论问题）的过程中，不能除，于是利用与原数积在模意义下为1的数，来替代除法

$a_n=n$的逆元求法

常用，但是常数有点大

推导：
$$\begin{array}{rl}& 模p意义下,求i的逆元\\ & p=ki+r(k=\lfloor \frac{p}{i}\rfloor ,r=p\mathrm{m}\mathrm{o}\mathrm{d}i)\\ & ki+r\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)\\ & i^{-1}{r}^{-1}(ki+r)\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)\\ & i^{-1}+k{r}^{-1}\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)\\ & i^{-1}\equiv -k{r}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)\\ & 由于r=p\mathrm{m}\mathrm{o}\mathrm{d}i,则r<i,求i时r已求出,可以递归求解\\ & 边界1^{-1}=1\end{array}$$

Code

a[1]=1;
for(int i=2;i<=n;++i) a[i]=(mod-(mod/i))%mod*a[mod%i]%mod;

任意数列的逆元求法

常数又小，适用范围又广
不学他学谁？

推导
对 于 { a n } ( 通 常 是 a n = n 的 数 列 ) , 记 p r i = ∏ j = 1 i a j , 第 一 遍 , i n v n = p r n p − 2 ( 即 p r n 的 逆 元 ) ≡ 1 ∏ i = 1 n a i i n v i = i n v i + 1 × a i + 1 , 递 推 下 去 第 二 遍 , i n v i = i n v i × p r i − 1 时 间 复 杂 度 O ( N ) \begin{aligned} &对于\{a_n\}(通常是a_n=n的数列),记pr_i=\prod_{j=1}^ia_j,\\ &第一遍,inv_n=pr_n^{p-2}(即pr_n的逆元)\equiv\frac{1}{\prod_{i=1}^na_i}\\ &inv_i=inv_{i+1}\times a_{i+1},递推下去\\ &第二遍,inv_i=inv_i\times pr_{i-1} \\\\&时间复杂度O(N) \end{aligned} ​对于{an​}(通常是an​=n的数列),记pri​=j=1∏i​aj​,第一遍,invn​=prnp−2​(即prn​的逆元)≡∏i=1n​ai​1​invi​=invi+1​×ai+1​,递推下去第二遍,invi​=invi​×pri−1​时间复杂度O(N)​

Code

pr[1]=1;
for(int i=2;i<=n;++i) pr[i]=pr[i-1]*i%mod;
b[n]=qpow(pr[n],mod-2);
for(int i=n-1;i;--i) b[i]=b[i+1]*(i+1)%mod;
b[1]=1;
for(int i=2;i<=n;++i) b[i]=b[i]*pr[i-1]%mod;

Tips

1. 注意边界：
   法一中inv[1]要赋初值
   法二中最后一层循环由于pr[0]=0，无法更新inv[1]，那么手动赋1
2. 注意求逆元的范围
   想清楚就行了