逆元模板、求1~n每个数的逆元模板

• 欧拉定理：若 $a$ 和 $n$ 互质，则 $a^{\phi (n)}\equiv 1(modn)$

• 欧拉定理的推论（费马定理）：在 $a$ 和 $p$ 互质的前提下，当 $p$ 为质数时，$a^{p-1}\equiv 1(modp)$

• 乘法逆元定义：若整数 $b$、$m$ 互质，并且 $b|a$，则存在一个整数 $x$，使得 $a/b\equiv a\ast x(modm)$，称 $x$ 为 $b$ 的模 $m$ 乘法逆元，记为 $b^{-1}(modm)$。

  $b$ 存在乘法逆元的充要条件是 $b$ 与模数 $m$ 互质。当模数 $m$ 为质数时，$b^{m-2}$ 即为 $b$ 的乘法逆元。

问题转换为快速幂的计算：
快速幂模板

---

逆元模板

#include<bits/stdc++.h>
using namespace std;

typedef long long LL;

LL ksm (LL a, LL k, LL p) {
	LL res = 1;
	while (k) {
		if (k & 1) res = res * a % p;
		a = a * a % p;
		k >>= 1;
	}
	return res;
}

int main()
{
	LL n, p;
	cin >> n >> p;
	
	cout << ksm (n, p-2, p) << endl;	

	return 0;
}

---

计算 $1\sim n$ 的乘法逆元模板：

题目链接

如果直接遍历 $1n$ 每个数都使用快速幂计算逆元会超时 $O(nlogp)$，需要转换为线性计算。

因为 $p$ 是质数，所以可以将 $p$ 表示为 $k\ast i+r$ ，$(1<r<i<p)$，其含义为 $k$ 是 $p/i$ 的商，$r$ 是余数。

$$k\ast i+r\equiv 0(modp)$$
等式两侧同时乘上 $i^{-1}$ 和 $r^{-1}$，得：

$$\begin{gathered} k\ast r^{-1}+i^{-1}\equiv 0(modp) \\ {i}^{-1}\equiv -k\ast r^{-1}(modp) \\ {i}^{-1}\equiv -\lfloor \frac{p}{i}\rfloor \ast (pmodi{)}^{-1}(modp) \end{gathered}$$

可以看出 $i$ 的逆元可以由 $pmodi$ 的逆元递推得到。

需要将 $1$ 的逆元初始化为 $1$，因为 $1^{-1}=1^{p-2}\%p=1$。

注意负数取模，$-\lfloor \frac{p}{i}\rfloor$ 对 $p$ 取模需要先将负数 $+p$ 使其变为整数再进行取模操作。

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 3e6+10;

LL inv[N], n, p;

int main()
{
	scanf ("%lld%lld", &n, &p); // cin 会TLE 
	
	inv[1] = 1; puts ("1");
	for (int i = 2;i <= n;i ++) 
		inv[i] = (p - p/i) * inv[p % i] % p,
		printf ("%lld\n", inv[i]); // cout 会TLE 

	return 0;
}

---

注意，数论题一般都要开LL。