线性求1~n逆元

转载自$\mathrm{j}\mathrm{z}\_\mathrm{t}\mathrm{e}\mathrm{r}\mathrm{r}\mathrm{y}$的博客

---

简述

本文要讲的是一种 $\mathcal{O}(n)$ 算法，用于求 $1\sim n$ 所有数在膜 $p$ 意义下的逆元。

算法

下面给出求 $x$ 在膜 $p$ 意义下逆元的过程：

假设 $p=ax+b$ ，则有 $ax+b\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$
式子两边同时除以 $xb$ ，得 $\frac{ax+b}{xb}\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$
整理得 $a\cdot b^{-1}+x^{-1}\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$
即 $x^{-1}\equiv -a\cdot b^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$
$\because p=ax+b$
$\therefore a=\lfloor \frac{p}{x}\rfloor ,b=p\mathrm{m}\mathrm{o}\mathrm{d}x$
将其带入即得 $x^{-1}\equiv -\lfloor \frac{p}{x}\rfloor \times (p\mathrm{m}\mathrm{o}\mathrm{d}x{)}^{-1}$
注意，为了避免出现负数，递推式中的 $-p/x$ 实际应是 $p-p/x$ 。

设 $f_i$ 表示 $i$ 在膜 $p$ 意义下的逆元，则有：
$f_i=(p-p/x)\times f_{p\mathrm{m}\mathrm{o}\mathrm{d}i}\mathrm{m}\mathrm{o}\mathrm{d}p$ （$/$为整除）

边界条件是 $f_1=1$

这样就可以求 $1\sim n$ 所有数在膜 $p$ 意义下的逆元了

---

$\mathrm{U}\mathrm{P}\mathrm{D}\mathrm{o}\mathrm{n}2019/7/28$ 补充一种自己口胡的证法

$\because f_{p\mathrm{m}\mathrm{o}\mathrm{d}i}\times p\mathrm{m}\mathrm{o}\mathrm{d}i\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$
$\therefore f_{p\mathrm{m}\mathrm{o}\mathrm{d}i}\times (p-\lfloor \frac{p}{i}\rfloor \times i)\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$
根据乘法分配律化简，得 $-\lfloor \frac{p}{i}\rfloor \times i\times f_{p\mathrm{m}\mathrm{o}\mathrm{d}i}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$
两边同乘 $f_i$ 得 $-\lfloor \frac{p}{i}\rfloor \times f_{p\mathrm{m}\mathrm{o}\mathrm{d}i}\equiv f_i(\mathrm{m}\mathrm{o}\mathrm{d}p)$
即 $f_i\equiv -p/i\times f_{p\mathrm{m}\mathrm{o}\mathrm{d}i}(\mathrm{m}\mathrm{o}\mathrm{d}p)$ （$/$为整除）
为避免出现负数， $-p/i$ 应变为 $p-p/i$ ，该变形对上式的正确性没有影响。

---

$\mathrm{U}\mathrm{P}\mathrm{D}\mathrm{o}\mathrm{n}2019/7/30$ 补充例题及程序

例题：$洛谷\mathsf{P}3811【模板】乘法逆元$

就是个模板题，下面是程序

#include<cstdio>
int n,p,inv[3000010];

char cw[1<<21],t[20]; int p1,p2=-1;
inline void flush() { fwrite(cw,1,p2+1,stdout),p2=-1; }
inline void print(int x) { //快写
    if (p2>(1<<20)) flush();
    do { t[++p1]=x%10^48; }while (x/=10);
    do { cw[++p2]=t[p1]; }while (--p1);
    cw[++p2]='\n';
}

int main() {
	scanf("%d%d",&n,&p);
	inv[1]=1; puts("1");
	for (int i=2; i<=n; ++i) {
		inv[i]=1ll*inv[p%i]*(p-p/i)%p; //递推式
		print(inv[i]);
	}
	flush();
	return 0;
}