C++ 线性筛逆元

分析

$O(n)$筛$1$到$n$所有数模$p$的逆元：

• 初值Inv[1]=1；
• 设$p=ki+r$，则$ki+r\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}p$；
• 两边同乘$i^{-1}{r}^{-1}$得$k{r}^{-1}+i^{-1}\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}p$；
• 移项：$i^{-1}\equiv -k{r}^{-1}$；
• $k=\lfloor \frac{p}{i}\rfloor$，所以$-k\equiv p-\lfloor \frac{p}{i}\rfloor \mathrm{m}\mathrm{o}\mathrm{d}p$；
• 所以Inv[i]=(p-p/i)*Inv[p%i]%p。

代码

#define MOD 998244353
#define MAXC 10000000
int Inv[MAXC+5];
void GetInv(int n){
    Inv[1]=1;
    for(int i=2;i<=n;i++)
        Inv[i]=(LL)Inv[MOD%i]*(MOD-MOD/i)%MOD;
}