【ybt高效进阶6-3-3】线性求逆元（数学）

线性求逆元

题目链接：ybt高效进阶6-3-3

题目大意

给你 n,p，要你求 1~n 每个数在模 p 意义下的逆元。

思路

线性求逆元的公式就是：
$in{v}_i=((p-p/i\ast in{v}_{p\mathrm{m}\mathrm{o}\mathrm{d}i})\mathrm{m}\mathrm{o}\mathrm{d}p+p)\mathrm{m}\mathrm{o}\mathrm{d}p$

这里给一下证明：
因为是线性，所以假设我们求 $i$，那 $1\sim i-1$ 的我们都求出来了。
那 $p$ 可以表示成 $\lfloor \frac{p}{i}\rfloor \ast i+r(0⩽r<i)$，所以就有
$\lfloor \frac{p}{i}\rfloor \ast i+r\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$

然后两边都乘 $in{v}_i\ast in{v}_r$，就有 $\lfloor \frac{p}{i}\rfloor \ast in{v}_r+in{v}_i\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$
$in{v}_i\equiv -\lfloor \frac{p}{i}\rfloor \ast in{v}_r(\mathrm{m}\mathrm{o}\mathrm{d}p)$

然后由于 $(0⩽r<i)$，$in{v}_r$ 是已知的，所以就得出了 $in{v}_i$。

代码

#include<cstdio>

using namespace std;

int n, p, inv[3000001];

int main() {
	scanf("%d %d", &n, &p);
	
	inv[1] = 1;
	for (int i = 2; i <= n; i++)//不开 long long 见祖宗
		inv[i] = ((1ll * p - 1ll * p / i * inv[p % i]) % p + p) % p;
	
	for (int i = 1; i <= n; i++) printf("%d\n", inv[i]);
	
	return 0;
}